Changeset 6024


Ignore:
Timestamp:
Nov 28, 2008, 11:43:16 AM (15 years ago)
Author:
kristy
Message:

added more data

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/carnarvon/build_carnarvon.py

    r5804 r6024  
    5454# input elevation directory filenames
    5555onshore_in_dir_name = project.onshore_in_dir_name
    56 coast_in_dir_name = project.coast_in_dir_name
    5756offshore_in_dir_name = project.offshore_in_dir_name
    5857offshore_in_dir_name1 = project.offshore_in_dir_name1
    59 offshore_in_dir_name2 = project.offshore_in_dir_name2
     58offshore_in_dir_aoi = project.offshore_in_dir_aoi
     59offshore_in_dir_aos = project.offshore_in_dir_aos
    6060
    6161# output elevation directory filenames
    6262onshore_dir_name = project.onshore_dir_name
    63 coast_dir_name = project.coast_dir_name
    6463offshore_dir_name = project.offshore_dir_name
    6564offshore_dir_name1 = project.offshore_dir_name1
    66 offshore_dir_name2 = project.offshore_dir_name2
     65offshore_dir_aoi = project.offshore_dir_aoi
     66offshore_dir_aos = project.offshore_dir_aos
    6767
    6868
     
    7070print "creates DEMs from asc data"
    7171convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True)
     72convert_dem_from_ascii2netcdf(offshore_in_dir_aoi, basename_out=offshore_dir_aoi, use_cache=True, verbose=True)
     73convert_dem_from_ascii2netcdf(offshore_in_dir_aos, basename_out=offshore_dir_aos, use_cache=True, verbose=True)
    7274
    7375# creates pts file for onshore DEM
    7476print "creates pts file for onshore DEM"
    7577dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
     78dem2pts(offshore_dir_aoi ,use_cache=True,verbose=True)
     79dem2pts(offshore_dir_aos ,use_cache=True,verbose=True)
     80
     81
    7682
    7783# create onshore pts files
    78 print'create Geospatial data1 objects from topographies'
     84print'create Geospatial onshore objects from topographies'
    7985G1 = Geospatial_data(file_name = onshore_dir_name + '.pts')
    8086
    8187# 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'
     88print'create Geospatial offshore objects from topographies'
    8589G_off1 = Geospatial_data(file_name = offshore_in_dir_name)
    86 print'create Geospatial data6 objects from topographies'
     90print'create Geospatial makeup objects from topographies'
    8791G_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)
     92print'create Geospatial aoi objects from topographies'
     93G_offaoi = Geospatial_data(file_name = offshore_dir_aoi+ '.pts')
     94print'create Geospatial aos objects from topographies'
     95G_offaos = Geospatial_data(file_name = offshore_dir_aos+ '.pts')
     96
    9097
    9198#-------------------------------------------------------------------------------
     
    94101
    95102print 'add all bathymetry files'
    96 G_off = G_off1 + G_off2 + G_off3
     103G_off = G_off1 + G_off2
     104
    97105print 'clipping data to coast'
    98106G_ocean = G_off.clip(project.poly_ocean, verbose=True) #mainland
     
    105113
    106114print'add all geospatial objects'
    107 G = G1_clip_m+ G1_clip_b+G1_clip_d+ G_coast + G_ocean_bd
     115G_aos = G1_clip_m+ G1_clip_b+G1_clip_d+ G_coast + G_ocean_bd
     116G_aosclip = G_aos.clip_outside(project.poly_aos_extract, verbose=True)
     117G_aoi = G_aosclip.clip_outside(project.poly_aoi1, verbose=True)
     118G = G_aoi +Goffaoi+Goffaos
    108119
    109120print'export combined DEM file'
Note: See TracChangeset for help on using the changeset viewer.