Changeset 4033 for trunk/LMDZ.GENERIC/libf/phystd/inifis_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/inifis_mod.F90 (modified) (1 diff)

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

@@ -292 +292 @@
      if(callthermos) then
        if (is_master) write(*,*) trim(rname)//&
+        ": call thermospheric conduction ?"
+       callconduc=.true. ! default value
+       call getin_p("conduction",callconduc)
+       if (is_master) write(*,*) trim(rname)//&
+         ": conduction = ",callconduc
+       if (is_master) write(*,*) trim(rname)//&
+         ": call thermospheric viscosity ?"
+         callmolvis=.true. ! default value
+       call getin_p("molvis",callmolvis)
+       if (is_master) write(*,*) trim(rname)//&
+          ": molvis = ",callmolvis
+       if (is_master) write(*,*) trim(rname)//&
          ": flux from thermosphere ? W/m2"
-       phitop = 0.0 ! default value
-       call getin_p("phitop",phitop)
+       phitop_conduc = 0.0 ! default value
+       call getin_p("phitop_conduc",phitop_conduc)
        if (is_master) write(*,*) trim(rname)//&
-         ": phitop = ",phitop
+         ": phitop_conduc = ",phitop_conduc
+        if (is_master) write(*,*) trim(rname)//&
+         ": momentum flux from thermosphere ?"
+       phitop_molvis = 0.0 ! default value
+       call getin_p("phitop_molvis",phitop_molvis)
+        if (is_master) write(*,*) trim(rname)//&
+         ": phitop_molvis = ",phitop_molvis
        if (is_master) write(*,*) trim(rname)//&
          ": Thickness of conduction ? m"