- Timestamp:
- Jun 13, 2024, 2:52:09 PM (6 months ago)
- Location:
- trunk/LMDZ.PLUTO/libf/phypluto
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.PLUTO/libf/phypluto/callkeys_mod.F90
r3361 r3372 202 202 real,save :: szangle 203 203 !$OMP THREADPRIVATE(global1d,szangle) 204 real,save :: alpha_top ! cooling constant at top of atmosphere 205 real,save :: pref ! pressure at mid transition fo alpha_top (Pa) 206 real,save :: deltap ! width of transition to alpha_top (Pa) 207 !$OMP THREADPRIVATE(alpha_top,pref,deltap) 208 204 209 205 210 integer,save :: iddist -
trunk/LMDZ.PLUTO/libf/phypluto/cooling_hcn_c2h2.F90
r3329 r3372 1 1 subroutine cooling_hcn_c2h2(ngrid,nlayer,pplay,pt,dtlw) 2 3 use callkeys_mod, only: alpha_top, deltap, pref 4 2 5 3 6 implicit none … … 33 36 INTEGER l,ig 34 37 REAL lonw 35 REAL alpha, alpha_top 36 REAL pref, deltap 38 REAL alpha !, alpha_top 37 39 REAL transition 38 40 REAL BB … … 40 42 41 43 lonw = 14.e-6 ! 14um 42 alpha_top=5.e-11 ! 1.e-13 ! cooling constant at top of atmosphere 43 !pref = 0.12 ! pressure at mid transition fo alpha_top (Pa) 44 pref = 0.02 ! pressure at mid transition fo alpha_top (Pa) 45 deltap = 0.1 ! width of transition to alpha_top (Pa) 44 ! alpha_top=5.e-11 ! 1.e-13 ! cooling constant at top of atmosphere 45 ! pref = 0.02 ! pressure at mid transition fo alpha_top (Pa) 46 ! deltap = 0.1 ! width of transition to alpha_top (Pa) 46 47 !lonw = 11.e-6 ! 14um 47 !alpha_top=1.e-9 ! 1.e-13 ! cooling constant at top of atmosphere 48 !pref = 0.06 ! pressure at mid transition fo alpha_top (Pa) 49 !deltap = 0.1 ! width of transition to alpha_top (Pa) 50 48 51 49 ! transition = 0 if p>pref+deltap/2 and 1 if p< pref-deltap/2 52 50 DO l = 1, nlayer -
trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90
r3353 r3372 855 855 fdch4_ampl=200. 856 856 fdch4_maxice=100. 857 call getin_p(" 858 call getin_p(" 859 call getin_p(" 860 call getin_p(" 861 call getin_p(" 862 call getin_p(" 863 call getin_p(" 864 call getin_p(" 865 call getin_p(" 857 call getin_p("fdch4_latn",fdch4_latn) 858 call getin_p("fdch4_lats",fdch4_lats) 859 call getin_p("fdch4_lone",fdch4_lone) 860 call getin_p("fdch4_lonw",fdch4_lonw) 861 call getin_p("fdch4_depalb",fdch4_depalb) 862 call getin_p("fdch4_finalb",fdch4_finalb) 863 call getin_p("fdch4_maxalb",fdch4_maxalb) 864 call getin_p("fdch4_maxice",fdch4_maxice) 865 call getin_p("fdch4_ampl",fdch4_ampl) 866 866 if (is_master)write(*,*)trim(rname)//& 867 867 " Values for albedo feedback = ",fdch4_latn,& … … 1099 1099 if (is_master)write(*,*)trim(rname)//& 1100 1100 " ITH2Od = ",ITH2Od 1101 1102 !************** COOLING *************** 1103 1104 alpha_top=5e-11 ! default value 1105 call getin_p("alpha_top",alpha_top) 1106 if (is_master)write(*,*)trim(rname)//& 1107 " alpha_top = ",alpha_top 1108 pref=0.02 ! default value 1109 call getin_p("pref",pref) 1110 if (is_master)write(*,*)trim(rname)//& 1111 " pref = ",pref 1112 deltap=0.1 ! default value 1113 call getin_p("deltap",deltap) 1114 if (is_master)write(*,*)trim(rname)//& 1115 " deltap = ",deltap 1101 1116 1102 1117 !=================================
Note: See TracChangeset
for help on using the changeset viewer.