commit abc2b96ee0c5daca9e2f7789177121c4c7b09ac2
parent ee096824646daeb8bdebc0b6fd51bc356b345595
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 3 Dec 2021 10:37:30 +0100
Makefile: dump all possible vars in ex file
Diffstat:
1 file changed, 7 insertions(+), 5 deletions(-)
diff --git a/Makefile b/Makefile
@@ -2,6 +2,7 @@ NPROC != awk '/^cpu cores/ {print $$4; exit}' /proc/cpuinfo
SLSERIES = sealvl.nc
INIT_T_END = 1e4
T_END = 1e3
+PISM_RES = -Mx 301 -My 3 -Mz 15
PISM_COMMON = \
-stress_balance ssa+sia\
-gradient eta\
@@ -14,7 +15,8 @@ PISM_COMMON = \
-sia_e 4.5 -ssa_e 0.512\
-bed_def lc\
-stress_balance.sia.max_diffusivity 1e4\
- -extra_vars thk,topg,effbwp,velbar_mag,flux_mag,mask,dHdt,usurf,hardav,velbase_mag,nuH,tauc,taud,taub,flux_divergence,cell_grounded_fraction\
+ -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\
-backup_size big\
-o_order zyx\
@@ -42,17 +44,17 @@ deltaSL-linear.nc: init-linear.nc ${SLSERIES} Makefile
${SLSERIES}: sealvl.py Makefile
printf '0\t0\n${T_END}\t200\n' | ./sealvl.py
-init-linear.nc: MISMIP_boot_1a_M1_A1.nc Makefile
+init-linear.nc: MISMIP_boot_1a_M1_A1.nc
mpiexec -n ${NPROC} \
pismr -i MISMIP_boot_1a_M1_A1.nc\
-bootstrap\
- -Mx 301 -My 3 -Mz 15\
+ ${PISM_RES}\
-Lz 5000\
${PISM_COMMON}\
-front_retreat_file MISMIP_boot_1a_M1_A1.nc\
- -ys 0 -ye ${INIT_T_END}\
+ -ys -${INIT_T_END} -ye 0\
-ts_file ts_$@\
- -ts_times 0:50:${INIT_T_END}\
+ -ts_times -${INIT_T_END}:50:0\
-o $@\
MISMIP_boot_1a_M1_A1.nc: prepare.py MISMIP.py