Workshop

Matplotlib_common/latex_labels.py

@@ -1,10 +1,16 @@
 var_labels = {
+
+  "lwp" : 'LWP [kg m$^{-2}$]',
+  "rwp" : 'RWP [kg m$^{-2}$]',
+  "cwp" : 'CWP [kg m$^{-2}$]',
+  "lwm" : 'LWM [kg]',
+  "cwm" : 'CWM [kg]',
+  "rwm" : 'RWM [kg]',
   "sgs_tke" : 'SGS TKE (Smagorinsky) [m$^{2}$ s$^{-2}$]',
   "sgs_tke_sd" : 'SGS TKE (Super-droplets) [m$^{2}$ s$^{-2}$]',
-  "lwp" : 'LWP [g m$^{-2}$]',
-  "rwp" : 'RWP [g m$^{-2}$]',
   "surf_precip" : 'Surface precip. [mm/day]',
   "acc_precip" : 'Accumulated precip. [mm]',
+  "acc_precip_vol" : 'Accumulated precip. [m^3]',
   "cl_nc" : '$N_c$ [cm$^{-3}$] (cloudy cells)',
   "cl_nr" : '$N_{r>25\mu m}$ [cm$^{-3}$]',
   "cl_gccn_conc" : '$N_{GCCN}$ [cm$^{-3}$] (cloudy cells)',
@@ -23,6 +29,8 @@
   "er" : 'Entrainment rate [cm s$^{-1}$]',
   "wvarmax" : 'Max. $w$ variance [m$^{2}$ s$^{-2}$]',
   "cloud_base" : 'Cloud base height [m]',
+  "cloud_top_height" :  'Cloud top height [m]',
+  "total_droplets_number" : 'Total Number of Cloud and Rain Droplets [#]',
   "gccn_rw_cl" : '$<r>$ of GCCN droplets (cloudy cells) [um]',
   "non_gccn_rw_cl" : '$<r>$ of CCN droplets (cloudy cells) [um]',
   "clb_bigrain_mean_rd" : '$<r_d>$ of (r$>$40um) @ clbase [m]',

Matplotlib_common/plot_profs.py

@@ -19,7 +19,8 @@ def plot_profiles(var_list, plot_iter, nplotx, nploty, axarr, xscaledict, yscale
       print file_name
     #  try:
       profiles_file = open(file_name, "r")
-      my_pos = read_my_var(profiles_file, "position")
+      ## my_pos = read_my_var(profiles_file, "position")
+      my_pos = read_my_var(profiles_file, "time")
       my_res = read_my_var(profiles_file, var)
       profiles_file.close()
 

Matplotlib_common/plot_series.py

@@ -17,11 +17,12 @@ def plot_series(var_list, plot_iter, nplotx, nploty, axarr, xscaledict, yscaledi
     for file_name in file_names:
       print file_name, var
       series_file = open(file_name, "r")
-      my_times = read_my_var(series_file, "position")
+      #my_times = read_my_var(series_file, "position")
+      my_time = read_my_var(series_file, "time")
       my_res = read_my_var(series_file, var)
 
       # rescale time to hours
-      my_times = my_times / 3600.
+      my_times = my_times / 3600. ##H
 
       series_file.close()
 

Matplotlib_common/read_UWLCM_arrays.py

@@ -7,7 +7,7 @@ def read_my_array(file_obj):
   line = line.split(" ")
   del line[0]
   del line[len(line)-1]
-  arr = map(float,line)
+  arr =list( map(float,line))
   return np.array(arr), arr_name
 
 def read_my_var(file_obj, var_name):

cases/ICMW2020_cumulus_congestus/ICMW2020_cumulus_congestus_series_comparison.py

@@ -0,0 +1,74 @@
+from matplotlib import rc
+import matplotlib.pyplot as plt
+from matplotlib.ticker import AutoMinorLocator, MultipleLocator
+
+import os
+import sys
+sys.path.append(os.path.dirname(os.path.abspath(__file__)) + "/../Matplotlib_common/")
+
+from plot_ranges import xscaledict, yscaledict, xlimdict_series, ylimdict_series
+from plot_series import *
+
+# activate latex text rendering
+rc('text', usetex=True)
+
+ICMW2020_vars = ["RH_max", "cloud_top_height", "surf_precip", "acc_precip", "acc_vol_precip"]
+
+# init the plot
+nplotx = 2
+nploty= 3
+fig, axarr = plt.subplots(nplotx,nploty)
+
+plot_series( ICMW2020_vars, 0, nplotx, nploty, axarr, xscaledict, yscaledict, xlimdict_series, ylimdict_series, xlabel='Time [h]')
+
+# legend font size
+plt.rcParams.update({'font.size': 8})
+
+# hide axes on empty plots
+if len( ICMW2020_vars) % nploty == 0:
+  nemptyplots = 0
+else:
+  nemptyplots = nploty - len( ICMW2020_vars) % nploty
+  emptyplots = np.arange(nploty - nemptyplots, nploty)
+for empty in emptyplots:
+  axarr[nplotx-1, empty].axis('off')
+
+# hide hrzntl tic labels
+x_empty_label = np.arange(0, nplotx-1)
+y_empty_label = np.arange(nploty)
+for x in x_empty_label:
+  for y in y_empty_label:
+    axarr[x,y].set_xticklabels([])
+
+#axes = plt.gca()
+#axes.tick_params(direction='in')
+x_arr = np.arange(nplotx)
+y_arr = np.arange(nploty)
+for x in x_arr:
+  for y in y_arr:
+    #tics inside
+    axarr[x,y].tick_params(direction='in', which='both', top=1, right=1)
+    #minor tics
+    axarr[x,y].xaxis.set_minor_locator(AutoMinorLocator())
+    axarr[x,y].yaxis.set_minor_locator(AutoMinorLocator())
+    #labels and tics font size
+    for item in ([axarr[x,y].xaxis.label, axarr[x,y].yaxis.label] + axarr[x,y].get_xticklabels() + axarr[x,y].get_yticklabels()):
+      item.set_fontsize(8)
+      # subplot numbering
+    if y < nploty - nemptyplots or x < (nplotx - 1):
+      axarr[x,y].text(0.2, 0.9, labeldict[y + x*nploty], fontsize=8, transform=axarr[x,y].transAxes)
+
+#single legend for the whole figure
+handles, labels = axarr[0,0].get_legend_handles_labels()
+
+lgd = fig.legend(handles, labels, handlelength=4, loc='lower center', bbox_to_anchor=(0.45,0))
+
+#figure size
+fig.set_size_inches(7.874, 5. + (len(labels) - 2) * 0.2)# 5.214)#20.75,13.74)
+#distances between subplots and from bottom of the plot
+fig.subplots_adjust(bottom=0.18 + (len(labels) - 2) * 0.03, hspace=0.1, wspace=0.4)
+
+#fig.tight_layout(pad=0, w_pad=0, h_pad=0)
+
+#plt.show()
+fig.savefig(argv[len(sys.argv)-1], bbox_inches='tight', dpi=300)#, bbox_extra_artists=(lgd,))

cases/ICMW2020_cumulus_congestus/plot_ranges.py

@@ -0,0 +1,124 @@
+xscaledict = {
+  "thl" : "linear",
+  "00rtot" : "linear",
+  "rliq" : "linear",
+  "prflux" : "linear",
+  "cl_nc" : "linear",
+  "clfrac" : "linear",
+  "wvar" : "linear",
+  "w3rd" : "linear",
+  "sat_RH" : "linear",
+  "rad_flx" : "linear",
+  "lwp" : "linear",
+  "RH_max" : "linear",
+  "cloud_top_height" : "linear",
+  "total_droplets_number" : "linear",
+  "acc_precip_vol" : "linear",
+  "rwp" : "linear",
+  "er" : "linear",
+  "wvarmax" : "linear",
+  "surf_precip" : "linear",
+  "acc_precip" : "linear",
+  "cloud_base" : "linear",
+  "gccn_rw_cl" : "linear",
+  "non_gccn_rw_cl" : "linear",
+  "base_prflux_vs_clhght" : "log",
+  "cl_gccn_conc" : "linear"
+  }
+
+yscaledict = {
+  "thl" : "linear",
+  "00rtot" : "linear",
+  "rliq" : "linear",
+  "prflux" : "linear",
+  "cl_nc" : "linear",
+  "clfrac" : "linear",
+  "wvar" : "linear",
+  "w3rd" : "linear",
+  "sat_RH" : "linear",
+  "rad_flx" : "linear",
+  "lwp" : "linear",
+  "RH_max" : "linear",
+  "cloud_top_height" : "linear",
+  "total_droplets_number" : "linear",
+  "ac_precip_vol" : "linear",
+  "rwp" : "linear",
+  "er" : "linear",
+  "wvarmax" : "linear",
+  "surf_precip" : "linear",
+  "acc_precip" : "linear",
+  "cloud_base" : "linear",
+  "gccn_rw_cl" : "linear",
+  "non_gccn_rw_cl" : "linear",
+  "base_prflux_vs_clhght" : "linear",
+  "cl_gccn_conc" : "log"
+                                                  }
+
+xlimdict_profs = {
+  "thl" : None,
+  "00rtot" : None,
+  "rliq" : None,
+  "prflux" : (0,20),
+  "cl_nc" : (0,90),
+  "clfrac" : None,
+  "wvar" : None,
+  "w3rd" : None,
+  "sat_RH" : None,
+  ##"RH_max" : None,
+  ##"cloud_top_height" : None,    
+  "rad_flx" : None,
+  "gccn_rw_cl" : (0,90),
+  "non_gccn_rw_cl" : (0,12),
+  "base_prflux_vs_clhght" : (1,10000)
+  }
+ylimdict_profs = {
+  "thl" : (0,3000),
+  "00rtot" : (0,3000),
+  "rliq" : (0,3000),
+  "prflux" : (0,3000),
+  "cl_nc" : (0,3000),
+  "clfrac" : (0,3000),
+  "wvar" : (0,3000),
+  "w3rd" : (0,3000),
+  "sat_RH" : (0,3000),
+  "rad_flx" : (0,3000),
+  "RH_max" : (0,3000),
+  "cloud_top_height" : (0,3000),
+  "gccn_rw_cl" : (0,3000),
+  "non_gccn_rw_cl" : (0,3000),
+  "base_prflux_vs_clhght" : (0,2500)
+  }
+
+xlimdict_series = {
+  "clfrac" :  None,
+  "cl_nc" :  None,
+  "lwp" :  None,
+  "RH_max" : None,
+  "cloud_top_height" : None,
+  "total_droplets_number" : None,
+  "rwp" :  None,
+  "er" :  None,
+  "wvarmax" :  None,
+  "surf_precip" :  None,
+  "acc_precip" :  None,
+  "acc_precip_vol" : None,
+  "cl_gccn_conc" :  None,
+  "cloud_base" :  None
+  }
+
+ylimdict_series = {
+  "clfrac" : None,
+  "cl_nc" : None,
+  "lwp" : None,
+  "RH_max" : None,
+  "cloud_top_height" : None,    
+  "total_droplets_number" : None,
+  "acc_precip_vol" : None,
+  "rwp" : None,
+  "er" : None,
+  "wvarmax" : None,
+  "surf_precip" : None,
+  "acc_precip" : None,#(0,0.07),
+  "cl_gccn_conc" : (1e-6, 1),
+  "cloud_base" : None
+  }

cases/Lasher_Trapp/plot_ranges.py

@@ -10,6 +10,10 @@
   "sat_RH" : "linear",
   "rad_flx" : "linear",
   "lwp" : "linear",
+  "RH_max" : "linear",
+  "cloud_top_height" : "linear",
+  "total_droplets_number" : "linear",
+  "acc_precip_vol" : "linear",
   "rwp" : "linear",
   "er" : "linear",
   "wvarmax" : "linear",
@@ -34,6 +38,10 @@
   "sat_RH" : "linear",
   "rad_flx" : "linear",
   "lwp" : "linear",
+  "RH_max" : "linear",
+  "cloud_top_height" : "linear",
+  "total_droplets_number" : "linear",
+  "acc_precip_vol" : "linear",
   "rwp" : "linear",
   "er" : "linear",
   "wvarmax" : "linear",
@@ -56,6 +64,8 @@
   "wvar" : None,
   "w3rd" : None,
   "sat_RH" : None,
+  "RH_max" : None,
+  "cloud_top_height" : None,    
   "rad_flx" : None,
   "gccn_rw_cl" : (0,90),
   "non_gccn_rw_cl" : (0,12),
@@ -72,6 +82,8 @@
   "w3rd" : (0,3000),
   "sat_RH" : (0,3000),
   "rad_flx" : (0,3000),
+  "RH_max" : (0,3000),
+  "cloud_top_height" : (0,3000),
   "gccn_rw_cl" : (0,3000),
   "non_gccn_rw_cl" : (0,3000),
   "base_prflux_vs_clhght" : (0,2500)
@@ -81,11 +93,15 @@
   "clfrac" :  None,
   "cl_nc" :  None,
   "lwp" :  None,
+  "RH_max" : None,
+  "cloud_top_height" : None,
+  "total_droplets_number" : None,
   "rwp" :  None,
   "er" :  None,
   "wvarmax" :  None,
   "surf_precip" :  None,
   "acc_precip" :  None,
+  "acc_precip_vol" : None,
   "cl_gccn_conc" :  None,
   "cloud_base" :  None
   }
@@ -94,6 +110,10 @@
   "clfrac" : None,
   "cl_nc" : None,
   "lwp" : None,
+  "RH_max" : None,
+  "cloud_top_height" : None,    
+  "total_droplets_number" : None,
+  "acc_precip_vol" : None,
   "rwp" : None,
   "er" : None,
   "wvarmax" : None,

drawbicyc/CMakeLists.txt

@@ -46,9 +46,9 @@ foreach(exe drawbicyc plot_series_compare average)
       $<BUILD_INTERFACE:${CMAKE_SOURCE_DIR}/include>
   )
   # enabling c++14
-  target_compile_features(${exe} PUBLIC cxx_std_14)
-  # enabling c++11
-  # target_compile_features(${exe} PUBLIC cxx_std_11)
+  #target_compile_features(${exe} PUBLIC cxx_std_14)
+  #enabling c++11
+  target_compile_features(${exe} PUBLIC cxx_std_11)
   target_link_libraries(${exe} PRIVATE uwlcm::uwlcm_includes)
   target_link_libraries(${exe} PRIVATE ${HDF5_LIBRARIES})
   target_include_directories(${exe} PUBLIC ${HDF5_INCLUDE_DIRS})

drawbicyc/drawbicyc.cpp

@@ -19,7 +19,7 @@ int main(int argc, char** argv)
     ("qv_qc_2_6_10_min", po::value<bool>()->default_value(false) , "plot comparison of qv and qc fields at 2, 6 and 10 min?")
     ("dir", po::value<std::string>()->required() , "directory containing out_lgrngn")
     ("micro", po::value<std::string>()->required(), "one of: blk_1m, blk_2m, lgrngn")
-    ("type", po::value<std::string>()->required(), "one of: dycoms, moist_thermal, rico, Lasher_Trapp")//, base_prflux_vs_clhght")
+    ("type", po::value<std::string>()->required(), "one of: dycoms, moist_thermal, rico, Lasher_Trapp, ICMW2020")//, base_prflux_vs_clhght")
   ;
 
   po::variables_map vm;
@@ -33,8 +33,8 @@ int main(int argc, char** argv)
 
   // handling the "type" option
   std::string type = vm["type"].as<std::string>();
-  if(type != "dycoms" && type != "moist_thermal" && type != "rico" && type != "pi_chamber" && type != "Lasher_Trapp")
-    throw std::runtime_error("Unrecognized 'type' option, only dycoms, rico, moist_thermal, pi_chamber, Lasher_Trapp  available now");//, base_prflux_vs_clhght available now");
+  if(type != "dycoms" && type != "moist_thermal" && type != "rico" && type != "pi_chamber" && type != "Lasher_Trapp" && type != "ICMW2020")
+    throw std::runtime_error("Unrecognized 'type' option, only dycoms, rico, moist_thermal, pi_chamber, Lasher_Trapp, ICMW2020  available now");//, base_prflux_vs_clhght available now");
 
   // should profiles be normalized by inversion height
   const bool normalize_prof = type == "dycoms";

drawbicyc/include/drawbicyc/Plotter3d.hpp

@@ -25,7 +25,7 @@ class Plotter_t<3> : public PlotterCommon
   auto h5load(
     const string &file, 
     const string &dataset,
-    bool srfc = false // load a horiznotal slice 
+    bool srfc = false // load a horizontal slice 
   ) -> decltype(blitz::safeToReturn(arr_t() + 0))
   {
     parent_t::h5load(file, dataset);
@@ -189,7 +189,7 @@ class Plotter_t<3> : public PlotterCommon
     this->map["dz"] = tmp(0,0,1) - tmp(0,0,0);
     this->CellVol = this->map["dx"] * this->map["dy"] * this->map["dz"];
     this->DomainSurf = this->map["dx"] * this->map["dy"] * this->map["x"] * this->map["y"];
-
+    this->DomainVolume = this->map["dx"] * this->map["dy"] * this->map["dz"] * this->map["x"]  * this->map["y"] * this->map["z"];
 
     // other dataset are of the size x*z, resize tmp
     tmp.resize(n[0]-1, n[1]-1, n[2]-1);

drawbicyc/include/drawbicyc/PlotterCommon.hpp

@@ -10,7 +10,7 @@ class PlotterCommon
   const string file;
   std::map<std::string, double> map;
   blitz::Array<float, 1> timesteps, p_e;
-  double CellVol, DomainSurf;
+  double CellVol, DomainSurf, DomainVolume;
 
   protected:
   H5::H5File h5f;