tMakefile - pism-exp-gsw - ice stream and sediment transport experiments
(HTM) git clone git://src.adamsgaard.dk/pism-exp-gsw
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) README
(DIR) LICENSE
---
tMakefile (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