Changeset 5245 for anuga_work/production/west_tas_2008

Timestamp:

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

Author:

sexton

Message:

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

Location:

anuga_work/production/west_tas_2008

Files:

• convert2eastnorth.py (modified) (1 diff)
• project_smf.py (modified) (4 diffs)
• run_west_tas_smf.py (modified) (1 diff)

anuga_work/production/west_tas_2008/convert2eastnorth.py

@@ -31 +31 @@
 #filename = project_smf.polygondir + 'domain.txt'
 #filename = project_smf.polygondir + 'local.txt'
-filename = project_smf.polygondir + 'region.txt'
+#filename = project_smf.polygondir + 'region.txt'
+#filename = project_smf.datadir + 'Sorell_Basin_150m_DD_Interp.txt'
+#filename = project_smf.datadir + 'S_Tasman_Rise_150m_DD_Interp.txt'
+filename = project_smf.polygondir + 'origin.txt'
 #filename_out = project_smf.polygondir + 'domain.csv'
 #filename_out = project_smf.polygondir + 'local.csv'
-filename_out = project_smf.polygondir + 'region.csv'
+#filename_out = project_smf.polygondir + 'region.csv'
+#filename_out = project_smf.datadir + 'Sorell_Basin_150m_EN_Interp.txt'
+#filename_out = project_smf.datadir + 'S_Tasman_Rise_150m_EN_Interp.txt'
+filename_out = project_smf.polygondir + 'origin_EN.csv'
 alter_file(filename,filename_out)

anuga_work/production/west_tas_2008/project_smf.py

@@ -24 +24 @@
 state = 'tasmania'
 scenario_dir_name = 'smf_scenario_2008'
-'''
-# onshore data provided by NSW LPI, Krishna merged with existing
-# PMD data to make 100m DEM
-on_offshore_name100 = 'bathyland100'
+
+# onshore data sourced from 250m grid
+onshore_name250 = 'bathyland250'
+
+# coastline from GA
+coast_line_name = 'aus_cst'
 
 # survey data
-offshore_name1 = 'surveyAreaA'
-offshore_name2 = 'surveyAreaB'
-offshore_name3 = 'surveyAreaC'
+offshore_name1 = 'S_Tasman_Rise_150m_EN_Interp'
+offshore_name2 = 'Sorell_Basin_150m_EN_Interp'
+offshore_name3 = 'West_Tas_EastNorth_z'
 
-# AHO data
-offshore_name4 = '1000003611export'
-offshore_name5 = '1000003613export'
-offshore_name6 = '1000003614export'
-offshore_name7 = '1000003627export'
-offshore_name8 = '1000003628export'
-offshore_name9 = 'AHDexport'
-'''
 #swollen/ all data output
 basename = 'source'
@@ -61 +55 @@
 codedirname = codedir + 'project_smf.py'
 meshname = outputtimedir + 'mesh_' + basename
-'''
+
 # Necessary if using point datasets, rather than grid
-on_offshore100_dem_name = datadir + on_offshore_name100
+onshore_250_dem_name = datadir + onshore_name250
 offshore_dem_name1 = datadir + offshore_name1
 offshore_dem_name2 = datadir + offshore_name2
 offshore_dem_name3 = datadir + offshore_name3
+coast_line = datadir + coast_line_name
 
 combined_dem_name   = datadir + 'west_tas_combined_elevation'
-'''
+
 ###############################
 # Domain definitions
@@ -94 +89 @@
 ###################################################################
 
-
 # exporting asc grid
 eastingmin = 333000
@@ -107 +101 @@
 
 # historical slides
-#slide_origin = [,]
+slide_origin = [338270.5407, 5217273.89]
 
 depth = 1750.0

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'