Changeset 3851 for anuga_work/production/dampier_2006/build_dampier.py

Timestamp:

Oct 25, 2006, 3:19:29 PM (17 years ago)

Author:

nick

Message:

build_dampier.py is working
run_dampier.py should work... it is still running
project.py is working

File:

• anuga_work/production/dampier_2006/build_dampier.py (modified) (6 diffs)

anuga_work/production/dampier_2006/build_dampier.py

@@ -30 +30 @@
 from anuga.shallow_water import Reflective_boundary
 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
-
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
-
 from anuga.geospatial_data.geospatial_data import *
 
@@ -38 +36 @@
 import project   # Definition of file names and polygons
 
-
-
 #------------------------------------------------------------------------------
 # Copy scripts to time stamped output directory and capture screen
 # output to file
 #------------------------------------------------------------------------------
+#print 'project.build_time', project.build_time
+#print 'project.output_build_time_dir',project.output_build_time_dir
+#print '__file__, project.output_time_dir + basename(__file__)', __file__, project.output_build_time_dir + basename(__file__)
+#print 'dirname(project.__file__) +sep+ project.__name__', dirname(project.__file__) +sep+ project.__name__+'.py'
+#print'project.topographies_time_dir',project.topographies_time_dir, project.combined_time_dir_name + '.pts'
+#print project.boundaries_time_dir, project.boundaries_time_dir_name
+#import sys; sys.exit()
+if access(project.output_build_time_dir,F_OK) == 0:
+    mkdir (project.output_build_time_dir)
+copy (dirname(project.__file__) +sep+ project.__name__+'.py',
+      project.output_build_time_dir + project.__name__+'.py')     #copies project.py
+copy (__file__, project.output_build_time_dir + basename(__file__))
+print 'files '+ project.__name__+'.py and '+ basename(__file__)+' copied to '+ project.output_build_time_dir       #copies this file
+#import sys; sys.exit()
 
-
-
-# creates copy of code in output dir if dir doesn't exist
-if access(project.output_time_dir,F_OK) == 0:
-    mkdir (project.output_time_dir)
-copy (dirname(project.__file__) +sep+ project.__name__+'.py',
-      project.output_time_dir + project.__name__+'.py')     #copies project.py
-copy (__file__, project.output_time_dir + basename(__file__))
-print 'project.output_time_dir',project.output_time_dir       #copies this file
-
-
-
-
-boundary_dir_name = project.boundary_dir_name
-
-
-from anuga.shallow_water.data_manager import urs2sww
-
-print 'minlat=project.north_boundary, maxlat=project.south_boundary',project.north_boundary, project.south_boundary
-print 'minlon= project.west_boundary, maxlon=project.east_boundary',project.west_boundary, project.east_boundary
-
-#import sys; sys.exit()
-urs2sww(boundary_dir_name,
-#        minlat=project.south, maxlat=project.north,
-#        minlon= project.west, maxlon=project.east,
-#        mint=0, maxt= 35000,
-        verbose='true')
-        
-import sys; sys.exit()
 #-------------------------------------------------------------------------------
 # Preparation of topographic data
@@ -102 +83 @@
 offshore_dir_name14 = project.offshore_dir_name14
 
-
-
-# files to be used
-#file_used = []
-
 # creates DEM from asc data
 convert_dem_from_ascii2netcdf(onshore_dir_name, use_cache=True, verbose=True)
@@ -120 +96 @@
         verbose=True)
 
-
 #creates pts file for islands DEM
 dem2pts(islands_dir_name, use_cache=True, verbose=True)
 
-
-print'create G1'
+print'create Geospatial data objects from topographies'
 G1 = Geospatial_data(file_name = project.onshore_dir_name + '.pts')
 G2 = Geospatial_data(file_name = project.coast_dir_name + '.xya')
@@ -145 +119 @@
 G_off14 = Geospatial_data(file_name = project.offshore_dir_name14 + '.xya')
 
-print'add G1+G2+G3+all offshore data'
+print'add all geospatial objects'
 G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5 \
     + G_off6 + G_off7 + G_off8 + G_off9 + G_off10 + G_off11 + G_off12 \
     + G_off13 + G_off14
 
+print'clip combined geospatial object by bounding polygon'
 G.clip(project.bounding_polygon)
 #FIXME: add a clip function to pts
 
-print'export G'
-G.export_points_file(project.combined_dir_name + '.pts')
-import sys; sys.exit()
-
-
+print'export combined DEM file'
+if access(project.topographies_time_dir,F_OK) == 0:
+    mkdir (project.topographies_time_dir)
+G.export_points_file(project.combined_time_dir_name + '.pts')
 
 #-------------------------------------------------------------------------
 # Convert URS to SWW file for boundary conditions 
 #-------------------------------------------------------------------------
-
-# filenames
-meshname = project.mesh_name+'.msh'
-source_dir = project.boundarydir
-boundary_file = project.boundaryname
-
-print 'Available boundary tags', domain.get_boundary_tags()
+print 'starting to create boundary conditions'
+boundaries_dir_name = project.boundaries_dir_name
 
 from anuga.shallow_water.data_manager import urs2sww
 
-urs2sww(boundary_file,
-        minlat=project.north_boundary, maxlat=project.south_boundary,
+print 'minlat=project.north_boundary, maxlat=project.south_boundary',project.north_boundary, project.south_boundary
+print 'minlon= project.west_boundary, maxlon=project.east_boundary',project.west_boundary, project.east_boundary
+
+if access(project.boundaries_time_dir,F_OK) == 0:
+    mkdir (project.boundaries_time_dir)
+urs2sww(boundaries_dir_name,basename_out= project.boundaries_time_dir_name,
+        minlat=project.south_boundary, maxlat=project.north_boundary,
         minlon= project.west_boundary, maxlon=project.east_boundary,
+        mint=0, maxt= 35000,
         verbose='true')
+        
 
 
@@ -182 +158 @@
 
 
+
+