Commit 96143c45 authored by Lennart Ramme's avatar Lennart Ramme
Browse files

[hamocc_quickplots]: qp-driver.py: add a few more HAMOCC plots

parent c987ebf2
......@@ -362,7 +362,8 @@ if do_hamocc_plots:
fig_names += ['ts_global_OMexp90', 'ts_global_calcexp90', 'ts_global_opalexp90']
fig_names += ['sec:Hamocc surface maps']
fig_names += ['srf_phyp', 'srf_zoop', 'srf_cya', 'srf_silicate', 'srf_nitrate', 'srf_phosphate']
fig_names += ['srf_alk', 'srf_dic', 'srf_pH', 'srf_co2flux'] # srf_hion
fig_names += ['srf_iron', 'srf_alk', 'srf_dic', 'srf_pH', 'srf_co2flux'] # srf_hion
fig_names += ['srf_NP_ratio', 'age_100m', 'age_500m', 'age_2000m']
fig_names += ['sec:Hamocc bias plots']
fig_names += ['po4_srf_bias_ham', 'po4_1000_bias_ham', 'no3_srf_bias_ham', 'no3_1000_bias_ham']
fig_names += ['si_srf_bias_ham', 'si_1000_bias_ham', 'o2_srf_bias_ham', 'o2_1000_bias_ham']
......@@ -3084,6 +3085,23 @@ for tave_int in tave_ints:
)
save_fig('Surface pH', path_pics, fig_name)
# ---
fig_name = 'srf_NP_ratio'
if fig_name in fig_names:
Nitrate, it_ave = pyic.time_average(IcD_ham_inv, 'no3', t1=t1, t2=t2, iz=0)
Phosphate, it_ave = pyic.time_average(IcD_ham_inv, 'po4', t1=t1, t2=t2, iz=0)
ntop= Nitrate/Phosphate
IaV = pyic.IconVariable('N/P-ratio', '', 'N/P-ratio')
IaV.data = ntop
IaV.interp_to_rectgrid(fpath_ckdtree)
pyic.hplot_base(IcD_ham_inv, IaV, clim=[0,16], cincr=1.0, cmap='RdYlBu_r',
projection=projection, xlim=[-180.,180.], ylim=[-90.,90.],
do_write_data_range=True,
title='surface N/P ratio',
save_data=save_data, fpath_nc=path_nc+fig_name+'.nc',
)
save_fig('Surface N/P ratio', path_pics, fig_name)
# ---
fig_name = 'srf_nitrate'
if fig_name in fig_names:
......@@ -3120,6 +3138,18 @@ for tave_int in tave_ints:
**Ddict_global)
save_fig('Surface silicate', path_pics, fig_name, FigInf)
# ---
fig_name = 'srf_iron'
if fig_name in fig_names:
FigInf = pyicqp.qp_hplot(fpath=path_data+fname_ham_inv, var='dfe', depth=0, it=0,
t1=t1, t2=t2,
var_fac=1e9, units='$ umol Fe m$^{-3}$',
clim=[0,6], cincr=0.3, cmap='cmo.matter',
IcD=IcD_ham_inv,
save_data=save_data, fpath_nc=path_nc+fig_name+'.nc',
**Ddict_global)
save_fig('Surface iron', path_pics, fig_name, FigInf)
# ---
fig_name = 'srf_co2flux'
if fig_name in fig_names:
......@@ -3132,6 +3162,42 @@ for tave_int in tave_ints:
**Ddict_global)
save_fig('CO2 flux', path_pics, fig_name, FigInf)
# ---
fig_name = 'age_100m'
if fig_name in fig_names:
FigInf = pyicqp.qp_hplot(fpath=path_data+fname_ham_inv, var='agesc', iz=k100, it=0,
t1=t1, t2=t2,
units='years',
clim=[0,2000], cincr=100, cmap='cubehelix_r',
IcD=IcD_ham_inv,
save_data=save_data, fpath_nc=path_nc+fig_name+'.nc',
**Ddict_global)
save_fig('Age in 100m', path_pics, fig_name, FigInf)
# ---
fig_name = 'age_500m'
if fig_name in fig_names:
FigInf = pyicqp.qp_hplot(fpath=path_data+fname_ham_inv, var='agesc', iz=k500, it=0,
t1=t1, t2=t2,
units='years',
clim=[0,2000], cincr=100, cmap='cubehelix_r',
IcD=IcD_ham_inv,
save_data=save_data, fpath_nc=path_nc+fig_name+'.nc',
**Ddict_global)
save_fig('Age in 500m', path_pics, fig_name, FigInf)
# ---
fig_name = 'age_2000m'
if fig_name in fig_names:
FigInf = pyicqp.qp_hplot(fpath=path_data+fname_ham_inv, var='agesc', iz=k2000, it=0,
t1=t1, t2=t2,
units='years',
clim=[0,2000], cincr=100, cmap='cubehelix_r',
IcD=IcD_ham_inv,
save_data=save_data, fpath_nc=path_nc+fig_name+'.nc',
**Ddict_global)
save_fig('Age in 2000m', path_pics, fig_name, FigInf)
# ---
fig_name = 'sedflux_prcaca'
if fig_name in fig_names:
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment