Changeset 7606

Timestamp:

Jan 12, 2010, 1:30:41 PM (13 years ago)

Author:

jgriffin

Message:

Added extra events for probabilistic run

Location:

anuga_work/production/new_south_wales/batemans_bay

Files:

• build_urs_boundary.py (modified) (1 diff)
• export_results_max.py (modified) (3 diffs)
• project.py (modified) (1 diff)
• run_multiple_events.py (modified) (1 diff)

anuga_work/production/new_south_wales/batemans_bay/build_urs_boundary.py

@@ -15 +15 @@
 from anuga.shallow_water.data_manager import urs2sts
 
-import numpy as num
+import Numeric as num
 
 

anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

@@ -23 +23 @@
 
 #Specify output directories
-time_dir1 = '20090601_172248_run_final_0.0_58284__250m_jgriffin'
+time_dir1 = '20090528_153025_run_final_0.0_58346_jgriffin'
 ##time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin'
 
@@ -30 +30 @@
 # sww filename extensions ie. if batemans_bay_time_37860_0.sww, input into list 37860
 # make sure numbers are in sequential order
-times = []
+times = [37860]
 
 #Modify the cellsize value to set the size of the raster you require
 #Take into account mesh size when aplying this paramater
-cellsize = 250   #dependent on data resolution in area of interest.
+cellsize = 20   #dependent on data resolution in area of interest.
 
 #Now set the timestep at which you want the raster generated.
@@ -57 +57 @@
                  'depth':     'stage-elevation',
                  'speed':     '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6)',
+                 'energy':  '(((xmomentum/(stage-elevation+1.e-6))**2 + (ymomentum/(stage-elevation+1.e-6))**2)*0.5*1000*(stage-elevation+1.e-6))+(9.81*stage*1000)',
                  'elevation': 'elevation' }
 
 # one or more key strings from var_equations above
 #var = ['depth', 'speed','stage']
-var = ['stage','speed','elevation']
+var = ['energy']
 ######
 # Start script, running through variables, area, folder, sww file (determine by times)

anuga_work/production/new_south_wales/batemans_bay/project.py

@@ -53 +53 @@
 else:    
     event_number = 58284
-#event_number_list = [31583,31602,31714,31913] #7871,#1.5-2.0m
-##event_number_list = [51268,51375,51436,51452] #31997, #1.5-2.0m
-##event_number_list =  [58286,58318,58337,58350] #58176,#1.5-2.0m
-
-##event_number_list = [31831,31897,31961,31973]#31708,
-##event_number_list = [51270,51328,51415,51460]#31981,
-##event_number_list = [58274,58331,58348,58357,58360]#51470,
+##event_number_list = [31602,31714,31913] #7871,31583,#1.5-2.0m
+##event_number_list = [51375,51436,51452] #31997,51268, #1.5-2.0m
+##event_number_list =  [58318,58337,58350] #58176,58286,#1.5-2.0m
+##event_number_list = [31897,31961,31973]#31708,31831,#2.1-3.0m
+
+##event_number_list = [51328,51415,51460]#31981,51270,#2.1-3.0m
+##event_number_list = [58331,58348,58357,58360]#51470,58274,#2.1-3.0m
+    
+##event_number_list = [2009,2045,2060,5716] #1964,#0.3-1.0m
+##event_number_list = [5890,6017,6052,6143] #5867,#0.3-1.0m
+    
+##event_number_list = [6262,6360,6358,6464] #6157,#0.3-1.0m
+##event_number_list = [6512,7741,7771,7813] #6475,#0.3-1.0m
+##event_number_list = [15442,15491,15505,16964]#7865,#0.3-1.0m
+##event_number_list = [16990,17006,17020,20766]#16968,#0.3-1.0m
+
+##event_number_list = [20993,22691,22784,22847] #20902,#0.3-1.0m
+##event_number_list = [24300,24356,24382,28082] #24247,#0.3-1.0m
+##event_number_list = [30207,30379,30532,30541] #28084,#0.3-1.0m
+##event_number_list = [31012,31173,31177,31361] #30858,#0.3-1.0m
+##event_number_list = [31536,31593,31964,31803]#31441,#0.3-1.0m
+
+# Re-run those that didn't work first time
+event_number_list = [15442,5890,6017,6052]
 
 

anuga_work/production/new_south_wales/batemans_bay/run_multiple_events.py

@@ -9 +9 @@
 for event_number in project.event_number_list:
     s = 'python2.5 run_model.py %d' % event_number
-    print s
+    #print s
     os.system(s)