CALL SHTORH (-1,NGP,T0,PSA,1. _d 0,Q0,WORK,QSAT0,myThid)

        DENVV(J,3)=(PRD*PSA(J)/T0(J,1))*SPEED0(J)

      CALL SHTORH (0,NGP,TLAND,PSA,1. _d 0,QDUMMY,RDUMMY,QSAT0(1,1),myThid)

      CALL SHTORH (0,NGP,TSEA ,PSA,1. _d 0,QDUMMY,RDUMMY,QSAT0(1,2),myThid)

      _RL PSA(NGP), UA(NGP,NLEV), VA(NGP,NLEV), TA(NGP,NLEV),

      SUBROUTINE SUFLUX (PSA,UA,VA,TA,QA,RH,QSAT,Vsurfsq,PHI,

      _RL PSA(NGP), QA(NGP,NLEV), QSAT(NGP,NLEV)

      SUBROUTINE LSCOND (PSA,QA,QSAT,

        PRECLS(J) = PRECLS(J)*PSA(J)

        FMASS=FM0(J)*PSA(J)*(QA(J,K)-RHBL*QSAT(J,K))/DQSAT

          ENMASS=ENTR(J,K)*PSA(J)*FMASS

      _RL PSA(NGP), Ta(NGP,NLEV), QA(NGP,NLEV), QSAT(NGP,NLEV)

      SUBROUTINE CONVMF (PSA,Ta,QA,QSAT,

      _RL PSA(NGP), QA(NGP,NLEV), RH(NGP,NLEV)

        TAU(J,1)=EXP(-ABSSW*PSA(J)*DSIG(1))

        TAUOZ(J)=EXP(-EPSSW*OZONE(J)*PSA(J))

     *           +ABCSW*CLOUDC(J)*QA(J,NL1(J)))*PSA(J)*DSIG(K))

     &                           *PSA(J)*DSIG(NLEVxy(J,myThid)))

        TAU(J,1)=EXP(-ABSLW*PSA(J)*DSIG(1))

     *           +ABCLW*CLOUDC(J)*QA(J,NL1(J)))*PSA(J)*DSIG(K))

     &   EXP(-(ABSLW+ABWLW*QA(J,NLEVxy(J,myThid)))*PSA(J)

      SUBROUTINE RADSW (PSA,QA,RH,

     I                   PSA,TA,QA,RH,ThA,Vsurf2,WVS,CLAT,FOROG,

        T1(J) = ThA(J,Ktmp)*(PSA(J)**kappa)

      CALL SHTORH (-1,NGP,T0, PSA, 1. _d 0, Q0, tmpRH, QSAT0, myThid)

        DENVV(J)=(PRD*PSA(J)/T0(J))*SPEED0(J)

        dTskin(J)=CTDAY*CLAT(J)*SSR(J)*PSA(J)

      _RL  PSA(NGP), TA(NGP,NLEV), QA(NGP,NLEV), RH(NGP,NLEV)

     I                   PSA, FMASK, EMISloc,

      CALL SHTORH (2, NGP, TSKIN, PSA, 1. _d 0, QDUMMY, dEvp,

      CALL SHTORH (0, NGP, TSFC, PSA, 1. _d 0, QDUMMY, RDUMMY,

      _RL  PSA(NGP), FMASK(NGP), EMISloc

      CALL SHTORH (0, NGP, Tsurf, PSA, 1. _d 0, QDUMMY, RDUMMY,

     I                   PSA, FMASK,

      _RL  PSA(NGP)

     I                   PSA, FMASK, EMISloc,

      CALL SHTORH (2, NGP, TSKIN, PSA, 1. _d 0, QDUMMY, dEvp,

      CALL SHTORH (0, NGP, TSFC, PSA, 1. _d 0, QDUMMY, RDUMMY,

      _RL  PSA(NGP), FMASK(NGP), EMISloc

          WVSurf(I2,myThid) = (LOG(PSA(I2))-SIGL(K))*WVI(K-1,1)

      _RL Vsurf2(NGP), PSA(NGP), dpFac(NGP,NLEV)

     O           TA, QA, ThA, Vsurf2, PSA, dpFac,

         PSA(I2) = Ro_surf(i,j,bi,bj)/atm_po

         PSA(I2) = 1.

      _RL PSA(NGP), dpFac(NGP,NLEV), QA(NGP,NLEV), QSAT(NGP,NLEV)

      SUBROUTINE LSCOND (PSA,dpFac,QA,QSAT,

        PSA2(J)    = PSA(J)*PSA(J)

           ENTR(K)= ( MAX( 0. _d 0, SIG(K)/PSA(J) - 0.5 _d 0) )**2

           ENTR_PS(J,K) = ENTR(K)*SENTR*PSA(J)

        IF (QA(J,Ktmp).LT.QATHR(J).OR.PSA(J).LT.PSMIN)

      _RL PSA(NGP), SE(NGP,NLEV), QA(NGP,NLEV), QSAT(NGP,NLEV)

      SUBROUTINE CONVMF (PSA,dpFac,SE,QA,QSAT,

        FPSA = MIN(1. _d 0 ,(PSA(J)-PSMIN)*RDPS)

      SUBROUTINE RADSW (PSA,dpFac,QA,RH,ALB,

      _RL PSA(NGP),dpFac(NGP,NLEV),QA(NGP,NLEV),RH(NGP,NLEV)

         ABS1=ABSDRY+ABSAER*(SIG(K)/PSA(J))**2

       ABS1=ABSDRY+ABSAER*(SIG(K)/PSA(J))**2

        STRATC(J)=STRATZ(J)*PSA(J)

        FLUX (J,1)=TAU2(J,1,1)*(FLUX(J,1)-OZONE(J)*PSA(J))

     I                   PSA, dpFac, ThA,

      _RL PSA(NGP), dpFac(NGP,NLEV), ThA(NGP,NLEV)

        T1(J) = ThA(J,kGrd(J))*(PSA(J)**kappa)