Changeset 5797 for LMDZ6/branches/contrails/libf/phylmd/physiq_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/physiq_mod.F90 (modified) (12 diffs)

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

@@ -39 +39 @@
     USE ioipsl_getin_p_mod, ONLY : getin_p
     USE indice_sol_mod
-    USE infotrac_phy, ONLY: nqtot, nbtr, nqo, nqtke, tracers, type_trac, addPhase, ivap, iliq, isol, ibs, icf, iqvc, itke, icfc, iqtc, inic
+    USE infotrac_phy, ONLY: nqtot, nbtr, nqo, nqtke, tracers, type_trac, addPhase, &
+        ivap, iliq, isol, ibs, icf, iqvc, itke, icfc, iqtc, iqic, inic
     USE strings_mod,  ONLY: strIdx
     USE iophy
@@ -341 +342 @@
        dqvc_adj, dqvc_sub, dqvc_con, dqvc_mix, qsatliq, qsatice, &
        !-- LSCP - aviation and contrails variables
-       cfc_seri, qtc_seri, nic_seri, d_cfc_dyn, d_qtc_dyn, d_nic_dyn, &
+       cfc_seri, qtc_seri, qic_seri, nic_seri, d_cfc_dyn, d_qtc_dyn, d_qic_dyn, d_nic_dyn, &
        d_q_avi, flight_dist, flight_fuel, &
        qice_cont, contfra, qradice_cont, &
@@ -2491 +2492 @@
           cfc_seri(i,k)= 0.
           qtc_seri(i,k)= 0.
+          qic_seri(i,k)= 0.
           nic_seri(i,k)= 0.
           !CR: ATTENTION, on rajoute la variable glace
@@ -2508 +2510 @@
              cfc_seri(i,k) = qx(i,k,icfc)
              qtc_seri(i,k) = qx(i,k,iqtc)
+             qic_seri(i,k) = qx(i,k,iqic)
              nic_seri(i,k) = qx(i,k,inic)
              !--DYNAMICO can return NaNs for children tracers
              IF (ISNAN(cfc_seri(i,k))) cfc_seri(i,k) = 0.
              IF (ISNAN(qtc_seri(i,k))) qtc_seri(i,k) = 0.
+             IF (ISNAN(qic_seri(i,k))) qic_seri(i,k) = 0.
              IF (ISNAN(nic_seri(i,k))) nic_seri(i,k) = 0.
           ENDIF
@@ -2603 +2607 @@
        d_cfc_dyn(:,:)= (cfc_seri(:,:)-cfc_ancien(:,:))/phys_tstep
        d_qtc_dyn(:,:)= (qtc_seri(:,:)-qtc_ancien(:,:))/phys_tstep
+       d_qic_dyn(:,:)= (qic_seri(:,:)-qic_ancien(:,:))/phys_tstep
        d_nic_dyn(:,:)= (nic_seri(:,:)-nic_ancien(:,:))/phys_tstep
        d_tke_dyn(:,:)= (pbl_tke(:,:,is_ave)-tke_ancien(:,:))/phys_tstep
@@ -2628 +2633 @@
        d_cfc_dyn(:,:)= 0.0
        d_qtc_dyn(:,:)= 0.0
+       d_qic_dyn(:,:)= 0.0
        d_nic_dyn(:,:)= 0.0
        d_tke_dyn(:,:)= 0.0
@@ -2755 +2761 @@
             qvc_seri(i,k) = qvc_seri(i,k) / ( q_seri(i,k) + ql_seri(i,k) + qs_seri(i,k) )
             qtc_seri(i,k) = qtc_seri(i,k) / ( q_seri(i,k) + ql_seri(i,k) + qs_seri(i,k) )
+            qic_seri(i,k) = qic_seri(i,k) / ( q_seri(i,k) + ql_seri(i,k) + qs_seri(i,k) )
           ELSE
             ql_seri_lscp(i,k) = 0.
@@ -2760 +2767 @@
             qvc_seri(i,k) = 0.
             qtc_seri(i,k) = 0.
+            qic_seri(i,k) = 0.
           ENDIF
         ENDDO
@@ -4033 +4041 @@
         DO i = 1, klon
           qtc_seri(i,k) = qtc_seri(i,k) * q_seri(i,k)
+          qic_seri(i,k) = qic_seri(i,k) * q_seri(i,k)
         ENDDO
       ENDDO
@@ -4063 +4072 @@
          dqi_adj, dqi_sub, dqi_con, dqi_mix, &
          dqvc_adj, dqvc_sub, dqvc_con, dqvc_mix, qsatliq, qsatice, &
-         cfc_seri, qtc_seri, nic_seri, qice_cont, &
+         cfc_seri, qtc_seri, qic_seri, nic_seri, qice_cont, &
          flight_dist, flight_fuel, contfra, qradice_cont, &
          Tcritcont, qcritcont, potcontfraP, potcontfraNP, &
@@ -5823 +5832 @@
              d_qx(i,k,icfc) = ( cfc_seri(i,k) - qx(i,k,icfc) ) / phys_tstep
              d_qx(i,k,iqtc) = ( qtc_seri(i,k) - qx(i,k,iqtc) ) / phys_tstep
+             d_qx(i,k,iqic) = ( qic_seri(i,k) - qx(i,k,iqic) ) / phys_tstep
              d_qx(i,k,inic) = ( nic_seri(i,k) - qx(i,k,inic) ) / phys_tstep
           ENDIF
@@ -5866 +5876 @@
     cfc_ancien(:,:)= cfc_seri(:,:)
     qtc_ancien(:,:)= qtc_seri(:,:)
+    qic_ancien(:,:)= qic_seri(:,:)
     nic_ancien(:,:)= nic_seri(:,:)
     tke_ancien(:,:)= pbl_tke(:,:,is_ave)