Changeset 2613 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Aug 23, 2016, 5:21:06 PM (8 years ago)

Author:

fhourdin

Message:

Intronduction d'une limite sur l'eau liquide nuageuse, oliqmax
Typiquement oliqmax>1e-3 (1g/kg)
Peut permettre d'éviter certains plantages.

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -578 +578 @@
     LOGICAL,SAVE :: ok_adjwk=.FALSE.
     !$OMP THREADPRIVATE(ok_adjwk) 
+    REAL,SAVE :: oliqmax=999.
+    !$OMP THREADPRIVATE(oliqmax) 
     REAL, SAVE :: alp_offset
     !$OMP THREADPRIVATE(alp_offset)
@@ -1197 +1199 @@
        CALL getin_p('random_notrig_max',random_notrig_max)
        CALL getin_p('ok_adjwk',ok_adjwk)
+       CALL getin_p('oliqmax',oliqmax)
        CALL getin_p('ratqsp0',ratqsp0)
        CALL getin_p('ratqsdp',ratqsdp)
@@ -3056 +3059 @@
     CALL add_phys_tend(du0,dv0,d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc,paprs, &
          'lsc',abortphy,flag_inhib_tend)
+    rain_num(:)=0.
+    DO k = 1, klev  ! re-evaporation de l'eau liquide nuageuse
+       DO i = 1, klon
+          IF (ql_seri(i,k)>oliqmax) THEN
+             rain_num(i)=rain_num(i)+(ql_seri(i,k)-oliqmax)*zmasse(i,k)/pdtphys
+             ql_seri(i,k)=oliqmax
+          ENDIF
+       ENDDO
+    ENDDO
+
     !---------------------------------------------------------------------------
     DO k = 1, klev