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

Timestamp:

Aug 5, 2025, 2:22:14 PM (10 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/lmdz_aviation.f90 (modified) (1 diff)

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

@@ -238 +238 @@
         dist_contrails = flight_dist(i) * dtime * potcontfraP(i)
         dcfc_ini(i) = dist_contrails * section_contrails(i)
-        dqtc_ini(i) = icesat_ratio * qsat(i) * dcfc_ini(i)
-        dqic_ini(i) = dqtc_ini(i) - qsat(i) * dcfc_ini(i)
+        dqtc_ini(i) = qsat(i) * dcfc_ini(i)
+        dqic_ini(i) = ( icesat_ratio - 1. ) * qsat(i) * dcfc_ini(i)
         dnic_ini(i) = Nice_per_m_init_contrails * dist_contrails / rho(i)
       ENDIF