def test_find_t0(self):
gdir = init_hef(border=DOM_BORDER, invert_with_sliding=False)
flowline.init_present_time_glacier(gdir)
glacier = gdir.read_pickle('model_flowlines')
df = pd.read_csv(utils.get_demo_file('hef_lengths.csv'), index_col=0)
df.columns = ['Leclercq']
df = df.loc[1950:]
vol_ref = flowline.FlowlineModel(glacier).volume_km3
init_bias = 100. # 100 so that "went too far" comes once on travis
rtol = 0.005
flowline.find_inital_glacier(gdir,
y0=df.index[0],
init_bias=init_bias,
rtol=rtol,
write_steps=False)
past_model = gdir.read_pickle('past_model')
vol_start = past_model.volume_km3
bef_fls = copy.deepcopy(past_model.fls)
mylen = []
for y in df.index:
past_model.run_until(y)
mylen.append(past_model.fls[-1].length_m)
df['oggm'] = mylen
df = df - df.iloc[-1]
vol_end = past_model.volume_km3
np.testing.assert_allclose(vol_ref, vol_end, rtol=0.05)
rmsd = utils.rmsd(df.Leclercq, df.oggm)
self.assertTrue(rmsd < 1000.)
if do_plot: # pragma: no cover
df.plot()
plt.ylabel('Glacier length (relative to 2003)')
plt.show()
fig = plt.figure()
lab = 'ref (vol={:.2f}km3)'.format(vol_ref)
plt.plot(glacier[-1].surface_h, 'k', label=lab)
lab = 'oggm start (vol={:.2f}km3)'.format(vol_start)
plt.plot(bef_fls[-1].surface_h, 'b', label=lab)
lab = 'oggm end (vol={:.2f}km3)'.format(vol_end)
plt.plot(past_model.fls[-1].surface_h, 'r', label=lab)
plt.plot(glacier[-1].bed_h, 'gray', linewidth=2)
plt.legend(loc='best')
plt.show()
def test_grow(self):
gdirs = up_to_inversion()
for gd in gdirs[0:2]:
if gd.rgi_id in ['RGI40-11.00719']:
# Bad bad glacier
continue
flowline.init_present_time_glacier(gd)
flowline.find_inital_glacier(gd)
def test_grow(self):
gdirs = up_to_inversion()
for gd in gdirs[0:2]:
if gd.rgi_id in ['RGI40-11.00719']:
# Bad bad glacier
continue
flowline.init_present_time_glacier(gd)
flowline.find_inital_glacier(gd)
def test_find_t0(self):
gdir = init_hef(border=DOM_BORDER, invert_with_sliding=False)
flowline.init_present_time_glacier(gdir)
glacier = gdir.read_pickle('model_flowlines')
df = pd.read_csv(utils.get_demo_file('hef_lengths.csv'), index_col=0)
df.columns = ['Leclercq']
df = df.loc[1950:]
vol_ref = flowline.FlowlineModel(glacier).volume_km3
init_bias = 100. # 100 so that "went too far" comes once on travis
rtol = 0.005
flowline.find_inital_glacier(gdir, y0=df.index[0], init_bias=init_bias,
rtol=rtol, write_steps=False)
past_model = gdir.read_pickle('past_model')
vol_start = past_model.volume_km3
bef_fls = copy.deepcopy(past_model.fls)
mylen = []
for y in df.index:
past_model.run_until(y)
mylen.append(past_model.fls[-1].length_m)
df['oggm'] = mylen
df = df-df.iloc[-1]
vol_end = past_model.volume_km3
np.testing.assert_allclose(vol_ref, vol_end, rtol=0.05)
rmsd = utils.rmsd(df.Leclercq, df.oggm)
self.assertTrue(rmsd < 1000.)
if do_plot: # pragma: no cover
df.plot()
plt.ylabel('Glacier length (relative to 2003)')
plt.show()
fig = plt.figure()
lab = 'ref (vol={:.2f}km3)'.format(vol_ref)
plt.plot(glacier[-1].surface_h, 'k', label=lab)
lab = 'oggm start (vol={:.2f}km3)'.format(vol_start)
plt.plot(bef_fls[-1].surface_h, 'b', label=lab)
lab = 'oggm end (vol={:.2f}km3)'.format(vol_end)
plt.plot(past_model.fls[-1].surface_h, 'r', label=lab)
plt.plot(glacier[-1].bed_h, 'gray', linewidth=2)
plt.legend(loc='best')
plt.show()