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":