#
stitle = "dstat.adv.late (K/hPa) lenday=%02d\n"%(len(ldday))
stitle = stitle + "%04d-%02d-%02d"%(year,mon,day)
#---
if (day == iday):
cbarname = figdir + "/dsta.cbar.png"
else:
cbarname = False
#---
bnd = [-0.05, -0.03, -0.01, 0.01, 0.03, 0.05]
#--
if len(ldday)==1:
print "AAAA"
bnd = [-0.09,-0.07,-0.05, -0.03, -0.01, 0.01, 0.03, 0.05, 0.07,0.09]
#--
ctrack_fig.mk_pict_saone_reg_symm(a2dsta_adv_lateral, stitle=stitle, soname=figname, bnd=bnd, mycm="RdBu",cbarname="./dstat.late.cbar.png")
print figname
##--- figure sta ----
#bnd = [-0.025, -0.015, -0.005, 0.005, 0.015, 0.025]
##bnd = [-0.05,-0.03,-0.01,0.0, 0.01, 0.03, 0.05]
#ctrack_fig.mk_pict_saone_reg_symm(a2sta, soname="./sta.png", bnd=bnd, mycm="RdBu",cbarname="./sta.cbar.png")
##--- figure dsta vertical ----
#bnd = [-0.05, -0.03, -0.01, 0.01, 0.03, 0.05]
#ctrack_fig.mk_pict_saone_reg_symm(a2dsta_adv_vertical, soname="./dsta.vert.png", bnd=bnd, mycm="RdBu",cbarname="./dstat.vert.cbar.png")
##--- figure dsta ----
#bnd = [-0.05, -0.03, -0.01, 0.01, 0.03, 0.05]
#ctrack_fig.mk_pict_saone_reg_symm(a2dsta_adv, soname="./dsta.png", bnd=bnd, mycm="RdBu",cbarname="./dstat.cbar.png")
##--- figure adv850 ----
#bnd = [-20,-15,-10,-5,0,5,10,15,20]
sodir = idir_fut soname = sodir + "/dif.freq.t.front.%s.%s.rad%04dkm.M1_%4.2f_M2_%3.1f.%04d-%04d.%s.sa.one"%(model,ens,countrad,thfmask1t,thfmask2t,iyear_fut,eyear_fut,season) a2freq_dif.tofile(soname) #******************************* # Figure #----------- bnd = [-21,-15,-9,-3,3,9,15,21] #bnd = [-1.0, -0.5, -0.1, 0.1, 0.5, 1.0] #bnd = [0.01,0.25,0.5,1.0,2.0,4.0,8.0] #bnd = [0.25,0.5,1.0,2.0,4.0,8.0] #bnd = [5,10,15,20,25,30,35,40,45] #bnd = [10,20,30,40,50,60,70,80] figdir = sodir figname = soname[:-7] + ".png" cbarname = soname[:-7] + ".cbar.png" #---------- stitle = "freq. %s.ExC season:%s %04d-%04d %s"%(model,season,iyear_fut, eyear_fut,ens) #mycm = "Spectral" mycm = "RdBu_r" datname = soname a2figdat = fromfile(datname, float32).reshape(ny,nx) #------------------------------- a2figdat = ma.masked_equal(a2figdat, miss).filled(0.0) * 100.0 ctrack_fig.mk_pict_saone_reg_symm(a2figdat, lllat=lllat, lllon=lllon, urlat=urlat, urlon=urlon, bnd=bnd, mycm=mycm, soname=figname, stitle=stitle, miss=miss, a2shade=a2shade, cbarname=cbarname) print figname
##--- figure gradqmean ------
#bnd = [0.5, 1.0, 1.5, 2.0, 2.5, 3.0]
##bnd = False
#ctrack_fig.mk_pict_saone_reg(a2gradqmean*1.e+8, soname="./grad.qmean.png", bnd=bnd, mycm="jet",cbarname="./grad.q.cbar.png")
##--- figure gradtmean ------
##bnd = [0.0, 0.2, 0.4, 0.6, 0.8, 1.0]
#bnd = [0.5, 1.0, 1.5, 2.0, 2.5, 3.0]
#ctrack_fig.mk_pict_saone_reg(a2gradtmean*1.e+5, soname="./grad.tmean.png", bnd=bnd, mycm="jet",cbarname="./grad.t.cbar.png")
#--- figure dsta lateral----
bnd = [-0.05, -0.03, -0.01, 0.01, 0.03, 0.05]
if len(ldday)==1:
print "AAAA"
bnd = [-0.09,-0.07,-0.05, -0.03, -0.01, 0.01, 0.03, 0.05, 0.07,0.09]
ctrack_fig.mk_pict_saone_reg_symm(a2dsta_adv_lateral, soname="./dsta.late.png", bnd=bnd, mycm="RdBu",cbarname="./dstat.late.cbar.png")
#--- figure sta ----
bnd = [-0.025, -0.015, -0.005, 0.005, 0.015, 0.025]
#bnd = [-0.05,-0.03,-0.01,0.0, 0.01, 0.03, 0.05]
ctrack_fig.mk_pict_saone_reg_symm(a2sta, soname="./sta.png", bnd=bnd, mycm="RdBu",cbarname="./sta.cbar.png")
##--- figure dsta vertical ----
#bnd = [-0.05, -0.03, -0.01, 0.01, 0.03, 0.05]
#ctrack_fig.mk_pict_saone_reg_symm(a2dsta_adv_vertical, soname="./dsta.vert.png", bnd=bnd, mycm="RdBu",cbarname="./dstat.vert.cbar.png")
##--- figure dsta ----
#bnd = [-0.05, -0.03, -0.01, 0.01, 0.03, 0.05]
#ctrack_fig.mk_pict_saone_reg_symm(a2dsta_adv, soname="./dsta.png", bnd=bnd, mycm="RdBu",cbarname="./dstat.cbar.png")
#--- figure adv850 ----
bnd = [-20,-15,-10,-5,0,5,10,15,20]
