Repository for the linear response function (LRF) calculations for BISICLES and analysis, based on Levermann et al. 2020. This project takes a directory of BISICLES plot files for a LRF simulation run. It then runs the stats file tool to create a csv of summary statistics and subsequently calculate the LRF and rolling mean for that simulation.
- You need to have the BISICLES filetools compiled on your system, instructions on how to do this can be found here
- Clone the repository to your desired location
- Create and activate a new virtual environment
- Install the required Python environment using
pip install -r requirements.txt
To calculate the LRFs you run the calculate_lrf.py file. Change the paths at the top to correspond to yours and change the EXPERIMENT variable to the amount of addition basal melt for that simulation. The inputs required are:
PATHto input data directory with plot filesMASK_PATHto directory with Antarctica region masksCSV_PATHoutput path for csv filesSTATS_TOOLpath to compiled stats filetool driverEXPERIMENTnumber by which basal melt was increased for the run
If you do not have Antarctica region masks, you can create some from the BISICLES/applications/bedmachine/maskeMasksfromVel directory.
Additionally, if you want to create a time series with a rolling 10-year mean, you can run create_timeseries.py, here only the path to the csv directory and an plot output directory are needed.