Changeset 4033 for trunk/LMDZ.GENERIC/libf/phystd/conduction.F90

Timestamp:

Jan 28, 2026, 3:06:52 PM (40 hours ago)

Author:

mmaurice

Message:

Generic PCM

Introduce molecular viscosity module. Add callconduc and callmolvis
flags to call conduction and molecular viscosity, respectively. Change
phitop in conduction to phitop_conduc (there is also phitop_molvis now).
Move thermal conductivity lambda to conc_mod.F90.

MM

File:

• trunk/LMDZ.GENERIC/libf/phystd/conduction.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/conduction.F90

@@ -9 +9 @@
     
       use comcstfi_mod, only: r, cpp, mugaz
-      use callkeys_mod, only: phitop,zztop,a_coeff,s_coeff,force_conduction
+      use callkeys_mod, only: phitop_conduc,zztop,a_coeff,s_coeff,force_conduction
+      use conc_mod,     only: lambda
       use gases_h
 
@@ -48 +49 @@
       INTEGER :: i,ig,l,igas,kgas
       REAL :: alpha(ngrid,nlayer)
-      REAL :: lambda(ngrid,nlayer)
       REAL :: muvol(ngrid,nlayer)   ! kg.m-3
       REAL :: C(ngrid,nlayer)
@@ -111 +111 @@
               molar_frac(:,:,igas) = gfrac(igas)
               here(igas) = .true.
+            elseif(igas.eq.igas_CO2) then
+              !Interpolated from Hurly et al., (2007)
+              !valid between 20 and 1000 K (max 3 percent of error)
+              akk(igas) = 3.072e-4
+              skk(igas) = 0.69
+              molar_mass(igas) = 44.010e-3
+              akk_visc(igas) = 4.5054e-7
+              skk_visc(igas) = 0.6658
+              molar_frac(:,:,igas) = gfrac(igas)
+              here(igas) = .true.
               ! Add more molecules here and the reference PLEASE!
             else 
@@ -205 +215 @@
                    * (1-D(:,nlayer-1))
       C(:,nlayer)  =C(:,nlayer-1)+zt(:,nlayer-1)-zt(:,nlayer) 
-      C(:,nlayer)  =(C(:,nlayer)*lambda(:,nlayer)+phitop) &
+      C(:,nlayer)  =(C(:,nlayer)*lambda(:,nlayer)+phitop_conduc) &
                    / den(:,nlayer)