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)