commit 25e1f8797bdf47bc4b6ac637295ba4d3f1f70123
parent 4617d9994263057846e20bbb18d4891bf633aa93
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Tue, 7 Dec 2021 15:01:25 +0100
Makefile: add constant sea-level experiment
Diffstat:
| M | Makefile | | | 26 | +++++++++++++++++++++++--- |
1 file changed, 23 insertions(+), 3 deletions(-)
diff --git a/Makefile b/Makefile
@@ -1,7 +1,7 @@
NPROC != awk '/^cpu cores/ {print $$4; exit}' /proc/cpuinfo
SLSERIES = sealvl.nc
INIT_T_END = 1e4
-T_END = 1e3
+T_END = 1e4
PISM_RES = -Mx 301 -My 3 -Mz 15
PISM_COMMON = \
-stress_balance ssa+sia\
@@ -17,13 +17,14 @@ PISM_COMMON = \
-sia_e 4.5 -ssa_e 0.512\
-stress_balance.sia.max_diffusivity 1e4\
-extra_vars \
- basal_mass_flux_floating,basal_mass_flux_grounded,basal_melt_rate_grounded,bedtoptemp,beta,bfrict,bheatflx,bmelt,bwat,bwatvel,bwp,bwprel,cell_grounded_fraction,climatic_mass_balance,cts,dHdt,dbdt,deviatoric_stresses,diffusivity,diffusivity_staggered,effbwp,effective_viscosity,elastic_bed_displacement,enthalpy,enthalpybase,enthalpysurf,flux,flux_divergence,flux_mag,flux_staggered,grounding_line_flux,h_x,h_y,hardav,hardness,heat_flux_from_bedrock,height_above_flotation,hydraulic_potential,ice_area_specific_volume,ice_base_elevation,ice_margin_pressure_difference,ice_mass,ice_surface_liquid_water_fraction,ice_surface_temp,lat,liqfrac,lon,mask,nuH,pressure,rank,schoofs_theta,sea_level,sftflf,sftgif,sftgrf,shelfbmassflux,shelfbtemp,strain_rates,strainheat,subglacial_water_flux_mag,subglacial_water_input_rate,surface_layer_mass,surface_layer_thickness,taub,taub_mag,tauc,taud,taud_mag,tauxz,tauyz,temp,temp_pa,tempbase,tempicethk,tempicethk_basal,temppabase,tempsurf,tendency_of_ice_amount,tendency_of_ice_amount_due_to_basal_mass_flux,tendency_of_ice_amount_due_to_calving,tendency_of_ice_amount_due_to_conservation_error,tendency_of_ice_amount_due_to_discharge,tendency_of_ice_amount_due_to_flow,tendency_of_ice_amount_due_to_surface_mass_flux,tendency_of_ice_mass,tendency_of_ice_mass_due_to_basal_mass_flux,tendency_of_ice_mass_due_to_calving,tendency_of_ice_mass_due_to_conservation_error,tendency_of_ice_mass_due_to_discharge,tendency_of_ice_mass_due_to_flow,tendency_of_ice_mass_due_to_surface_mass_flux,tendency_of_subglacial_water_mass,tendency_of_subglacial_water_mass_at_domain_boundary,tendency_of_subglacial_water_mass_at_grounded_margins,tendency_of_subglacial_water_mass_at_grounding_line,tendency_of_subglacial_water_mass_due_to_conservation_error,tendency_of_subglacial_water_mass_due_to_flow,tendency_of_subglacial_water_mass_due_to_input,thk,thk_bc_mask,thksmooth,tillphi,tillwat,topg,topg_sl_adjusted,topgsmooth,usurf,uvel,vel_bc_mask,vel_bc_values,velbar,velbar_mag,velbase,velbase_mag,velsurf,velsurf_mag,viscous_bed_displacement,vonmises_stress,vvel,wallmelt,wvel,wvel_rel,wvelbase,wvelsurf,tillflux_mag,tillflux\
+ basal_mass_flux_floating,basal_mass_flux_grounded,basal_melt_rate_grounded,bedtoptemp,beta,bfrict,bheatflx,bmelt,bwat,bwatvel,bwp,bwprel,cell_grounded_fraction,climatic_mass_balance,cts,dHdt,dbdt,deviatoric_stresses,diffusivity,diffusivity_staggered,effbwp,effective_viscosity,elastic_bed_displacement,enthalpy,enthalpybase,enthalpysurf,flux,flux_divergence,flux_mag,flux_staggered,grounding_line_flux,h_x,h_y,hardav,hardness,heat_flux_from_bedrock,height_above_flotation,hydraulic_potential,ice_area_specific_volume,ice_base_elevation,ice_margin_pressure_difference,ice_mass,ice_surface_liquid_water_fraction,ice_surface_temp,lat,liqfrac,lon,mask,nuH,pressure,rank,schoofs_theta,sea_level,sftflf,sftgif,sftgrf,shelfbmassflux,shelfbtemp,strain_rates,strainheat,subglacial_water_flux_mag,subglacial_water_input_rate,surface_layer_mass,surface_layer_thickness,taub,taub_mag,tauc,taud,taud_mag,tauxz,tauyz,temp,temp_pa,tempbase,tempicethk,tempicethk_basal,temppabase,tempsurf,tendency_of_ice_amount,tendency_of_ice_amount_due_to_basal_mass_flux,tendency_of_ice_amount_due_to_calving,tendency_of_ice_amount_due_to_conservation_error,tendency_of_ice_amount_due_to_discharge,tendency_of_ice_amount_due_to_flow,tendency_of_ice_amount_due_to_surface_mass_flux,tendency_of_ice_mass,tendency_of_ice_mass_due_to_basal_mass_flux,tendency_of_ice_mass_due_to_calving,tendency_of_ice_mass_due_to_conservation_error,tendency_of_ice_mass_due_to_discharge,tendency_of_ice_mass_due_to_flow,tendency_of_ice_mass_due_to_surface_mass_flux,tendency_of_subglacial_water_mass,tendency_of_subglacial_water_mass_at_domain_boundary,tendency_of_subglacial_water_mass_at_grounded_margins,tendency_of_subglacial_water_mass_at_grounding_line,tendency_of_subglacial_water_mass_due_to_conservation_error,tendency_of_subglacial_water_mass_due_to_flow,tendency_of_subglacial_water_mass_due_to_input,thk,thk_bc_mask,thksmooth,tillphi,tillwat,topg,topg_sl_adjusted,topgsmooth,usurf,uvel,vel_bc_mask,vel_bc_values,velbar,velbar_mag,velbase,velbase_mag,velsurf,velsurf_mag,viscous_bed_displacement,vonmises_stress,vvel,wallmelt,wvel,wvel_rel,wvelbase,wvelsurf,tilldeposit,tillflux_mag,tillflux\
-backup_size big\
-o_order zyx\
-log_view\
all: \
ex_deltaSL-linear_1d-profile.pdf\
+ ex_constant-linear_1d-profile.pdf\
ex_deltaSL-linear_1d-profile.pdf: ex_deltaSL-linear_1d.nc plot-time-series.py
./plot-time-series.py ex_deltaSL-linear_1d.nc
@@ -31,12 +32,30 @@ ex_deltaSL-linear_1d-profile.pdf: ex_deltaSL-linear_1d.nc plot-time-series.py
ex_deltaSL-linear_1d.nc: deltaSL-linear.nc
flowline.py -o $@ --collapse -d y ex_deltaSL-linear.nc
+ex_constant-linear_1d-profile.pdf: ex_constant-linear_1d.nc plot-time-series.py
+ ./plot-time-series.py ex_constant-linear_1d.nc
+
+ex_constant-linear_1d.nc: constant-linear.nc
+ flowline.py -o $@ --collapse -d y ex_constant-linear.nc
+
deltaSL-linear.nc: init-linear.nc ${SLSERIES}
mpiexec -n ${NPROC} \
pismr -i init-linear.nc\
${PISM_COMMON}\
- -till_flux\
-sea_level constant,delta_sl -ocean_delta_sl_file ${SLSERIES}\
+ -till_flux\
+ -extra_file ex_$@\
+ -extra_times 0:50:${T_END}\
+ -ts_file ts_$@\
+ -ts_times 0:50:${T_END}\
+ -ys 0 -ye ${T_END}\
+ -o $@
+
+constant-linear.nc: init-linear.nc ${SLSERIES}
+ mpiexec -n ${NPROC} \
+ pismr -i init-linear.nc\
+ ${PISM_COMMON}\
+ -till_flux\
-extra_file ex_$@\
-extra_times 0:50:${T_END}\
-ts_file ts_$@\
@@ -44,6 +63,7 @@ deltaSL-linear.nc: init-linear.nc ${SLSERIES}
-ys 0 -ye ${T_END}\
-o $@
+
${SLSERIES}: sealvl.py
printf '0\t0\n${T_END}\t200\n' | ./sealvl.py
