FM[0, 0, :],
FY[0, 0, :],
cfignm=cdiag + '_' + CRUN,
dt_year=ittic,
cyunit=cyu,
ctitle=CRUN + ': ' + clnm,
ymin=ym,
ymax=yp)
# Global for different depth:
vlab = ['All', '0m-100m', '100-1000m', '1000m-bottom']
bp.plot_1d_multi(XM[0, :],
XM[1:, :],
vlab,
cfignm=cdiag + '_lev_' + CRUN,
dt_year=ittic,
cyunit=cyu,
ctitle=CRUN + ': ' + clnm,
ymin=ym0,
ymax=yp0)
# Show each ocean (All depth):
bp.plot_1d_multi(XM[0, :],
FM[:, 0, :],
bo.voce2treat,
cfignm=cdiag + '_basins_' + CRUN,
dt_year=ittic,
cyunit=cyu,
ctitle=CRUN + ': ' + clnm,
ymin=ym0,
ymax=yp0)
nby = nbr/12
vice_c = nmp.zeros(nby*2) ; vice_c.shape = [2, nby]
vt_y = nmp.zeros(nby)
for jy in range(nby):
j03 = jy*12+2 ; j09 = jy*12+8
vt_y[jy] = nmp.round(vtime[jy*12],0)+0.5
vice_c[0,jy] = X1d[jv1d,j03]
vice_c[1,jy] = X1d[jv1d,j09]
ittic = bt.iaxe_tick(nby)
bp.plot_1d_multi(vt_y, vice_c, vlabels=['March', 'September'], cfignm=cfn,
dt_year=ittic, cyunit=list_units[jv], ctitle = CRUN+': Sea-Ice extent ('+cv[13:]+')',
cfig_type='svg', l_tranparent_bg=False)
elif cv == 'votemper' or cv == 'vosaline':
if not l3d:
print ' ERROR: '+sys.argv[0]+' => variable '+cv+' are here but '+cv_depth+' was not there!'
sys.exit(0)
VY, ZY = bt.monthly_2_annual(vtime, nmp.flipud(nmp.rot90(X2d[jv2d,:,:])))