PSM for lake varve thickness

In this tutorial, we introduce the PSM for lake varve thickness in cfr.

%load_ext autoreload
%autoreload 2

import cfr
print(cfr.__version__)

Data preparation

Proxy

pdb = cfr.ProxyDatabase().fetch('PAGES2kv2')

pobj = pdb['Arc_001']
fig, ax = pobj.plot()

Model

model_tas = cfr.ClimateField().fetch('iCESM_past1000historical/tas')

>>> The target file seems existed at: ./data/tas_sfc_Amon_iCESM_past1000historical_085001-200512.nc . Loading from it instead of downloading ...

Get climate data for a specific ProxyRecord

%%time

pobj.del_clim()
pobj.get_clim(model_tas, tag='model')

CPU times: user 2.59 ms, sys: 39.6 ms, total: 42.1 ms
Wall time: 41.2 ms

Create a PSM object

mdl = cfr.psm.Lake_VarveThickness(pobj)

%%time
pp = mdl.forward()

CPU times: user 127 ms, sys: 36.2 ms, total: 164 ms
Wall time: 165 ms

fig, ax = pp.plot()