def mk_a2chart_gradtv_highside_saone(sresol, year, mon,day, hour): stime = "%04d%02d%02d%02d"%(year, mon, day, hour) plev = 850*100. nx,ny = 360,180 miss_out= -9999. highsidedist = ctrack_para.ret_highsidedist() #-- chart front loc -------------------------- frontname = "/media/disk2/out/chart/ASAS/front/%04d%02d/front.ASAS.%04d.%02d.%02d.%02d.saone"%(year,mon,year,mon,day,hour) a2loc = fromfile(frontname, float32).reshape(ny,nx) #-- q: mixing ratio -------------------------- qname = "/media/disk2/data/JRA25/sa.one.%s/6hr/SPFH/%04d%02d/anl_p.SPFH.%04dhPa.%04d%02d%02d%02d.sa.one"%(sresol, year, mon, plev*0.01, year, mon, day, hour) a2q = fromfile(qname, float32).reshape(ny,nx) #-- t: --------------------------------------- tname = "/media/disk2/data/JRA25/sa.one.%s/6hr/TMP/%04d%02d/anl_p.TMP.%04dhPa.%04d%02d%02d%02d.sa.one"%(sresol, year, mon, plev*0.01, year, mon, day, hour) a2t = fromfile(tname, float32).reshape(ny,nx) #-- tv: -------------------------------------- a2tv = a2t * (1.0+0.61*a2q) #-- a2gradtv --------------------------------- a2gradtv = dtanl_fsub.mk_a2grad_abs_saone(a2tv.T).T #-- theta_e ----------------------------------- a2theta_e = dtanl_fsub.mk_a2theta_e(plev, a2t.T, a2q.T).T #-- grad.theta_e ------------------------------------ a2thermo = a2theta_e a2gradtheta_e = dtanl_fsub.mk_a2grad_abs_saone(a2thermo.T).T a2out = ctrack_fsub.find_highsidevalue_saone(a2gradtheta_e.T, a2loc.T, a2gradtv.T, highsidedist, miss_out).T return a2out
eday = calendar.monthrange(year,mon)[1]
#-------------
for day in range(iday,eday+1):
for ftype in lftype:
if ftype == "t":
thfmask1 = thfmask1t # sample for "t"
thfmask2 = thfmask2t # sample for "t"
elif ftype == "q":
thfmask1 = thfmask1q # sample for "q"
thfmask2 = thfmask2q # sample for "q"
#-----------------------------
#-----------------------------
miss_out = -9999.0
ny = 180
nx = 360
highsidedist = ctrack_para.ret_highsidedist()
stime = "%04d%02d%02d%02d"%(year, mon, day, hour)
sodir = "/home/utsumi/temp/tenkizu"
cbardir = sodir
ctrack_func.mk_dir(sodir)
dlat = 1.0
dlon = 1.0
a1lat = arange(-89.5, 89.5 + dlat*0.1, dlat)
a1lon = arange(0.5, 359.5 + dlon*0.1, dlon)
meridians = 10.0
resol = "anl_p"
