def wrap_front_t_saone(year,mon,day,hour, miss=-9999.0): #-- ny,nx = 180,360 thorog = ctrack_para.ret_thorog() thgradorog = ctrack_para.ret_thgradorog() thgrids = front_para.ret_thgrids() thfmask1t, thfmask2t, thfmask1q, thfmask2q = front_para.ret_thfmasktq() # orogname = "/media/disk2/data/JRA25/sa.one.125/const/topo/topo.sa.one" gradname = "/media/disk2/data/JRA25/sa.one.125/const/topo/maxgrad.0200km.sa.one" a2orog = fromfile(orogname, float32).reshape(ny,nx) a2gradorog = fromfile(gradname, float32).reshape(ny,nx) iname1t = "/media/disk2/out/JRA25/sa.one.anl_p/6hr/front.t/%04d%02d/front.t.M1.%04d.%02d.%02d.%02d.sa.one"%(year,mon,year,mon,day,hour) iname2t = "/media/disk2/out/JRA25/sa.one.anl_p/6hr/front.t/%04d%02d/front.t.M2.%04d.%02d.%02d.%02d.sa.one"%(year,mon,year,mon,day,hour) a2pott1 = fromfile(iname1t,float32).reshape(ny,nx) a2pott2 = fromfile(iname2t,float32).reshape(ny,nx) a2loct = complete_front_t_saone(a2pott1, a2pott2, thfmask1t, thfmask2t, a2orog, a2gradorog, thorog, thgradorog, thgrids, miss ) return a2loct
def wrap_front_q_saone(year,mon,day,hour, miss=-9999.0): #-- para --- ny,nx = 180,360 thorog = ctrack_para.ret_thorog() thgradorog = ctrack_para.ret_thgradorog() thgrids = front_para.ret_thgrids() thfmask1t, thfmask2t, thfmask1q, thfmask2q = front_para.ret_thfmasktq() #-- names --- orogname = "/media/disk2/data/JRA25/sa.one.125/const/topo/topo.sa.one" gradname = "/media/disk2/data/JRA25/sa.one.125/const/topo/maxgrad.0200km.sa.one" iname1q = "/media/disk2/out/JRA25/sa.one.anl_p/6hr/front.q/%04d%02d/front.q.M1.%04d.%02d.%02d.%02d.sa.one"%(year,mon,year,mon,day,hour) iname2q = "/media/disk2/out/JRA25/sa.one.anl_p/6hr/front.q/%04d%02d/front.q.M2.%04d.%02d.%02d.%02d.sa.one"%(year,mon,year,mon,day,hour) #-- load ---- a2orog = fromfile(orogname, float32).reshape(ny,nx) a2gradorog = fromfile(gradname, float32).reshape(ny,nx) a2potq1 = fromfile(iname1q,float32).reshape(ny,nx) a2potq2 = fromfile(iname2q,float32).reshape(ny,nx) a2loct = wrap_front_t_saone(year,mon,day,hour, miss=-9999.0) #-- make locq -- a2locq = complete_front_q_saone(a2loct, a2potq1, a2potq2, thfmask1q, thfmask2q, a2orog, a2gradorog, thorog, thgradorog, thgrids, miss) # return a2locq
dlat = 1.0
dlon = 1.0
a1lat = arange(-89.5, 89.5 +dlat*0.5, dlat)
a1lon = arange(0.5, 359.5 +dlon*0.5, dlon)
#--------------
a2one = ones([ny,nx],float32)
#--------------
llkey = [[expr,model,season] for expr in lexpr for model in lmodel for season in lseason]
for expr, model, season in llkey:
lmon = ctrack_para.ret_lmon(season)
#----
ens = cmip_para.ret_ens(model, expr, "psl")
sunit, scalendar = cmip_para.ret_unit_calendar(model,expr)
iyear,eyear = dyrange[expr]
lyear = range(iyear,eyear+1)
thfmask1t, thfmask2t, thfmask1q, thfmask2q = front_para.ret_thfmasktq(model)
#-- orog ------------------------
orogname = "/media/disk2/data/CMIP5/sa.one.%s.%s/orog/orog.%s.sa.one"%(model,expr,model)
#--- init ------------------
a2count = zeros([ny,nx],float32)
#-----------------------------------------
a1dtime,a1tnum = cmip_func.ret_times(iyear,eyear,lmon,sunit,scalendar,stepday)
for year in lyear:
for mon in lmon:
#----------------------------
print "agg.exist.front.cmip","rad",countrad,year,mon
#-- load ----------------
#idir_root = "/media/disk2/out/CMIP5/sa.one.MIROC5.historical/6hr/front.t/freq"
idir_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/front.t/freq"%(model,expr)
region= "ASAS"
ny = 180
nx = 360
miss = -9999.0
thdist = front_para.ret_thdistkm() # (km)
#
sresol = "anl_p"
lftype = ["t","q"]
#-- para for objective locator -------------
plev = 850*100.0 # (Pa)
thorog = ctrack_para.ret_thorog()
thorog_front = ctrack_para.ret_thorog_front()
#thgradorog = ctrack_para.ret_thgradorog()
thgradorog = 9999e+5
thgrids = front_para.ret_thgrids()
thfmask1t, thfmask2t, thfmask1q, thfmask2q = front_para.ret_thfmasktq(sresol)
orogname = "/media/disk2/data/JRA25/sa.one.125/const/topo/topo.sa.one"
maxorogname = "/media/disk2/data/JRA25/sa.one.125/const/topo/maxtopo.0300km.sa.one"
gradname = "/media/disk2/data/JRA25/sa.one.125/const/topo/maxgrad.0200km.sa.one"
a2maxorog = fromfile(maxorogname, float32).reshape(ny,nx)
a2gradorog = fromfile(gradname, float32).reshape(ny,nx)
#-------------------------------------------
a2one = ones([ny,nx],float32)
#******************************
#-- out dir -----------------
odir_root = "/media/disk2/out/JRA25/sa.one.%s/6hr/tenkizu/front/agg"%(sresol)
#-----------------------
for ftype in lftype:
for year in range(iyear, eyear+1):
#---------------------
