Changeset 6577

Timestamp:

Mar 23, 2009, 10:52:41 AM (14 years ago)

Author:

jgriffin

Message:

different model runs; changed side boundary to transmissive

Location:

anuga_work/production/new_south_wales/batemans_bay

Files:

• Arc_asc2raster_GDA94z56.py (added)
• export_results_max.py (modified) (2 diffs)
• project.py (modified) (2 diffs)
• run_model.py (modified) (3 diffs)

anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

@@ -29 +29 @@
 ##time_dirs = [time_dir1, time_dir2]#, time_dir3]
 
-time_dir1 = '20090306_165253_run_final_0_58368_jgriffin'
+time_dir1 = '20090317_112148_run_final_0_58368_jgriffin'
 ##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
@@ -63 +63 @@
 
 # one or more key strings from var_equations above
-var = ['elevation', 'stage', 'speed','depth']
+var = ['stage','depth','speed','elevation']#, 'stage', 'speed','depth']
 
 ######

anuga_work/production/new_south_wales/batemans_bay/project.py

@@ -27 +27 @@
 # One or all can be changed each time the run_model script is executed
 tide = 0                # difference between MSL and HAT
-event_number = 51414    # the event number or the mux file name
+event_number = 58368    # the event number or the mux file name
 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=6000         # final time for simulation
 
 setup = 'final'         # This can be one of three values
@@ -78 +78 @@
 point_filenames = ['SD100031996_jgriffin_clip.csv',
                    'tomaga_offshore_AHD_MGA_1997.csv',
-                   'Batemans_BBHD_MGA_1995.csv']
+                   'Batemans_BBHD_MGA_1995.csv',
+                   'moruya_AHD_MGA_2000.csv']
           
 

anuga_work/production/new_south_wales/batemans_bay/run_model.py

@@ -141 +141 @@
                     domain, mean_stage=project.tide,
                     time_thinning=1,
-                    default_boundary=Bd,
+                    default_boundary=Dirichlet_boundary([0, 0, 0]),
                     boundary_polygon=bounding_polygon_sts,                    
                     use_cache=True,
@@ -147 +147 @@
 
 domain.set_boundary({'back': Br,
-                     'side': Bd,
+                     'side': Bt,
                      'ocean': Bf}) 
 
@@ -156 +156 @@
 t0 = time.time()
 for t in domain.evolve(yieldstep=project.yieldstep, 
-                       finaltime=project.finaltime,
+                       finaltime=1440,
                        skip_initial_step=False): 
     print domain.timestepping_statistics()
     print domain.boundary_statistics(tags='ocean')
 
+for t in domain.evolve(yieldstep=5, 
+                       finaltime=project.finaltime,
+                       skip_initial_step=True): 
+    print domain.timestepping_statistics()
+    print domain.boundary_statistics(tags='ocean')
+
 print 'Simulation took %.2f seconds' % (time.time()-t0)