C&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
      SUBROUTINE DIFCOF(LMHK,GM,GH,EL,Q2,Z,AKM,AKH)
C     ******************************************************************
C     *                                                                *
C     *                LEVEL 2.5 DIFFUSION COEFFICIENTS                *
C     *                                                                *
C     ******************************************************************
C-----------------------------------------------------------------------
      INCLUDE "parmeta"
      INCLUDE "mpp.h"
#include "sp.h"
C-----------------------------------------------------------------------
                             P A R A M E T E R
     &(LP1=LM+1,LM1=LM-1)
C-----------------------------------------------------------------------
                             P A R A M E T E R
C-----------------------------------------------------------------------
     &(G=9.8,BETA=1./270.,BTG=BETA*G
     &,PRT=1.0,GAM1=0.2222222222222222222
     &,A1=0.659888514560862645,A2=0.6574209922667784586
     &,B1=11.87799326209552761,B2=7.226971804046074028
     &,C1=0.000830955950095854396)
C-----------------------------------------------------------------------
CUN     &(G=9.8,BETA=1./270.,BTG=BETA*G
CUN     &,PRT=1.0,GAM1=0.2222222222222222222
CUN     &,A1=0.3310949523016403346,A2=0.8273378704055731278
CUN     &,B1=5.959709141429526024,B2=3.626088092074591135
CUN     &,C1=-0.3330651924968952113)
C-----------------------------------------------------------------------
CMY     &(G=9.8,BETA=1./270.,BTG=BETA*G
CMY     &,PRT=0.8,GAM1=0.2222222222222222222
CMY     &,A1=0.9222222350809054114,A2=0.7350190142719400952
CMY     &,B1=16.60000023145629741,B2=10.10000014082581951
CMY     &,C1=0.0805318118080613468)
C-----------------------------------------------------------------------
                             P A R A M E T E R
C--------------COEFFICIENTS FOR THE SM AND SH DETERMINANTS--------------
     &(BSMH=-3.*A1*A2*(3.*A2+3.*B2*C1+12.*A1*C1-B2)*BTG
     &,CESM=    A1*(1.-3.*C1)
C
     &,BSHM=18.*A1*A1*A2*C1
     &,BSHH= 9.*A1*A2*A2*BTG
     &,CESH=    A2
C--------------COEFFICIENTS OF THE TERMS IN THE DENOMINATOR-------------
     &,ADNM=18.*A1*A1*A2*(B2-3.*A2)*BTG
     &,ADNH= 9.*A1*A2*A2*(12.*A1+3.*B2)*BTG*BTG
     &,BDNM= 6.*A1*A1
     &,BDNH= 3.*A2*(7.*A1+B2)*BTG)
C-----------------------------------------------------------------------
                             D I M E N S I O N
     & Q2    (LM)
                             D I M E N S I O N
     & GM    (LM1),GH    (LM1),EL    (LM1)
     &,AKM   (LM1),AKH   (LM1)
     &,Z     (LP1)
C-----------------------------------------------------------------------
C***********************************************************************
      LMHM=LMHK-1
      LMHP=LMHK+1
          DO 100 L=1,LMHM
      ELL=EL(L)
C-----------------------------------------------------------------------
      ELOQ2=ELL*ELL/Q2(L)
      ELOQ4=ELOQ2*ELOQ2
C-----------------------------------------------------------------------
      GML=GM(L)
      GHL=GH(L)
C--------------COEFFICIENTS OF THE TERMS IN THE DENOMINATOR-------------
      ADEN=(ADNM*GML+ADNH*GHL)*GHL
      BDEN= BDNM*GML+BDNH*GHL
      CDEN= 1.
C--------------COEFFICIENTS FOR THE SM DETERMINANT----------------------
      BESM=BSMH*GHL
C--------------COEFFICIENTS FOR THE SH DETERMINANT----------------------
      BESH=BSHM*GML+BSHH*GHL
C--------------1./DENOMINATOR-------------------------------------------
      RDEN=1./(ADEN*ELOQ4+BDEN*ELOQ2+CDEN)
C--------------SM AND SH------------------------------------------------
      ESM=(BESM*ELOQ2+CESM)*RDEN
      ESH=(BESH*ELOQ2+CESH)*RDEN
C--------------DIFFUSION COEFFICIENTS-----------------------------------
      RDZ=2./(Z(L)-Z(L+2))
      Q1L=SQRT(Q2(L))
      ELQDZ=ELL*Q1L*RDZ
      AKM(L)=ELQDZ*ESM
      AKH(L)=ELQDZ*ESH
C-----------------------------------------------------------------------
  100 CONTINUE
C-----------------------------------------------------------------------
                           RETURN
			   END