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