case = '20130815' #ensembletype = 'STCH' ensembletype = 'MXMP' IC = 'GEFSR2' ICens = [ 'p09', ] #MP = 'Morrison_Hail' MP = 'ICBC' # NUMBER 3: SET THE TIME(s) YOU WANT PLOT # Loop over times... itime = (2013, 8, 15, 18, 0, 0) ftime = (2013, 8, 16, 12, 0, 0) times = utils.generate_times(itime, ftime, 60 * 60) # ... or set just one time skewT_time = (2013, 8, 16, 3, 0, 0) # NUMBER 4: SET THE LAT/LON skewT_latlon = (35.2435, -97.4708) # NUMBER 5: # This is a loop over different folders # For one plot, just use the (last line) plot_skewT command """ for en in ensnames: for ex in experiments: # Reload settings p.C = Settings()
#ensnames = ['c00'] + ['p'+"%02d" %n for n in range(1,11)] ensnames = ['p'+"%02d" %n for n in range(8,11)] #ensnames = ['p04'] #ensnames = ['anl'] experiment = {'ICBC':'CTRL'} itime = (2006,5,26,22,0,0) #itime = (2009,9,10,23,0,0) #itime = (2011,4,19,18,0,0) #itime = (2013,8,15,18,0,0) ftime = (2006,5,27,11,0,0) #ftime = (2009,9,11,14,0,0) #ftime = (2011,4,20,10,30,0) #ftime = (2013,8,16,11,30,0) times = utils.generate_times(itime,ftime,60*60) shear_times = utils.generate_times(itime,ftime,3*60*60) thresh = 10 variables = {'cref':{}, 'wind10':{}, 'buoyancy':{},'shear':{}, 'strongestwind':{}} variables['cref'] = {'lv':2000,'pt':times} variables['wind10'] = {'lv':2000,'pt':times} variables['buoyancy'] = {'lv':2000,'pt':times} variables['shear'] = {'top':3,'bottom':0,'pt':shear_times} variables['strongestwind'] = {'lv':2000, 'itime':itime, 'ftime':ftime, 'range':(thresh,27.5,1.25)} # variables = {'shear':{'pt':shear_times, 'top':3, 'bottom':0}} # variables = {'thetae':{'lv':2000,'pt':times}, 'CAPE':{'pt':times}} # variables = {'cref':{'lv':2000,'pt':times}, 'shear':{'pt':shear_times, 'top':3, 'bottom':0}} # variables = {'wind10':{'lv':2000,'pt':times}} #shear06 = {'shear':{'top':6,'bottom':0,'pt':shear_times}}
if wy == 2006 or wy == 2013: radar_view = "cent_plains" elif wy == 2009: radar_view = "north_plains" elif wy == 2011: radar_view = "upper_missvly" elif plot == "shear03": va_str = "shear" lv_str = "0to3" elif plot == "shear06": va_str = "shear" lv_str = "0to6" times = utils.generate_times(itime, ftime, 60 * 60 * 1) elif plot == "maxwind": va_str = "strongestwind" lv_str = "sfc" elif plot == "streamlines": va_str = "streamlines" lv_str = "sfc" if plot != "maxwind": time_strs = get_times(itime, ftime, int_hr) else: time_strs = ["range"] if plot == "streamlines":
p = WRFEnviron(config) if case[:4] == '2006': itime = (2006, 5, 26, 0, 0, 0) ftime = (2006, 5, 27, 15, 0, 0) elif case[:4] == '2009': itime = (2009, 9, 10, 0, 0, 0) ftime = (2009, 9, 11, 18, 0, 0) elif case[:4] == '2011': itime = (2011, 4, 19, 0, 0, 0) ftime = (2011, 4, 20, 15, 0, 0) else: itime = (2013, 8, 15, 0, 0, 0) ftime = (2013, 8, 16, 15, 0, 0) times = utils.generate_times(itime, ftime, 3 * 3600) if experiment == 'STCH': picklefolder = os.path.join(config.wrfout_root, case, IC, ens, MP) p.C.output_root = os.path.join(config.output_root, case, IC, ens, MP) path_to_wrfouts = [] for s in stoch_names: fpath = os.path.join(config.wrfout_root, case, IC, ens, MP, s) path_to_wrfouts.append(p.wrfout_files_in(fpath, dom=1)[0]) sensitivity = 0 elif experiment == 'ICBC': picklefolder = os.path.join(config.wrfout_root, case, IC) p.C.output_root = os.path.join(config.output_root, case, IC) sensitivity = ens_names
#ensnames = ['c00'] + ['p'+"%02d" %n for n in range(1,11)] ensnames = ['p' + "%02d" % n for n in range(8, 11)] #ensnames = ['p04'] #ensnames = ['anl'] experiment = {'ICBC': 'CTRL'} itime = (2006, 5, 26, 22, 0, 0) #itime = (2009,9,10,23,0,0) #itime = (2011,4,19,18,0,0) #itime = (2013,8,15,18,0,0) ftime = (2006, 5, 27, 11, 0, 0) #ftime = (2009,9,11,14,0,0) #ftime = (2011,4,20,10,30,0) #ftime = (2013,8,16,11,30,0) times = utils.generate_times(itime, ftime, 60 * 60) shear_times = utils.generate_times(itime, ftime, 3 * 60 * 60) thresh = 10 variables = { 'cref': {}, 'wind10': {}, 'buoyancy': {}, 'shear': {}, 'strongestwind': {} } variables['cref'] = {'lv': 2000, 'pt': times} variables['wind10'] = {'lv': 2000, 'pt': times} variables['buoyancy'] = {'lv': 2000, 'pt': times} variables['shear'] = {'top': 3, 'bottom': 0, 'pt': shear_times} variables['strongestwind'] = {
p = WRFEnviron(config) if case[:4] == '2006': itime = (2006,5,26,0,0,0) ftime = (2006,5,27,15,0,0) elif case[:4] == '2009': itime = (2009,9,10,0,0,0) ftime = (2009,9,11,18,0,0) elif case[:4] == '2011': itime = (2011,4,19,0,0,0) ftime = (2011,4,20,15,0,0) else: itime = (2013,8,15,0,0,0) ftime = (2013,8,16,15,0,0) times = utils.generate_times(itime,ftime,3*3600) if experiment=='STCH': picklefolder = os.path.join(config.wrfout_root,case,IC,ens,MP) p.C.output_root = os.path.join(config.output_root,case,IC,ens,MP) path_to_wrfouts = [] for s in stoch_names: fpath = os.path.join(config.wrfout_root,case,IC,ens,MP,s) path_to_wrfouts.append(p.wrfout_files_in(fpath,dom=1)[0]) sensitivity=0 elif experiment=='ICBC': picklefolder = os.path.join(config.wrfout_root,case,IC) p.C.output_root = os.path.join(config.output_root,case,IC) sensitivity=ens_names