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

Timestamp:

Jan 28, 2026, 3:06:52 PM (4 weeks 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/physiq_mod.F90 (modified) (5 diffs)

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

@@ -71 +71 @@
                               ok_slab_ocean, photochem, rings_shadow, rmean, &
                               season, sedimentation,generic_condensation, &
-                              specOLR, callthermos, callvolcano, &
+                              specOLR, callvolcano, &
+                              callthermos, callconduc, callmolvis, &
                               startphy_file, testradtimes, tlocked, &
                               tracer, UseTurbDiff, water, watercond, &
@@ -83 +84 @@
       use callcorrk_mod, only: callcorrk
       use conduction_mod, only: conduction
+      use molvis_mod, only: molvis
       use volcano_mod, only: volcano
       use pindex_mod, only: pindex
@@ -382 +384 @@
       real zdhdif(ngrid,nlayer)                             ! Turbdiff/vdifc routines.
       real zdhadj(ngrid,nlayer)                             ! Convadj routine.
+      real zdumolvis(ngrid,nlayer)                          ! Molecular viscosity (thermosphere)
+      real zdvmolvis(ngrid,nlayer)                          ! Molecular viscosity (thermosphere)
 
       ! For Pressure and Mass :
@@ -541 +545 @@
          zdtlw(:,:) = 0.0
          zdtconduc(:,:) = 0.0
+         zdumolvis(:,:) = 0.0
+         zdvmolvis(:,:) = 0.0
 
 !        Initialize fixed variable mu
@@ -2384 +2390 @@
 ! -------------------------
       if (callthermos) then
-        call conduction(ngrid,nlayer,nq,ptimestep,pplay,pplev,  &
-                          pt,tsurf,zzlev,zzlay,muvar,pq,firstcall,zdtconduc)
-        pdt(1:ngrid,1:nlayer)=pdt(1:ngrid,1:nlayer)+zdtconduc(1:ngrid,1:nlayer)
+
+        if (callconduc) then
+          call conduction(ngrid,nlayer,nq,ptimestep,pplay,pplev,  &
+                            pt,tsurf,zzlev,zzlay,muvar,pq,firstcall,zdtconduc)
+          pdt(1:ngrid,1:nlayer)=pdt(1:ngrid,1:nlayer)+zdtconduc(1:ngrid,1:nlayer)
+        endif
+      
+        if (callmolvis) then
+          ! u
+          call molvis(ngrid,nlayer,ptimestep,pplay,pplev,pt,pdt, &
+                      pu,zzlev,zzlay,zdumolvis)
+          pdu(1:ngrid,1:nlayer) = pdu(1:ngrid,1:nlayer) + zdumolvis(1:ngrid,1:nlayer)
+          ! v
+          call molvis(ngrid,nlayer,ptimestep,pplay,pplev,pt,pdt, &
+                      pv,zzlev,zzlay,zdvmolvis)
+          pdv(1:ngrid,1:nlayer) = pdv(1:ngrid,1:nlayer) + zdvmolvis(1:ngrid,1:nlayer)
+        endif
+
       endif ! of if (callthermos)