Changeset 5065 for anuga_validation

Timestamp:

Feb 20, 2008, 3:12:21 PM (17 years ago)

Author:

nick

Message:

update run_circular.py with exporting gauges

File:

• anuga_validation/circular_island/run_circular.py (modified) (8 diffs)

anuga_validation/circular_island/run_circular.py

@@ -12 +12 @@
 from anuga.shallow_water import Transmissive_Momentum_Set_Stage_boundary
 from anuga.shallow_water import Dirichlet_boundary, Time_boundary
-from anuga.shallow_water.data_manager import get_maximum_inundation_data, start_screen_catcher
-from anuga.abstract_2d_finite_volumes.util import file_function, sww2csv_gauges
+from anuga.shallow_water.data_manager import get_maximum_inundation_data, start_screen_catcher, copy_code_files
+from anuga.abstract_2d_finite_volumes.util import file_function, sww2csv_gauges,csv2timeseries_graphs
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 from anuga.utilities.polygon import read_polygon#, plot_polygons
@@ -19 +19 @@
 from Numeric import array, zeros, Float, allclose,resize
 from anuga.shallow_water.data_manager import csv2dict
-
+from time import localtime, strftime, gmtime
+from anuga.utilities.system_tools import get_user_name, get_host_name
+from os import sep, environ, getenv
 #-------------------------
 # Create Domain from mesh
 #-------------------------
 
-#start_screen_catcher()
+user = get_user_name()
+home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir  
+time = strftime('%Y%m%d_%H%M%S',localtime())
+
+res=0.01
+
+dir_comment = time+'_res_'+str(res)+'_'+str(user)
+
+anuga_dir = home+'anuga_validation'+sep+'circular_island_tsunami_benchmark'+sep+'anuga'+sep
+output_dir = anuga_dir+'outputs'+sep+dir_comment+sep
+copy_code_files(output_dir,__file__)
+
+start_screen_catcher(output_dir)
 
 def get_xy(x=0,y=0,r=1,angle=1):
@@ -91 +105 @@
 
 
-
+# angle to make circular interior polygons
 angles1 = (0.0, 22.5, 45.0, 67.5, 90.0, 112.5, 135.0, 157.5, 180.0, 202.5, 225.0,
           247.5, 270.0, 292.5, 315.0, 337.5)
@@ -98 +112 @@
 center_y = 13.80
 
-d=0.75
-res=0.05
-res=1.
+#d=0.75
 #L=3.6 and therefore the gauges 1,2,3 and 4 are located 1.8 from the toe of the island
 # 13.8 - 3.6(radii of island) - 1.8 = 8.4
@@ -202 +214 @@
 model_name='wavetank_model'
 domain.set_name(model_name)  # Name of output sww file
+domain.set_datadir(output_dir)
 domain.set_default_order(1)               # Apply second order scheme 
 #domain.set_all_limiters(0.9)              # Max second order scheme (old lim)
@@ -208 +221 @@
 domain.tight_slope_limiters = 1
 domain.beta_h = 0.0
+sww_file=output_dir+sep+model_name+'.sww'
 
 
@@ -261 +275 @@
 r2 = 2.6
 #degrees to check
-d=2
+d=1
 
 poly_segment = []
@@ -276 +290 @@
     #print i,poly[i]
     
-    run_up, location = get_maximum_inundation_data(filename=model_name+'.sww',polygon=poly_segment, verbose=False)
+    run_up, location = get_maximum_inundation_data(filename=sww_file,polygon=poly_segment, verbose=False)
     
     print 'maxd %.3f, %.6f, %.6f, %.5s, %.5s '%(((angle/pi)*180)+180, run_up*100, (run_up-water_height)*100, location[0],location[1]) 
 
-
-
-
-
-
-
+sww_file=anuga_dir+'outputs'+sep+'20080220_010628_res_0.05_nbartzis'+sep+'wavetank_model.sww'
+sww_file=anuga_dir+'outputs'+sep+'20080219_233611_res_0.05_nbartzis'+sep+'wavetank_model.sww'
+sww2csv_gauges(sww_file=sww_file,
+                   gauge_file='gauges.csv',
+                   quantities = ['stage','depth', 'elevation', 'xmomentum', 'ymomentum'],
+                   verbose=True,
+                   use_cache = True)
+
+csv2timeseries_graphs(directories_dic={anuga_dir+sep+'boundaries'+sep:['wavetank',0, 0],
+                                       anuga_dir+'outputs'+sep+'20080219_233611_res_0.05_nbartzis'+sep:['Fixed Wave',0,0]},
+                            output_dir=anuga_dir+'outputs'+sep+'20080219_233611_res_0.05_nbartzis'+sep,
+                            base_name='gauge_',
+                            plot_numbers='',
+                            quantities=['stage'],
+                            extra_plot_name='',
+                            assess_all_csv_files=True,                            
+                            create_latex=False,
+                            verbose=True)
+
+
+
+
+