The python script in this folder will read the csv-file containing concetration of black carbon as measured at the Upper Freeman Glacier. It will also read a selection of model files and calculate the yearly concentration of BC from the model CanESM5.
import pandas as pd
import numpy as np
import sys,os,glob
import xarray as xr
It is expected that you run the code distributed in this project from the project_root folder (where you find this file). To have an exact copy of the relevant environment you need to install poetry. You can then run the visualization.ipynb or examplescript.py:
$[~/PROJECT_ROOT] poetry shell
$[~/PROJECT_ROOT] jupyter notebook
$[~/PROJECT_ROOT] python src/example.py "examplescript.py" contains functions which is called in the end of the script
ice_conc, iceyear, icelat, icelon = icecore_info()
model_conc = make_concentration(icelat,icelon)
print(ice_conc,model_conc)The concentrations are currently just printed, but plotting them together is encouraged for comparing model concentration of black carbon to what is measured in an ice core.
Pull requests are welcome. For major changes, please open an issue first to discuss what you would like to change.
Please make sure to update tests as appropriate.