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

Timestamp:

May 1, 2025, 6:00:03 PM (2 months ago)

Author:

aborella

Message:

Changed the contrails classes that are advected

File:

• LMDZ6/branches/contrails/libf/phylmd/phyetat0_mod.f90 (modified) (3 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, qvc_ancien, qvcon, qccon, cfa_ancien, pcf_ancien, qva_ancien, qia_ancien, &
+       cf_ancien, qvc_ancien, qvcon, qccon, cfl_ancien, cfc_ancien, qtl_ancien, qtc_ancien, &
        radpas, radsol, rain_fall, &
        ratqs, rnebcon, rugoro, sig1, snow_fall, bs_fall, solaire_etat0, sollw, sollwdown, &
@@ -425 +425 @@
   ! cas specifique des variables de l'aviation
   IF ( ok_plane_contrail ) THEN
-    ancien_ok=ancien_ok.AND.phyetat0_get(cfa_ancien,"CFAANCIEN","CFAANCIEN",0.)
-    ancien_ok=ancien_ok.AND.phyetat0_get(pcf_ancien,"PCFANCIEN","PCFANCIEN",0.)
-    ancien_ok=ancien_ok.AND.phyetat0_get(qva_ancien,"QVAANCIEN","QVAANCIEN",0.)
-    ancien_ok=ancien_ok.AND.phyetat0_get(qia_ancien,"QIAANCIEN","QIAANCIEN",0.)
-  ELSE
-    cfa_ancien(:,:)=0.
-    pcf_ancien(:,:)=0.
-    qva_ancien(:,:)=0.
-    qia_ancien(:,:)=0.
+    ancien_ok=ancien_ok.AND.phyetat0_get(cfl_ancien,"CFLANCIEN","CFLANCIEN",0.)
+    ancien_ok=ancien_ok.AND.phyetat0_get(cfc_ancien,"CFCANCIEN","CFCANCIEN",0.)
+    ancien_ok=ancien_ok.AND.phyetat0_get(qtl_ancien,"QTLANCIEN","QTLANCIEN",0.)
+    ancien_ok=ancien_ok.AND.phyetat0_get(qtc_ancien,"QTCANCIEN","QTCANCIEN",0.)
+  ELSE
+    cfl_ancien(:,:)=0.
+    cfc_ancien(:,:)=0.
+    qtl_ancien(:,:)=0.
+    qtc_ancien(:,:)=0.
   ENDIF
 
@@ -464 +464 @@
 
   IF ( ok_plane_contrail ) THEN
-    IF ( ( maxval(cfa_ancien).EQ.minval(cfa_ancien) ) .OR. &
-         ( maxval(pcf_ancien).EQ.minval(pcf_ancien) ) .OR. &
-         ( maxval(qva_ancien).EQ.minval(qva_ancien) ) .OR. &
-         ( maxval(qia_ancien).EQ.minval(qia_ancien) ) ) THEN
+    IF ( ( maxval(cfl_ancien).EQ.minval(cfl_ancien) ) .OR. &
+         ( maxval(cfc_ancien).EQ.minval(cfc_ancien) ) .OR. &
+         ( maxval(qtl_ancien).EQ.minval(qtl_ancien) ) .OR. &
+         ( maxval(qtc_ancien).EQ.minval(qtc_ancien) ) ) THEN
        ancien_ok=.false.
      ENDIF