Changeset 4045 for trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

Timestamp:

Feb 5, 2026, 5:50:32 PM (4 weeks ago)

Author:

aslmd

Message:

Titan CRM

• add new output variable rsat_ch4 (methane relative humidity)

EMo

File:

• trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (modified) (3 diffs)

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

@@ -49 +49 @@
                      comm_FLUXTOP_LW, comm_FLUXSURF_SW,&
                      comm_FLUXSURF_LW, comm_FLXGRD,&
-                     comm_zqfi_omp, comm_zdtlc
+                     comm_zqfi_omp, comm_zdtlc,&
+                     comm_rsat_ch4
 #endif
 #ifdef CPP_XIOS      
@@ -2039 +2040 @@
 #endif
 #else 
-      !MESOSCALE outputs      
+      !! =================
+      !! MESOSCALE outputs
+      !! =================
       comm_HR_SW(1:ngrid,1:nlayer) = zdtsw(1:ngrid,1:nlayer)
       comm_HR_LW(1:ngrid,1:nlayer) = zdtlw(1:ngrid,1:nlayer)
@@ -2051 +2054 @@
       comm_zqfi_omp(1:ngrid,1:nlayer,1:nq) = zq(1:ngrid,1:nlayer,1:nq)
       comm_zdtlc(1:ngrid,1:nlayer) = zdtlc(1:ngrid,1:nlayer)
+
+      !CH4 saturation ratio [kg/kg]
+      do iq=1, nice
+         if(trim(nameOfTracer(gazs_indx(iq))) .eq. "CH4") then
+            comm_rsat_ch4(1:ngrid,1:nlayer) = mmd_gazs_sat(1:ngrid,1:nlayer,iq)
+         endif
+      enddo
 #endif