Changeset 6359


Ignore:
Timestamp:
Feb 18, 2009, 9:14:00 AM (16 years ago)
Author:
myall
Message:

phase 2 scripts for ceduna, eucla and esperance still not working...

Location:
anuga_work/production/australia_ph2
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/ceduna/build_elevation.py

    r6293 r6359  
    5050# Create Geospatial data from ASCII files
    5151geospatial_data = {}
    52 for filename in project.ascii_grid_filenames:
     52##for filename in project.ascii_grid_filenames:
     53##    absolute_filename = join(project.topographies_folder, filename)
     54##    convert_dem_from_ascii2netcdf(absolute_filename,
     55##                                  basename_out=absolute_filename,
     56##                                  use_cache=True,
     57##                                  verbose=True)
     58##    dem2pts(absolute_filename, use_cache=True, verbose=True)
     59
     60##    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
     61##                                                verbose=True)
     62
     63# Create Geospatial data from TXT files
     64for filename in project.point_filenames:
    5365    absolute_filename = join(project.topographies_folder, filename)
    54     convert_dem_from_ascii2netcdf(absolute_filename,
    55                                   basename_out=absolute_filename,
    56                                   use_cache=True,
    57                                   verbose=True)
    58     dem2pts(absolute_filename, use_cache=True, verbose=True)
    59 
    60     geospatial_data[filename] = Geospatial_data(file_name=absolute_filename+'.pts',
     66    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
    6167                                                verbose=True)
    62 
    63 ### Create Geospatial data from TXT files
    64 ##for filename in project.point_filenames:
    65 ##    absolute_filename = join(project.topographies_folder, filename)
    66 ##    geospatial_data[filename] = Geospatial_data(file_name=absolute_filename,
    67 ##                                                verbose=True)
    6868
    6969
  • anuga_work/production/australia_ph2/ceduna/project.py

    r6342 r6359  
    4646finaltime = 1000 #80000         # final time for simulation
    4747
    48 setup = 'final'  # Final can be replaced with trial or basic.
     48setup = 'trial'  # Final can be replaced with trial or basic.
    4949               # Either will result in a coarser mesh that will allow a
    5050               # faster, but less accurate, simulation.
     
    8181# Used in build_elevation.py
    8282# Format for ascii grids, as produced in ArcGIS + a projection file
    83 ascii_grid_filenames = ['grid250m'] # 250m grid 2005
     83#ascii_grid_filenames = ['grid250m'] # 250m grid 2005
    8484
    8585# Format for point is x,y,elevation (with header)
    86 ##point_filenames = ['Busselton_Contour0.txt',     # Coastline
    87 ##                   'Busselton_BeachSurvey.txt',    # Beach survey
    88 ##                   'Busselton_NavyFinal.txt',  # Bathymetry
    89 ##                   'Busselton_Chart.txt', # Bathymetry Charts
    90 ##                   'Busselton_Digitised.txt', # Digitised Fairsheet
    91 ##                   'Busselton_250m.txt', # 250m
    92 ##                   'Bunbury_TIN.txt', # Bunbury aoi TIN'd in ArcGIS
    93 ##                   'Busselton_TIN.txt', # Busselton aoi TIN'd in ArcGIS
    94 ##                   'XYAHD_clip.txt'] # To extend boundary
    95 
     86point_filenames = ['grid250m.txt']
     87                 
    9688# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
    9789# Used in build_elevation.py
  • anuga_work/production/australia_ph2/eucla_motel/project.py

    r6316 r6359  
    4646finaltime = 1000 #80000         # final time for simulation
    4747
    48 setup = 'final'  # Final can be replaced with trial or basic.
     48setup = 'trial'  # Final can be replaced with trial or basic.
    4949               # Either will result in a coarser mesh that will allow a
    5050               # faster, but less accurate, simulation.
Note: See TracChangeset for help on using the changeset viewer.