/[MITgcm]/MITgcm/pkg/aim_v23/phy_lscond.F
ViewVC logotype

Diff of /MITgcm/pkg/aim_v23/phy_lscond.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jmc, Fri Nov 22 17:17:03 2002 UTC revision 1.2 by jmc, Thu Jan 26 01:14:11 2006 UTC
# Line 76  c_FM  QSMAX = 50. Line 76  c_FM  QSMAX = 50.
76  C--   2. Tendencies of temperature and moisture  C--   2. Tendencies of temperature and moisture
77  C        NB. A maximum heating rate is imposed to avoid  C        NB. A maximum heating rate is imposed to avoid
78  C            grid-point-storm instability  C            grid-point-storm instability
 C-jmc: this breaks heat conservation !!  
 C   for checking the heat budget, set a very large QSMAX (e.g. 1.e+9)  
 C   or use the simplyfied form below.  
79    
80        DO K=2,NLEV        DO K=2,NLEV
81          SIG2=SIG(K)*SIG(K)          SIG2=SIG(K)*SIG(K)
# Line 89  c_FM    DQMAX = (1.1-RHREF)*QSMAX*SIG2*R Line 86  c_FM    DQMAX = (1.1-RHREF)*QSMAX*SIG2*R
86            DQMAX = (1.1 _d 0-RHREF)*QSMAX*SIG2*RTLSC            DQMAX = (1.1 _d 0-RHREF)*QSMAX*SIG2*RTLSC
87            DQLSC(J,K) = MIN(0. _d 0,(RHREF*QSAT(J,K)-QA(J,K)))*RTLSC            DQLSC(J,K) = MIN(0. _d 0,(RHREF*QSAT(J,K)-QA(J,K)))*RTLSC
88  c_FM      DTLSC(J,K) = TFACT*MIN(-DQLSC(J,K),DQMAX*PSA2(J))  c_FM      DTLSC(J,K) = TFACT*MIN(-DQLSC(J,K),DQMAX*PSA2(J))
89            DTLSC(J,K) = TFACT*MIN(-DQLSC(J,K),DQMAX)  c         DTLSC(J,K) = TFACT*MIN(-DQLSC(J,K),DQMAX)
90  C-jmc: Heat_Conserve:    C-jmc: the threshold on latent heat (above) breaks the energy conservation;
91  c         DTLSC(J,K) = -TFACT*DQLSC(J,K)  C      to fix this, apply the threshold directly to the condensation
92              DQLSC(J,K) = MAX(-DQMAX, DQLSC(J,K) )
93              DTLSC(J,K) = -TFACT*DQLSC(J,K)
94          ENDDO          ENDDO
95        ENDDO        ENDDO
96    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22