lmon = range(1,12+1) iday = 1 lhour = [0,6,12,18] miss = -9999.0 ny = 180 nx = 360 model = "org" dpgradrange = ctrack_para.ret_dpgradrange() thpgrad = dpgradrange[1][0] thdura = 36 #thwind = 15 #m/s #thrvort = 3.5e-5 #thwcore = 1.5 # (K) thwind = tc_para.ret_thwind() thrvort = tc_para.ret_thrvort() thwcore = tc_para.ret_thwcore() thsst = tc_para.ret_thsst() # #thdura = -9999.0 #thwind = -9999.0 #m/s #thrvort = -9999.0 #thwcore = -9999.0 # (K) plev_low = 850*100.0 # (Pa) plev_mid = 500*100.0 # (Pa) plev_up = 250*100.0 # (Pa) psldir_root = "/media/disk2/data/JRA25/sa.one.%s/6hr/PRMSL"%(model) pgraddir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/pgrad"%(model)
urlat = 80.0 lllon = 0.0 urlon = 360.0 # ---------------------------- for expr, model in [[expr, model] for expr in lexpr for model in lmodel]: # ---- ens = cmip_para.ret_ens(model, expr, "psl") sunit, scalendar = cmip_para.ret_unit_calendar(model, expr) # ---- # thsst = 273.15 + 25.0 # thwind = 0.0 #m/s # thrvort = 7.0e-5 # thwcore = 0.5 # (K) thsst = tc_para.ret_thsst() thwind = tc_para.ret_thwind() thrvort = tc_para.ret_thrvort(model) thwcore = tc_para.ret_thwcore(model) plev_low = 850 * 100.0 # (Pa) plev_mid = 500 * 100.0 # (Pa) plev_up = 250 * 100.0 # (Pa) tplev_low = 850 * 100.0 psldir_root = "/media/disk2/data/CMIP5/sa.one.%s.%s/psl" % (model, expr) pgraddir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/pgrad" % (model, expr) lifedir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/life" % (model, expr) lastposdir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/lastpos" % (model, expr) iposdir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/ipos" % (model, expr) nextposdir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/nextpos" % (model, expr) #
pgradmin = ctrack_para.ret_dpgradrange()[2][0] # Pa/1000km #-- orog ------------------------ orogname = "/media/disk2/data/JRA25/sa.one.125/const/topo/topo.sa.one" a2orog = fromfile(orogname, float32).reshape(ny,nx) a2shade = ma.masked_where(a2orog >thorog, a2orog).filled(miss) #---------------------------- for sresol in lsresol: #******************* # tc params #------------------- thsst = tc_para.ret_thsst() thwind = tc_para.ret_thwind() thrvort = tc_para.ret_thrvort(sresol) thwcore = tc_para.ret_thwcore(sresol) #------------------- psldir_root = "/media/disk2/data/JRA25/sa.one.%s/6hr/PRMSL"%(sresol) pgraddir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/pgrad"%(sresol) lifedir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/life"%(sresol) #************************************ # init #-------------------- a2count = zeros([ny,nx],float32) #------------------------------------ if calcflag ==True: for year in range(iyear, eyear+1): # TC ---- dbstxy = tc_func.ret_ibtracs_dpyxy_saone(year)