WVEL always updraught

yoctoyotta1024 · Jan 16, 2024 · 643ecca · 643ecca
1 parent 47a6af4
commit 643ecca
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/pySD/thermobinary_src/thermogen.py b/pySD/thermobinary_src/thermogen.py
@@ -559,6 +559,8 @@ def wvel_profile(self, gbxbounds, ndims, ntime):
     zfaces = rgrid.coords_forgridboxfaces(gbxbounds, ndims, "z")[0]
     WVEL = self.WMAX * np.sin(np.pi * zfaces/(2*self.Wlength))
 
+    WVEL[WVEL < 0.0] = 0.0
+
     return np.tile(WVEL, ntime)
 
   def generate_winds(self, gbxbounds, ndims, ntime, THERMODATA):

diff --git a/scripts/create_thermobinaries_script.py b/scripts/create_thermobinaries_script.py
@@ -70,8 +70,8 @@
 TEMPlapses  = [9.8, 6.5]            # -dT/dz [K/km]
 qvaplapses  = [2.97, "saturated"]   # -dvap/dz [g/Kg km^-1]
 qcond       = 0.0                   # [Kg/Kg]
-WMAX        = 1.0                   # [m/s]
-Wlength     = Zbase                 # [m] use constant W (Wlength=0.0), or sinusoidal 1-D profile below cloud base
+WMAX        = 3.0                   # [m/s]
+Wlength     = 800                   # [m] use constant W (Wlength=0.0), or sinusoidal 1-D profile below cloud base
 
 thermodyngen = thermogen.ConstHydrostaticLapseRates(configfile, constsfile,
                                                     PRESS0, TEMP0, qvap0,