Changeset 3235 for trunk

Timestamp:

Feb 26, 2024, 10:53:39 AM (9 months ago)

Author:

gmilcareck

Message:

Add thermal conduction from thermosphere.

Location:

trunk/LMDZ.GENERIC

Files:

• changelog.txt (modified) (1 diff)
• libf/phystd/callkeys_mod.F90 (modified) (1 diff)
• libf/phystd/conduction.F90 (added)
• libf/phystd/inifis_mod.F90 (modified) (1 diff)
• libf/phystd/physiq_mod.F90 (modified) (5 diffs)

trunk/LMDZ.GENERIC/changelog.txt

@@ -1858 +1858 @@
 Add the possibility to use a fixed vertical molar fraction profile for the 
 collision-induced absorption like the correlated-k.
+
+== 26/02/2024 == GM
+
+Add thermal conduction from thermosphere. The conductivity for each gas is 
+calculated with the interpolated akk*T**skk formula and the total conductivity 
+of a mix of gases is calculated according to Mason & Saxena (1958) semi-empirical
+formulation.

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

@@ -14 +14 @@
       logical,save :: strictboundcia                                     
 !$OMP THREADPRIVATE(strictboundcia)
+      logical,save :: callthermos
+      logical,save :: force_conduction
+      real,save    :: phitop
+      real,save    :: zztop
+      real,save    :: a_coeff
+      real,save    :: s_coeff
+!$OMP THREADPRIVATE(callthermos,phitop,zztop,force_conduction,a_coeff,s_coeff)
 
       logical,save :: enertest

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

@@ -260 +260 @@
      if (is_master) write(*,*) trim(rname)//&
        ": strictboundcorrk = ",strictboundcorrk
+       
+     if (is_master) write(*,*) trim(rname)//&
+       ": call thermosphere ?"
+     callthermos=.false. ! default value
+     call getin_p("callthermos",callthermos)
+     if (is_master) write(*,*) trim(rname)//&
+       ": callthermos = ",callthermos
+     if(callthermos) then
+       if (is_master) write(*,*) trim(rname)//&
+         ": flux from thermosphere ? W/m2"
+       phitop = 0.0 ! default value
+       call getin_p("phitop",phitop)
+       if (is_master) write(*,*) trim(rname)//&
+         ": phitop = ",phitop
+       if (is_master) write(*,*) trim(rname)//&
+         ": Thickness of conduction ? m"
+       zztop = 1000.
+       call getin_p("zztop",zztop)
+       if (is_master) write(*,*) trim(rname)//&
+         ": zztop = ",zztop
+       if (is_master) write(*,*) trim(rname)//&
+         ": Force conduction ? "
+       force_conduction=.false. ! default value
+       call getin_p("force_conduction",force_conduction)
+       if (is_master) write(*,*) trim(rname)//&
+         ": force_conduction = ",force_conduction
+       if(force_conduction) then
+         if (is_master) write(*,*) trim(rname)//&
+           ": a_coeff ?"
+         a_coeff = 0.0 ! default value
+         call getin_p("a_coeff",a_coeff)
+         if (is_master) write(*,*) trim(rname)//&
+           ": a_coeff = ",a_coeff
+         if (is_master) write(*,*) trim(rname)//&
+           ": s_coeff ?"
+         s_coeff = 0.0 ! default value
+         call getin_p("s_coeff",s_coeff)
+         if (is_master) write(*,*) trim(rname)//&
+           ": s_coeff = ",s_coeff
+       endif
+     endif
 
      if (is_master) write(*,*) trim(rname)//&

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

@@ -66 +66 @@
                               ok_slab_ocean, photochem, rings_shadow, rmean, &
                               season, sedimentation,generic_condensation, &
-                              specOLR, &
+                              specOLR, callthermos, &
                               startphy_file, testradtimes, tlocked, &
                               tracer, UseTurbDiff, water, watercond, &
@@ -77 +77 @@
       use phys_state_var_mod
       use callcorrk_mod, only: callcorrk
+      use conduction_mod
       use pindex_mod, only: pindex
       use vdifc_mod, only: vdifc
@@ -308 +309 @@
       real zdtsdif(ngrid)     ! Turbdiff/vdifc routines.
       real zdtsurf_hyd(ngrid) ! Hydrol routine.
+      
+      ! For Thermosphere : (K/s)
+      real zdtconduc(ngrid,nlayer)
 
       ! For Atmospheric Temperatures : (K/s)
@@ -518 +522 @@
          zdtsw(:,:) = 0.0
          zdtlw(:,:) = 0.0
+         zdtconduc(:,:)= 0.0
          
 !        Initialize fixed variable mu
@@ -2311 +2316 @@
 
 
+! -------------------------
+!   IX. Thermosphere
+! -------------------------
+      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)
+      endif ! of if (callthermos)
+
+
 !----------------------------------------------------------------------
 !        Writing NetCDF file  "RESTARTFI" at the end of the run