#Prepros log date and time stamp series#

Here we computed the AAR from the “model perspective”, where the glacier is located in bins. if the glacier flowlines would have different temperature sensitivities, you need to take

For a different glacier experiment you might need Keep in mind that this is just an example for computing the AAR. Over time you can expect that this difference becomes larger. You can see in the plot that the difference between the AAR calculated from the dynamic glacier and the one of the time ) dyna_aar = tot_area_above / tot_area dyna_aar = fixed_aar dyna_aar. sum ( dim = 'dis_along_flowline' ) # Write the output dyna_aar = pd. area_m2 tot_area_above += ( area * is_above ). bed_h # In xarray things are much easier - the broadcasting is automatic is_above = surface_h > ela_ds area = is_above * ds. sum ( dim = 'dis_along_flowline' ) # Compute the surface elevation surface_h = ds. get_filepath ( 'fl_diagnostics' ), group = 'fl_' + str ( fn )) as ds : ds = ds. DataArray ( ela_df, dims = ( 'time' )) for fn in np.

# We start by reading the flowline geometry diagnostic files tot_area = 0 tot_area_above = 0 ela_ds = xr. Numerical differences in the dynamical model (these are assumed to be small as well) Melt at the tongue which might be larger than the glacier thickness is not accounted for in the fixed geometry data (these are assumed to be small) The changes in geometry during the simulation time (i.e. The difference between OGGM dynamics and fixed geometry (so small that they are barely visible for this short time period) is due to: The systematic difference between WGMS and OGGM is due to the calibration of the model, which is done with fixed geometry (the bias is zero over the entire calibration period but changes sign over the 50 year period, see plot above) OGGM, label = 'OGGM (fixed geom)' ) plt. time, smb, label = 'OGGM (dynamics)' ) # The SMB from WGMS and fixed geometry we already have plt. with changing geometry feedbacks) smb = ( ds_diag. format ( url )) 496 # if reset delete also the log directory 497 if reset : InvalidParamsError: base url seems unreachable with these parameters: url_exists ( url ): -> 493 raise InvalidParamsError ( 'base url seems unreachable with these ' 494 'parameters: '. get_prepro_base_url ( base_url = prepro_base_url, 489 border = prepro_border, 490 prepro_level = from_prepro_level, 491 rgi_version = prepro_rgi_version ) 492 if cfg. init_present_time_glacier ( gdir ) File /usr/local/pyenv/versions/3.10.7/lib/python3.10/site-packages/oggm/workflow.py:493, in init_glacier_directories (rgidf, reset, force, from_prepro_level, prepro_border, prepro_rgi_version, prepro_base_url, from_tar, delete_tar) 488 url = utils. init_glacier_directories (, from_prepro_level = 3 ) 3 gdir = gdirs 4 tasks. InvalidParamsError Traceback (most recent call last)Ĭell In, line 2 1 # (The second glacier here we will only use in the second part of this notebook.) -> 2 gdirs = workflow. Small overview of HoloViz capability of data explorationĬompare different DEMs for individual glaciers Merge, analyse and visualize OGGM GCM runs Ingest gridded products such as ice velocity into OGGM Using your our own glacier inventory with OGGMĪ quick look into the mass-balance calibration procedureĬompute smoother centerlines for shapefile output

Merging glaciers for past climate simulations The Oerlemans & Nick frontal ablation parameterization in OGGM Ice thickness inversion with frontal ablation: a case study OGGM-Shop and Glacier Directories in OGGMĮrror analysis of the global pre-processing workflowĬreate local topography maps from different DEM sources with OGGM

#Prepros log date and time stamp series#

Plotting the OGGM surface mass-balance, the ELA and AARįilter the glacier length and area time series OGGM centerlines versus elevation-band flowlines: a quick comparison Storing glacier directories for later use Getting started with OGGM: a real case study, step by step