oggm.tasks.process_gcm_data#

oggm.tasks.process_gcm_data(gdir, filesuffix='', prcp=None, temp=None, year_range=('1961', '1990'), scale_stddev=True, time_unit=None, calendar=None, source='', apply_bias_correction=True)[source]#

Applies the anomaly method to GCM climate data

This function can be applied to any GCM data, if it is provided in a suitable xarray.DataArray. See Parameter description for format details.

For CESM-LME a specific function tasks.process_cesm_data() is available which does the preprocessing of the data and subsequently calls this function.

Parameters
gdiroggm.GlacierDirectory

where to write the data

filesuffixstr

append a suffix to the filename (useful for ensemble experiments).

prcpxarray.DataArray
monthly total precipitation [mm month-1]
Coordinates:
lat float64
lon float64
time: cftime object
tempxarray.DataArray
monthly temperature [K]
Coordinates:
lat float64
lon float64
time cftime object
year_rangetuple of str

the year range for which you want to compute the anomalies. Default is (‘1961’, ‘1990’)

scale_stddevbool

whether or not to scale the temperature standard deviation as well

time_unitstr

The unit conversion for NetCDF files. It must be adapted to the length of the time series. The default is to choose it ourselves based on the starting year. For example: ‘days since 0850-01-01 00:00:00’

calendarstr

If you use an exotic calendar (e.g. ‘noleap’)

sourcestr

For metadata: the source of the climate data

apply_bias_correctionboolean

if a bias-correction should be applied. Default is True, only set it to False if the GCM has already been externally bias-corrected to the applied observational calibration dataset (true for ISIMIP 3b that is bias-corrected to W5E5). !!! We assume that temp is in Kelvin and convert to CELSIUS !!!

Notes

Files written to the glacier directory:
gcm_data.nc

The monthly GCM climate timeseries stored in a netCDF file.