Changeset 3062
- Timestamp:
- Sep 28, 2023, 3:32:59 PM (15 months ago)
- Location:
- trunk/LMDZ.MARS
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/changelog.txt
r3061 r3062 4246 4246 - Explicit comments in several sections of the code with these infos 4247 4247 4248 == 28/09/2023 == JN 4249 Following previous commit, correction of a spelling mistake "temp_dependant_m" 4250 is actually "temp_dependent_m". Sorry about that. -
trunk/LMDZ.MARS/libf/phymars/callkeys.h
r2966 r3062 17 17 & ,CLFvarying,satindexco2,rdstorm,topflows,calllott_nonoro & 18 18 & ,latentheat_surfwater,gwd_convective_source,startphy_file & 19 & ,hdo,hdofrac,cst_cap_albedo,temp_depend ant_m,refill_watercap &19 & ,hdo,hdofrac,cst_cap_albedo,temp_dependent_m,refill_watercap & 20 20 & ,cloud_adapt_ts 21 21 !$OMP THREADPRIVATE(/callkeys_l/) … … 75 75 logical latentheat_surfwater ! latent heat release from ground water ice sublimation/condensation 76 76 logical cst_cap_albedo ! polar cap albedo remains unchanged by water frost deposition 77 logical temp_depend ant_m ! temperature-dependant water contact parameter77 logical temp_dependent_m ! temperature-dependent water contact parameter 78 78 logical refill_watercap ! h2o_ice_s is converted to watercap when above threshold 79 79 logical cloud_adapt_ts ! adaptative timestep for cloud microphysics -
trunk/LMDZ.MARS/libf/phymars/conf_phys.F
r3061 r3062 746 746 write(*,*) "water contact parameter ?" 747 747 mteta = 0.95 ! default value 748 temp_depend ant_m = .false. ! default value749 call getin_p("temp_depend ant_m",temp_dependant_m)750 if (temp_depend ant_m) then !(JN 2023)751 print*,'You have chosen a temperature-depend ant water'748 temp_dependent_m = .false. ! default value 749 call getin_p("temp_dependent_m",temp_dependent_m) 750 if (temp_dependent_m) then !(JN 2023) 751 print*,'You have chosen a temperature-dependent water' 752 752 print*,'contact parameter ! From Maattanen et al. 2014' 753 else if (.not.temp_depend ant_m) then753 else if (.not.temp_dependent_m) then 754 754 print*,'Water contact parameter is constant' 755 755 call getin_p("mteta",mteta) … … 769 769 ! global water cycle is nuts !) 770 770 771 if ((cloud_adapt_ts.and..not.temp_depend ant_m).or.772 & (temp_depend ant_m.and..not.cloud_adapt_ts)) then771 if ((cloud_adapt_ts.and..not.temp_dependent_m).or. 772 & (temp_dependent_m.and..not.cloud_adapt_ts)) then 773 773 print*,'Water cycle v6 : cloud_adapt_ts and ' 774 print*,'temp_depend ant_m must be used together !'774 print*,'temp_dependent_m must be used together !' 775 775 print*,'Otherwise water cycle is off. At least one ' 776 776 print*,'of these flags is turned off. ' -
trunk/LMDZ.MARS/libf/phymars/nuclea.F
r3061 r3062 59 59 mtetalocal = mteta !! use mtetalocal for better performance 60 60 61 IF (temp_depend ant_m) THEN61 IF (temp_dependent_m) THEN 62 62 c J.Naar - sep 2023 : 63 63 c Simple linear parametrisation from Maattaanen 2014 … … 67 67 c With adapt_ts need to use this relation (same paper, tanh): 68 68 mtetalocal = 0.469+((0.972-0.469)*tanh((temp/158.282)**4.244)) 69 ENDIF ! (temp_depend ant_m) THEN69 ENDIF ! (temp_dependent_m) THEN 70 70 cccccccccccccccccccccccccccccccccccccccccccccccccc 71 71 ccccccccccc ESSAIS TN MTETA = F (T) cccccccccccccc … … 83 83 cccccccccccccccccccccccccccccccccccccccccccccccccc 84 84 cccccccccccccccccccccccccccccccccccccccccccccccccc 85 IF (firstcall.and.temp_depend ant_m) THEN85 IF (firstcall.and.temp_dependent_m) THEN 86 86 print*, ' ' 87 87 print*, 'dear user, please keep in mind that' … … 92 92 print*, ' ' 93 93 firstcall=.false. 94 ELSE IF (firstcall.and.(.not.(temp_depend ant_m))) THEN94 ELSE IF (firstcall.and.(.not.(temp_dependent_m))) THEN 95 95 print*, ' ' 96 96 print*, 'dear user, please keep in mind that'
Note: See TracChangeset
for help on using the changeset viewer.