Changeset 6027 for anuga_work/production/geraldton/project.py

Timestamp:

Nov 28, 2008, 11:48:45 AM (15 years ago)

Author:

kristy

Message:

general maintenance of scripts, added boundary to mesh in run_geraldton.py

File:

• anuga_work/production/geraldton/project.py (modified) (7 diffs)

anuga_work/production/geraldton/project.py

@@ -43 +43 @@
 # Model specific parameters. One or all can be changed each time the
 # run_scenario script is executed
-tide = 0.6                #0.6
+tide = 0.6         #0.6
 #event_number = 27255 # Java 9.3 worst case for Perth
 event_number = 68693 # Sumatra 9.2
@@ -50 +50 @@
 friction=0.01           # manning's friction coefficient
 starttime=0             
-finaltime=80000         # final time for simulation
+finaltime=80000       # final time for simulation
 
 setup='final'  # Final can be replaced with trial or basic.
@@ -56 +56 @@
                # faster, but less accurate, simulation.
 
+if setup =='coarse':
+    print'coarse'
+    res_factor=10000000
+    time_thinning=48
+    yieldstep=240
 if setup =='trial':
     print'trial'
@@ -100 +105 @@
 gauge_name2 = 'thinned_MGA50.csv'
 
+#buildings - used in run_building_inundation.py
+building = 'geraldton_res_clip'
+
+# barrier
+barrier = 'CBD.csv'
+barrier1 = 'wall.csv'
+
 # BOUNDING POLYGON - used in build_boundary.py and run_geraldton.py respectively
 # NOTE: when files are put together the points must be in sequence - for ease go clockwise!
@@ -156 +168 @@
 # where the mesh sits (this is created during the run_geraldton.py)
 meshes_dir_name = meshes_dir + scenario_name+'.msh'
+barrier_dir_name = meshes_dir + barrier
+barrier_dir_name1 = meshes_dir + barrier1
 
 # where the boundary ordering files sit (this is used within build_boundary.py)
@@ -172 +186 @@
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
 
-#w here the directory of the gauges sit
+# where the directory of the gauges sit
 gauges_dir_name = gauges_dir + gauge_name       #used for get_timeseries.py
 gauges_dir_name2 = gauges_dir + gauge_name2     #used for get_timeseries.py
+building_in_dir_name = gauges_dir + building + '.csv'    #used for run_building_inundation.py
 
 #------------------------------------------------------------------------------
@@ -184 +199 @@
 poly_mainland=read_polygon(polygons_dir +'land_initial_condition.csv')
 poly_ocean=read_polygon(polygons_dir +'ocean_initial_condition.csv')
+poly_buffer_20m=read_polygon(polygons_dir +'buffer_20m.csv')
 
 # Initial bounding polygon for data clipping