pism/pism — Issue #23: Initialization of the Lingle-Clark bed deformation model is broken
Author: ckhroulev
Created: 2011-11-23 17:08:08+00:00
State: open
Labels: bug, urgent

The Lingle-Clark (`-bed_def lc`) initialization-from-the-uplift code seems to be broken. This makes `-bed_def lc` runs not re-startable and results in rather nasty artifacts in `wvelbase` and `wvelsurf`.

Run the following script, then compare `ex-1.nc` and `ex-2.nc`. Look for jumps in `topg` and `dbdt` at around year 100.

```bash
#!/bin/bash -x

pisms -y 1000 -Mx 31 -My 31 -o in.nc

EX1_FILE="ex-1.nc"
EX2_FILE="ex-2.nc"
OPTS="-extra_vars topg,dbdt -extra_times 0:1:200 -bed_def lc"

rm -f $EX1_FILE $EX2_FILE ex-*.nc out-*.nc

# Run for 200 years straight
pismr -i in.nc $OPTS -extra_file $EX1_FILE -ys 0 -ye 200 -o out-200.nc

# Run for 100 years, then for 100 more years:
pismr -i in.nc $OPTS -extra_file $EX2_FILE -ys 0 -ye 100 -o out-100.nc

pismr -i out-100.nc $OPTS -extra_file $EX2_FILE -extra_append -ye 200 -o out-200-2.nc
```