Changeset 6254


Ignore:
Timestamp:
Jan 31, 2009, 4:43:13 AM (10 years ago)
Author:
ole
Message:

Established work area for the development of a modelling demo

Location:
anuga_work/production/busselton/standardised_version
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/standardised_version/build_busselton.py

    r6253 r6254  
    4949# Fine pts file to be clipped to area of interest
    5050#-------------------------------------------------------------------------------
     51
     52# FIXME(Ole): Here's the Shark Bay way.
     53
     54## FIXME (Ole): Clip data by interior regions if possible
     55#
     56#print 'creating geospatial data objects from asc data (via dem and pts formats)'
     57#for filename in project.ascii_grid_filenames:
     58#    convert_dem_from_ascii2netcdf(filename,
     59#                                  basename_out=filename,
     60#                                  use_cache=True, verbose=True)
     61#    dem2pts(filename, use_cache=True, verbose=True)#
     62#
     63#    geospatial_data += Geospatial_data(file_name=filename + '.pts',
     64#                                       verbose=True)
     65##
     66#
     67#print 'creating geospatial data objects from txt data'
     68#for filename in project.point_filenames:
     69#    geospatial_data += Geospatial_data(file_name=filename + '.txt',
     70#                                       verbose=True)
     71#
     72#
     73#print 'clip combined geospatial object by bounding polygon'
     74#G = geospatial_data.clip(project.bounding_polygon)
     75#
     76#
     77#print 'export combined geospatial data'
     78#if access(project.topographies_dir, F_OK) == 0:
     79#    mkdir (project.topographies_dir)
     80#G.export_points_file(project.combined_dir_name + '.pts')
     81
     82
     83
    5184print "project.poly_all", project.poly_all
    5285print "project.combined_dir_name", project.combined_dir_name
  • anuga_work/production/busselton/standardised_version/project.py

    r6253 r6254  
    102102#------------------------------------------------------------------------------
    103103
     104# FIXME(Ole): Use lists of elevation data files as in Shark Bay study:
     105## elevation data filenames
     106#ascii_grid_filenames = ['10m_dem_without_survey', '50m_dem_without_10m_dem',
     107#                        'bathysteeppt', 'bathyleft', 'bathyright']
     108#point_filenames = ['field_survey_north', 'field_survey_south',
     109#                   'clipped_bathymetry_final', 'coast_points_final']
     110
     111
     112
    104113# elevation data used in build_busselton.py
    105114# onshore data: format ascii grid with accompanying projection file
  • anuga_work/production/busselton/standardised_version/run_busselton.py

    r6253 r6254  
    5454def run_model(**kwargs):
    5555   
    56     #------------------------------------------------------------------------------
     56    #-----------------------------------------------------------------------
    5757    # Copy scripts to time stamped output directory and capture screen
    5858    # output to file
    59     #------------------------------------------------------------------------------
     59    #-----------------------------------------------------------------------
    6060    print "Processor Name:",get_processor_name()
    6161
     
    102102    # causes problems with the ability to cache set quantity which takes alot of times
    103103       
     104
     105    # FIXME(Ole): Introduce create_domain_from_regions: Simpler and caches well
    104106    print 'start create mesh from regions'
    105107
Note: See TracChangeset for help on using the changeset viewer.