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 )