Changeset 5645 for anuga_work/production/onslow_2008/run_onslow.py

Timestamp:

Aug 13, 2008, 8:21:21 AM (16 years ago)

Author:

kristy

Message:

Updated evolve in run file so that a tide wave takes over when the sts file finishes

File:

• anuga_work/production/onslow_2008/run_onslow.py (modified) (4 diffs)

anuga_work/production/onslow_2008/run_onslow.py

@@ -42 +42 @@
 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
 from anuga.geospatial_data.geospatial_data import find_optimal_smoothing_parameter
-
+from Scientific.IO.NetCDF import NetCDFFile
 # Application specific imports
 import project                 # Definition of file names and polygons
@@ -83 +83 @@
     # Reading the landward defined points, this incorporates the original clipping
     # polygon minus the 100m contour
-    landward_bounding_polygon = read_polygon(project.polygons_dir+'landward_bounding_polygon.txt')
+    landward_bounding_polygon = read_polygon(project.boundaries_dir+'landward_bounding_polygon.txt')
 
     # Combine sts polyline with landward points
@@ -216 +216 @@
     Bd = Dirichlet_boundary([kwargs['tide'],0,0])
 
-    print dir(Bf)
-    print 'start reading boundary file'
-    
-
+    fid = NetCDFFile(boundary_urs_out+'.sts', 'r')    #Open existing file for read
+    sts_time=fid.variables['time'][:]+fid.starttime
+    tmin=min(sts_time)
+    tmax=max(sts_time)
+    fid.close()
+
+    print 'Boundary end time ', tmax-tmin
+     
 ##    Bf = Field_boundary(kwargs['boundary_file'],
 ##                domain, time_thinning=kwargs['time_thinning'], mean_stage=kwargs['tide'], 
@@ -241 +245 @@
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean')
+
+        if t >= tmax-tmin:
+            print 'changed to tide boundary condition at ocean'
+            domain.set_boundary({'ocean': Bd})
             
     x, y = domain.get_maximum_inundation_location()