Changeset 3963 for trunk/LMDZ.TITAN

Timestamp:

Nov 18, 2025, 1:41:14 PM (4 days ago)

Author:

debatzbr

Message:

Titan PCM: Corrects the calculation of latent heat + add HC3N calculation.
BBT

File:

• trunk/LMDZ.TITAN/libf/phytitan/calc_condlh.F90 (modified) (7 diffs)

trunk/LMDZ.TITAN/libf/phytitan/calc_condlh.F90

@@ -24 +24 @@
 !   ------
 !   B. de Batz de Trenquelléon (07/2023)
-!
-!   !!! WARNING !!!
-!   ---------------
-!   We suppose a given order of tracers !
-!   Tracers come from microphysics (ices_indx) :
-!   CH4  --> 1
-!   C2H2 --> 2
-!   C2H6 --> 3
-!   HCN  --> 4
-!   AC6H6--> 5
 !====================================================================    
 
@@ -56 +46 @@
 real, parameter :: mmolC2H2 = 26.e-3 ! Molar mass of C2H2 [kg.mol-1]
 real, parameter :: mmolC2H6 = 30.e-3 ! Molar mass of C2H6 [kg.mol-1]
+real, parameter :: mmolC6H6 = 78.e-3 ! Molar mass of AC6H6 [kg.mol-1]
 real, parameter :: mmolHCN  = 27.e-3 ! Molar mass of HCN [kg.mol-1]
-real, parameter :: mmolAC6H6 = 78.e-3 ! Molar mass of AC6H6 [kg.mol-1]
+real, parameter :: mmolHC3N = 51.e-3 ! Molar mass of HC3N [kg.mol-1]
 
 real, parameter :: TcCH4  = 190.56 ! Critical point of CH4  [K]
 real, parameter :: TcC2H2 = 308.35 ! Critical point of C2H2 [K]
 real, parameter :: TcC2H6 = 305.35 ! Critical point of C2H6 [K]
+real, parameter :: TcC6H6 = 562.10 ! Critical point of AC6H6 [K]
 real, parameter :: TcHCN  = 456.70 ! Critical point of HCN [K]
-real, parameter :: TcAC6H6 = 562.10 ! Critical point of AC6H6 [K]
+real, parameter :: TcHC3N = 527.00 ! Critical point of HC3N [K]
 
 ! Local variables :
@@ -83 +75 @@
         do ig = 1, ngrid ! Main loop on the horizontal grid
             do l = 1, nlayer ! Loop on vertical layers
-                ftx = (1. - pt(ig,l)) / TcCH4
+                ftx = (1. - pt(ig,l) / TcCH4)
                 if (ftx .le. 1.e-3) ftx = 1.e-3
                 Lcx(ig,l,iq) = 8.314 * TcCH4 * (7.08 * ftx**0.354 + 10.95 * 1.1e-2 * ftx**0.456) / mmolCH4
@@ -96 +88 @@
         do ig = 1, ngrid ! Main loop on the horizontal grid
             do l = 1, nlayer ! Loop on vertical layers       
-                ftx = (1. - pt(ig,l)) / TcC2H2
+                ftx = (1. - pt(ig,l) / TcC2H2)
                 if (ftx .le. 1.e-3) ftx = 1.e-3
                 Lcx(ig,l,iq) = 8.314 * TcC2H2 * (7.08 * ftx**0.354 + 10.95 * 1.9e-1 * ftx**0.456) / mmolC2H2
@@ -109 +101 @@
         do ig = 1, ngrid ! Main loop on the horizontal grid
             do l = 1, nlayer ! Loop on vertical layers
-                ftx = (1. - pt(ig,l)) / TcC2H6
+                ftx = (1. - pt(ig,l) / TcC2H6)
                 if (ftx .le. 1.e-3) ftx = 1.e-3
                 Lcx(ig,l,iq) = 8.314 * TcC2H6 * (7.08 * ftx**0.354 + 10.95 * 9.9e-2 * ftx**0.456) / mmolC2H6
@@ -115 +107 @@
         enddo ! End of main loop on the horizontal grid
     
-    ! AC6H6 :
-    !--------
+    ! C6H6 :
+    !-------
     else if(trim(nameOfTracer(gazs_indx(iq))) .eq. "AC6H6") then
         
         do ig = 1, ngrid ! Main loop on the horizontal grid
             do l = 1, nlayer ! Loop on vertical layers
-                ftx = (1. - pt(ig,l)) / TcAC6H6
+                ftx = (1. - pt(ig,l) / TcC6H6)
                 if (ftx .le. 1.e-3) ftx = 1.e-3
-                Lcx(ig,l,iq) = 8.314 * TcAC6H6 * (7.08 * ftx**0.354 + 10.95 * 9.9e-2 * ftx**0.456) / mmolAC6H6
+                Lcx(ig,l,iq) = 8.314 * TcC6H6 * (7.08 * ftx**0.354 + 10.95 * 9.9e-2 * ftx**0.456) / mmolC6H6
             enddo ! End of loop on vertical layers
         enddo ! End of main loop on the horizontal grid
@@ -133 +125 @@
         do ig = 1, ngrid ! Main loop on the horizontal grid
             do l = 1, nlayer ! Loop on vertical layers
-                ftx = (1. - pt(ig,l)) / TcHCN
+                ftx = (1. - pt(ig,l) / TcHCN)
                 if (ftx .le. 1.e-3) ftx = 1.e-3
                 Lcx(ig,l,iq) = 8.314 * TcHCN * (7.08 * ftx**0.354 + 10.95 * 3.88e-1 * ftx**0.456) / mmolHCN
+            enddo ! End of loop on vertical layers
+        enddo ! End of main loop on the horizontal grid
+
+    ! HC3N :
+    !-------
+    else if(trim(nameOfTracer(gazs_indx(iq))) .eq. "HC3N") then
+
+        do ig = 1, ngrid ! Main loop on the horizontal grid
+            do l = 1, nlayer ! Loop on vertical layers
+                ftx = (1. - pt(ig,l) / TcHC3N)
+                if (ftx .le. 1.e-3) ftx = 1.e-3
+                Lcx(ig,l,iq) = 8.314 * TcHC3N * (7.08 * ftx**0.354 + 10.95 * 3.88e-1 * ftx**0.456) / mmolHC3N
             enddo ! End of loop on vertical layers
         enddo ! End of main loop on the horizontal grid