Makefile (6388B)
1 NPROC != awk '/^cpu cores/ {print $$4; exit}' /proc/cpuinfo 2 SLSERIES = sealvl.nc 3 INIT_T_END = 1e4 4 T_END = 1e4 5 #PISM_RES = -Mx 301 -My 3 -Mz 15 6 PISM_RES = -Mx 601 -My 3 -Mz 15 7 PISM_COMMON = \ 8 -stress_balance ssa+sia\ 9 -gradient eta\ 10 -cfbc -part_grid -kill_icebergs -subgl\ 11 -front_retreat_file MISMIP_boot_1a_M1_A1.nc\ 12 -options_left\ 13 -plastic_phi 24\ 14 -till_cohesion 0\ 15 -hydrology routing\ 16 -bootstrapping.defaults.geothermal_flux 70e-3\ 17 -bed_def lc\ 18 -sia_e 4.5 -ssa_e 0.512\ 19 -stress_balance.sia.max_diffusivity 1e4\ 20 -backup_size big\ 21 -log_view\ 22 23 OUT_VARS =\ 24 -extra_vars \ 25 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 26 27 OUT_COMMON =\ 28 -extra_times 0:50:${T_END}\ 29 -ts_times 0:50:${T_END}\ 30 -ys 0 -ye ${T_END}\ 31 32 TILL_FLUX_OPTS =\ 33 -till_flux\ 34 -till_flux_empirical\ 35 -till_sea_floor_diffusivity 6.341958396752918e-07 \ 36 37 all: \ 38 ex_deltaSL-linear_1d-till-evol.pdf\ 39 ex_deltaSL-linear_1d-profile.pdf\ 40 ex_constant-linear_1d-till-evol.pdf\ 41 ex_constant-linear_1d-profile.pdf\ 42 ex_nt-deltaSL-linear_1d-profile.pdf\ 43 ex_nt-constant-linear_1d-profile.pdf\ 44 45 ex_deltaSL-linear_1d-till-evol.pdf: ex_deltaSL-linear_1d.nc plot-till-evolution.py 46 ./plot-till-evolution.py ex_deltaSL-linear_1d.nc 47 48 ex_deltaSL-linear_1d-profile.pdf: ex_deltaSL-linear_1d.nc plot-time-series.py 49 ./plot-time-series.py ex_deltaSL-linear_1d.nc 50 51 ex_constant-linear_1d-till-evol.pdf: ex_constant-linear_1d.nc plot-till-evolution.py 52 ./plot-till-evolution.py ex_constant-linear_1d.nc 53 54 ex_constant-linear_1d-profile.pdf: ex_constant-linear_1d.nc plot-time-series.py 55 ./plot-time-series.py ex_constant-linear_1d.nc 56 57 ex_nt-deltaSL-linear_1d-profile.pdf: ex_nt-deltaSL-linear_1d.nc plot-time-series.py 58 ./plot-time-series.py ex_nt-deltaSL-linear_1d.nc 59 60 ex_nt-constant-linear_1d-profile.pdf: ex_nt-constant-linear_1d.nc plot-time-series.py 61 ./plot-time-series.py ex_nt-constant-linear_1d.nc 62 63 ex_deltaSL-linear_1d.nc: deltaSL-linear.nc 64 flowline.py -o $@ --collapse -d y ex_deltaSL-linear.nc 65 66 ex_test_1d.nc: test.nc 67 flowline.py -o $@ --collapse -d y ex_test.nc 68 69 ex_constant-linear_1d.nc: constant-linear.nc 70 flowline.py -o $@ --collapse -d y ex_constant-linear.nc 71 72 ex_nt-deltaSL-linear_1d.nc: nt-deltaSL-linear.nc 73 flowline.py -o $@ --collapse -d y ex_nt-deltaSL-linear.nc 74 75 ex_nt-constant-linear_1d.nc: nt-constant-linear.nc 76 flowline.py -o $@ --collapse -d y ex_nt-constant-linear.nc 77 78 deltaSL-linear.nc: init-linear.nc ${SLSERIES} 79 mpiexec -n ${NPROC} \ 80 pismr -i init-linear.nc\ 81 ${PISM_COMMON}\ 82 ${OUT_COMMON}\ 83 ${OUT_VARS},tilldeposit,tillseafloorchange,tillseafloordeposit,tillflux_mag,tillflux\ 84 -sea_level constant,delta_sl -ocean_delta_sl_file ${SLSERIES}\ 85 ${TILL_FLUX_OPTS}\ 86 -extra_file ex_$@\ 87 -ts_file ts_$@\ 88 -o $@ 89 90 # from constant-linear.nc 91 test.nc: init-linear.nc 92 mpiexec -n ${NPROC} \ 93 pismr -i init-linear.nc\ 94 ${PISM_COMMON}\ 95 -extra_times 0:10:${T_END}\ 96 -ts_times 0:10:${T_END}\ 97 -ys 0 -ye 1e2\ 98 ${OUT_VARS},tilldeposit,tillseafloorchange,tillseafloordeposit,tillflux_mag,tillflux\ 99 ${TILL_FLUX_OPTS}\ 100 -extra_file ex_$@\ 101 -ts_file ts_$@\ 102 -o $@ 103 104 constant-linear.nc: init-linear.nc 105 mpiexec -n ${NPROC} \ 106 pismr -i init-linear.nc\ 107 ${PISM_COMMON}\ 108 ${OUT_COMMON}\ 109 ${OUT_VARS},tilldeposit,tillseafloorchange,tillseafloordeposit,tillflux_mag,tillflux\ 110 ${TILL_FLUX_OPTS}\ 111 -extra_file ex_$@\ 112 -ts_file ts_$@\ 113 -o $@ 114 115 nt-deltaSL-linear.nc: init-linear.nc ${SLSERIES} 116 mpiexec -n ${NPROC} \ 117 pismr -i init-linear.nc\ 118 ${PISM_COMMON}\ 119 ${OUT_COMMON}\ 120 ${OUT_VARS}\ 121 -sea_level constant,delta_sl -ocean_delta_sl_file ${SLSERIES}\ 122 -extra_file ex_$@\ 123 -ts_file ts_$@\ 124 -o $@ 125 126 nt-constant-linear.nc: init-linear.nc 127 mpiexec -n ${NPROC} \ 128 pismr -i init-linear.nc\ 129 ${PISM_COMMON}\ 130 ${OUT_COMMON}\ 131 ${OUT_VARS}\ 132 -extra_file ex_$@\ 133 -ts_file ts_$@\ 134 -o $@ 135 136 ${SLSERIES}: sealvl.py 137 printf '0\t0\n${T_END}\t100\n' | ./sealvl.py 138 139 init-linear.nc: MISMIP_boot_1a_M1_A1.nc 140 mpiexec -n ${NPROC} \ 141 pismr -i MISMIP_boot_1a_M1_A1.nc\ 142 -bootstrap\ 143 ${PISM_RES}\ 144 -Lz 5000\ 145 ${PISM_COMMON}\ 146 -ys -${INIT_T_END} -ye 0\ 147 -ts_file ts_$@\ 148 -ts_times -${INIT_T_END}:50:0\ 149 -o $@\ 150 151 MISMIP_boot_1a_M1_A1.nc: prepare.py MISMIP.py 152 ./prepare.py -o $@ -e 1a -m 1 153 154 clean: 155 rm -f *.nc *.nc~ *.pdf ${SLSERIES} SSAFD_kspdivergederror.petsc* 156 157 .PHONY: all clean