Changeset 4521

Timestamp:

May 30, 2007, 2:35:04 PM (18 years ago)

Author:

sexton

Message:

update to script for looking at "cannot find triangle for point"

Location:

anuga_work/production/onslow_2006

Files:

• export_results.py (modified) (3 diffs)
• test_export_results.py (modified) (1 diff)

anuga_work/production/onslow_2006/export_results.py

@@ -5 +5 @@
 from os import sep
 
-#time_dir = '20070507_215717_run_final_1.5_nbartzis' # HAT 1 in Dampier 10000 yr
+time_dir = '20070507_215717_run_final_1.5_nbartzis' # HAT 1 in Dampier 10000 yr
 #time_dir = '20070517_061654_run_final_0.0_nbartzis' # MSL 1 in Dampier 10000 yr
 #time_dir = '' # HAT 1 in Broome 10000 yr
@@ -15 +15 @@
 #time_dir = '' # HAT 1 in Exmouth 10000 yr
 #time_dir = '' # MSL 1 in Exmouth 10000 yr
-time_dir = '20060704_063005' # 2006 simulation HAT Mw9
-#cellsize = 30 
-cellsize = 250
+#time_dir = '20060704_063005' # 2006 simulation HAT Mw9
+cellsize = 25
+#cellsize = 150
 timestep = None
 directory = project.outputdir
 #name = directory + time_dir + sep + project.scenario_name
-#name = directory+time_dir+sep+project.scenario_name
-name = directory+time_dir+sep+'source'
+name = directory+time_dir+sep+project.scenario_name
+#name = directory+time_dir+sep+'source'
 
-is_parallel = False
-#is_parallel = True
+#is_parallel = False
+is_parallel = True
 if is_parallel == True: nodes = 4
 print 'output dir:', name
 
-#var = [2,3,4] # depth and speed
-var = [3] # depth and speed
+var = [2,3,4] # depth and speed
+#var = [3] # depth and speed
 
 for which_var in var:
@@ -46 +46 @@
 
     if which_var == 3:  # Speed
-        outname = name + '_speed'
+        outname = name + '_speed_150'
         quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
 

anuga_work/production/onslow_2006/test_export_results.py

@@ -39 +39 @@
             timestep = timestep,
             cellsize = cellsize,      
-            easting_min = project.e_min_area,
-            easting_max = project.e_max_area,
-            northing_min = project.n_min_area,
-            northing_max = project.n_max_area,        
+            #easting_min = project.e_min_area,
+            #easting_max = project.e_max_area,
+            #northing_min = project.n_min_area,
+            #northing_max = project.n_max_area,        
             reduction = max, 
             verbose = True,