Changeset 5645

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

Location:

anuga_work/production

Files:

• busselton/run_busselton.py (modified) (3 diffs)
• onslow_2008/run_onslow.py (modified) (4 diffs)
• perth/run_perth.py (modified) (3 diffs)

anuga_work/production/busselton/run_busselton.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
@@ -199 +199 @@
     Bd = Dirichlet_boundary([kwargs['tide'],0,0])
 
-    print dir(Bf)
-    print 'finished 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
 
     domain.set_boundary({'back': Bd,
@@ -220 +225 @@
         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()

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()

anuga_work/production/perth/run_perth.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
@@ -213 +214 @@
                    boundary_polygon=bounding_polygon)
 
+    
     Br = Reflective_boundary(domain)
     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 +247 @@
         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()