$        phis_var(i,j)) )/grav )

     .            plstr,plestr,dpresstr,grav,rgas,cp,

        tmpdiag(i,j) = dragx(i,j) + sumu(i,j)*pz(i,j)/grav*100

        tmpdiag(i,j) = dragy(i,j) + sumu(i,j)*pz(i,j)/grav*100

     .                  grav,rgas,cp,irun,Lm,nthin,nbase,lstar )

      _RL grav, rgas, cp

      agrav  = 1.0/grav

      gocp   = grav/cp

        vai2    = 2.0*grav/( t(i,Lm+1-L)+t(i,Lm+2-L) )

        vai2  = 2.0*grav/(t(i,Lm+1-L)+t(i,Lm+2-L))

       coef = -grav*ps(i)/dpres(i,Lm+1-L)

       coef = grav*ps(i)/(ple(i,Lm+1-nthin(i))-ple(i,Lm+1))

      _RL grav, rgas, cp, cpinv, lstar

      grav  = getcon('GRAVITY')

     *,      cp,grav,rkappa,alhl,rhfrac,rasmax )

      _RL cp,grav,rkappa,alhl,rhfrac(lng),rasmax

       tx1(i) = cmb2pa *  (prs(i,ib+1) - prs(i,ib))/(grav*dt)

         tx1(i)   = 100*(prs(i,L+1)-prs(i,L))/grav

      _RL tmstp,tminv,sday,grav,alhl,cp,elocp,gamfac

      _RL frac, cp,  alhl, rkap, grav, p00, crtmsf

     8, cp,grav,rkappa,alhl,rhfrac(num),rasmax )

     .  CLSBTH,TMSTP,one,cp,grav,alhl,gamfac,cldlz,rhcrit,offset,alpha)

      grav  = getcon('GRAVITY')

      _RL      getcon,grav,cp,undef

      alf = grav*(ple(i,Lm+1)-ptop)/(cp*dpstrip(i,L)*100)

     4 irun,nlev,nltop,NYMD,NHMS,grav,cp,rgas,faceps,virtcon,undef,

      _RL grav,cp,rgas,faceps,virtcon,undef

      _RL faceps,grav,caltoj,virtcon,getcon

       grav     = getcon('GRAVITY')

     4 stq10m,istrip,nlay,nltop,nymd,nhms,grav,cp,rgas,faceps,virtcon,

         tempdiag(i,j) = topoZ(i,j,bi,bj)*grav

       _RL grav, getcon

       grav = getcon('GRAVITY')