Changeset 2657

Timestamp:

Oct 10, 2016, 11:41:47 AM (8 years ago)

Author:

fhourdin

Message:

Seuil max sur eau glace nuageuse oicemax (en plus de oliqmax)
Truc pour éviter des plantages dans des atmosphère super saturées.

File:

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

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -584 +584 @@
     LOGICAL,SAVE :: ok_adjwk=.FALSE.
     !$OMP THREADPRIVATE(ok_adjwk) 
-    REAL,SAVE :: oliqmax=999.
-    !$OMP THREADPRIVATE(oliqmax) 
+    REAL,SAVE :: oliqmax=999.,oicemax=999.
+    !$OMP THREADPRIVATE(oliqmax,oicemax) 
     REAL, SAVE :: alp_offset
     !$OMP THREADPRIVATE(alp_offset)
@@ -1201 +1201 @@
        CALL getin_p('ok_adjwk',ok_adjwk)
        CALL getin_p('oliqmax',oliqmax)
+       CALL getin_p('oicemax',oicemax)
        CALL getin_p('ratqsp0',ratqsp0)
        CALL getin_p('ratqsdp',ratqsdp)
@@ -3040 +3041 @@
          'lsc',abortphy,flag_inhib_tend)
     rain_num(:)=0.
-    DO k = 1, klev  ! re-evaporation de l'eau liquide nuageuse
+    DO k = 1, klev
        DO i = 1, klon
           IF (ql_seri(i,k)>oliqmax) THEN
@@ -3048 +3049 @@
        ENDDO
     ENDDO
+    IF (nqo==3) THEN
+    DO k = 1, klev
+       DO i = 1, klon
+          IF (qs_seri(i,k)>oicemax) THEN
+             rain_num(i)=rain_num(i)+(qs_seri(i,k)-oicemax)*zmasse(i,k)/pdtphys
+             qs_seri(i,k)=oicemax
+          ENDIF
+       ENDDO
+    ENDDO
+    ENDIF
 
     !---------------------------------------------------------------------------