Changeset 6024
- Timestamp:
- Nov 28, 2008, 11:43:16 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/carnarvon/build_carnarvon.py
r5804 r6024 54 54 # input elevation directory filenames 55 55 onshore_in_dir_name = project.onshore_in_dir_name 56 coast_in_dir_name = project.coast_in_dir_name57 56 offshore_in_dir_name = project.offshore_in_dir_name 58 57 offshore_in_dir_name1 = project.offshore_in_dir_name1 59 offshore_in_dir_name2 = project.offshore_in_dir_name2 58 offshore_in_dir_aoi = project.offshore_in_dir_aoi 59 offshore_in_dir_aos = project.offshore_in_dir_aos 60 60 61 61 # output elevation directory filenames 62 62 onshore_dir_name = project.onshore_dir_name 63 coast_dir_name = project.coast_dir_name64 63 offshore_dir_name = project.offshore_dir_name 65 64 offshore_dir_name1 = project.offshore_dir_name1 66 offshore_dir_name2 = project.offshore_dir_name2 65 offshore_dir_aoi = project.offshore_dir_aoi 66 offshore_dir_aos = project.offshore_dir_aos 67 67 68 68 … … 70 70 print "creates DEMs from asc data" 71 71 convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True) 72 convert_dem_from_ascii2netcdf(offshore_in_dir_aoi, basename_out=offshore_dir_aoi, use_cache=True, verbose=True) 73 convert_dem_from_ascii2netcdf(offshore_in_dir_aos, basename_out=offshore_dir_aos, use_cache=True, verbose=True) 72 74 73 75 # creates pts file for onshore DEM 74 76 print "creates pts file for onshore DEM" 75 77 dem2pts(onshore_dir_name ,use_cache=True,verbose=True) 78 dem2pts(offshore_dir_aoi ,use_cache=True,verbose=True) 79 dem2pts(offshore_dir_aos ,use_cache=True,verbose=True) 80 81 76 82 77 83 # create onshore pts files 78 print'create Geospatial data1objects from topographies'84 print'create Geospatial onshore objects from topographies' 79 85 G1 = Geospatial_data(file_name = onshore_dir_name + '.pts') 80 86 81 87 # create coastal and offshore txt files 82 print'create Geospatial data4 objects from topographies' 83 G_coast = Geospatial_data(file_name = coast_in_dir_name) 84 print'create Geospatial data5 objects from topographies' 88 print'create Geospatial offshore objects from topographies' 85 89 G_off1 = Geospatial_data(file_name = offshore_in_dir_name) 86 print'create Geospatial data6objects from topographies'90 print'create Geospatial makeup objects from topographies' 87 91 G_off2 = Geospatial_data(file_name = offshore_in_dir_name1) 88 print'create Geospatial data7 objects from topographies' 89 G_off3 = Geospatial_data(file_name = offshore_in_dir_name2) 92 print'create Geospatial aoi objects from topographies' 93 G_offaoi = Geospatial_data(file_name = offshore_dir_aoi+ '.pts') 94 print'create Geospatial aos objects from topographies' 95 G_offaos = Geospatial_data(file_name = offshore_dir_aos+ '.pts') 96 90 97 91 98 #------------------------------------------------------------------------------- … … 94 101 95 102 print 'add all bathymetry files' 96 G_off = G_off1 + G_off2 + G_off3 103 G_off = G_off1 + G_off2 104 97 105 print 'clipping data to coast' 98 106 G_ocean = G_off.clip(project.poly_ocean, verbose=True) #mainland … … 105 113 106 114 print'add all geospatial objects' 107 G = G1_clip_m+ G1_clip_b+G1_clip_d+ G_coast + G_ocean_bd 115 G_aos = G1_clip_m+ G1_clip_b+G1_clip_d+ G_coast + G_ocean_bd 116 G_aosclip = G_aos.clip_outside(project.poly_aos_extract, verbose=True) 117 G_aoi = G_aosclip.clip_outside(project.poly_aoi1, verbose=True) 118 G = G_aoi +Goffaoi+Goffaos 108 119 109 120 print'export combined DEM file'
Note: See TracChangeset
for help on using the changeset viewer.