Changeset 6895

Timestamp:

Apr 24, 2009, 10:58:04 AM (16 years ago)

Author:

jgriffin

Message:

 

Location:

anuga_work/production/australia_ph2/bamaga

Files:

• build_elevation.py (modified) (3 diffs)
• build_urs_boundary.py (modified) (1 diff)
• project.py (modified) (5 diffs)
• run_model.py (modified) (1 diff)

anuga_work/production/australia_ph2/bamaga/build_elevation.py

@@ -59 +59 @@
         dem2pts(absolute_filename, use_cache=True, verbose=True)
 
-        geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
+        G_grid = Geospatial_data(file_name=absolute_filename+'.pts',
                                                     verbose=True)
+        print 'Clip geospatial object'
+        geospatial_data[filename] = G_grid.clip(project.bounding_polygon)
 
 # Create Geospatial data from TXT files
@@ -66 +68 @@
     for filename in project.point_filenames:
         absolute_filename = join(project.topographies_folder, filename)
-        geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
+        G_points = Geospatial_data(file_name=absolute_filename,
                                                     verbose=True)
+        print 'Clip geospatial object'
+        geospatial_data[filename] = G_points.clip(project.bounding_polygon)
 
 #-------------------------------------------------------------------------------
@@ -82 +86 @@
     G += geospatial_data[key]
         
-print 'Clip combined geospatial data'
+#print 'Clip combined geospatial data'
 ##G_clip = G.clip_outside(project.poly_aoi1)
 ##G_all = G_clip + geospatial_data[project.point_filenames[4]] 
 #G_clipped = G_all.clip(project.poly_all)
-G_clip = G.clip(project.bounding_polygon)
+#G_clip = G.clip(project.bounding_polygon)
 
 
 print 'Export combined DEM file'
-G_clip.export_points_file(project.combined_elevation + '.pts')
+G.export_points_file(project.combined_elevation + '.pts')
 print 'Do txt version too'
 # Use for comparision in ARC
-G_clip.export_points_file(project.combined_elevation + '.txt')
+G.export_points_file(project.combined_elevation + '.txt')
 

anuga_work/production/australia_ph2/bamaga/build_urs_boundary.py

@@ -139 +139 @@
                 ordering_filename=project.urs_order,
                 weights=mux_weights,
-                central_meridian=project.central_meridian,
+                #central_meridian=project.central_meridian,
                 zone=project.zone,
                 verbose=True)

anuga_work/production/australia_ph2/bamaga/project.py

@@ -17 +17 @@
 # Note, the user needs to set up the directory system accordingly
 state = 'australia_ph2'
-scenario_name = 'bamaga'
+scenario_name = 'gulf'
 
 #-------------------------------------------------------------------------------
@@ -26 +26 @@
 # One or all can be changed each time the run_model script is executed
 
-central_meridian = 143.5 # Central meridian for projection (optional)
-zone = None
-import sys
-if len(sys.argv) > 1:
-    event_number = int(sys.argv[1])
-else:    
-    event_number = 31975    # the event number or the mux file name
-
-event_number_list = [31975, 51390, 63735] # To piggy back multiple events
+#central_meridian = 143.5 # Central meridian for projection (optional)
+zone = 54
+##import sys
+##if len(sys.argv) > 1:
+##    event_number = int(sys.argv[1])
+##else:
+
+event_number = 64977    # the event number or the mux file name
+
+#event_number_list = [31975, 51390, 63735] # To piggy back multiple events
 
 tide = 0                # difference between MSL and HAT
@@ -40 +41 @@
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=60000 # final time for simulation
+finaltime=80000 # final time for simulation
 setup = 'final'         # This can be one of three values
                         #    trial - coarsest mesh, fast
@@ -47 +48 @@
 
 # index is only used when wave = Tb
-index = 2000            # index from the PTHA - Y2000 0.257m
+index = 1914            # index from the PTHA - Y2000 0.257m
 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
 
@@ -78 +79 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = ['grid250m_pro'] # 250m grid 2005
+ascii_grid_filenames = ['gulf_bg'] # 250m grid 2005
 
 # Format for points is x,y,elevation (with header)

anuga_work/production/australia_ph2/bamaga/run_model.py

@@ -92 +92 @@
 # Boundary tags refer to project.landward_boundary
 # 4 points equals 5 segments start at N
-boundary_tags={'back': range(num_ocean_segments+1,
-                             num_ocean_segments+num_land_points),
-               'side': [num_ocean_segments,
-                        num_ocean_segments+num_land_points],
-               'ocean': range(num_ocean_segments)}
+##boundary_tags={'back': range(num_ocean_segments+1,
+##                             num_ocean_segments+num_land_points),
+##               'side': [num_ocean_segments,
+##                        num_ocean_segments+num_land_points],
+##               'ocean': range(num_ocean_segments)}
+boundary_tags = {'back':[3],
+                 'side':[2,4],
+                 'ocean':[0,1,5]}
+
+print 'boundary tags', boundary_tags
 
 # Build mesh and domain