Changeset 6443

Timestamp:

Mar 4, 2009, 12:48:51 PM (14 years ago)

Author:

kristy

Message:

New boundary created through one csv time series file

Location:

anuga_work/production/australia_ph2/dampier

Files:

• export_results_max.py (modified) (1 diff)
• prepare_timeboundary.py (added)
• project.py (modified) (4 diffs)
• run_model.py (modified) (4 diffs)

anuga_work/production/australia_ph2/dampier/export_results_max.py

@@ -23 +23 @@
 
 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt'
-##time_dir2 = ''
-##time_dir3 = ''
-##time_dir4 = ''
-##time_dir5 = ''
-##time_dir6 = ''
-##
-time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]
-  
+time_dirs = [time_dir1]  
 
 cellsize = 250

anuga_work/production/australia_ph2/dampier/project.py

@@ -8 +8 @@
 from time import localtime, strftime, gmtime
 from os.path import join, exists
+from os import sep
 
 
@@ -29 +30 @@
 zone = 50               # specify zone of model
 event_number = 70028    # the event number or the mux file name
+index = 1203
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=1000          # final time for simulation
+finaltime=80000          # final time for simulation
 
 setup = 'trial'         # This can be one of three values
@@ -115 +117 @@
 mux_input_filename = 'event.list'
 multi_mux = True
+
+
 
 
@@ -215 +219 @@
 mux_input = join(event_folder, mux_input_filename)
 
+boundary_csv = event_folder + sep + 'sts_gauge_' + str(index) +'.csv'

anuga_work/production/australia_ph2/dampier/run_model.py

@@ -35 +35 @@
 from anuga.interface import Reflective_boundary
 from anuga.interface import Field_boundary
+from anuga.interface import Time_boundary
+from anuga.interface import file_function
+
 from anuga.interface import create_sts_boundary
 from anuga.interface import csv2building_polygons
@@ -47 +50 @@
 from setup_model import project
 import build_urs_boundary as bub
+import prepare_timeboundary as TB
 
 #-------------------------------------------------------------------------------
@@ -135 +139 @@
 print 'Set boundary - available tags:', domain.get_boundary_tags()
 
+# Prepare time boundary
+TB.prepare_timeboundary(project.boundary_csv)
+f = file_function(project.boundary_csv[:-4] + '.tms')
+
 Br = Reflective_boundary(domain)
 Bt = Transmissive_stage_zero_momentum_boundary(domain)
@@ -145 +153 @@
                     use_cache=True,
                     verbose=True)
+Tb = Time_boundary(domain,f)
 
 domain.set_boundary({'back': Br,
                      'side': Bd,
-                     'ocean': Bf}) 
+                     'ocean': Tb}) 
 
 #-------------------------------------------------------------------------------