PSM for coral Sr/Ca

In this tutorial, we introduce the PSM for coral Sr/Ca in cfr.

%load_ext autoreload
%autoreload 2

import cfr
print(cfr.__version__)

Data preparation

Proxy

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

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

Model

# note here the nc file has been interpolated to the atmos grid

%%time
dirpath = './data/icesm_last_millennium_historical'
model_tos = cfr.ClimateField().load_nc(os.path.join(dirpath, 'tos_sfc_Omon_iCESM_past1000historical_085001-200512_interp.nc'))

CPU times: user 980 ms, sys: 7.35 ms, total: 987 ms
Wall time: 1.05 s

Get climate data for a specific ProxyRecord

%%time

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

CPU times: user 208 ms, sys: 3.1 s, total: 3.31 s
Wall time: 1min 24s

Create a PSM object

mdl = cfr.psm.Coral_SrCa(pobj)

%%time
pp = mdl.forward()

CPU times: user 498 µs, sys: 537 µs, total: 1.03 ms
Wall time: 1.35 ms

fig, ax = pp.plot()