Changeset 2705

Timestamp:

Nov 20, 2016, 10:54:34 PM (8 years ago)

Author:

fhourdin

Message:

Extraction du calcul de la reevaporation dans un sous programme.
Ajout des tendances par add_phys_tend.
Travail preparatoire a la conservation de l'energie.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• physiq_mod.F90 (modified) (2 diffs)
• reevap.F90 (added)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -54 +54 @@
        d_t_ajs_x,d_q_ajs_x, &
        !
-       d_t_eva,d_q_eva, &
+       d_t_eva,d_q_eva,d_ql_eva,d_qi_eva, &
        d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc, &
        d_t_lscst,d_q_lscst, &
@@ -1958 +1958 @@
     ! Re-evaporer l'eau liquide nuageuse
     !
-    DO k = 1, klev  ! re-evaporation de l'eau liquide nuageuse
-       DO i = 1, klon
-          zlvdcp=RLVTT/RCPD/(1.0+RVTMP2*q_seri(i,k))
-          !jyg<
-          !  Attention : Arnaud a propose des formules completement differentes
-          !                  A verifier !!!
-          zlsdcp=RLSTT/RCPD/(1.0+RVTMP2*q_seri(i,k))
-          IF (iflag_ice_thermo .EQ. 0) THEN
-             zlsdcp=zlvdcp
-          ENDIF
-          !>jyg
-
-          IF (iflag_ice_thermo.eq.0) THEN   
-             !pas necessaire a priori
-
-             zdelta = MAX(0.,SIGN(1.,RTT-t_seri(i,k)))
-             zb = MAX(0.0,ql_seri(i,k))
-             za = - MAX(0.0,ql_seri(i,k)) &
-                  * (zlvdcp*(1.-zdelta)+zlsdcp*zdelta)
-             t_seri(i,k) = t_seri(i,k) + za
-             q_seri(i,k) = q_seri(i,k) + zb
-             ql_seri(i,k) = 0.0
-             d_t_eva(i,k) = za
-             d_q_eva(i,k) = zb
-
-          ELSE
-
-             !CR: on r\'e-\'evapore eau liquide et glace
-
-             !        zdelta = MAX(0.,SIGN(1.,RTT-t_seri(i,k)))
-             !        zb = MAX(0.0,ql_seri(i,k))
-             !        za = - MAX(0.0,ql_seri(i,k)) &
-             !             * (zlvdcp*(1.-zdelta)+zlsdcp*zdelta)
-             zb = MAX(0.0,ql_seri(i,k)+qs_seri(i,k))
-             za = - MAX(0.0,ql_seri(i,k))*zlvdcp & 
-                  - MAX(0.0,qs_seri(i,k))*zlsdcp
-             t_seri(i,k) = t_seri(i,k) + za
-             q_seri(i,k) = q_seri(i,k) + zb
-             ql_seri(i,k) = 0.0
-             !on \'evapore la glace
-             qs_seri(i,k) = 0.0
-             d_t_eva(i,k) = za
-             d_q_eva(i,k) = zb
-          ENDIF
-
-       ENDDO
-    ENDDO
-    !IM
+     CALL reevap (klon,klev,iflag_ice_thermo,t_seri,q_seri,ql_seri,qs_seri, &
+   &         d_t_eva,d_q_eva,d_ql_eva,d_qi_eva)
+
+     CALL add_phys_tend &
+            (du0,dv0,d_t_eva,d_q_eva,d_ql_eva,d_qi_eva,paprs,&
+               'eva',abortphy,flag_inhib_tend)
+
     IF (ip_ebil_phy.ge.2) THEN 
        ztit='after reevap'