Changeset 2223 for LMDZ5/trunk/libf/phylmd

Timestamp:

Mar 10, 2015, 8:36:12 PM (10 years ago)

Author:

crio

Message:

Variation de qsat avec la température lors de la condensation grande-échelle:
Recalcul de la température après condensation au lieu d'utiliser la température Tbef calculée dans la boucle pour assurer la conservation de l'énergie malgré l'imprécision du calcul issu de la boucle de convergence.

File:

• LMDZ5/trunk/libf/phylmd/fisrtilp.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/fisrtilp.F90

@@ -783 +783 @@
         else if (iflag_fisrtilp_qsat.gt.0) then
            DO i= 1, klon
-             if (lognormale(i)) then
-             zt(i)=Tbef(i)
-             else
              zt(i) = zt(i) + zcond(i) * RLVTT/RCPD/(1.0+RVTMP2*(zq(i)+zcond(i)))
-             endif
            ENDDO
         endif
@@ -800 +796 @@
 !                                     t_glace_max, exposant_glace)
               if (iflag_t_glace.eq.0) then
-                    zfice(i) = 1.0 - (Tbef(i)-t_glace_min_old) / (RTT-t_glace_min_old)
+                    zfice(i) = 1.0 - (zt(i)-t_glace_min_old) / (RTT-t_glace_min_old)
                     zfice(i) = MIN(MAX(zfice(i),0.0),1.0)
                     zfice(i) = zfice(i)**exposant_glace_old
@@ -809 +805 @@
          else
            DO i=1, klon
-              if (lognormale(i)) then
-                 zt(i)=Tbef(i)
-              else
 ! JBM: icefrac_lsc is now a function contained in icefrac_lsc_mod
 !              zfice(i) = icefrac_lsc(zt(i), t_glace_min, &
 !                                     t_glace_max, exposant_glace)
               if (iflag_t_glace.eq.0) then
-                    zfice(i) = 1.0 - (Tbef(i)-t_glace_min_old) / (RTT-t_glace_min_old)
+                    zfice(i) = 1.0 - (zt(i)-t_glace_min_old) / (RTT-t_glace_min_old)
                     zfice(i) = MIN(MAX(zfice(i),0.0),1.0)
                     zfice(i) = zfice(i)**exposant_glace_old
@@ -822 +815 @@
               zt(i) = zt(i) + (1.-zfice(i))*zcond(i) * RLVTT/RCPD/(1.0+RVTMP2*(zq(i)+zcond(i))) &
                        +zfice(i)*zcond(i) * RLSTT/RCPD/(1.0+RVTMP2*(zq(i)+zcond(i)))
-              endif
            ENDDO
          endif