Changeset 5601 for LMDZ6/branches/contrails/libf/phylmd/phyetat0_mod.f90

Timestamp:

Apr 2, 2025, 4:04:40 PM (3 months ago)

Author:

aborella

Message:

Multiple changes:

• tracers which were ratios are now absolute quantities. This is needed because when the ratio

is not defined, some aberrations may occur

• added a new tracer for total specific humidity in contrails
• rework of the mixing process for cirrus clouds (and contrails)
• changed the numerical integration of ice crystals' sublimation
• subroutines do not take real inputs anymore (at least klon tables)
• added more radiative diagnostics for contrails

File:

• LMDZ6/branches/contrails/libf/phylmd/phyetat0_mod.f90 (modified) (4 diffs)

LMDZ6/branches/contrails/libf/phylmd/phyetat0_mod.f90 ¶

@@ -23 +23 @@
        falb_dir, falb_dif, prw_ancien, prlw_ancien, prsw_ancien, prbsw_ancien, &
        ftsol, pbl_tke, pctsrf, q_ancien, ql_ancien, qs_ancien, qbs_ancien, &
-       cf_ancien, rvc_ancien, rcont_ancien, radpas, radsol, rain_fall, ratqs, &
+       cf_ancien, qvc_ancien, cfa_ancien, qta_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, &
@@ -414 +414 @@
   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.)
+    ancien_ok=ancien_ok.AND.phyetat0_get(qvc_ancien,"QVCANCIEN","QVCANCIEN",0.)
   ELSE
     cf_ancien(:,:)=0.
-    rvc_ancien(:,:)=0.
+    qvc_ancien(:,:)=0.
   ENDIF
 
   ! cas specifique des variables de l'aviation
   IF ( ok_plane_contrail ) THEN
-    ancien_ok=ancien_ok.AND.phyetat0_get(rcont_ancien,"RCONTANCIEN","RCONTANCIEN",0.)
-  ELSE
-    rcont_ancien(:,:)=0.
+    ancien_ok=ancien_ok.AND.phyetat0_get(cfa_ancien,"CFAANCIEN","CFAANCIEN",0.)
+    ancien_ok=ancien_ok.AND.phyetat0_get(cfa_ancien,"QTAANCIEN","QTAANCIEN",0.)
+  ELSE
+    cfa_ancien(:,:)=0.
+    qta_ancien(:,:)=0.
   ENDIF
 
@@ -449 +451 @@
   IF ( ok_ice_supersat ) THEN
     IF ( (maxval(cf_ancien).EQ.minval(cf_ancien))     .OR. &
-         (maxval(rvc_ancien).EQ.minval(rvc_ancien)) ) THEN
+         (maxval(qvc_ancien).EQ.minval(qvc_ancien)) ) THEN
        ancien_ok=.false.
      ENDIF
@@ -455 +457 @@
 
   IF ( ok_plane_contrail ) THEN
-    IF ( maxval(rcont_ancien).EQ.minval(rcont_ancien) ) THEN
+    IF ( ( maxval(cfa_ancien).EQ.minval(cfa_ancien) ) .OR. &
+         ( maxval(qta_ancien).EQ.minval(qta_ancien) ) ) THEN
        ancien_ok=.false.
      ENDIF