Changeset 6021 for anuga_work/production/carnarvon/project.py

Timestamp:

Nov 28, 2008, 11:38:37 AM (15 years ago)

Author:

kristy

Message:

addition of more elevation data

File:

• anuga_work/production/carnarvon/project.py (modified) (4 diffs)

anuga_work/production/carnarvon/project.py

@@ -45 +45 @@
 tide = 0                #0.6
 #event_number = 27255 # Java 9.3 worst case for Perth
-event_number = 68693 # Sumatra 9.2
-#event_number = 27283  # Java 9.3 original
+#event_number = 68693 # Sumatra 9.2
+event_number = 27283  # Java 9.3 original
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
 starttime=0             
-finaltime=1000         # final time for simulation
+finaltime=80000         # final time for simulation
 
 setup='final'  # Final can be replaced with trial or basic.
@@ -86 +86 @@
 # onshore data: format ascii grid with accompanying projection file
 onshore_name = 'dem_onshore' 
-# coastline: format x,y,elevation (with title)
-coast_name = 'coastline.txt'
 # bathymetry: format x,y,elevation (with title)
-offshore_name = 'Shark_Bay_Clip.txt'
-offshore_name1 = 'XYAHD_Clip.txt'
-offshore_name2 = 'DPI_data.txt'
+offshore_name = 'carnarvon_data.txt'
+offshore_name1 = 'makeup.txt'
+#grids made for aoi and aos as anuga interpretation was bad
+topo_aoi = 'topo_20m_aoi'
+topo_aos = 'topo_50m_aos'
 
 # gauges - used in get_timeseries.py
@@ -131 +131 @@
 # where the input data sits
 onshore_in_dir_name = topographies_in_dir + onshore_name
-coast_in_dir_name = topographies_in_dir + coast_name
 offshore_in_dir_name = topographies_in_dir + offshore_name
 offshore_in_dir_name1 = topographies_in_dir + offshore_name1
-offshore_in_dir_name2 = topographies_in_dir + offshore_name2
+offshore_in_dir_aoi = topographies_in_dir + topo_aoi
+offshore_in_dir_aos = topographies_in_dir + topo_aos
 
 # where the output data sits
 onshore_dir_name = topographies_dir + onshore_name
-coast_dir_name = topographies_dir + coast_name
 offshore_dir_name = topographies_dir + offshore_name
 offshore_dir_name1 = topographies_dir + offshore_name1
-offshore_dir_name2 = topographies_dir + offshore_name2
+offshore_dir_aoi = topographies_dir + topo_aoi
+offshore_dir_aos = topographies_dir + topo_aos
 
 # where the combined elevation file sits
@@ -179 +179 @@
 poly_island2=read_polygon(polygons_dir +'land_initial_condition_dorne.csv')
 poly_ocean=read_polygon(polygons_dir +'ocean_initial_condition.csv')
+poly_aos_extract = read_polygon(polygons_dir+'extract_aos.csv')
 
 # Initial bounding polygon for data clipping