Changeset 5797 for LMDZ6/branches/contrails/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Aug 5, 2025, 2:22:14 PM (4 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/phys_state_var_mod.F90 (modified) (3 diffs)

LMDZ6/branches/contrails/libf/phylmd/phys_state_var_mod.F90

@@ -142 +142 @@
       REAL, ALLOCATABLE, SAVE :: qvcon(:,:), qccon(:,:)
 !$OMP THREADPRIVATE(qvcon, qccon)
-      REAL, ALLOCATABLE, SAVE :: cfc_ancien(:,:), qtc_ancien(:,:), nic_ancien(:,:)
-!$OMP THREADPRIVATE(cfc_ancien, qtc_ancien, nic_ancien)
+      REAL, ALLOCATABLE, SAVE :: cfc_ancien(:,:), qtc_ancien(:,:)
+!$OMP THREADPRIVATE(cfc_ancien, qtc_ancien)
+      REAL, ALLOCATABLE, SAVE :: qic_ancien(:,:), nic_ancien(:,:)
+!$OMP THREADPRIVATE(qic_ancien, nic_ancien)
        REAL, ALLOCATABLE, SAVE :: tke_ancien(:,:)
 !$OMP THREADPRIVATE(tke_ancien)
@@ -680 +682 @@
       ALLOCATE(u_ancien(klon,klev), v_ancien(klon,klev))
       ALLOCATE(cf_ancien(klon,klev), qvc_ancien(klon,klev))
-      ALLOCATE(cfc_ancien(klon,klev), qtc_ancien(klon,klev), nic_ancien(klon,klev))
+      ALLOCATE(cfc_ancien(klon,klev), qtc_ancien(klon,klev))
+      ALLOCATE(qic_ancien(klon,klev), nic_ancien(klon,klev))
       ALLOCATE(qvcon(klon,klev), qccon(klon,klev))
       ALLOCATE(tke_ancien(klon,klev+1))
@@ -951 +954 @@
       DEALLOCATE(u_ancien, v_ancien)
       DEALLOCATE(cf_ancien, qvc_ancien)
-      DEALLOCATE(cfc_ancien, qtc_ancien, nic_ancien)
+      DEALLOCATE(cfc_ancien, qtc_ancien, qic_ancien, nic_ancien)
       DEALLOCATE(qvcon, qccon)
       DEALLOCATE(tke_ancien)