Changeset 3281

Timestamp:

Jul 6, 2006, 12:22:16 PM (19 years ago)

Author:

sexton

Message:

updated scripts for pt hedland - data received

Location:

production/pt_hedland_2006

Files:

• project.py (modified) (3 diffs)
• run_pt_hedland.py (modified) (3 diffs)

production/pt_hedland_2006/project.py

@@ -26 +26 @@
 offshore_name2 = 'pt_hedland_offshore_points_fairsheet'
 
-# coastline developed from 30m DTED and aerial photography
-coast_name = 'pt_hedland_coastline_points_dted2_new'
+# coastline developed from aerial photography and 1.5m DLI contour
+#coast_name = 'pt_hedland_coastline_points_dted2_new'
+coast_name = 'coast_with_extracted_z'
+
 
 boundary_basename = 'SU-AU_clip'
@@ -39 +41 @@
     home = getenv('INUNDATIONHOME') #Sandpit's parent dir
 #    python_home = getenv('PWD')     
-    #user = basename(getenv('USERPROFILE'))
+    user = getenv('USERPROFILE')
     #print 'USER:', user
 else:
@@ -76 +78 @@
 
 meshname = meshdir + basename
-onshore_dem_name = datadir + onshore_name_dted
+#onshore_dem_name = datadir + onshore_name_dted
+onshore_dem_name = datadir + onshore_name_dli
 offshore_dem_name1 = datadir + offshore_name1
 offshore_dem_name2 = datadir + offshore_name2

production/pt_hedland_2006/run_pt_hedland.py

@@ -63 +63 @@
 
 # filenames
-onshore_dem_name = project.onshore_dem_name
-coast_dem_name = project.coast_dem_name
-offshore_points1 = project.offshore_dem_name1
-offshore_points2 = project.offshore_dem_name2
 meshname = project.meshname+'.msh'
 source_dir = project.boundarydir
@@ -72 +68 @@
 # fine data (clipping the points file to smaller area)
 # creates DEM from asc data 
-convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True)
+convert_dem_from_ascii2netcdf(project.onshore_dem_name, use_cache=True, verbose=True)
 
 #creates pts file from DEM
-dem2pts(onshore_dem_name,
+dem2pts(project.onshore_dem_name,
         easting_min=project.eastingmin,
         easting_max=project.eastingmax,
@@ -105 +101 @@
 
 region_res = 100000
-coast_res = 25000
-pt_hedland_res = 5000
+coast_res = 2500
+pt_hedland_res = 1500
 interior_regions = [[project.poly_pt_hedland, pt_hedland_res],
                     [project.poly_region, region_res]]