Commit 6404928f authored by Lennart Ramme's avatar Lennart Ramme
Browse files

[hamocc_quickplots]: qp_driver.py: improvements for HAMOCC quickplots

parent 2fc4b5e2
...@@ -358,7 +358,7 @@ if do_hamocc_plots: ...@@ -358,7 +358,7 @@ if do_hamocc_plots:
fig_names += ['sec:Hamocc time series'] fig_names += ['sec:Hamocc time series']
fig_names += ['ts_global_npp', 'ts_global_nppcya', 'ts_global_zoograzing', 'ts_global_netco2flux'] fig_names += ['ts_global_npp', 'ts_global_nppcya', 'ts_global_zoograzing', 'ts_global_netco2flux']
fig_names += ['ts_global_surface_alk', 'ts_global_surface_dic', 'ts_global_surface_sil', 'ts_global_surface_phos'] fig_names += ['ts_global_surface_alk', 'ts_global_surface_dic', 'ts_global_surface_sil', 'ts_global_surface_phos']
fig_names += ['ts_WC_denit', 'ts_sed_denit', 'ts_n2fix', 'ts_global_opal_prod', 'ts_global_caco3_prod'] fig_names += ['ts_global_surface_nitrate', 'ts_WC_denit', 'ts_sed_denit', 'ts_n2fix', 'ts_global_opal_prod', 'ts_global_caco3_prod']
fig_names += ['ts_global_OMexp90', 'ts_global_calcexp90', 'ts_global_opalexp90'] fig_names += ['ts_global_OMexp90', 'ts_global_calcexp90', 'ts_global_opalexp90']
fig_names += ['sec:Hamocc surface maps'] fig_names += ['sec:Hamocc surface maps']
fig_names += ['srf_phyp', 'srf_zoop', 'srf_cya', 'srf_silicate', 'srf_nitrate', 'srf_phosphate'] fig_names += ['srf_phyp', 'srf_zoop', 'srf_cya', 'srf_silicate', 'srf_nitrate', 'srf_phosphate']
...@@ -2869,6 +2869,10 @@ for tave_int in tave_ints: ...@@ -2869,6 +2869,10 @@ for tave_int in tave_ints:
if fig_name in fig_names: if fig_name in fig_names:
FigInf, Dhandles = pyicqp.qp_timeseries(IcD_ham_mon, fname_ham_mon, ['global_surface_silicate'], t1=t1, t2=t2, ave_freq=ave_freq, var_fac=1e6, units='mmol Si m$^{-3}$', omit_last_file=omit_last_file, save_data=save_data, fpath_nc=path_nc+fig_name+'.nc') FigInf, Dhandles = pyicqp.qp_timeseries(IcD_ham_mon, fname_ham_mon, ['global_surface_silicate'], t1=t1, t2=t2, ave_freq=ave_freq, var_fac=1e6, units='mmol Si m$^{-3}$', omit_last_file=omit_last_file, save_data=save_data, fpath_nc=path_nc+fig_name+'.nc')
save_fig(fig_name, path_pics, fig_name) save_fig(fig_name, path_pics, fig_name)
fig_name = 'ts_global_surface_nitrate'
if fig_name in fig_names:
FigInf, Dhandles = pyicqp.qp_timeseries(IcD_ham_mon, fname_ham_mon, ['global_surface_nitrate'], t1=t1, t2=t2, ave_freq=ave_freq, var_fac=1e6, units='mmol N m$^{-3}$', omit_last_file=omit_last_file, save_data=save_data, fpath_nc=path_nc+fig_name+'.nc')
save_fig(fig_name, path_pics, fig_name)
fig_name = 'ts_global_opal_prod' fig_name = 'ts_global_opal_prod'
if fig_name in fig_names: if fig_name in fig_names:
FigInf, Dhandles = pyicqp.qp_timeseries(IcD_ham_mon, fname_ham_mon, ['global_opal_production'], t1=t1, t2=t2, ave_freq=ave_freq, omit_last_file=omit_last_file, save_data=save_data, fpath_nc=path_nc+fig_name+'.nc') FigInf, Dhandles = pyicqp.qp_timeseries(IcD_ham_mon, fname_ham_mon, ['global_opal_production'], t1=t1, t2=t2, ave_freq=ave_freq, omit_last_file=omit_last_file, save_data=save_data, fpath_nc=path_nc+fig_name+'.nc')
...@@ -2887,7 +2891,7 @@ for tave_int in tave_ints: ...@@ -2887,7 +2891,7 @@ for tave_int in tave_ints:
save_fig(fig_name, path_pics, fig_name) save_fig(fig_name, path_pics, fig_name)
fig_name = 'ts_global_opalexp90' fig_name = 'ts_global_opalexp90'
if fig_name in fig_names: if fig_name in fig_names:
FigInf, Dhandles = pyicqp.qp_timeseries(IcD_ham_mon, fname_ham_mon, ['global_opal_export_at_90m'], t1=t1, t2=t2, ave_freq=ave_freq, var_fac=365.25*24*60*60*1e-9, units='Tmol Si yr$^{-1}$', omit_last_file=omit_last_file, save_data=save_data, fpath_nc=path_nc+fig_name+'.nc') FigInf, Dhandles = pyicqp.qp_timeseries(IcD_ham_mon, fname_ham_mon, ['global_opal_export_at_90m'], t1=t1, t2=t2, ave_freq=ave_freq, var_fac=365.25*24*60*60*1e-9, units=' [Tmol Si yr$^{-1}$]', omit_last_file=omit_last_file, save_data=save_data, fpath_nc=path_nc+fig_name+'.nc')
save_fig(fig_name, path_pics, fig_name) save_fig(fig_name, path_pics, fig_name)
# ------------------------------------------------------------------------------ # ------------------------------------------------------------------------------
...@@ -3025,9 +3029,9 @@ for tave_int in tave_ints: ...@@ -3025,9 +3029,9 @@ for tave_int in tave_ints:
fig_name = 'srf_cya' fig_name = 'srf_cya'
if fig_name in fig_names: if fig_name in fig_names:
FigInf = pyicqp.qp_hplot(fpath=path_data+fname_ham_inv, var='phydiaz', depth=0, it=0, FigInf = pyicqp.qp_hplot(fpath=path_data+fname_ham_inv, var='phydiaz', depth=0, it=0,
t1=t1, t2=t2, t1=t1, t2=t2, logplot=True,
clim=[0,0.5], cincr=0.025, cmap='cmo.amp', clim=[-5.0,0.0], cincr=0.5, cmap='cmo.amp',
var_fac=1e6, units='mmol P m$^{-3}$', var_fac=1e6, units='log$_{10}$ mmol P m$^{-3}$',
IcD=IcD_ham_inv, IcD=IcD_ham_inv,
save_data=save_data, fpath_nc=path_nc+fig_name+'.nc', save_data=save_data, fpath_nc=path_nc+fig_name+'.nc',
**Ddict_global) **Ddict_global)
...@@ -3143,7 +3147,7 @@ for tave_int in tave_ints: ...@@ -3143,7 +3147,7 @@ for tave_int in tave_ints:
if fig_name in fig_names: if fig_name in fig_names:
FigInf = pyicqp.qp_hplot(fpath=path_data+fname_ham_inv, var='dfe', depth=0, it=0, FigInf = pyicqp.qp_hplot(fpath=path_data+fname_ham_inv, var='dfe', depth=0, it=0,
t1=t1, t2=t2, t1=t1, t2=t2,
var_fac=1e9, units='$ umol Fe m$^{-3}$', var_fac=1e9, units='umol Fe m$^{-3}$',
clim=[0,6], cincr=0.3, cmap='cmo.matter', clim=[0,6], cincr=0.3, cmap='cmo.matter',
IcD=IcD_ham_inv, IcD=IcD_ham_inv,
save_data=save_data, fpath_nc=path_nc+fig_name+'.nc', save_data=save_data, fpath_nc=path_nc+fig_name+'.nc',
......
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