Changeset 3951


Ignore:
Timestamp:
Nov 8, 2006, 6:41:19 PM (16 years ago)
Author:
sexton
Message:

update broome script to incorporate additional data delivered today

Location:
anuga_work/production/broome_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/broome_2006/project.py

    r3944 r3951  
    5656coast_name = 'coastline'
    5757
     58# TIN model to fill in data gap
     59bathy_interp = 'interpolate'
     60
    5861boundary_basename = 'SU-AU' # Mw ?
    5962
     
    8285# Necessary if using point datasets, rather than grid
    8386onshore_dem_name = datadir + onshore_name
     87offshore_interp_dem_name = datadir + bathy_interp
    8488offshore_dem_name1 = datadir + offshore_name1
    8589offshore_dem_name2 = datadir + offshore_name2
  • anuga_work/production/broome_2006/run_broome.py

    r3897 r3951  
    6363# filenames
    6464onshore_dem_name = project.onshore_dem_name
     65offshore_interp_dem_name = project.offshore_interp_dem_name
    6566coast_points = project.coast_dem_name
    6667meshname = project.meshname+'.msh'
     
    7172#creates pts file for onshore DEM
    7273dem2pts(onshore_dem_name, use_cache=True, verbose=True)
     74
     75# creates DEM from asc data
     76convert_dem_from_ascii2netcdf(offshore_interp_dem_name, use_cache=True, verbose=True)
     77
     78#creates pts file for offshore interpolated DEM
     79dem2pts(offshore_interp_dem_name, use_cache=True, verbose=True)
    7380
    7481print 'create offshore'
     
    94101     Geospatial_data(file_name = project.offshore_dem_name20 + '.xya')+\
    95102     Geospatial_data(file_name = project.offshore_dem_name21 + '.xya')+\
    96      Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')
     103     Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')+\
     104     Geospatial_data(file_name = project.offshore_interp_dem_name + '.pts')
    97105print 'create onshore'
    98106G2 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
Note: See TracChangeset for help on using the changeset viewer.