예제 #1
0
#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)
예제 #2
0
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()