#include "dims.h"
      SUBROUTINE CMPAR
      use cons_module,only: len1,len3
      implicit none
C     ****     ADVANCE AR COMPOSITION BY ONE TIME STEP
#include "params.h"
#include "fcom.h"
#include "vscr.h"
#include "index.h"
#include "buff.h"
#include "mwt.h"
#include "diffk.h"
#include "phys.h"
#include "cmpdat.h"
!
! Local:
      real,parameter :: ALFAAR=0.17, PSARB=0.1146E-1
      integer :: ibnd,ibndb,i
      real :: xyar
      real :: PHIAR(3) = (/1.042,1.509,1.176/)
!
!     ****     UPPER BOUNDARY - ZERO FLUX
!     ****     LOWER BOUNDARY - CONSTANT VALUE
      IBND = 1
      IBNDB = 0
      DO 1 I=1,LEN1
      T4(I)=0.
      T1(I)=0.
      T2(I)=1.
      T3(I)=-SQRT(WAR(1)*WAR(2))
    1 CONTINUE
C     ****     SOURCES ARE ZERO
      DO 2 I=1,LEN3
      S1(I,1)=0
      S2(I,1)=0.
    2 CONTINUE
      XYAR = 1.E-10
      CALL MINOR(NPAR,NPARNM,RMAR,PHIAR,ALFAAR,IBND,IBNDB,WAR,XYAR,
     1 NPDHAR,difkk)
C*************************************************
C     NPARK = NJNP + NPAR - 1
C     DO K = 1,KMAXP1,4
C       DO I = 1,IMAX
C         PLOT(I,J,K) = F(I+2,NPARK + K)
C       ENDDO
C     ENDDO
C     IF(J.EQ.JMAX)THEN
C       DO K = 1,KMAXP1,4
C         CALL EZCNTR(PLOT(1,1,K),IMAX,JMAX)
C         WRITE(6,*)'K = ', K
C       ENDDO
C     ENDIF
C*************************************************
      RETURN
      END