#####################################################
#####################################################
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)
