|
35 | 35 | "\n", |
36 | 36 | "# Locals\n", |
37 | 37 | "import oggm.cfg as cfg\n", |
38 | | - "from oggm import utils, workflow, tasks, DEFAULT_BASE_URL\n", |
39 | | - "from oggm.shop import gcm_climate" |
| 38 | + "from oggm import utils, workflow, DEFAULT_BASE_URL" |
| 39 | + ] |
| 40 | + }, |
| 41 | + { |
| 42 | + "cell_type": "code", |
| 43 | + "execution_count": null, |
| 44 | + "metadata": {}, |
| 45 | + "outputs": [], |
| 46 | + "source": [ |
| 47 | + "DEFAULT_BASE_URL" |
40 | 48 | ] |
41 | 49 | }, |
42 | 50 | { |
|
250 | 258 | "cell_type": "markdown", |
251 | 259 | "metadata": {}, |
252 | 260 | "source": [ |
253 | | - "This is not really visible in the plots above, but the \"old\" method of initialisation in OGGM had another issue. It assumed dynamical steady state at the begining of the simulation (the RGI date), which was required by the bed inversion process. This could lead to artifacts (mainly in the glacier length and area, as well as velocities) during the first few years of the simulation. The dynamical spinup addresses this issue by starting the simulation in 1980. \n", |
| 261 | + "This is not really visible in the plots above, but the \"old\" method of initialisation in OGGM had another issue. It assumed a dynamical steady state at the beginning of the simulation (the RGI date), which was required by the bed inversion process. This could lead to artifacts (mainly in the glacier length and area, as well as velocities) during the first few years of the simulation. The dynamical spinup addresses this issue by starting the simulation in 1980.\n", |
254 | 262 | "\n", |
255 | | - "One of the way to see the importance of the spinup is to have a look at glacier velocities. Let's plot glacier volocities along the flowline in the year 2005 (the first year we have velocities from both the dynamical spinup, and without the spinup (\"cold start\" from an equilibrium):" |
| 263 | + "One of the way to see the importance of the spinup is to have a look at glacier velocities. Let's plot glacier velocities along the flowline in the year 2005 (the first year we have velocities from both the dynamical spinup, and without the spinup (\"cold start\" from an equilibrium):" |
256 | 264 | ] |
257 | 265 | }, |
258 | 266 | { |
|
262 | 270 | "outputs": [], |
263 | 271 | "source": [ |
264 | 272 | "f = gdir.get_filepath('fl_diagnostics', filesuffix='_historical')\n", |
265 | | - "with xr.open_dataset(f, group=f'fl_0') as dg:\n", |
| 273 | + "with xr.open_dataset(f, group='fl_0') as dg:\n", |
266 | 274 | " dgno = dg.load()\n", |
267 | 275 | "f = gdir.get_filepath('fl_diagnostics', filesuffix='_spinup_historical')\n", |
268 | | - "with xr.open_dataset(f, group=f'fl_0') as dg:\n", |
| 276 | + "with xr.open_dataset(f, group='fl_0') as dg:\n", |
269 | 277 | " dgspin = dg.load()\n", |
270 | 278 | "\n", |
271 | 279 | "year = 2005\n", |
272 | | - "dgno.ice_velocity_myr.sel(time=year).plot(label='No spinup');\n", |
273 | | - "dgspin.ice_velocity_myr.sel(time=year).plot(label='With spinup');\n", |
| 280 | + "dgno.ice_velocity_myr.sel(time=year).plot(label='No spinup')\n", |
| 281 | + "dgspin.ice_velocity_myr.sel(time=year).plot(label='With spinup')\n", |
274 | 282 | "plt.title(f'Velocity along the flowline at year {year}'); plt.legend();" |
275 | 283 | ] |
276 | 284 | }, |
|
305 | 313 | "- return to the [OGGM documentation](https://docs.oggm.org)\n", |
306 | 314 | "- back to the [table of contents](../welcome.ipynb)" |
307 | 315 | ] |
| 316 | + }, |
| 317 | + { |
| 318 | + "cell_type": "code", |
| 319 | + "execution_count": null, |
| 320 | + "metadata": {}, |
| 321 | + "outputs": [], |
| 322 | + "source": [] |
308 | 323 | } |
309 | 324 | ], |
310 | 325 | "metadata": { |
|
0 commit comments