# 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]