##################################################### ##################################################### dlat_one = 1.0 dlon_one = 1.0 a1lat_one = arange(-89.5, 89.5+dlat_one*0.1, dlat_one) a1lon_one = arange(0.5, 359.5+dlon_one*0.1, dlat_one) ##################################################### for model in lmodel: for expr in lexpr: for var in lvar: #------------------ odir_root ="/media/disk2/data/CMIP5/sa.one.%s.%s/%s"%(model, expr, var) odir_dump = odir_root #------------------ ens = cmip_para.ret_ens(model, expr, var) ######################### # set nc dir ######################### #incdir = "/home/utsumi/mnt/export/nas_d/data/CMIP5/fx" incdir = "/home/utsumi/mnt/iis.data2/CMIP5/cmip5.working/%s.%s"%(model,expr) ######################### ihead = var + "_" + tstp + "_" +model + "_" + expr +"_"\ +ens ohead = "%s.%s"%(var,model) ##################################################### odir_dump = odir_root try: os.makedirs(odir_dump) except OSError:
a1lon = arange(0.5, 359.5 + dlon * 0.5, dlon) # ---------------------------- dpgradrange = ctrack_para.ret_dpgradrange() lclass = dpgradrange.keys() nclass = len(lclass) thpgrad = dpgradrange[1][0] # ---------------------------- lllat = -80.0 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)
dyrange = {"historical":[1980,1999], "rcp85":[2080,2099]} #dyrange = {"historical":[1980,1981], "rcp85":[2080,2081]} lseason = ["ALL"] ny,nx = 180,360 miss = -9999. #--------------------- #region = "GLOB" region = "JPN" lllon, lllat, urlon, urlat = chart_para.ret_domain_corner_rect_forfig(region) #********************************************** llkey = [[season,exprfut,model] for season in lseason for exprfut in lexprfut for model in lmodel] for season, exprfut, model in llkey: #------ ens = cmip_para.ret_ens(model, exprfut, "psl") sunit, scalendar = cmip_para.ret_unit_calendar(model,exprfut) iyear_fut,eyear_fut = dyrange[exprfut] iyear_his,eyear_his = dyrange[exprhis] lyear = range(iyear_fut,eyear_fut+1) lmon = ctrack_para.ret_lmon(season) #------ futdir = "/media/disk2/data/CMIP5/sa.one.%s.%s/pr/%04d-%04d.%s"%(model,exprfut,iyear_fut,eyear_fut,season) hisdir = "/media/disk2/data/CMIP5/sa.one.%s.%s/pr/%04d-%04d.%s"%(model,exprhis,iyear_his,eyear_his,season) futname = futdir + "/pr.%s.%s.%s.%04d-%04d.%s.sa.one"%(model,exprfut,ens,iyear_fut,eyear_fut,season) hisname = hisdir + "/pr.%s.%s.%s.%04d-%04d.%s.sa.one"%(model,exprhis,ens,iyear_his,eyear_his,season) a2fut = fromfile(futname, float32).reshape(ny,nx) a2his = fromfile(hisname, float32).reshape(ny,nx)