- Timestamp:
- Nov 14, 2024, 12:57:31 PM (10 days ago)
- Location:
- trunk/LMDZ.COMMON/libf/evolution
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/evolution/changelog.txt
r3514 r3516 483 483 == 13/11/2024 == JBC 484 484 Small correction for the launching script following the r3498 (time step from integer to real). 485 486 == 14/11/2024 == JBC 487 Small corrections related to r3498 (time step from integer to real) for the launching script + making a proper initilization of 'year_bp_ini' in case of there is no orbital evolution. -
trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90
r3498 r3516 62 62 call getin('evol_orbit_pem',evol_orbit_pem) 63 63 64 year_ bp_ini = 0.64 year_earth_bp_ini = 0. 65 65 call getin('year_earth_bp_ini',year_earth_bp_ini) 66 66 year_bp_ini = year_earth_bp_ini/convert_years … … 68 68 var_obl = .true. 69 69 call getin('var_obl',var_obl) 70 write(*,*) 'Does obliquity vary 70 write(*,*) 'Does obliquity vary?',var_obl 71 71 72 72 var_ecc = .true. 73 73 call getin('var_ecc',var_ecc) 74 write(*,*) 'Does eccentricity vary 74 write(*,*) 'Does eccentricity vary?',var_ecc 75 75 76 76 var_lsp = .true. 77 77 call getin('var_lsp',var_lsp) 78 write(*,*) 'Does Ls peri vary 78 write(*,*) 'Does Ls peri vary?',var_lsp 79 79 80 80 !#---------- Stopping criteria parameters ----------# 81 81 Max_iter_pem = 100000000 82 82 call getin('Max_iter_pem',Max_iter_pem) 83 write(*,*) 'Max_iter_pem =',Max_iter_pem 83 84 84 85 h2o_ice_crit = 0.2 85 86 call getin('h2o_ice_crit',h2o_ice_crit) 87 write(*,*) 'h2o_ice_crit =',h2o_ice_crit 86 88 87 89 co2_ice_crit = 0.2 88 90 call getin('co2_ice_crit',co2_ice_crit) 91 write(*,*) 'co2_ice_crit =',co2_ice_crit 89 92 90 93 ps_criterion = 0.15 91 94 call getin('ps_criterion',ps_criterion) 95 write(*,*) 'ps_criterion =',ps_criterion 92 96 93 97 dt = 1. … … 128 132 call abort_physic(modname,"Ice table must be used when soil_pem = T",1) 129 133 endif 130 if (icetable_equilibrium . or. icetable_dynamic) then134 if (icetable_equilibrium .and. icetable_dynamic) then 131 135 write(*,*) 'Ice table is asked to be computed both by the equilibrium and dynamic method.' 132 136 write(*,*) 'The dynamic method is then chosen.' 137 icetable_dynamic = .false. 133 138 endif 134 139 -
trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh
r3514 r3516 52 52 53 53 if [ -v n_mars_years ] && [ ! -z "$n_mars_years" ]; then 54 if [ $ n_mars_years -lt 0.]; then54 if [ $(echo "$n_mars_years < 0." | bc -l) ]; then 55 55 echo "Error: the value of 'n_mars_years' must be >0!" 56 56 errlaunch 57 57 fi 58 58 elif [ -v n_earth_years ] && [ ! -z "$n_earth_years" ]; then 59 if [ $ n_earth_years -lt 0.]; then59 if [ $(echo "$n_earth_years < 0." | bc -l) ]; then 60 60 echo "Error: the value of 'n_earth_years' must be >0!" 61 61 errlaunch -
trunk/LMDZ.COMMON/libf/evolution/pem.F90
r3512 r3516 1012 1012 ! II_f Update the tendencies 1013 1013 !------------------------ 1014 call recomp_tend_co2_slope(ngrid,nslope,timelen,d_co2ice,d_co2ice_ini,co2_ice,emis,vmr_co2_PCM,vmr_co2_PEM_phys,ps_timeseries, & 1015 global_avg_press_PCM,global_avg_press_new) 1014 call recomp_tend_co2_slope(ngrid,nslope,timelen,d_co2ice,d_co2ice_ini,co2_ice,emis,vmr_co2_PCM,vmr_co2_PEM_phys,ps_timeseries,global_avg_press_PCM,global_avg_press_new) 1016 1015 1017 1016 !------------------------
Note: See TracChangeset
for help on using the changeset viewer.