Changeset 5245 for anuga_work/production/west_tas_2008/run_west_tas_smf.py

Timestamp:

Apr 28, 2008, 8:00:14 AM (15 years ago)

Author:

sexton

Message:

(1) update of production processes document (2) test for generating points for URS output (3) update of graduate proposal

File:

• anuga_work/production/west_tas_2008/run_west_tas_smf.py (modified) (1 diff)

anuga_work/production/west_tas_2008/run_west_tas_smf.py

@@ -54 +54 @@
 
 # filenames
-on_offshore25_dem_name = project_smf.on_offshore25_dem_name
+onshore_250_dem_name = project_smf.onshore_250_dem_name
 meshname = project_smf.meshname+'.msh'
 
 # creates DEM from asc data
-convert_dem_from_ascii2netcdf(on_offshore25_dem_name, use_cache=True, verbose=True)
+convert_dem_from_ascii2netcdf(onshore_250_dem_name, use_cache=True, verbose=True)
 
 #creates pts file for onshore DEM
-dem2pts(on_offshore25_dem_name,
-        easting_min=project_smf.eastingmin25,
-        easting_max=project_smf.eastingmax25,
-        northing_min=project_smf.northingmin25,
-        northing_max= project_smf.northingmax25,
-        use_cache=True, verbose=True)
+dem2pts(onshore_250_dem_name, use_cache=True, verbose=True)
 
 print 'create offshore'
-G1 = Geospatial_data(file_name = project_smf.offshore_dem_name1 + '.xya')+\
-     Geospatial_data(file_name = project_smf.offshore_dem_name2 + '.xya')+\
-     Geospatial_data(file_name = project_smf.offshore_dem_name3 + '.xya')
+G1 = Geospatial_data(file_name = project_smf.offshore_dem_name1 + '.txt')+\
+     Geospatial_data(file_name = project_smf.offshore_dem_name2 + '.txt')+\
+     Geospatial_data(file_name = project_smf.offshore_dem_name3 + '.txt')
 
 print 'create onshore'
-G2 = Geospatial_data(file_name = project_smf.on_offshore25_dem_name + '.pts')
+G2 = Geospatial_data(file_name = project_smf.onshore_250_dem_name + '.pts')
+
+print 'create coastline'
+G3 = Geospatial_data(file_name = project_smf.coast_line + '.txt')
 
 print 'add'
-G = G1 + G2
+G = G1.clip(Geospatial_data(project_smf.polyAll)) + G2 + G3
 
 print 'export points'