bestdir = "/media/disk2/data/ibtracs/v03r04"
odir = "/media/disk2/out/obj.valid/exc.pgrad"
#------------------
def positive_lon(lon):
if lon <0.0:
positive_lon = 360 + lon
else:
positive_lon = lon
#
return positive_lon
#------------------
for model in lmodel:
for region in lregion:
#-- region mask ---
lllat, lllon, urlat, urlon = ctrack_para.ret_tcregionlatlon(region)
a2domain = ctrack_func.mk_region_mask(lllat, urlat, lllon, urlon, nx, ny, lat_first, lon_first, dlat, dlon)
#-- init ----------
lvort_rad = []
lvort = []
ldt = []
#------------------
for year_loop in range(iyear, eyear+1):
bestname = bestdir + "/Year.%04d.ibtracs_all.v03r04.csv"%(year_loop)
#-- open -------
f = open(bestname, "r")
lines = f.readlines()
f.close()
#-- init dictionary --
#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() thwind = tc_para.ret_thwind() thrvort = tc_para.ret_thrvort(sresol) thwcore = tc_para.ret_thwcore(sresol) #--------------------- nx_eqgrid = nradeqgrid*2 + 1 ny_eqgrid = nradeqgrid*2 + 1 #--------------------- (latmin, lonmin, latmax, lonmax) = ctrack_para.ret_tcregionlatlon(region) #------------------------------------------- miss_int = -9999 miss_out = -9999.0 miss_gpcp = -99999. #------------------------------------------ def readlatlon(fname): f = open(fname, "r") lines = f.readlines() f.close() lines = map(float, lines) return lines #-------- def latlon2yx(lat, lon, lat_first, lon_first, dlat, dlon): iy = int( (lat + 0.5*dlat - lat_first)/dlat )
