Changeset 2393

Timestamp:

Feb 14, 2006, 12:52:07 PM (18 years ago)

Author:

sexton

Message:

Allowing domain to be passed in slump/slide function to allow georeferencing.

File:

• inundation/pyvolution/smf.py (modified) (2 diffs)

inundation/pyvolution/smf.py

@@ -51 +51 @@
                   gravity=9.8, gamma=1.85, \
                   massco=1, dragco=1, frictionco=0, psi=0, \
-                  dx=None, kappa=3.0, kappad=0.8, zsmall=0.01,
+                  dx=None, kappa=3.0, kappad=0.8, zsmall=0.01, \
+                  domain=None,
                   verbose=False):
 
     from math import sin, tan, radians, pi, sqrt, exp
 
+    if domain is not None:
+        xllcorner = domain.geo_reference.get_xllcorner()
+        yllcorner = domain.geo_reference.get_yllcorner()
+        x0 = x0 - xllcorner  # slump origin (relative)
+        y0 = y0 - yllcorner
+        
     #if width not provided, set to typical value
     if width is None:
@@ -177 +184 @@
                   gravity=9.8, gamma=1.85, \
                   massco=1, dragco=1, frictionco=0, \
-                  dx=None, kappa=3.0, kappad=0.8, zsmall=0.01,
+                  dx=None, kappa=3.0, kappad=0.8, zsmall=0.01, \
+                  domain=None,
                   verbose=False):
 
     from math import sin, radians, sqrt
+
+    if domain is not None:
+        xllcorner = domain.geo_reference.get_xllcorner()
+        yllcorner = domain.geo_reference.get_yllcorner()
+        x0 = x0 - xllcorner  # slump origin (relative)
+        y0 = y0 - yllcorner
 
     #if width not provided, set to typical value