Changeset 4045 for trunk

Timestamp:

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

Author:

aslmd

Message:

Titan CRM

• add new output variable rsat_ch4 (methane relative humidity)

EMo

Location:

trunk

Files:

• LMDZ.TITAN/libf/phytitan/comm_wrf.F90 (modified) (3 diffs)
• LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (modified) (3 diffs)
• WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/update_outputs_physiq_mod.F (modified) (3 diffs)
• WRF.COMMON/INTERFACES_V4/module_lmd_driver.F (modified) (3 diffs)

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

@@ -20 +20 @@
   REAL,SAVE,ALLOCATABLE :: comm_zqfi_omp(:,:,:)
   REAL,SAVE,ALLOCATABLE :: comm_zdtlc(:,:)
+  REAL,SAVE,ALLOCATABLE :: comm_rsat_ch4(:,:)
 
 contains
@@ -43 +44 @@
   allocate(comm_zqfi_omp(ngrid,nlayer,nq))
   allocate(comm_zdtlc(ngrid,nlayer))
+  allocate(comm_rsat_ch4(ngrid,nlayer))
 
   end subroutine allocate_comm_wrf
@@ -63 +65 @@
   deallocate(comm_zqfi_omp)
   deallocate(comm_zdtlc)
+  deallocate(comm_rsat_ch4)
 
   end subroutine deallocate_comm_wrf

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      
 

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/update_outputs_physiq_mod.F

@@ -165 +165 @@
             DTLSC,DTRAIN,DT_MOIST,&
             H2OICE_REFF,LATENT_HF,&
-            DT_COND)
+            DT_COND,RSAT_CH4)
 
    USE comm_wrf !! to get fields to be written from physiq
@@ -176 +176 @@
    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: &
      HR_SW,HR_LW,CLOUDFRAC,HR_DYN,DT_RAD,RH,DQICE,DQVAP,&
-     DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,DT_COND
+     DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,DT_COND,RSAT_CH4
    INTEGER :: i,j,subs
 
@@ -205 +205 @@
      FLXGRD(i,j) = comm_FLXGRD(subs)
      DT_COND(i,kps:kpe,j) = comm_ZDTLC(subs,kps:kpe)
+     RSAT_CH4(i,kps:kpe,j) = comm_rsat_ch4(subs,kps:kpe)
 !    DTLSC(i,kps:kpe,j) = comm_DTLSC(subs,kps:kpe)
 !    DTRAIN(i,kps:kpe,j) = comm_DTRAIN(subs,kps:kpe)

trunk/WRF.COMMON/INTERFACES_V4/module_lmd_driver.F

@@ -49 +49 @@
         USTM,HFX,&
         SLPX,SLPY,RESTART,&
-        DT_COND)
+        DT_COND,RSAT_CH4)
 ! NB: module_lmd_driver_output1.inc : output arguments generated from Registry
 
@@ -116 +116 @@
      HR_SW,HR_LW,HR_DYN,DT_RAD,&
      CLOUDFRAC,RH,DQICE,DQVAP,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,&
-     DT_COND
+     DT_COND, RSAT_CH4
 REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: &
      P_Q2
@@ -664 +664 @@
             ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,&
             FLUXSURF_LW,FLXGRD,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,LATENT_HF,&
-            DT_COND)
+            DT_COND,RSAT_CH4)
 !!!
 !print *, '** ',planet_type,'** OUTPUT PHYSICS DONE'