trun.sh - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
(HTM) git clone git://src.adamsgaard.dk/pism
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) LICENSE
---
trun.sh (901B)
---
1 #!/bin/bash
2
3 xx=151
4 yy=$xx
5 length=400
6
7 input="-i input.nc -bootstrap"
8
9 grid="-Mx $xx -My $yy -Mz 11 -Mbz 1 -Lz 1500 -Lbz 0 -y $length"
10
11 physics="-stress_balance ssa+sia -ssa_dirichlet_bc -cfbc -part_grid"
12
13 extra="-extra_vars thk,mask,velbar_mag,Href,velbar"
14
15 extra="$extra -extra_times 5 -extra_file ex.nc"
16 ts="-ts_file ts.nc -ts_times 1"
17
18 output="-o o.nc -o_order zyx $extra $ts"
19
20 ocean="-ssa_method fd -ocean constant,delta_MBP -ocean_delta_MBP_file delta_MBP.nc"
21 mpiexec -n 4 pismr $input $grid $physics $ocean $output
22
23 # Cut out a slice with x == 0 (through the center of the domain).
24 ncks -O -d x,$(( $xx / 2 )) ex.nc center-mbp.nc
25
26 ocean="-ssa_method fd -ocean constant"
27 mpiexec -n 4 pismr $input $grid $physics $ocean $output
28
29 # Cut out a slice with x == 0 (through the center of the domain).
30 ncks -O -d x,$(( $xx / 2 )) ex.nc center-no-mbp.nc
31
32 ncdiff -O center-mbp.nc center-no-mbp.nc diff.nc