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

Timestamp:

Aug 5, 2025, 2:22:14 PM (9 days ago)

Author:

aborella

Message:

Bugfix for saturation adjustment in cirrus mixing + bugfix for contrails sedimentation + new diagnostics + support for unadjusted contrails

File:

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

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

@@ -26 +26 @@
        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, cfc_ancien, qtc_ancien, nic_ancien, &
+       cf_ancien, qvc_ancien, qvcon, qccon, cfc_ancien, qtc_ancien, qic_ancien, nic_ancien, &
        tke_ancien, radpas, radsol, rain_fall, ratqs, &
        rnebcon, rugoro, sig1, snow_fall, bs_fall, solaire_etat0, sollw, sollwdown, &
@@ -498 +498 @@
     ancien_ok=ancien_ok.AND.phyetat0_get(cfc_ancien,"CFCANCIEN","CFCANCIEN",0.)
     ancien_ok=ancien_ok.AND.phyetat0_get(qtc_ancien,"QTCANCIEN","QTCANCIEN",0.)
+    ancien_ok=ancien_ok.AND.phyetat0_get(qic_ancien,"QICANCIEN","QICANCIEN",0.)
     ancien_ok=ancien_ok.AND.phyetat0_get(nic_ancien,"NICANCIEN","NICANCIEN",0.)
   ELSE
     cfc_ancien(:,:)=0.
     qtc_ancien(:,:)=0.
+    qic_ancien(:,:)=0.
     nic_ancien(:,:)=0.
   ENDIF
@@ -535 +537 @@
     IF ( ( maxval(cfc_ancien).EQ.minval(cfc_ancien) ) .OR. &
          ( maxval(qtc_ancien).EQ.minval(qtc_ancien) ) .OR. &
+         ( maxval(qic_ancien).EQ.minval(qic_ancien) ) .OR. &
          ( maxval(nic_ancien).EQ.minval(nic_ancien) ) ) THEN
        ancien_ok=.false.