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_vol_precip" : 'Accumulated volume precipitation [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,13 @@ 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
+      my_times = my__times / 60.  ##Minutes
 
       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/Lasher_Trapp/Lasher_Trapp_series_comparison.py

@@ -12,7 +12,8 @@
 # activate latex text rendering
 rc('text', usetex=True)
 
-Lasher_Trapp_vars = ["lwp", "rwp", "surf_precip", "acc_precip", "cl_nc"]
+##Lasher_Trapp_vars = ["lwp", "rwp", "surf_precip", "acc_precip", "cl_nc"]
+Lasher_Trapp_vars = ["RH_max", "cloud_top_height", "surf_precip", "acc_precip", "acc_vol_precip"]
 
 # init the plot
 nplotx = 2

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_vol_precip" : "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_vol_precip" : "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_vol_precip" : 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_vol_precip" : 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/include/drawbicyc/Plotter3d.hpp

@@ -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;

drawbicyc/include/drawbicyc/PlotterMicro.hpp

@@ -307,7 +307,16 @@ class PlotterMicro_t : public Plotter_t<NDims>
                      / 1e3         // to m^3 of water
                      * 1e3;        // to mm
   }
-
+  // accumulated volume precipitation [m^3]
+  double calc_acc_volume_precip(double prec_vol)
+  {
+    if(this->micro == "lgrngn")
+      return prec_vol / this->DomainVolume;
+    if(this->micro == "blk_1m")
+      return prec_vol * this->map["dt"]
+                    / (this->map["x"] * this->map["y"] * this->map["z"])
+                    / 1e3;         // to m^3 of water
+  }
   //ctor
   PlotterMicro_t(const string &file, const string &micro):
     parent_t(file),

drawbicyc/include/drawbicyc/cases/Lasher_Trapp/plots.hpp

@@ -2,26 +2,34 @@
 
 
 const std::vector<std::string> series_Lasher_Trapp({
-  "clfrac",
+  //"clfrac",
   "lwp",
   "rwp",
+  "cwp",
+  "lwm",
+  "cwm",
+  "rwm",
   "surf_precip",
   "acc_precip",
-  "cl_nc",
+  "acc_vol_precip",
+  //"cl_nc",
   "cloud_base",
-  "surf_flux_latent",
-  "surf_flux_sensible",
-  "sd_conc_avg",
+  "RH_max",
+  "cloud_top_height",
+  "total_droplets_number",
+  //"surf_flux_latent",
+  //"surf_flux_sensible",
+  //"sd_conc_avg",
   //"mass_dry",
-  "cl_gccn_conc",
-  "gccn_conc",
-  "cl_non_gccn_conc",
-  "non_gccn_conc",
-  "cl_gccn_to_non_gccn_conc_ratio"
+  //"cl_gccn_conc",
+  //"gccn_conc",
+  //"cl_non_gccn_conc",
+  //"non_gccn_conc",
+  //"cl_gccn_to_non_gccn_conc_ratio"
   //, "cl_gccn_meanr"
   //,"cl_avg_cloud_rad"
   // "sd_conc_std_dev",
-  // // "tot_water"
+ // "tot_water"
 });
 
 std::vector<std::string> profs_Lasher_Trapp({

drawbicyc/include/drawbicyc/plot_prof.hpp

@@ -785,7 +785,8 @@ void plot_profiles(Plotter_t plotter, Plots plots, std::string type, const bool
 //    res_pos = i * n["dz"] / z_i; 
     if(!res_pos_out_done)
     {
-      oprof_file << "position" << endl;
+      //oprof_file << "position" << endl;
+      oprof_file << "time" << endl;
       oprof_file << res_pos;
       res_pos_out_done = true;
     }