if( q(i).lt.minval ) minval = q(i)

      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 )

      function minval (q,im)

      _RL q(im), minval

      subroutine SRCLOUDS (th,q,plk,pl,plke,cloud,cldwat,irun,irise,

      _RL    q(irun,irise)

      rh    = q(i,L)/qs(i,L)

      call CTEI ( th,q,cloud,cldwat,pl,plk,plke,irun,irise )

      subroutine CTEI ( th,q,cldfrc,cldwat,pl,plk,plke,im,lm )

      _RL  th(im,lm),q(im,lm),plke(im,lm+1),cldwat(im,lm)

          dq = q(i,L)+cldwat(i,L)-q(i,L-1)-cldwat(i,L-1)

      _RL q(sNx,sNy,Nrphys,ntracer)

       q(i,j,L,1) = sphy(i,j,L,bi,bj)

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

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

      subroutine MAXVAL1 (q,im,qval,qmax,imax)

      _RL   q(im), qmax, qval

      if( q(i).gt.qmax ) then

      qmax = q(i)

      if(2.eq.1 )then

      if(2.eq.1 )then

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

      call RAYLEIGH(myid,pl,pkl,pkht,surfz,u,v,t,q,im,jm,lm,

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

      _RL q(im,jm,lm,ntracer)

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

      _RL q(sNx,sNy,Nrphys,ntracer)

      _RL q(sNx,sNy,Nrphys,ntracer)

       q(i,j,L,1) = sphy(i,j,L,bi,bj)

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

      _RL del_P ,del_th, th, q

      q      = del_th

      th     = th + (1 - 1/sqrt(two))*(del_th - q)

      q      = (two-sqrt(two))*del_th + (-two+three/sqrt(two))*q

      th     = th + (1 + 1/sqrt(two))*(del_th - q)

      q      = (two + sqrt(two))*del_th + (-two-three/sqrt(two))*q

      SW_PTMP     = th + (del_th - two*q)/(two*three)