Changeset 4048 for trunk

Timestamp:

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

Author:

aslmd

Message:

Titan PCM:
Some comments
EMo

Location:

trunk

Files:

• LMDZ.TITAN/libf/phytitan/calc_condlh.F90 (modified) (1 diff)
• LMDZ.TITAN/libf/phytitan/calmufi.F90 (modified) (2 diffs)
• WRF.COMMON/INTERFACES_V4/dynphy_wrf_titan_lmd/callphysiq_mod.F (modified) (1 diff)
• WRF.COMMON/INTERFACES_V4/module_lmd_driver.F (modified) (1 diff)

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

@@ -19 +19 @@
 !   OUTPUT
 !   ------
-!   Lcx = Condensation latente heat [J.kg-1]
+!   Lcx = Condensation latent heat [J.kg-1]
 !
 !   Author

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

@@ -35 +35 @@
   REAL(kind=8), DIMENSION(:,:,:), INTENT(IN)  :: pq    !! Tracers (\(X.kg^{-1}}\)).
   REAL(kind=8), DIMENSION(:,:,:), INTENT(IN)  :: zdqfi !! Tendency from former processes for tracers (\(X.kg^{-1}}\)).
-  REAL(kind=8), DIMENSION(:,:,:), INTENT(OUT) :: zdq   !! Microphysical tendency for tracers (\(X.kg^{-1}}\)).
+  REAL(kind=8), DIMENSION(:,:,:), INTENT(OUT) :: zdq   !! Microphysical tendency for tracers (X.kg_air-1.s-1)
   
   REAL(kind=8), DIMENSION(:,:,:), ALLOCATABLE :: zq !! Local tracers updated from former processes (\(X.kg^{-1}}\)).
@@ -58 +58 @@
   REAL(kind=8), DIMENSION(:,:), ALLOCATABLE :: dgazs !! Tendencies of each condensible gaz species !(\(mol.mol^{-1}\)). 
 
-  REAL(kind=8), DIMENSION(:,:), ALLOCATABLE ::  int2ext !! (\(m^{-2}\)).
+  REAL(kind=8), DIMENSION(:,:), ALLOCATABLE ::  int2ext !! (kg.m-2)
   REAL(kind=8), DIMENSION(:), ALLOCATABLE   :: tmp 
   TYPE(error) :: err

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

@@ -40 +40 @@
   ! these two quantities do not even have the same units (kg/kg and kg/m^2)!!!
   ! seems it is the same for some other physics. It is probably a bug.
-  !EMo25 zqfi_omp in kg/kghumidair?
+  !EMo25 zqfi_omp in kg/kghumidair for gases, m3/kghumidair for ice tracers or in (nb of particles)/kghumidair for m0 tracers
 
   call allocate_comm_wrf(klon,llm,nqtot)

trunk/WRF.COMMON/INTERFACES_V4/module_lmd_driver.F

@@ -752 +752 @@
       SCALAR(i,kps:kpe,j,P_MARKER) = SCALAR(i,kps:kpe,j,P_MARKER)*exp(-dt/tau_decay)
       SCALAR(i,1,j,P_MARKER) = 1. !! this tracer is emitted in the surface layer
-    CASE(61) !emoisan tobechecked
+    CASE(61) !emoisan
       scalar(i,kps:kpe,j,P_mu_m0as)=scalar(i,kps:kpe,j,P_mu_m0as) &
            +zdqfi_omp(subs,kps:kpe,1)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))