#urlat = 80.0 #urlon = 210.0 #----- miss = -9999.0 miss_int = -9999 chartdir_root = "/media/disk2/out/chart/ASAS/front" tdir_root = "/media/disk2/data/JRA25/sa.one.%s/6hr/TMP"%(sresol) qdir_root = "/media/disk2/data/JRA25/sa.one.%s/6hr/SPFH"%(sresol) #------------------------ # for cyclone center #---------- pgraddir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/pgrad"%(sresol) lifedir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/life"%(sresol) thdura_c = 48 pgradmin = ctrack_para.ret_dpgradrange()[2][0] # Pa/1000km #************************ lat_first = -89.5 dlat = 1.0 dlon = 1.0 #******************************************** # Function #----------- def load_a2var(year,mon,day,hour, plev, var): tdir_root = "/media/disk2/data/JRA25/sa.one.anl_p/6hr/TMP" qdir_root = "/media/disk2/data/JRA25/sa.one.anl_p/6hr/SPFH" wdir_root = "//media/disk2/data/JRA25/sa.one.anl_chipsi/6hr/VVEL" tdir = tdir_root + "/%04d%02d"%(year,mon) qdir = qdir_root + "/%04d%02d"%(year,mon)
mnum_min = 1.0
#***************************************
diyear = {"his": iyear_his, "fut": iyear_fut}
deyear = {"his": eyear_his, "fut": eyear_fut}
#***************************************
(imon, emon) = ctrack_para.ret_im_em(season)
#***************************************
dexpr = {"his": "historical", "fut": "rcp85", "dif": "dif"}
#***************************************
lera = ["fut", "his"]
lvar = ["num", "sp", "mnum"]
ldifvar = ["drnum", "drtotnum", "dmp"]
#***************************************
# class
#-----------------------------
dpgradrange = ctrack_para.ret_dpgradrange()
cmin = dpgradrange[0][0]
lclass = dpgradrange.keys()
nclass = len(lclass) -1
#***************************************
# wbin
#-----------------------------
dlwbin = ctrack_para.ret_dlwbin()
liw = dlwbin.keys()
nwbin = len(liw)
#***************************************
oekakidir = "/home/utsumi/bin/dtanl/ctrack/oekaki"
#----------------------
# colormap
iyear = 2001 eyear = 2009 lmon = range(1,12+1) tstp = "6hr" nx_org = 360 ny_org = 180 sresol = "anl_p" var = "pr" thdura = 48 iday = 1 lhour = [0, 6, 12, 18] #region = "PNW" region = "GLB" #thgrad_min = 500.0 # Pa/1000km #thgrad_max = 1000.0 # Pa/1000km dpgradrange = ctrack_para.ret_dpgradrange() lclass = dpgradrange.keys() #vtype = "GSMaP" # "JRA", "GPCP1DD", "GSMaP" #vtype = "JRA" # "JRA", "GPCP1DD", "GSMaP" #vtype = "GPCP1DD" # "JRA", "GPCP1DD", "GSMaP" #lvtype = ["GSMaP", "JRA", "GPCP1DD"] lvtype = ["GSMaP"] #lvtype = ["GPCP1DD"] #lvtype = ["APHRO_MA"] #--------------------- dkm = 100.0 # equal area grid resolution [km] nradeqgrid = 30 #-- TC para ----- thpgrad = ctrack_para.ret_dpgradrange(sresol)[2][0] thsst = tc_para.ret_thsst()