Changeset 4359

Timestamp:

Apr 4, 2007, 12:44:33 PM (18 years ago)

Author:

nick

Message:

update perth

File:

• anuga_work/production/perth_2006/run_perth.py (modified) (6 diffs)

anuga_work/production/perth_2006/run_perth.py

@@ -53 +53 @@
 
 tide = project.tide
-
+print 'hello'
 # creates copy of code in output dir
 if myid == 0:
@@ -63 +63 @@
 print 'min triangles', project.trigs_min,
 print 'Note: This is generally about 20% less than the final amount'
-
+import sys
+sys.exit()
 #--------------------------------------------------------------------------
 # Create the triangular mesh based on overall clipping polygon with a
@@ -72 +73 @@
 
 if myid == 0:
-
-    print 'start create mesh from regions'
-#    interior_regions = [[project.poly_pos20_neg20, project.res_pos20_neg20],
-#                    [project.poly_cbd, project.res_cbd],
-#                    [project.poly_penguin, project.res_penguin]]    
-#                    [project.poly_interior, 1000]]       
-
-#    from anuga.utilities.polygon import plot_polygons
-#    figname = project.output_run_time_dir + 'poly_pic'
-#    plot_polygons([project.poly_coast,project.poly_pipeline,project.poly_facility,
-#               project.bounding_polygon],
-#               figname,
-#               verbose = True)
-#    import sys; sys.exit()
 
     print 'start create mesh from regions'
@@ -108 +95 @@
 domain = Domain(meshes_dir_name, use_cache=True, 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()
-'''
 
 #-------------------------------------------------------------------------
@@ -154 +109 @@
 
     domain.set_quantity('elevation', 
-#                    filename = project.combined_dir_name+'.txt',
-                    filename = project.combined_smaller_name_dir+'.txt',
+                    filename = project.combined_dir_name+'.txt',
+#                    filename = project.combined_smaller_name_dir+'.txt',
                     use_cache = True,
                     verbose = True,
@@ -176 +131 @@
 domain.set_datadir(project.output_run_time_dir)
 domain.set_default_order(2) # Apply second order scheme
-domain.set_minimum_storable_height(0.01) # Don't store anything less than 1cm
+domain.set_minimum_storable_height(0.001) # Don't store anything less than 0.1cm
+#domain.set_minimum_storable_height(0.01) # Don't store anything less than 1cm
+#domain.set_minimum_storable_height(0.1) # Don't store anything less than 10cm
 domain.set_store_vertices_uniquely(False)
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])