Changeset 2965 for inundation/pyvolution/smf.py

Timestamp:

May 25, 2006, 9:25:03 AM (19 years ago)

Author:

sexton

Message:

update to velocity calculation

File:

• inundation/pyvolution/smf.py (modified) (1 diff)

inundation/pyvolution/smf.py

@@ -360 +360 @@
         z = zeros(N, Float)
         maxz = 0.0
+        minz = 0.0
         print 'here i am about to calculate z with dx = ', dx
         for i in range(N):
             try:
-                z[i] =  -am / ((cosh(kappa*(yr[i]-y0)/(wi+wa)))**2) \
+                z[i] =  -am / (amin*(cosh(kappa*(yr[i]-y0)/(wi+wa)))**2) \
                             * (exp(-((xr[i]-x0)/wa)**2) - \
                                 kappad*exp(-((xr[i]-dx-x0)/wa)**2))
                 if z[i] > maxz: maxz = z[i]
+                if z[i] < minz: minz = z[i]
             except OverflowError:
                 pass
         
         print 'maximum of Double Gaussian function', maxz
+        print 'minimum of Double Gaussian function', minz
         return z