def ret_regionmask(region,year,mon):
if region == "ASAS":
xydatadir = "/media/disk2/out/chart/ASAS/const"
paradate = datetime.date(year,mon,1)
if (paradate < datetime.date(2006,1,1)):
name_domain_mask = xydatadir + "/domainmask_saone.%s.2000.01.bn"%(region)
if ( datetime.date(2006,1,1)<=paradate<datetime.date(2006,3,1)):
name_domain_mask = xydatadir + "/domainmask_saone.%s.2006.01.bn"%(region)
if ( datetime.date(2006,3,1)<=paradate):
name_domain_mask = xydatadir + "/domainmask_saone.%s.2006.03.bn"%(region)
#
a2regionmask = fromfile(name_domain_mask, float32).reshape(ny,nx)
a2regionmask = ma.masked_equal(a2regionmask, 0.0).filled(miss)
else:
lllon, lllat, urlon, urlat = chart_para.ret_domain_corner_rect_forfig(region)
a2regionmask = ctrack_func.mk_region_mask(lllat,urlat,lllon,urlon, nx,ny, lat_first, lon_first, dlat, dlon)
a2regionmask = ma.masked_equal(a2regionmask, 0.0).filled(miss)
#--
return a2regionmask
#dist_f = 400 # [km]
nx,ny =[360,180]
thdura_c = 48
thdura_tc = thdura_c
thorog = 1500 # [m]
miss = -9999.0
miss_out = -9999.0
miss_gpcp = -99999.
#---------------------
region = "GLOB"
lllon, lllat, urlon, urlat = chart_para.ret_domain_corner_rect_forfig(region)
#---------------------
for prtype in lprtype:
for ldist in lldist:
dist_tc, dist_c, dist_f = ldist
for bstflag in lbstflag:
#-- TC type ----------
if bstflag ==True:
tctype = "bst"
else:
tctype = ""
#--- tag dir_root -------------------
tagdir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/tag/c%02dh.tc%02dh.%stc%04d.c%04d.f%04d"%(sresol, thdura_c, thdura_tc, tctype, dist_tc, dist_c, dist_f)
#--- precipitation directory & timestep-----
if prtype == "GSMaP":
