historical_delta_method(gdir, ref_filesuffix='', hist_filesuffix='', output_filesuffix='', ref_year_range=None, delete_input_files=True, scale_stddev=True, replace_with_ref_data=True)¶
Applies the anomaly method to historical climate data.
This function can be used to prolongate historical time series, for example by bias-correcting CERA-20C to ERA5 or ERA5-Land.
The timeseries must be already available in the glacier directory
- gdir :
where to write the data
- ref_filesuffix : str
the filesuffix of the historical climate data to take as reference
- hist_filesuffix : str
the filesuffix of the historical climate data to apply to the reference
- output_filesuffix : str
the filesuffix of the output file (usually left empty - i.e. this file will become the default)
- ref_year_range : tuple of str
the year range for which you want to compute the anomalies. The default is to take the entire reference data period, but you could also choose (‘1961’, ‘1990’) for example
- delete_input_files : bool
delete the input files after use - useful for operational runs where you don’t want to carry too many files
- scale_stddev : bool
whether or not to scale the temperature standard deviation as well (you probably want to do that)
- replace_with_ref_data : bool
the default is to paste the bias-corrected data where no reference data is available, i.e. creating timeseries which are not consistent in time but “better” for recent times (e.g. CERA-20C until 1980, then ERA5). Set this to False to present this and make a consistent time series of CERA-20C (but bias corrected to the reference data, so “better” than CERA-20C out of the box).
- Files writen to the glacier directory:
- The historical monthly climate timeseries stored in a netCDF file.
- gdir :