Changeset 6398

Timestamp:

Feb 24, 2009, 12:28:04 PM (16 years ago)

Author:

myall

Message:

ceduna working using older scripts - going to upgrade to the latest ones anyway, so this isn't very important

Location:

anuga_work/production/australia_ph2/ceduna

Files:

• build_boundary_64468.py (modified) (5 diffs)
• project.py (modified) (2 diffs)
• run_model.py (modified) (1 diff)

anuga_work/production/australia_ph2/ceduna/build_boundary_64468.py

@@ -63 +63 @@
             ordering_filename = urs_order,
             weights = weights,
+            zone=53,
             mean_stage = project.tide,
             verbose = True)
@@ -93 +94 @@
 
     maxname = 'max_sts_stage.csv'
-    fid_max = open(project.boundaries_dir_event+sep+maxname,'w')
+    fid_max = open(project.event_folder+sep+maxname,'w')
     s = 'index, x, y, max_stage \n'
     fid_max.write(s)    
@@ -110 +111 @@
 
     minname = 'min_sts_stage.csv'
-    fid_min = open(project.boundaries_dir_event+sep+minname,'w')
+    fid_min = open(project.event_folder+sep+minname,'w')
     s = 'index, x, y, max_stage \n'
     fid_min.write(s)    
@@ -124 +125 @@
 
 
-        fid_sts = open(project.boundaries_dir_event+sep+basename+'_'+ str(index)+'.csv', 'w')
+        fid_sts = open(project.event_folder+sep+basename+'_'+ str(index)+'.csv', 'w')
         s = 'time, stage, xmomentum, ymomentum \n'
         fid_sts.write(s)
@@ -142 +143 @@
     return quantities,elevation,time
 
-quantities,elevation,time=get_sts_gauge_data(os.path.join(project.boundaries_dir_event,project.scenario_name),verbose=False)
+quantities,elevation,time=get_sts_gauge_data(project.event_sts,verbose=False)
 
 print len(elevation), len(quantities['stage'][0,:])

anuga_work/production/australia_ph2/ceduna/project.py

@@ -40 +40 @@
 # run_scenario script is executed
 tide = 0 
-event_number = 27255 # Java 9.3 worst case for Perth
+#event_number = 27255 # Java 9.3 worst case for Perth
+event_number = 64468 # event for SA coast
 alpha = 0.1             # smoothing parameter for mesh
 friction = 0.01           # manning's friction coefficient
@@ -104 +105 @@
 # Thinned ordering file from Hazard Map (geographic)
 # Format is index,latitude,longitude (with header)
-urs_order_filename = 'urs_order.csv'
+urs_order_filename = 'urs_order_austtg.csv'
 
 # Landward bounding points

anuga_work/production/australia_ph2/ceduna/run_model.py

@@ -60 +60 @@
 # Read in boundary from ordered sts file
 event_sts = create_sts_boundary(project.event_sts)
+print 'HELLO event_sts', event_sts
 
 # Reading the landward defined points, this incorporates the original clipping
 # polygon minus the 100m contour
 landward_boundary = read_polygon(project.landward_boundary)
+print 'landward boundary', landward_boundary
 
 # Combine sts polyline with landward points
 bounding_polygon_sts = event_sts + landward_boundary
+print 'bounding polygon', bounding_polygon_sts
 
 # Number of boundary segments