Changeset 4193

Timestamp:

Jan 24, 2007, 3:03:26 PM (17 years ago)

Author:

nick

Message:

updates

Location:

anuga_work/production/dampier_2006

Files:

• project.py (modified) (7 diffs)
• run_dampier.py (modified) (8 diffs)

anuga_work/production/dampier_2006/project.py

@@ -76 +76 @@
 tide_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'tide_data'+sep
 
-boundaries_source = '????'
+boundaries_source = ''
+boundaries_name = 'o_velocity'
 #boundaries locations
 boundaries_in_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep+'urs'+sep+boundaries_source+sep
-boundaries_in_dir_name = boundaries_in_dir + scenario_name
+boundaries_in_dir_name = boundaries_in_dir + boundaries_name
 boundaries_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep
-boundaries_dir_name = boundaries_dir + scenario_name
+boundaries_dir_name = boundaries_dir + boundaries_name
 #boundaries_time_dir = home+sep+state+sep+scenario+sep+'anuga'+sep+'boundaries'+sep+build_time+sep
 #boundaries_time_dir_name = boundaries_time_dir + boundaries_name  #Used by post processing
@@ -102 +103 @@
 
 refzone = 50 
-##south = degminsec2decimal_degrees(-20,55,0)
-##north = degminsec2decimal_degrees(-20,15,0)
-##west = degminsec2decimal_degrees(116,17,0)
-##east = degminsec2decimal_degrees(117,10,0)
+south_boundary = degminsec2decimal_degrees(-20,58,0)
+north_boundary = degminsec2decimal_degrees(-20,13,0)
+west_boundary = degminsec2decimal_degrees(116,15,0)
+east_boundary = degminsec2decimal_degrees(117,11,0)
 ##
 ##p0 = [south, degminsec2decimal_degrees(116,32,0)]
@@ -122 +123 @@
 print 'Area of bounding polygon', polygon_area(poly_all)/1000000.0
 
-res_poly_all = 100000
+#res_poly_all = 100000
+res_poly_all = 500000
 
 ###############################
@@ -132 +134 @@
 
 poly_dampier = read_polygon(polygons_dir+'dampier_town.csv')
-res_dampier = 500
+#res_dampier = 500
+res_dampier = 5000
 
 poly_karratha = read_polygon(polygons_dir+'karrathav2.csv')
@@ -138 +141 @@
 
 poly_karratha_town = read_polygon(polygons_dir+'karratha_townv2.csv')
-res_karratha_town = 500
+#res_karratha_town = 500
+res_karratha_town = 5000
+
+poly_facility = read_polygon(polygons_dir+'facility.csv')
+res_facility = 1000
 
 poly_delambre = read_polygon(polygons_dir+'delambre.csv')
-res_delambre = 1000
+res_delambre = 5000
 
 poly_coast = read_polygon(polygons_dir+'coastpoly.csv')
-res_coast = 1000
+#res_coast = 1000
+res_coast = 10000
 
 poly_NWislands = read_polygon(polygons_dir+'nw_islands_area.csv')
@@ -161 +169 @@
 res_island0 = res_poly_all
 
-res_islands = 5000
+#res_islands = 5000
+res_islands = 15000
 
 poly_ref_nw4 = read_polygon(polygons_dir+'ref_nw4.csv')
@@ -203 +212 @@
                     [poly_karratha,res_karratha],[poly_karratha_town,res_karratha_town],
                     [poly_delambre,res_delambre],[poly_coast,res_coast],
+                    [poly_facility,res_facility],
                     #[poly_NWislands,res_NWislands],
                     [poly_island0,res_island0],[poly_island1,res_island0],

anuga_work/production/dampier_2006/run_dampier.py

@@ -67 +67 @@
 # resolutions (maximal area of per triangle) for each polygon
 #--------------------------------------------------------------------------
+
+poly = [[0,0],[0,100],[100,100],[100,0]]
+
+create_mesh_from_regions(poly,
+                             boundary_tags={'back': [0], 'side': [1,3],
+                                            'ocean': [2]},
+                             maximum_triangle_area=1,
+                             interior_regions=None,
+                             filename=meshes_dir_name,
+                             use_cache=True,
+                             verbose=True)
+
+sys.exit()
 
 if myid == 0:
@@ -90 +103 @@
 print domain.statistics()
 
-"""
+
 print 'starting to create boundary conditions'
 boundaries_in_dir_name = project.boundaries_in_dir_name
@@ -102 +115 @@
     cache(urs2sww,
           (boundaries_in_dir_name,
-    #       boundaries_time_dir_name), 
            boundaries_dir_name), 
           {'verbose': True,
@@ -109 +121 @@
            '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(),
@@ -121 +129 @@
            )
 barrier()
-"""
+
 
 #-------------------------------------------------------------------------
@@ -131 +139 @@
 
     from polygon import *
+    #following sets the stage/water to be offcoast only
     IC = Polygon_function( [(project.poly_bathy, tide)], default = 0.)
     domain.set_quantity('stage', IC)
@@ -170 +179 @@
 print 'domain id', id(domain)
 print 'Reading Boundary file'
-#Bf = File_boundary(boundaries_dir_name + '.sww',
-#                  domain, time_thinning=5, use_cache=True, verbose=True)
+Bf = File_boundary(boundaries_dir_name + '.sww',
+                  domain, time_thinning=10, use_cache=True, verbose=True)
 
 print 'finished reading boundary file'
@@ -184 +193 @@
 domain.set_boundary({'back': Br,
                      'side': Bd,
-                     'ocean': Bd}) 
+                     'ocean': Bf}) 
 print'finish set boundary'