Coupling of ice sheet model PISM with solid Earth model VILMA
Submit the heterogeneous slurm batch job sbatch slurm_vilma3d-pism.cmd, with selected time parameters and SBATCH options. Default is
- timeint=.1 # kyr, meaning 100 years
- ys=-246.0 # kyr before present
- ye=0 # kyr before present
- it=1 # number of iterations to decrease the misfit of preent-day bed topography
The slurm job calls the python coupler script python run_coupler.py, which prepares the input data and the calls the model steps bash run_pism_step.sh and tcsh run_vilma3d_step.cli and remaps the data from a global Gauss-Legendre grid to the projected PISM grid, using CDO and NCO tools, and vise versa.
- inp_ice6g_pism: densi.inp, load_hist.inp, tint.inp, visko.inp, wepochs.inp, see README
- data: pism_boot.nc, visc3d_new.nc
- data/ice6g: loadh.inp, ice6g_100years.nc, topo0.nc, Ice0.nc
- out:
- restart:
- pism/results: