Changeset 6254

Timestamp:

Jan 31, 2009, 4:43:13 AM (15 years ago)

Author:

ole

Message:

Established work area for the development of a modelling demo

Location:

anuga_work/production/busselton/standardised_version

Files:

• README.txt (added)
• build_busselton.py (modified) (1 diff)
• project.py (modified) (1 diff)
• run_busselton.py (modified) (2 diffs)

anuga_work/production/busselton/standardised_version/build_busselton.py

@@ -49 +49 @@
 # Fine pts file to be clipped to area of interest
 #-------------------------------------------------------------------------------
+
+# FIXME(Ole): Here's the Shark Bay way.
+
+## FIXME (Ole): Clip data by interior regions if possible
+#
+#print 'creating geospatial data objects from asc data (via dem and pts formats)'
+#for filename in project.ascii_grid_filenames:
+#    convert_dem_from_ascii2netcdf(filename,
+#                                  basename_out=filename,
+#                                  use_cache=True, verbose=True)
+#    dem2pts(filename, use_cache=True, verbose=True)#
+#
+#    geospatial_data += Geospatial_data(file_name=filename + '.pts',
+#                                       verbose=True)
+##
+#
+#print 'creating geospatial data objects from txt data'
+#for filename in project.point_filenames:
+#    geospatial_data += Geospatial_data(file_name=filename + '.txt',
+#                                       verbose=True)
+#
+#
+#print 'clip combined geospatial object by bounding polygon'
+#G = geospatial_data.clip(project.bounding_polygon)
+#
+#
+#print 'export combined geospatial data'
+#if access(project.topographies_dir, F_OK) == 0:
+#    mkdir (project.topographies_dir)
+#G.export_points_file(project.combined_dir_name + '.pts')
+
+
+
 print "project.poly_all", project.poly_all
 print "project.combined_dir_name", project.combined_dir_name

anuga_work/production/busselton/standardised_version/project.py

@@ -102 +102 @@
 #------------------------------------------------------------------------------
 
+# FIXME(Ole): Use lists of elevation data files as in Shark Bay study:
+## elevation data filenames
+#ascii_grid_filenames = ['10m_dem_without_survey', '50m_dem_without_10m_dem',
+#                        'bathysteeppt', 'bathyleft', 'bathyright']
+#point_filenames = ['field_survey_north', 'field_survey_south',
+#                   'clipped_bathymetry_final', 'coast_points_final']
+
+
+
 # elevation data used in build_busselton.py
 # onshore data: format ascii grid with accompanying projection file

anuga_work/production/busselton/standardised_version/run_busselton.py

@@ -54 +54 @@
 def run_model(**kwargs):
     
-    #------------------------------------------------------------------------------
+    #-----------------------------------------------------------------------
     # Copy scripts to time stamped output directory and capture screen
     # output to file
-    #------------------------------------------------------------------------------
+    #-----------------------------------------------------------------------
     print "Processor Name:",get_processor_name()
 
@@ -102 +102 @@
     # causes problems with the ability to cache set quantity which takes alot of times
         
+
+    # FIXME(Ole): Introduce create_domain_from_regions: Simpler and caches well
     print 'start create mesh from regions'