# "RCE_3km_2mom_db_rlongup", # "RCE_3km_2mom_dm", # "RCE_3km_2mom_dm_lrz", "RCE_CHAN_S300", ] Nsims = len(SimList) InFileTemplate = "DIAGS/profiles_<SIM>.h5" OutFileTemplate = "Plots.py/RceSkewT_<SIM>.png" TempVname = "/tempk_prof" RhVname = "/relhum_prof" Zname = "/z_coords" LabelScheme = conf.SetLabelScheme() for isim in range(Nsims): Sim = SimList[isim] print( "********** Creating sounding plots for simulation: {0:s}".format(Sim)) print("") InFname = InFileTemplate.replace("<SIM>", Sim) InFile = h5py.File(InFname, mode='r') print(" Reading: {0:s} ({1:s})".format(InFname, TempVname)) print(" Reading: {0:s} ({1:s})".format(InFname, RhVname)) print(" Reading: {0:s} ({1:s})".format(InFname, Zname)) print("")
#!/usr/bin/env python3 import os import sys sys.path.append("{0:s}/etc/python/common".format(os.environ['HOME'])) sys.path.append("{0:s}/etc/python/nasa_wisc".format(os.environ['HOME'])) import matplotlib.pyplot as plt import numpy as np import h5py import PlotUtils as plu import ConfigRce as conf # color and label schemes for plotting ColorScheme = conf.SetColorScheme() LabelScheme = conf.SetLabelScheme() PlotList = [ # [ [ 'RCE_1km', 'RCE_1km_SM', 'RCE_1km_DM', 'RCE_1km_DP' ], r'$IPLAWS = 0$', 'RceIpl0', 0, 50 ], # [ [ 'RCE_1km_IPL2', 'RCE_1km_SM_IPL2', 'RCE_1km_DM_IPL2', 'RCE_1km_DP_IPL2' ], r'$IPLAWS = 2$', 'RceIpl2', 0, 50 ], # [ [ 'RCE_3km_1mom', 'RCE_3km_2mom', 'RCE_3km_2mom_db', 'RCE_3km_2mom_dm' ], '', 'Rce3km', 0, 80 ], [['RCE_3km_1mom', 'RCE_3km_1mom_db', 'RCE_3km_1mom_dm'], '', 'Rce3km1mom', 0, 120], [['RCE_3km_2mom_db', 'RCE_3km_2mom_dm', 'RCE_3km_2mom_dm_lrz'], '', 'Rce3km2mom', 0, 120], ] Nplots = len(PlotList) VarList = [ #[ 'DIAGS/dom_avg_pcprr_<SIM>.h5', 'pcprr', [ 0, 12 ], r'$PR\ (mm\ day^{-1})$', 'Prate', 24 ], #[ 'DIAGS/dom_avg_vint_vapor_<SIM>.h5', 'vint_vapor', [ 0, 60 ], r'$PW\ (mm)$', 'Pwater', 24 ],
#!/usr/bin/env python3 import os import sys sys.path.append("{0:s}/etc/python/common".format(os.environ['HOME'])) sys.path.append("{0:s}/etc/python/nasa_wisc".format(os.environ['HOME'])) import h5py import numpy as np import numpy.ma as ma import ConfigRce as conf import Hdf5Utils as h5u Tstring = conf.SetTimeString() # One argument which is the name of the simulation. Sim = sys.argv[1] Xname = '/x_coords' Yname = '/y_coords' Zname = '/z_coords' Tname = '/t_coords' RtFileTemplate = "SIMS/<SIM>/HDF5/rtotal-<SIM>-LS-2012-01-01-000000-g1.h5" RtVname = "/rtotal_orig" RvFileTemplate = "SIMS/<SIM>/HDF5/vapor-<SIM>-LS-2012-01-01-000000-g1.h5" RvVname = "/vapor" OutFileTemplate = "DIAGS/tcond_xsections_<SIM>.h5" OutVnameBase = "/tcond"
#!/usr/bin/env python3 import os import sys sys.path.append("{0:s}/etc/python/common".format(os.environ['HOME'])) sys.path.append("{0:s}/etc/python/nasa_wisc".format(os.environ['HOME'])) import matplotlib.pyplot as plt import numpy as np import h5py import PlotUtils as plu import ConfigRce as conf # label schemes for plotting LabelScheme = conf.SetLabelScheme() PlotList = [ # [ [ 'RCE_1km', 'RCE_1km_SM', 'RCE_1km_DM', 'RCE_1km_DP' ], r'$IPLAWS = 0$', 'RceIpl0' ], # [ [ 'RCE_1km_IPL2', 'RCE_1km_SM_IPL2', 'RCE_1km_DM_IPL2', 'RCE_1km_DP_IPL2' ], r'$IPLAWS = 2$', 'RceIpl2' ], # [ [ 'RCE_3km_1mom', 'RCE_3km_2mom', 'RCE_3km_2mom_db', 'RCE_3km_2mom_dm' ], '', 'Rce3km' ], [['RCE_3km_1mom', 'RCE_3km_1mom_db', 'RCE_3km_1mom_dm'], '', 'Rce3km1mom'], [['RCE_3km_2mom_db', 'RCE_3km_2mom_dm', 'RCE_3km_2mom_dm_lrz'], '', 'Rce3km2mom'], [[ 'RCE_3km_1mom_db', 'RCE_3km_2mom_db', ], '', 'RceDb'], ] Nplots = len(PlotList) VarList = [