Changeset 6342

Timestamp:

Feb 16, 2009, 8:56:21 AM (16 years ago)

Author:

myall

Message:

 

Location:

anuga_work/production

Files:

• australia_ph2/ceduna/build_boundary_27255.py (modified) (5 diffs)
• australia_ph2/ceduna/build_boundary_64468.py (added)
• australia_ph2/ceduna/file_length.py (added)
• australia_ph2/ceduna/project.py (modified) (2 diffs)
• australia_ph2/ceduna/run_model.py (modified) (6 diffs)
• australia_ph2/eucla_motel/run_model.py (modified) (1 diff)
• busselton/get_timeseries.py (modified) (1 diff)
• busselton/project.py (modified) (1 diff)
• geraldton/get_timeseries.py (modified) (1 diff)
• geraldton/project.py (modified) (1 diff)
• perth/get_timeseries.py (modified) (1 diff)
• perth/project.py (modified) (1 diff)

anuga_work/production/australia_ph2/ceduna/build_boundary_27255.py

@@ -5 +5 @@
 Input: order_filename from project.py
        event_number needs to be reflected in the project file 
-Output: creates a sts file and csv files stored in project.boundaries_dir_event.
+Output: creates a sts file and csv files stored in project.event_folder.
 The run_busselton.py is reliant on the output of this script.
 """
@@ -95 +95 @@
 
     maxname = 'max_sts_stage.csv'
-    fid_max = open(project.boundaries_dir_event+sep+maxname,'w')
+    fid_max = open(project.event_folder+sep+maxname,'w')
     s = 'index, x, y, max_stage \n'
     fid_max.write(s)    
@@ -112 +112 @@
 
     minname = 'min_sts_stage.csv'
-    fid_min = open(project.boundaries_dir_event+sep+minname,'w')
+    fid_min = open(project.event_folder+sep+minname,'w')
     s = 'index, x, y, max_stage \n'
     fid_min.write(s)    
@@ -126 +126 @@
 
 
-        fid_sts = open(project.boundaries_dir_event+sep+basename+'_'+ str(index)+'.csv', 'w')
+        fid_sts = open(project.event_folder+sep+basename+'_'+ str(index)+'.csv', 'w')
         s = 'time, stage, xmomentum, ymomentum \n'
         fid_sts.write(s)
@@ -144 +144 @@
     return quantities,elevation,time
 
-quantities,elevation,time=get_sts_gauge_data(os.path.join(project.boundaries_dir_event,project.scenario_name),verbose=False)
+quantities,elevation,time=get_sts_gauge_data(project.event_sts,verbose=False)
 
 print len(elevation), len(quantities['stage'][0,:])

anuga_work/production/australia_ph2/ceduna/project.py

@@ -44 +44 @@
 friction = 0.01           # manning's friction coefficient
 starttime = 0             
-finaltime = 80000         # final time for simulation
+finaltime = 1000 #80000         # final time for simulation
 
-setup = 'trial'  # Final can be replaced with trial or basic.
+setup = 'final'  # Final can be replaced with trial or basic.
                # Either will result in a coarser mesh that will allow a
                # faster, but less accurate, simulation.
@@ -156 +156 @@
 event_sts = join(boundaries_folder, str(event_number), scenario_name)
 
+# The absolute pathname of the event folder
+event_folder = join(boundaries_folder, str(event_number))
+
 # The absolute pathname for the output folder names
 # Used for build_elevation.py

anuga_work/production/australia_ph2/ceduna/run_model.py

@@ -26 +26 @@
 # Related major packages
 from anuga.interface import create_domain_from_regions
-from anuga.interface import Transmissive_stage_zero_momentum_boundary
 from anuga.interface import Dirichlet_boundary
 from anuga.interface import Reflective_boundary
@@ -32 +31 @@
 from anuga.interface import create_sts_boundary
 from anuga.interface import csv2building_polygons
+from file_length import file_length
 
 from anuga.shallow_water.data_manager import start_screen_catcher
@@ -70 +70 @@
 # Number of boundary segments
 N = len(event_sts)-1
+# Number of landward_boundary points
+M = file_length(project.landward_boundary)
 
 # Boundary tags refer to project.landward_boundary
 # 4 points equals 5 segments start at N
-boundary_tags={'back': [N+1, N+2, N+3, N+4, N+5, N+6, N+7],
-               'side': [N, N+8],
+boundary_tags={'back': range(N+1,N+M),
+               'side': [N,N+M],
                'ocean': range(N)}
 
@@ -101 +103 @@
                       default=project.tide,
                       geo_reference=domain.geo_reference)
-domain.set_quantity('stage', IC, use_cache=True, verbose=True)
+domain.set_quantity('stage', 0, use_cache=True, verbose=True)
 domain.set_quantity('friction', project.friction) 
 domain.set_quantity('elevation', 
@@ -116 +118 @@
 
 Br = Reflective_boundary(domain)
-Bt = Transmissive_stage_zero_momentum_boundary(domain)
+Bd = Dirichlet_boundary([project.tide,0,0])
 Bf = Field_boundary(project.event_sts+'.sts',
                     domain, mean_stage=project.tide,
@@ -125 +127 @@
                     verbose=True)
 
+
 domain.set_boundary({'back': Br,
-                     'side': Bt,
+                     'side': Bd,
                      'ocean': Bf}) 
 

anuga_work/production/australia_ph2/eucla_motel/run_model.py

@@ -103 +103 @@
                       default=project.tide,
                       geo_reference=domain.geo_reference)
-domain.set_quantity('stage', IC, use_cache=True, verbose=True)
+domain.set_quantity('stage', 0, use_cache=True, verbose=True)
 domain.set_quantity('friction', project.friction) 
 domain.set_quantity('elevation', 

anuga_work/production/busselton/get_timeseries.py

@@ -27 +27 @@
 ##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
 
-time_dir1 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt'
+##time_dir1 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt'
+time_dir6 = '20081111_113832_run_final_0.6_27255_alpha0.1_kvanputt'
 
-time_dirs = [time_dir1]
+time_dirs = [time_dir6]
 
 for time_dir in time_dirs:

anuga_work/production/busselton/project.py

@@ -118 +118 @@
 # gauges - used in get_timeseries.py
 #gauge_name = scenario_name+'.txt'
-gauge_name = 'Gauges.csv'
+#gauge_name = 'Gauges.csv'
+gauge_name = 'MH_gauges.csv'
 # buildings - used in run_building_inundation.py
 building = 'busselton_res_clip'

anuga_work/production/geraldton/get_timeseries.py

@@ -18 +18 @@
 directory = project.output_dir
 
-time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
-time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
-time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
-time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'+sep+'gauge'
-time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
-time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
+time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'
+##time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
 
-time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir4] #[time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
 
 for time_dir in time_dirs:

anuga_work/production/geraldton/project.py

@@ -102 +102 @@
 
 # gauges - used in get_timeseries.py
-gauge_name = 'geraldton.csv'
+gauge_name = 'MH_gauges.csv'  #'geraldton.csv'
 gauge_name2 = 'thinned_MGA50.csv'
 

anuga_work/production/perth/get_timeseries.py

@@ -18 +18 @@
 directory = project.output_dir
 
-time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
-time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
-time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
-time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'+sep+'gauge'
-time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
-time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge'
+time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
+##time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge'
+##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge'
 
-time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir4] #[time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
 
 for time_dir in time_dirs:

anuga_work/production/perth/project.py

@@ -100 +100 @@
 
 # gauges - used in get_timeseries.py
-gauge_name = 'perth.csv'
+##gauge_name = 'perth.csv'
+gauge_name = 'MH_gauges.csv'
 
 #buildings - used in run_building_inundation.py