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

Timestamp:

Nov 8, 2006, 4:04:55 PM (17 years ago)

Author:

nick

Message:

update to dampier

File:

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

anuga_work/production/dampier_2006/build_dampier.py

@@ -32 +32 @@
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 from anuga.geospatial_data.geospatial_data import *
-from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
+from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
 
 # Application specific imports
@@ -42 +42 @@
 #------------------------------------------------------------------------------
 
-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()
-
-
-#normal screen output is stored in 
-screen_output_name = project.output_build_time_dir + "screen_output.txt"
-screen_error_name = project.output_build_time_dir + "screen_error.txt"
-
-#used to catch screen output to file
-sys.stdout = Screen_Catcher(screen_output_name)
-sys.stderr = Screen_Catcher(screen_error_name)
+copy_code_files(project.output_build_time_dir,__file__, 
+               dirname(project.__file__)+sep+ project.__name__+'.py' )
+
+start_screen_catcher(project.output_build_time_dir)
 
 print 'USER:    ', project.user
@@ -125 +113 @@
 G_off13 = Geospatial_data(file_name = project.offshore_dir_name13 + '.xya')
 G_off14 = Geospatial_data(file_name = project.offshore_dir_name14 + '.xya')
-
+'''
+
+
+print'clip nw', project.clip_poly_nw
+print'clip e', project.clip_poly_e
+
+print'reading combined_dir_name'
+G_offshore_data = Geospatial_data(file_name = project.topographies_dir+'dampier_combined_elevation_final.pts')
+
+print'reading offshore_dir_name_old'
+G_offshore_old = Geospatial_data(file_name = project.offshore_dir_name_old + '.pts')
+
+'''
+print 'G_offshore_data_old_nw', 
+G_nw_name = project.topographies_dir + 'nw_old_data'
+G_offshore_nw = G_offshore_old.clip(project.clip_poly_nw)
+G_offshore_nw.export_points_file(G_nw_name + '.pts')
+G_offshore_nw.export_points_file(G_nw_name + '.xya')
+G_nw = array(G_offshore_nw.get_data_points())
+print'shape of arr nw data', G_nw.shape
+print' max and min of array 0',max(G_nw[:,0]),min(G_nw[:,0])
+print' max and min of array 1',max(G_nw[:,1]),min(G_nw[:,1])
+
+print 'G_offshore_data_old_e'#, G_offshore_old_nw.get_data_points()
+G_e_name = project.topographies_dir+'e_old_data'
+G_offshore_e = G_offshore_old.clip(project.clip_poly_e)
+G_offshore_e.export_points_file(G_e_name + '.pts')
+G_offshore_e.export_points_file(G_e_name + '.xya')
+G_e = array(G_offshore_e.get_data_points())
+print'shape of arr e data', G_e.shape
+print' max and min of array 0',max(G_e[:,0]),min(G_e[:,0])
+print' max and min of array 1',max(G_e[:,1]),min(G_e[:,1])
+'''
+
+print 'G_offshore_data_mid_e'#, G_offshore_old_nw.get_data_points()
+G_mid_e_name = project.topographies_dir+'mid_e_old_data'
+G_offshore_mid_e = G_offshore_old.clip(project.clip_poly_mid_e)
+G_offshore_mid_e.export_points_file(G_mid_e_name + '.pts')
+G_offshore_mid_e.export_points_file(G_mid_e_name + '.xya')
+G_mid_e = array(G_offshore_mid_e.get_data_points())
+print'shape of arr e data', G_mid_e.shape
+print' max and min of array 0',max(G_mid_e[:,0]),min(G_mid_e[:,0])
+print' max and min of array 1',max(G_mid_e[:,1]),min(G_mid_e[:,1])
+
+print 'G_offshore_data_mid_w'#, G_offshore_old_nw.get_data_points()
+G_mid_w_name = project.topographies_dir+'mid_w_old_data'
+G_offshore_mid_w = G_offshore_old.clip(project.clip_poly_mid_w)
+G_offshore_mid_w.export_points_file(G_mid_w_name + '.pts')
+G_offshore_mid_w.export_points_file(G_mid_w_name + '.xya')
+G_mid_w = array(G_offshore_mid_w.get_data_points())
+print'shape of arr e data', G_mid_w.shape
+print' max and min of array 0',max(G_mid_w[:,0]),min(G_mid_w[:,0])
+print' max and min of array 1',max(G_mid_w[:,1]),min(G_mid_w[:,1])
+
+
+print 'G_offshore_data_old_e'#, G_offshore_old_e.get_data_points() 
 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)
+#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
+
+G = G_offshore_data + G_offshore_mid_w + G_offshore_mid_e
+print'shape of arr G data', G.get_data_points().shape
+
+
+#print'clip combined geospatial object by bounding polygon'
+G_clipped = G.clip(project.bounding_polygon)
 #FIXME: add a clip function to pts
+print'shape of clipped data', G_clipped.get_data_points().shape
 
 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')
+G_clipped.export_points_file(project.combined_time_dir_final_name + '.pts')
+G_clipped.export_points_file(project.combined_time_dir_final_name + '.xya')
 '''
 #-------------------------------------------------------------------------
@@ -156 +205 @@
 if access(project.boundaries_time_dir,F_OK) == 0:
     mkdir (project.boundaries_time_dir)
-'''
-urs2sww(boundaries_in_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')
-
-
-urs2sww(boundaries_in_dir_name, basename_out= project.boundaries_time_dir_name,
-#        minlat=project.south, maxlat=project.north,
-#        minlon= project.west, maxlon=project.east,
-        minlat=project.south_boundary, maxlat=project.north_boundary,
-        minlon= project.east_boundary, maxlon=project.west_boundary,
-        mint=0, maxt= 35000,
-        verbose='true')
-'''
+
 from caching import cache
 cache(urs2sww,
@@ -193 +227 @@
        )
 #       dependencies = source_dir + project.boundary_basename + '.sww')
-   
-
-
-
-
-
-
-
-
+'''   
+
+
+
+
+
+
+
+