Changeset 4429 for anuga_work/production/broome_2006/run_broome.py

Timestamp:

May 15, 2007, 4:13:58 PM (17 years ago)

Author:

nick

Message:

update broome

File:

• anuga_work/production/broome_2006/run_broome.py (modified) (6 diffs)

anuga_work/production/broome_2006/run_broome.py

@@ -81 +81 @@
                          interior_regions=project.interior_regions,
                          filename=meshes_dir_name,
-                         use_cache=True,
+                         use_cache=False,
                          verbose=True)
 
@@ -92 +92 @@
 print 'Setup computational domain'
 #domain = Domain(meshes_time_dir_name, use_cache=True, verbose=True)
-domain = Domain(meshes_dir_name, use_cache=True, verbose=True)
+domain = Domain(meshes_dir_name, use_cache=False, verbose=True)
 print domain.statistics()
-'''
-print 'starting to create boundary conditions'
-boundaries_in_dir_name = project.boundaries_in_dir_name
-
-from anuga.shallow_water.data_manager import urs2sww, ferret2sww
-
-# put above distribute
-print 'boundary file is: ',boundaries_dir_name
-from caching import cache
-if myid == 0:
-    cache(ferret2sww,
-          (boundaries_in_dir_name,
-    #       boundaries_time_dir_name), 
-           boundaries_dir_name), 
-          {'verbose': True,
-           'minlat': project.south_boundary,
-           'maxlat': project.north_boundary,
-           'minlon': project.west_boundary,
-           'maxlon': project.east_boundary,
-#           'minlat': project.south,
-#           'maxlat': project.north,
-#           'minlon': project.west,
-#           'maxlon': project.east,
-           'mint': 0, 'maxt': 35100,
-           'origin': domain.geo_reference.get_origin(),
-           'mean_stage': project.tide,
-#           'zscale': 1,                 #Enhance tsunami
-           'fail_on_NaN': False},
-           verbose = True,
-           )
-barrier()
-'''
 
 #-------------------------------------------------------------------------
@@ -134 +102 @@
     print 'Setup initial conditions'
 
-    domain.set_quantity('stage', tide)
-    domain.set_quantity('friction', 0.01) 
+    from polygon import Polygon_function
+    #following sets the stage/water to be offcoast only
+    IC = Polygon_function( [(project.poly_mainland, -1.0)], default = tide,
+                             geo_reference = domain.geo_reference)
+    domain.set_quantity('stage', IC)
     #combined_time_dir_name = project.topographies_dir+build_time+project.combined_name
     print 'Start Set quantity'
     print 'project.combined_dir_name_unclipped1',project.combined_dir_name_unclipped1+'.txt'
     domain.set_quantity('elevation', 
+#                    filename = project.combined_dir_name+'.txt',
                     filename = project.combined_dir_name_unclipped1+'.txt',
 #                    filename = project.combined_smaller_name_dir+'.xya',
@@ -175 +147 @@
 print 'domain id', id(domain)
 #Bf = File_boundary(boundaries_dir_name + '.sww',
-#                  domain, time_thinning=5, use_cache=True, verbose=True)
+#                  domain, time_thinning=12, use_cache=True, verbose=True)
+Bf = Field_boundary(boundaries_dir_name + '.sww',
+                  domain, time_thinning=1, mean_stage=tide, 
+                  use_cache=True, verbose=True)
 
 print 'finished reading boundary file'
@@ -185 +160 @@
 domain.set_boundary({'back': Br,
                      'side': Bd,
-                     'ocean': Bd}) 
+                     'ocean': Bf}) 
+#                     'ocean': Bd}) 
 print'finish set boundary'
 
@@ -194 +170 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 60, finaltime = 10000):
+for t in domain.evolve(yieldstep = 60, finaltime = 30000):
     domain.write_time()
     domain.write_boundary_statistics(tags = 'ocean')
-    if allclose(t, 120):
-        domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bo})
+#    if allclose(t, 120):
+#        domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bo})
 
-    if allclose(t, 1020):
-        domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})
+#    if allclose(t, 1020):
+#        domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})