from JRA55 import jra55
import Front
from datetime import datetime, timedelta
front = Front.front("JRA55","bn.anl_p125")
DTime = datetime(2004,3,5)
a2f = front.load_tfront(DTime)
print a2f
from numpy import * from detect_fsub import * import Front from datetime import datetime Year = 2004 Mon = 6 Day = 8 Hour = 0 DTime = datetime(Year,Mon,Day) new = Front.front(model="JRA25",res="sa.one") path_new = new.path_potloc(DTime) srcnew1 = path_new.srcPath1 srcnew2 = path_new.srcPath2 print srcnew1 print srcnew2 srcold1 = "/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) srcold2 = "/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) srcrpr1 = "/media/disk2/out/JRA25/sa.one.temp2/6hr/front.t/%04d%02d/front.t.M1.%04d.%02d.%02d.%02d.sa.one"%(Year,Mon,Year,Mon,Day,Hour) srcrpr2 = "/media/disk2/out/JRA25/sa.one.temp2/6hr/front.t/%04d%02d/front.t.M2.%04d.%02d.%02d.%02d.sa.one"%(Year,Mon,Year,Mon,Day,Hour) a2new1 =fromfile(srcnew1,float32).reshape(180,360) a2new2 =fromfile(srcnew2,float32).reshape(180,360) a2old1 =fromfile(srcold1,float32).reshape(180,360) a2old2 =fromfile(srcold2,float32).reshape(180,360)
Year = 2004
lMD = [[3,6]]
#lMD = [[6,7]]
#lMD = [[9,8]]
#lMD = [[12,6]]
Hour = 0
lM1 = [0.26, 0.3, 0.34]
lM2 = [1.0, 1.4, 1.8]
#lM1 = [0.34]
#lM2 = [1.4]
BBox = [[20,120],[60,170]]
front = Front.front(model="JRA55",res="bn")
a1lat = front.Lat
a1lon = front.Lon
miss = front.miss
for Mon, Day in lMD:
Dtime = datetime(Year,Mon,Day,Hour)
for M1,M2 in [[M1,M2] for M1 in lM1 for M2 in lM2]:
a2front = front.load_tfront(Dtime, M1=M1, M2=M2)
a2front = ma.masked_equal(a2front, miss)
oPath = "/home/utsumi/temp/tune/%02d.%02d.M1.%04.2f.M2.%03.1f.png"%(Mon,Day,M1,M2)
Fig.DrawMap(a2front, a1lat, a1lon, BBox, figname = oPath, stitle="%02d-%02d M1=%f M2=%f"%(Mon,Day, M1,M2))
#- join --
i = -1
da2dat = {}
from numpy import * from datetime import datetime, timedelta import Front import Fig import Image import ctrack_func model = "JRA25" res = "sa.one" #model = "JRA55" #res = "bn" front = Front.front(model=model, res=res) a1lat = front.Lat a1lon = front.Lon miss = front.miss Year = 2004 Mon = 6 Day = 1 Hour = 0 BBox =[[20,120],[60,170]] figdir = "/home/utsumi/temp/tenkizu/%04d%02d"%(Year,Mon) Dtime = datetime(Year,Mon,Day,Hour) a2tfront = ma.masked_not_equal( front.load_tfront(Dtime), miss).filled(2.0) a2qfront = ma.masked_not_equal( front.load_qfront(Dtime), miss).filled(1.0) print a2qfront.mean(),a2tfront.mean() a2front = ma.masked_where(a2tfront==2.0, a2qfront).filled(2.0)
