Changeset 6300

Timestamp:

Feb 10, 2009, 8:45:18 AM (15 years ago)

Author:

myall

Message:

 

Location:

anuga_work/production/australia_ph2/eucla_motel

Files:

• build_boundary_27255.py (modified) (5 diffs)
• build_elevation.py (modified) (1 diff)
• project.py (modified) (4 diffs)

anuga_work/production/australia_ph2/eucla_motel/build_boundary_27255.py

@@ -5 +5 @@
 Input: order_filename from project.py
        event_number needs to be reflected in the project file 
-Output: creates a sts file and csv files stored in project.boundaries_dir_event.
+Output: creates a sts file and csv files stored in project.event_folder.
 The run_busselton.py is reliant on the output of this script.
 """
@@ -95 +95 @@
 
     maxname = 'max_sts_stage.csv'
-    fid_max = open(project.boundaries_dir_event+sep+maxname,'w')
+    fid_max = open(project.event_folder+sep+maxname,'w')
     s = 'index, x, y, max_stage \n'
     fid_max.write(s)    
@@ -112 +112 @@
 
     minname = 'min_sts_stage.csv'
-    fid_min = open(project.boundaries_dir_event+sep+minname,'w')
+    fid_min = open(project.event_folder+sep+minname,'w')
     s = 'index, x, y, max_stage \n'
     fid_min.write(s)    
@@ -126 +126 @@
 
 
-        fid_sts = open(project.boundaries_dir_event+sep+basename+'_'+ str(index)+'.csv', 'w')
+        fid_sts = open(project.event_folder+sep+basename+'_'+ str(index)+'.csv', 'w')
         s = 'time, stage, xmomentum, ymomentum \n'
         fid_sts.write(s)
@@ -144 +144 @@
     return quantities,elevation,time
 
-quantities,elevation,time=get_sts_gauge_data(os.path.join(project.boundaries_dir_event,project.scenario_name),verbose=False)
+quantities,elevation,time=get_sts_gauge_data(project.event_sts,verbose=False)
 
 print len(elevation), len(quantities['stage'][0,:])

anuga_work/production/australia_ph2/eucla_motel/build_elevation.py

@@ -50 +50 @@
 # Create Geospatial data from ASCII files
 geospatial_data = {}
-for filename in project.ascii_grid_filenames:
-    absolute_filename = join(project.topographies_folder, filename)
-    convert_dem_from_ascii2netcdf(absolute_filename,
-                                  basename_out=absolute_filename,
-                                  use_cache=True,
-                                  verbose=True)
-    dem2pts(absolute_filename, use_cache=True, verbose=True)
-
-    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
-                                                verbose=True)
+##for filename in project.ascii_grid_filenames:
+##    absolute_filename = join(project.topographies_folder, filename)
+##    convert_dem_from_ascii2netcdf(absolute_filename,
+##                                  basename_out=absolute_filename,
+##                                  use_cache=True,
+##                                  verbose=True)
+##    dem2pts(absolute_filename, use_cache=True, verbose=True)
+##
+##    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
+##                                                verbose=True)
 
 ### Create Geospatial data from TXT files
-##for filename in project.point_filenames:
-##    absolute_filename = join(project.topographies_folder, filename)
-##    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
-##                                                verbose=True)
+for filename in project.point_filenames:
+    absolute_filename = join(project.topographies_folder, filename)
+    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
+                                                verbose=True)
 
 

anuga_work/production/australia_ph2/eucla_motel/project.py

@@ -32 +32 @@
 # Note, the user needs to set up the directory system accordingly
 state = 'australia_ph2'
-scenario_name = 'ceduna'
+scenario_name = 'eucla_motel'
 
 #------------------------------------------------------------------------------
@@ -44 +44 @@
 friction = 0.01           # manning's friction coefficient
 starttime = 0             
-finaltime = 80000         # final time for simulation
+finaltime = 1000 #80000         # final time for simulation
 
-setup = 'trial'  # Final can be replaced with trial or basic.
+setup = 'final'  # Final can be replaced with trial or basic.
                # Either will result in a coarser mesh that will allow a
                # faster, but less accurate, simulation.
@@ -81 +81 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = ['grid250m'] # 250m grid 2005
+##ascii_grid_filenames = ['grid250m'] # 250m grid 2005
 
 # Format for point is x,y,elevation (with header)
-##point_filenames = ['Busselton_Contour0.txt',     # Coastline
-##                   'Busselton_BeachSurvey.txt',    # Beach survey
-##                   'Busselton_NavyFinal.txt',  # Bathymetry
-##                   'Busselton_Chart.txt', # Bathymetry Charts
-##                   'Busselton_Digitised.txt', # Digitised Fairsheet
-##                   'Busselton_250m.txt', # 250m
-##                   'Bunbury_TIN.txt', # Bunbury aoi TIN'd in ArcGIS
-##                   'Busselton_TIN.txt', # Busselton aoi TIN'd in ArcGIS
-##                   'XYAHD_clip.txt'] # To extend boundary
+point_filenames = ['grid250m_pts.txt'] # To extend boundary
 
 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
@@ -156 +148 @@
 event_sts = join(boundaries_folder, str(event_number), scenario_name)
 
+# The absolute pathname of the event folder
+event_folder = join(boundaries_folder, str(event_number))
+
 # The absolute pathname for the output folder names
 # Used for build_elevation.py