Changeset 3372 for trunk/LMDZ.PLUTO

Timestamp:

Jun 13, 2024, 2:52:09 PM (5 months ago)

Author:

afalco

Message:

Pluto PCM:
Cooling parameters in .def.
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callkeys_mod.F90 (modified) (1 diff)
• cooling_hcn_c2h2.F90 (modified) (3 diffs)
• inifis_mod.F90 (modified) (2 diffs)

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

@@ -202 +202 @@
       real,save    :: szangle
 !$OMP THREADPRIVATE(global1d,szangle)
+      real,save    :: alpha_top ! cooling constant at top of atmosphere
+      real,save    :: pref      ! pressure at mid transition fo alpha_top (Pa)
+      real,save    :: deltap    ! width of transition to alpha_top (Pa)  
+!$OMP THREADPRIVATE(alpha_top,pref,deltap)
+      
 
       integer,save :: iddist

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

@@ -1 +1 @@
       subroutine cooling_hcn_c2h2(ngrid,nlayer,pplay,pt,dtlw)
+
+      use callkeys_mod, only: alpha_top, deltap, pref 
+      
           
       implicit none
@@ -33 +36 @@
       INTEGER l,ig
       REAL lonw
-      REAL alpha, alpha_top
-      REAL pref, deltap 
+      REAL alpha !, alpha_top
       REAL transition
       REAL BB 
@@ -40 +42 @@
 
       lonw = 14.e-6  ! 14um
-      alpha_top=5.e-11  ! 1.e-13  ! cooling constant at top of atmosphere
-      !pref = 0.12      ! pressure at mid transition fo alpha_top (Pa)
-      pref = 0.02      ! pressure at mid transition fo alpha_top (Pa)
-      deltap = 0.1    ! width of transition to alpha_top (Pa)  
+      ! alpha_top=5.e-11  ! 1.e-13  ! cooling constant at top of atmosphere
+      ! pref = 0.02      ! pressure at mid transition fo alpha_top (Pa)
+      ! deltap = 0.1    ! width of transition to alpha_top (Pa)  
       !lonw = 11.e-6  ! 14um
-      !alpha_top=1.e-9  ! 1.e-13  ! cooling constant at top of atmosphere
-      !pref = 0.06      ! pressure at mid transition fo alpha_top (Pa)
-      !deltap = 0.1    ! width of transition to alpha_top (Pa)  
-
+      
 !     transition = 0 if p>pref+deltap/2 and 1 if p< pref-deltap/2
       DO l = 1, nlayer

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

@@ -855 +855 @@
      fdch4_ampl=200.
      fdch4_maxice=100.
-     call getin_p("  fdch4_latn",fdch4_latn)
-     call getin_p("  fdch4_lats",fdch4_lats)
-     call getin_p("  fdch4_lone",fdch4_lone)
-     call getin_p("  fdch4_lonw",fdch4_lonw)
-     call getin_p("  fdch4_depalb",fdch4_depalb)
-     call getin_p("  fdch4_finalb",fdch4_finalb)
-     call getin_p("  fdch4_maxalb",fdch4_maxalb)
-     call getin_p("  fdch4_maxice",fdch4_maxice)
-     call getin_p("  fdch4_ampl",fdch4_ampl)
+     call getin_p("fdch4_latn",fdch4_latn)
+     call getin_p("fdch4_lats",fdch4_lats)
+     call getin_p("fdch4_lone",fdch4_lone)
+     call getin_p("fdch4_lonw",fdch4_lonw)
+     call getin_p("fdch4_depalb",fdch4_depalb)
+     call getin_p("fdch4_finalb",fdch4_finalb)
+     call getin_p("fdch4_maxalb",fdch4_maxalb)
+     call getin_p("fdch4_maxice",fdch4_maxice)
+     call getin_p("fdch4_ampl",fdch4_ampl)
      if (is_master)write(*,*)trim(rname)//&
       " Values for albedo feedback = ",fdch4_latn,&
@@ -1099 +1099 @@
      if (is_master)write(*,*)trim(rname)//&
       " ITH2Od = ",ITH2Od
+
+!************** COOLING ***************
+
+     alpha_top=5e-11 ! default value
+     call getin_p("alpha_top",alpha_top)
+     if (is_master)write(*,*)trim(rname)//&
+      " alpha_top = ",alpha_top
+     pref=0.02 ! default value
+     call getin_p("pref",pref)
+     if (is_master)write(*,*)trim(rname)//&
+      " pref = ",pref
+     deltap=0.1 ! default value
+     call getin_p("deltap",deltap)
+     if (is_master)write(*,*)trim(rname)//&
+      " deltap = ",deltap
 
 !=================================