_RL dp(im2,jm2,Nrphys,Nbi,Nbj)

     .             moistt(i,j,L,bi,bj)*pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)

     .             turbt(i,j,L,bi,bj)*pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)

     .             *pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)

     .             swdt(i,j,L,bi,bj)*pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)

      subroutine FIZHI_STEP_DIAG(myid,p,uphy,vphy,thphy,sphy,qq,pk,dp,

      qval     = pq(i  ,j,L)*d2p(i  ,j)*dp(i,j,L)

      array(1) = pq(ip1,j,L)*d2p(ip1,j)*dp(i,j,L)

      array(2) = pq(im1,j,L)*d2p(im1,j)*dp(i,j,L)

      array(3) = pq(i,j+1,L)*d2p(i,j+1)*dp(i,j,L)

      array(4) = pq(i,j-1,L)*d2p(i,j-1)*dp(i,j,L)

      array(5) = pq(i,j,L-1)*d2p(i,j)*dp(i,j,L)

      array(6) = pq(i,j,L+1)*d2p(i,j)*dp(i,j,L)

     .                                 pq(i,j,L)*dp(i,j,L) /dp(i,j,L-1)

     .                                 pq(i,j,L)*dp(i,j,L) /dp(i,j,L+1)

      _RL dp(idim1:idim2,jdim1:jdim2,ldim,Nsx,Nsy)

       ddsig = dp(i,j,LM1,bi,bj)/dp(i,j,L,bi,bj)                                     

      subroutine TRACER_FILL ( pq,im,jm,lm,dlam,dphi,dp)

     .     jm1,jm2,bi,bj,dp,qz)

      _RL    pq(im,jm,lm),dlam(im),dphi(jm),dp(im,jm,lm)

      _RL dp,xx,p1,dwe,dpe,a1,b1,fk1,a2,b2,fk2

         dp           = pl(i,j,k+1)-pl(i,j,k)

         dh2o (i,j,k) = 1.02*wa(i,j,k)*dp+1.e-10

         do3  (i,j,k) = 476.*oa(i,j,k)*dp+1.e-10

         dco2 (i,j,k) = 789.*co2*dp+1.e-10

         dch4 (i,j,k) = 789.*ch4(k)*dp+1.e-10

         dn2o (i,j,k) = 789.*n2o(k)*dp+1.e-10

         df11 (i,j,k) = 789.*cfc11*dp+1.e-10

         df12 (i,j,k) = 789.*cfc12*dp+1.e-10

         df22 (i,j,k) = 789.*cfc22*dp+1.e-10

         xx=pa(i,j,k)*0.001618*wa(i,j,k)*wa(i,j,k)*dp

             dp          =pl(i,j,k+1)-pl(i,j,k)

             cwp(i,j,k,1)=1.02*10000.0*cwc(i,j,k,1)*dp

             cwp(i,j,k,2)=1.02*10000.0*cwc(i,j,k,2)*dp

             cwp(i,j,k,3)=1.02*10000.0*cwc(i,j,k,3)*dp

          call CAL_CONVDATE( previousdate,yp,mp,dp,sp,lp,wp,mythid )

      integer di,df,dp,dn