datdir = odir
datname = oname
if region == "GLOB":
bnd = [-40,-20,-10,-5,-2,2,5,10,20,40]
figdir = datdir
dotstep = 5
if region == "JPN":
bnd = [-40,-20,-10,-5,-2,2,5,10,20,40]
figdir = datdir + ".JPN"
dotstep = 1
ctrack_func.mk_dir(figdir)
#----
figname = figdir + "/" + datname.split("/")[-1][:-7] + ".png"
cbarname = figdir + "/dpr.dec.cbar.png"
mons = len(ctrack_para.ret_lmon(season))
a2figdat = fromfile(datname,float32).reshape(ny,nx) *60*60*24.*totaltimes /mons # (mm/month)
##-------------------------------
##--- test ---
#a2shade = ones([ny,nx],float32)
##------------
ctrack_fig.mk_pict_saone_reg_symm_dotshade(a2figdat, lllat=lllat, lllon=lllon, urlat=urlat, urlon=urlon, bnd=bnd, mycm=mycm, soname=figname, stitle=stitle, miss=miss, a2shade=a2shade, cbarname=cbarname, dotcolor="k",dotstep=dotstep )
#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
# at the center of original grid box lllat = 25. urlat = 60. lllon = 110. urlon = 180. #lllat = -89.5 #urlat = 89.5 #lllon = 0.5 #urlon = 359.5 thdura = 72 #---------------------------- lmon = ctrack_para.ret_lmon(season) #lmon = [1] #---------------------------- dlat = 1.0 dlon = 1.0 a1lat = arange(-89.5, 89.5 +dlat*0.5, dlat) a1lon = arange(0.5, 359.5 +dlon*0.5, dlon) #---------------------------- #---------------------------- psldir_root = "/media/disk2/out/chart/ASAS/exc" lifedir_root = "/media/disk2/out/chart/ASAS/exc/life" nextposdir_root = "/media/disk2/out/chart/ASAS/exc/nextpos" #************************************ ltrack = []
if region == "GLOB":
bnd = [-40,-20,-10,-5,-2,2,5,10,20,40]
bnd = [-40,-20,-10,-5,-2,2,5,10,20,40]
figdir = datdir
dotstep = 5
if region == "JPN":
bnd = [-40,-20,-10,-5,-2,2,5,10,20,40]
figdir = datdir + ".JPN"
dotstep = 1
ctrack_func.mk_dir(figdir)
#----
figname = figdir + "/" + datname.split("/")[-1][:-7] + ".png"
cbarname = figdir + "/dpr.dec.cbar.png"
mons = len(lyear) * len(ctrack_para.ret_lmon(season))
a2figdat = fromfile(datname,float32).reshape(ny,nx) *60*60*24.*totaldays /mons # (mm/month)
##-------------------------------
##--- test ---
#a2shade = ones([ny,nx],float32)
##------------
ctrack_fig.mk_pict_saone_reg_symm_dotshade(a2figdat, lllat=lllat, lllon=lllon, urlat=urlat, urlon=urlon, bnd=bnd, mycm=mycm, soname=figname, stitle=stitle, miss=miss, a2shade=a2shade, cbarname=cbarname, dotcolor="k",dotstep=dotstep )
#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