Changeset 3437 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Sep 25, 2024, 1:50:44 PM (2 months ago)

Author:

emillour

Message:

Generic PCM Integrate update from newton branch by LT:

• remove call to old newtrelax routine and replace it with a new one, newtcool_MOCHA, tuned to the MOCHA intercomparison project

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• inifis_mod.F90 (modified) (2 diffs)
• newton_cooling_hotJ.F90 (added)
• physiq_mod.F90 (modified) (3 diffs)

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

@@ -27 +27 @@
   use ioipsl_getin_p_mod, only : getin_p
   use mod_phys_lmdz_para, only : is_parallel, is_master, bcast
+  use newton_cooling_hotJ, only: planetary_suffix
 
 !=======================================================================
@@ -406 +407 @@
      tau_relax=tau_relax*24*3600 ! convert Earth days --> seconds
 
+     if (is_master) write(*,*) trim(rname)//&
+     ": planetary suffix for loading file in Newtonian cooling HotJ ?"
+     planetary_suffix='default_value' ! default value
+     call getin_p("planetary_suffix",planetary_suffix)
+     if (is_master) write(*,*) trim(rname)//": planetary_suffix",planetary_suffix
+     
      if (is_master) write(*,*)trim(rname)//&
        ": call thermal conduction in the soil ?"

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

@@ -89 +89 @@
       use condensation_generic_mod, only: condensation_generic
       use datafile_mod, only: datadir
+      use newton_cooling_hotJ, only: newtcool_MOCHA ! LT, adding for MOCHA protocol
+      
 #ifndef MESOSCALE
       use vertical_layers_mod, only: presnivs, pseudoalt
@@ -1149 +1151 @@
 ! II.b Call Newtonian cooling scheme
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-               call newtrelax(ngrid,nlayer,mu0,sinlat,zpopsk,pt,pplay,pplev,dtrad,firstcall)
+               ! call newtrelax(ngrid,nlayer,mu0,sinlat,zpopsk,pt,pplay,pplev,dtrad,firstcall)
+               call newtcool_MOCHA(ngrid,nlayer,coslon,coslat,pt,pplay,firstcall,lastcall,dtrad)
 
                zdtsurf(1:ngrid) = +(pt(1:ngrid,1)-tsurf(1:ngrid))/ptimestep
                ! e.g. surface becomes proxy for 1st atmospheric layer ?
+               fluxrad_sky(1:ngrid) = 0. ! otherwise, not initialized!
 
             else
@@ -1181 +1185 @@
          ! Transformation of the radiative tendencies
          ! ------------------------------------------
+
          zplanck(1:ngrid)=tsurf(1:ngrid)*tsurf(1:ngrid)
          zplanck(1:ngrid)=emis(1:ngrid)*sigma*zplanck(1:ngrid)*zplanck(1:ngrid)