# ------------------------ # ************************ # for objective front locator # ----------------------------- # Name frontdir_t_root = "/media/disk2/out/CMIP5/sa.one.%s.%s/6hr/front.t" % (model, expr) frontdir_t = frontdir_t_root + "/%04d%02d" % (year, mon) # fronttname1 = frontdir_t + "/front.t.M1.%04d.%02d.%02d.%02d.sa.one" % (yeart, mont, dayt, hourt) fronttname2 = frontdir_t + "/front.t.M2.%04d.%02d.%02d.%02d.sa.one" % (yeart, mont, dayt, hourt) # -- front.t --- a2fbc1 = fromfile(fronttname1, float32).reshape(ny, nx) a2fbc2 = fromfile(fronttname2, float32).reshape(ny, nx) a2fbc = front_func.complete_front_t_saone( a2fbc1, a2fbc2, thfmask1t, thfmask2t, a2orog, a2gradorog, thorog, thgradorog, thgrids, miss ) # -- count baloclinic front loc -- a2loc = a2fbc a2fbc_terr = ctrack_fsub.mk_territory_saone(a2loc.T, countrad * 1000.0, miss, -89.5, 1.0, 1.0).T a2fbc_terr = ma.masked_where(a2fbc_terr == miss, a2one).filled(0.0) # ************************* # count ExC + Front # ------------------------- a2count_tmp = ma.masked_greater(a2c_terr + a2fbc_terr, 0.0).filled(1.0) a2count = a2count + a2count_tmp # **************************** # write
import front_func import ctrack_para, front_para thfmask1t, thfmask2t, thfmask1q, thfmask2q = front_para.ret_thfmasktq() ny,nx = 180,360 thorog = ctrack_para.ret_thorog() thgradorog = ctrack_para.ret_thgradorog() thgrids = front_para.ret_thgrids() miss = -9999.0 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/200404/front.t.M1.2004.04.06.00.sa.one" iname2t = "/media/disk2/out/JRA25/sa.one.anl_p/6hr/front.t/200404/front.t.M2.2004.04.06.00.sa.one" a2pott1 = fromfile(iname1t,float32).reshape(180,360) a2pott2 = fromfile(iname2t,float32).reshape(180,360) a2loct = front_func.complete_front_t_saone(a2pott1, a2pott2, thfmask1t, thfmask2t, a2orog, a2gradorog, thorog, thgradorog, thgrids, miss ) print a2loct.max() iname1q = "/media/disk2/out/JRA25/sa.one.anl_p/6hr/front.q/200404/front.q.M1.2004.04.06.00.sa.one" iname2q = "/media/disk2/out/JRA25/sa.one.anl_p/6hr/front.q/200404/front.q.M2.2004.04.06.00.sa.one" a2potq1 = fromfile(iname1q,float32).reshape(180,360) a2potq2 = fromfile(iname2q,float32).reshape(180,360) a2locq = front_func.complete_front_q_saone(a2loct, a2potq1, a2potq2, thfmask1q, thfmask2q, a2orog, a2gradorog, thorog, thgradorog, thgrids, miss)