Changeset 6596

Timestamp:

Mar 24, 2009, 2:14:11 PM (16 years ago)

Author:

jgriffin

Message:

changed event number

Location:

anuga_work/production/australia_ph2/sydney

Files:

• Arc_asc2raster_GDA94z56.py (modified) (2 diffs)
• export_results_max.py (modified) (5 diffs)
• get_timeseries.py (modified) (1 diff)
• project.py (modified) (2 diffs)

anuga_work/production/australia_ph2/sydney/Arc_asc2raster_GDA94z56.py

@@ -30 +30 @@
 
 time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
-time_dir2 = '20090311_145828_run_final_0_58152_2956_Bf_jgriffin'
-time_dir3 = '20090313_154712_run_final_0_58222_2956_Bf_jgriffin'
-time_dir4 = '20090313_154810_run_final_0_58326_2956_Bf_jgriffin'
-time_dir5 = '20090316_142137_run_final_0_58326_2956_Bf_jgriffin'
+##time_dir2 = '20090311_145828_run_final_0_58152_2956_Bf_jgriffin'
+##time_dir3 = '20090313_154712_run_final_0_58222_2956_Bf_jgriffin'
+##time_dir4 = '20090313_154810_run_final_0_58326_2956_Bf_jgriffin'
+##time_dir5 = '20090316_142137_run_final_0_58326_2956_Bf_jgriffin'
 ##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
 
 
-time_dirs = [time_dir3, time_dir4 , time_dir3, time_dir4, time_dir5]#, time_dir6]
+time_dirs = [time_dir1]#, time_dir4 , time_dir3, time_dir4, time_dir5]#, time_dir6]
 
 for time_dir in time_dirs:
@@ -48 +48 @@
     #ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_conditions_ocean1"
     
-    print 'Process: Create File GDB'
-    gp.CreateFileGDB_management(folder, "raster")
+##    print 'Process: Create File GDB'
+##    gp.CreateFileGDB_management(folder, "raster")
 
     gp.Workspace = raster_gbd

anuga_work/production/australia_ph2/sydney/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090316_142137_run_final_0_58326_2956_Bf_jgriffin'
+time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
 ##time_dir2 = '20090311_152929_run_final_0_70385_1203_Tb_jgriffin'
 ##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
@@ -28 +28 @@
 time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4]  
 
-cellsize = 250
+cellsize = 20
 ##cellsize = 5
 
-#timestep = None    # None means no timestep!
-timestep = 0
+timestep = None    # None means no timestep!
+#timestep = 0
 
 ######
@@ -52 +52 @@
 
 # one or more key strings from var_equations above
-var = ['elevation','stage']
+var = ['stage']
 
 ######
@@ -69 +69 @@
             northing_min = None
             northing_max = None
+##            easting_min = 290000
+##            easting_max = 316000
+##            northing_min = 6117000
+##            northing_max = 6170000
         else:
             try:
@@ -82 +86 @@
 
             name1 = join(directory, time_dir, project.scenario_name)
-##            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
+            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
 ##            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
 ##            name4 = join(directory, time_dir, project.scenario_name)+'_time_39960_0'

anuga_work/production/australia_ph2/sydney/get_timeseries.py

@@ -12 +12 @@
 """
 
-from os import sep
+from os import sep, rename, listdir, system
 from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
 import project
 
-timestamp1='20080815_103442_run_final_0.0_polyline_alpha0.1_kvanputt'
-timestamp2='20080815_103336_run_final_0.6_polyline_alpha0.1_kvanputt'
+directory = project.output_folder
 
-timestamps = [timestamp1, timestamp2]
-for timestamp in timestamps:
+##time_dir1 = '20081009_184721_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir2 = '20081031_120316_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir3 = '20081031_120515_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir4 = '20081031_120642_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir5 = '20081031_120753_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir6 = '20081111_113832_run_final_0.6_27255_alpha0.1_kvanputt'+sep+'gauge'
+##
+##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
 
-    filename1=project.output_dir+timestamp+sep+project.scenario_name+'.sww'
-    filename2=project.output_dir+timestamp+sep+'sww2'+sep+project.scenario_name+'_time_39900_0.sww'
+##time_dir1 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt'
+time_dir6 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
 
-    filenames = [filename1, filename2]
-    for filename in filenames:
+time_dirs = [time_dir6]
 
-        gauges = [project.gauges_dir_name, project.gauges_dir_name2]
-        for gauge in gauges:
-            
-            sww2csv_gauges(filename,gauge,
-                            quantities = ['stage','speed','depth','elevation'],
-                            verbose=True)
-                        
+for time_dir in time_dirs:
+    name = directory+sep+time_dir+sep+project.scenario_name 
+    gauge = project.gauges
+    sww2csv_gauges(name+'.sww',gauge,
+                   quantities = ['stage','speed','depth','elevation'],
+                   verbose=True)
+  
 
-    
-
-

anuga_work/production/australia_ph2/sydney/project.py

@@ -29 +29 @@
 tide = 0                # difference between MSL and HAT
 zone = 56              # specify zone of model
-event_number = 58326    # the event number or the mux file name
+event_number = 58152    # the event number or the mux file name
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
@@ -37 +37 @@
 # index is only used when wave = Tb
 index = 2956            # index from the PTHA
-wave = 'Bf'             # Bf (sts wave) Tb (index wave)
+wave = 'Tb'             # Bf (sts wave) Tb (index wave)
 
 setup = 'final'         # This can be one of three values