_RL p(sNx,sNy,nSx,nSy)

      COMMON /saver/ p,pl,ple,dpres,pkle,pkl

       p(i,j,bi,bj) = pephy(i,j,Nrphys+1,bi,bj)/100.

     & p(1,1,bi,bj),u,v,t,q,pl(1,1,1,bi,bj),ple(1,1,1,bi,bj),

     &   swdt(i,j,L,bi,bj)*radswt(i,j,bi,bj) )*p0kappa ) / p(i,j,bi,bj)

     &                                    /p(i,j,bi,bj)

      call FIZHI_STEP_DIAG(myid,p,u,v,t,q,qq,pkl,dpres,

     &   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,

     &   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,

     &  ptop,p,u,v,t,q,ntracer,ptracer,pl,ple,dpres,pkht,pkl,

      call GWDRAG (myid,p,pl,ple,dpres,pkl,u,v,t,q,phis_var,

     &    p,pl,ple,dpres,pkht,pkl,u,v,t,q,bi,bj,ntracer,ptracer,qq,

      _RL p(im,jm)

     & p,u,v,t,q,pl,ple,dpres,pkht,pkl,surfz,fracland,landtype,radswt,

      _RL p(lma)

      data p  / 93.81, 111.65, 133.46, 157.80, 186.51,

         if( pl(i).ge.p(L)   .and.

     .       pl(i).le.p(L+1)) then

         temp = ( pl(i)-p(L) )/( p(L+1)-p(L) )

      if( pl(i).lt.p(1)   ) acr(i) = a(1)

      if( pl(i).gt.p(lma) ) acr(i) = a(lma)

      _RL    getcon,cp,alhl,elocp,cpoel,t,p,s,qs,dqsdt,dq

      p = pl(i,L)

      call QSAT ( t,p,qs,dqsdt,.true. )

      _RL p(im2,jm2,Nbi,Nbj)

      pinv(i,j) = 1.0 / p(i,j,bi,bj)

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

      call DIAGNOSTICS_FILL(p(1,1,bi,bj),'PS      ',0,1,3,bi,bj,myid)

      _RL pk(istrip,nlay), qq(istrip,nlay),   p(istrip,nlay)

       call STRIP2TILE(plz,igrd,p,nchp,ijall,istrip,nlay,nn)

      call TRBFLX(nn,th,thv,sh,u,v,qq,p,pe,pk,pke,dpstr,stwatr,stz0,

      subroutine VQSAT (tt,p,q,dqdt,ldqdt,n)

      _RL tt(n), p(n), q(n), dqdt(n)

      call QSAT ( tt(i),p(i),q(i),dqdt(i),ldqdt )

         p       = locPres(i,j)*SItodBar

     &   + p*(     teos(13) + 2.*teos(14)*ct  + sa*2.*teos(16)

     &        + p*(teos(18) + 2.*teos(19)*ct))

     &   + p*(teos(38) + ct*(2.*teos(39) + 3.*teos(40)*ct)

     &   + p*(teos(44) + 2.*teos(45)*ct + teos(46)*sa

     &   + p*teos(48) ))

      _RL ct, sa, sqrtsa, p

      _RL ct, sa, sqrtsa, p

         p       = locPres(i,j)*SItodBar

     &   + p*(teos(15) + teos(16)*ct + p*teos(20))

     &   + p*(teos(41) + teos(42)*ct + p*teos(46))

       p       = pLoc*SItodBar

     &   + p*(teos(12) + ct*(teos(13) + teos(14)*ct)

     &   + p*(teos(17) + ct*(teos(18) + teos(19)*ct) + teos(20)*sa))

     %   + p*(teos(37) + ct*(teos(38) + ct*(teos(39) + teos(40)*ct))

     %   + p*(teos(43) + ct*(teos(44) + teos(45)*ct + teos(46)*sa)

     %   + p*(teos(47) + teos(48)*ct)))

      _RL t, t2, t3, t4, s, s3o2, p, p2

            p = locPres(i,j)*SItoBar

            p2 = p*p

            bMpres =  p * ( eosP1 +

     &             + p*s* ( eosP5 +

     &             + p*s3o2*eosP8

     &         p*( eosJMDCKP(1)

     &     + p*s*( eosJMDCKP(5)

     &      + p*s3o2*eosJMDCKP(8)

      _RL ct, sa, sqrtsa, p

        p       = locPres(i,j)*SItodBar

     &   + p*(teos(12) + ct*(teos(13) + teos(14)*ct)

     &   + p*(teos(17) + ct*(teos(18) + teos(19)*ct) + teos(20)*sa))

     &   + p*(teos(37) + ct*(teos(38) + ct*(teos(39) + teos(40)*ct))

     &   + p*(teos(43) + ct*(teos(44) + teos(45)*ct + teos(46)*sa)

     &   + p*(teos(47) + teos(48)*ct)))

      _RL ct, sa, sqrtsa, p

      INTEGER bi, bj, i, j, xnode, ynode, xq, yq, m, n, p, kx, ky

      _RL p(sNx,sNy,nSx,nSy)