_RL  ptop,p,undef,kappa,getcon

      pk    =    p**kappa

      subroutine PRESTOPRES (qprs,qinp,pkz,pksrf,ptop,p,im,jm,lm,myThid)

              p = plevs1(k)

              call PRESTOPRES(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,

              p = plevs2(k)

              call PRESTOPRES(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,

              p = plevs3(k)

              call PRESTOPRES(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,

      _RL p

      _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)

      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 )

      _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(sNx,sNy,Nsx,Nsy)

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

     . xlons,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 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,

      _RL p(sNx,sNy,Nsx,Nsy)

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

            p = locPres(i,j)*SItoBar

            p2 = p*p

     &         p*( eosJMDCKP(1)

     &     + p*s*( eosJMDCKP(5)

     &      + p*s3o2*eosJMDCKP(8)

      _RL p(sNx,sNy,Nsx,Nsy)

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

     . xlons,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,uphy,vphy,thphy,sphy,qq,pkl,dpres,

      _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,