Changeset 5224 for LMDZ6/branches/Amaury_dev/libf/phylmd/phyetat0_mod.F90

Timestamp:

Sep 24, 2024, 10:47:17 AM (4 weeks ago)

Author:

abarral

Message:

Merge r5204 r5205
Light lint
Correct missing IOIPSL includes

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/phyetat0_mod.F90 (modified) (6 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/phyetat0_mod.F90

@@ -1 @@
-! $Id$
-
 MODULE phyetat0_mod
   USE lmdz_abort_physic, ONLY: abort_physic
@@ -21 +19 @@
             du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
             falb_dir, falb_dif, prw_ancien, prlw_ancien, prsw_ancien, prbsw_ancien, &
-            ftsol, pbl_tke, pctsrf, q_ancien, ql_ancien, qs_ancien, qbs_ancien, rneb_ancien, radpas, radsol, rain_fall, ratqs, &
+            ftsol, pbl_tke, pctsrf, q_ancien, ql_ancien, qs_ancien, qbs_ancien, &
+            cf_ancien, rvc_ancien, radpas, radsol, rain_fall, ratqs, &
             rnebcon, rugoro, sig1, snow_fall, bs_fall, solaire_etat0, sollw, sollwdown, &
             solsw, solswfdiff, t_ancien, u_ancien, v_ancien, w01, wake_cstar, wake_deltaq, &
@@ -397 +396 @@
     ancien_ok = ancien_ok.AND.phyetat0_get(ql_ancien, "QLANCIEN", "QLANCIEN", 0.)
     ancien_ok = ancien_ok.AND.phyetat0_get(qs_ancien, "QSANCIEN", "QSANCIEN", 0.)
-    ancien_ok = ancien_ok.AND.phyetat0_get(rneb_ancien, "RNEBANCIEN", "RNEBANCIEN", 0.)
     ancien_ok = ancien_ok.AND.phyetat0_get(u_ancien, "UANCIEN", "UANCIEN", 0.)
     ancien_ok = ancien_ok.AND.phyetat0_get(v_ancien, "VANCIEN", "VANCIEN", 0.)
@@ -413 +411 @@
     ENDIF
 
+    ! cas specifique des variables de la sursaturation par rapport a la glace
+    IF (ok_ice_supersat) THEN
+      ancien_ok = ancien_ok.AND.phyetat0_get(cf_ancien, "CFANCIEN", "CFANCIEN", 0.)
+      ancien_ok = ancien_ok.AND.phyetat0_get(rvc_ancien, "RVCANCIEN", "RVCANCIEN", 0.)
+    ELSE
+      cf_ancien(:, :) = 0.
+      rvc_ancien(:, :) = 0.
+    ENDIF
+
     ! Ehouarn: addtional tests to check if t_ancien, q_ancien contain
     !          dummy values (as is the case when generated by ce0l,
@@ -419 +426 @@
             (maxval(ql_ancien)==minval(ql_ancien))     .OR. &
             (maxval(qs_ancien)==minval(qs_ancien))     .OR. &
-            (maxval(rneb_ancien)==minval(rneb_ancien)) .OR. &
             (maxval(prw_ancien)==minval(prw_ancien))   .OR. &
             (maxval(prlw_ancien)==minval(prlw_ancien)) .OR. &
@@ -431 +437 @@
               (maxval(prbsw_ancien)==minval(prbsw_ancien))) THEN
         ancien_ok = .FALSE.
+      ENDIF
+    ENDIF
+
+    IF (ok_ice_supersat) THEN
+      IF ((maxval(cf_ancien)==minval(cf_ancien))     .OR. &
+              (maxval(rvc_ancien)==minval(rvc_ancien))) THEN
+        ancien_ok = .false.
       ENDIF
     ENDIF