import matplotlib.pyplot as plt import matplotlib.dates from matplotlib.dates import DateFormatter, YearLocator, MonthLocator, DayLocator, HourLocator import os def epoch_to_datetime(i): return datetime.fromtimestamp(i), datetime.fromtimestamp(i).strftime('%H') def datetime_to_str(x): return x.strftime("%H") timer1 = datetime.now() model_start = datetime(2015,6,17,0) wrf_dir = '/uufs/chpc.utah.edu/common/home/horel-group4/model/bblaylock/WRF3.7_lake303_ember/WRFV3/test/em_real/' STN_NAME, STN_ID, LAT, LON = read_tslist.read_tslist(wrf_dir+"tslist") FIGDIR = '/uufs/chpc.utah.edu/common/home/u0553130/public_html/MS/timeheight/' if not os.path.exists(FIGDIR): os.makedirs(FIGDIR) for s in np.arange(0,len(STN_ID)): station = STN_ID[s] st_name = STN_NAME[s] print 'working on',station plt.cla() plt.clf() plt.close() # Open the geopotential height file dates, PH = WRF_timeseries.get_full_vert(wrf_dir+station+'.d02.PH',model_start)
from functions import wind_calcs, MesoWest_timeseries, WRF_timeseries, read_tslist, WRF_timeseries from datetime import datetime, timedelta import numpy as np import matplotlib.pyplot as plt from matplotlib.dates import DateFormatter, YearLocator, MonthLocator, DayLocator, HourLocator timer1 = datetime.now() start_date_early = datetime(2015, 6, 17, 0) start_date = datetime(2015, 6, 17, 0) end_date = datetime(2015, 6, 19, 7) FIGDIR = '/uufs/chpc.utah.edu/common/home/u0553130/public_html/MS/timeseries/early/303/' # Get list of TS station id tslist = '/uufs/chpc.utah.edu/common/home/horel-group4/model/bblaylock/WRF3.7_lake303_ember/WRFV3/test/em_real/tslist' ts = read_tslist.read_tslist(tslist) station_IDs = ts[1] station_names = ts[0] station = 'UKBKB' for i in np.arange(14, len(station_IDs)): plt.cla() plt.clf() plt.close() station = station_IDs[i] print "Working on", station station_name = station_names[i] if station == 'USDR1' or station == 'QHV': print "skipped" pass
from functions import wind_calcs, MesoWest_timeseries, WRF_timeseries, read_tslist, WRF_timeseries from datetime import datetime, timedelta import numpy as np import matplotlib.pyplot as plt from matplotlib.dates import DateFormatter, YearLocator, MonthLocator, DayLocator, HourLocator timer1 = datetime.now() start_date_early = datetime(2015,6,17,0) start_date = datetime(2015,6,17,0) end_date = datetime(2015,6,19,7) FIGDIR = '/uufs/chpc.utah.edu/common/home/u0553130/public_html/MS/timeseries/early/303/' # Get list of TS station id tslist = '/uufs/chpc.utah.edu/common/home/horel-group4/model/bblaylock/WRF3.7_lake303_ember/WRFV3/test/em_real/tslist' ts = read_tslist.read_tslist(tslist) station_IDs = ts[1] station_names = ts[0] station = 'UKBKB' for i in np.arange(14,len(station_IDs)): plt.cla() plt.clf() plt.close() station = station_IDs[i] print "Working on", station station_name = station_names[i] if station =='USDR1' or station =='QHV': print "skipped"