Changeset 4078 for trunk/LMDZ.PLUTO/libf/phypluto

Timestamp:

Feb 18, 2026, 8:45:27 PM (12 days ago)

Author:

tbertrand

Message:

Pluto PCM:
Allowing atmospheric temperatures to be fixed following an input profile (fixed_temp_prof option)
TB

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callkeys_mod.F90 (modified) (2 diffs)
• inifis_mod.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (6 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/callkeys_mod.F90

@@ -62 +62 @@
       logical,save :: kbo
 !$OMP THREADPRIVATE(kbo)
-      logical,save :: cooling
-!$OMP THREADPRIVATE(cooling)
+      logical,save :: cooling,fixed_temp_prof
+!$OMP THREADPRIVATE(cooling,fixed_temp_prof)
       logical,save :: source_haze,hazeconservch4,haze_ch4proffix
 !$OMP THREADPRIVATE(source_haze,hazeconservch4,haze_ch4proffix)
@@ -112 +112 @@
       real,save    :: tau_prechaze
 !$OMP THREADPRIVATE(tau_prechaze)
+      real,save    :: tau_temp
+!$OMP THREADPRIVATE(tau_temp)
       real,save    :: paleoyears
 !$OMP THREADPRIVATE(paleoyears)

trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

@@ -311 +311 @@
 
      if (is_master) write(*,*) trim(rname)//&
+       ": call simple fixed temperature profile"
+     fixed_temp_prof=.false. ! default value
+     call getin_p("fixed_temp_prof",fixed_temp_prof)
+     if (is_master) write(*,*) trim(rname)//": fixed_temp_prof = ",fixed_temp_prof
+     tau_temp=1. ! default value
+     call getin_p("tau_temp",tau_temp)
+
+     if (is_master) write(*,*) trim(rname)//&
        ": prohibit calculations outside corrk T grid?"
      strictboundcorrk=.true. ! default value

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -68 +68 @@
                               tracer, UseTurbDiff,                            &
                               global1d, szangle,                              &
-                              callmufi, callmuclouds, evol1d
+                              callmufi, callmuclouds, evol1d,fixed_temp_prof, &
+                              tau_temp
 
       use check_fields_mod, only: check_physics_fields
@@ -121 +122 @@
 !         II.2.a Option 1 : Call correlated-k radiative transfer scheme.
 !         II.2.b Option 2 : Atmosphere has no radiative effect.
+!         II.2.c Option 3 : Fixed temperature profile (using "Newtonian restoring force"-type equation).
 !
 !      II.3 Gravity wave and subgrid scale topography drag :
@@ -502 +504 @@
       real alpha,lay1 ! coefficients for building layers
       integer iloop
+
+      ! for fixed temperature profile (fixed_temp_prof option) 
+      real,save,allocatable :: tmean(:,:)
+      !!read altitudes and radius
+      integer ifine,Nfine
+      parameter(Nfine=1751)
+      character(len=100) :: file_path
+      real,save :: levdat(Nfine),tempdat(Nfine)
 
       ! flags to trigger extra sanity checks
@@ -1029 +1039 @@
 
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-! II.a Call correlated-k radiative transfer scheme
+! II.2.a Call correlated-k radiative transfer scheme
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                if(kastprof)then
@@ -1099 +1109 @@
             else
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-! II.b Atmosphere has no radiative effect
+! II.2.b Atmosphere has no radiative effect
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                fluxtop_dn(1:ngrid)  = fract(1:ngrid)*mu0(1:ngrid)*Fat1AU/dist_star**2
@@ -1119 +1129 @@
 
             endif ! end of corrk
+
+! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+! II.2.c Converging atm temperature toward a fixed profile with a characteristic timescale
+! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+            if (firstcall.and.fixed_temp_prof) then
+              allocate(tmean(ngrid,nlayer))
+              file_path=trim(datadir)//'/gas_prop/tempNH.txt'
+              open(323,file=file_path,form='formatted')
+              do ifine=1,Nfine
+                read(323,*) levdat(ifine), tempdat(ifine)
+              enddo
+              close(323)
+              print*, 'Read temp file from ',file_path
+              ! levs have been put in km
+              DO ig=1,ngrid
+                CALL interp_line(levdat,tempdat,Nfine,zzlay(ig,:)/1000.,tmean(ig,:),nlayer)
+              enddo
+            endif
+            if (fixed_temp_prof) then
+               DO ig=1,ngrid
+                 dtrad(ig,1:nlayer)=(tmean(ig,1:nlayer)-(pt(ig,1:nlayer)+    &
+                   (pdt(ig,1:nlayer)+dtrad(ig,1:nlayer))*ptimestep))         &
+                   *(1.-exp(-ptimestep/tau_temp))/ptimestep
+               ENDDO
+            endif
 
          endif ! of if(mod(icount-1,iradia).eq.0)