diff --git a/pySD/thermobinary_src/thermogen.py b/pySD/thermobinary_src/thermogen.py
index 08a73358..b1020454 100644
--- a/pySD/thermobinary_src/thermogen.py
+++ b/pySD/thermobinary_src/thermogen.py
@@ -553,17 +553,14 @@ def hydrostatic_lapserates_thermo(self, zfulls):
     return TEMP, PRESS, qvap
 
   def wvel_profile(self, gbxbounds, ndims, ntime):
-    ''' returns WVEL=0.0 above cloud base, and sinusoidal profile
-    with amplitude WMAX and wavelength 2*Wlength below cloud base'''
+    ''' returnssinusoidal profile with amplitude WMAX
+    and wavelength 2*Wlength '''
 
     zfaces = rgrid.coords_forgridboxfaces(gbxbounds, ndims, "z")[0]
     WVEL = self.WMAX * np.sin(np.pi * zfaces/(2*self.Wlength))
     
-    WVEL = np.where(zfaces < self.Zbase, WVEL, 0.0)
-    
     return np.tile(WVEL, ntime)
 
-
   def generate_winds(self, gbxbounds, ndims, ntime, THERMODATA):
 
     THERMODATA = constant_winds(ndims, ntime, THERMODATA,