thpgrad = dpgradrange[0][0]
#----------------------------
psldir_root = "/media/disk2/data/JRA25/sa.one.%s/6hr/PRMSL"%(model)
pgraddir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/pgrad"%(model)
lifedir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/life"%(model)
nextposdir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/nextpos"%(model)
#************************************
dtrack = {}
for iclass in lclass:
dtrack[iclass] = []
#------------------------------------
for year in range(iyear, eyear+1):
#************************
# load TCs
#------------------------
dTC = tc_func.ret_ibtracs_dpyxy_saone(year)
#------------------------
for mon in lmon:
eday = calendar.monthrange(year,mon)[1]
###############
## no leap
###############
#if (mon==2)&(calendar.isleap(year)):
# eday = calendar.monthrange(year,mon)[1] -1
#else:
# eday = calendar.monthrange(year,mon)[1]
#if singleday == True:
# eday = iday
##############
for iclass in lclass:
thgrad_min = dpgradrange[iclass][0]
thgrad_max = dpgradrange[iclass][1]
#---------------------
# for only figure, comment out from here
#---------------------
i = 0
for year in range(iyear, eyear+1):
#***********************
#-- TC best track
#-----------------------
if bsttc == True:
ddatbst = tc_func.ret_ibtracs_dpyxy_saone(year)
#***********************
for mon in lmon:
eday = calendar.monthrange(year, mon)[1]
if singleday == True:
eday = iday
#---- dummy ------------------------------
a2one = ones([ny_eqgrid, nx_eqgrid], float32)
a2sum_eqgrid = zeros([ny_eqgrid, nx_eqgrid], float32)
a2num_eqgrid = zeros([ny_eqgrid, nx_eqgrid], float32)
#---- dummy for cor. coef ----------------
a2cor_num = zeros([ny_eqgrid, nx_eqgrid], float32)
a2cor_SA = zeros([ny_eqgrid, nx_eqgrid], float32)
a2cor_SB = zeros([ny_eqgrid, nx_eqgrid], float32)
a2cor_SAB = zeros([ny_eqgrid, nx_eqgrid], float32)