Ignore:
Timestamp:
Apr 28, 2008, 8:00:14 AM (16 years ago)
Author:
sexton
Message:

(1) update of production processes document (2) test for generating points for URS output (3) update of graduate proposal

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/west_tas_2008/run_west_tas_smf.py

    r4968 r5245  
    5454
    5555# filenames
    56 on_offshore25_dem_name = project_smf.on_offshore25_dem_name
     56onshore_250_dem_name = project_smf.onshore_250_dem_name
    5757meshname = project_smf.meshname+'.msh'
    5858
    5959# creates DEM from asc data
    60 convert_dem_from_ascii2netcdf(on_offshore25_dem_name, use_cache=True, verbose=True)
     60convert_dem_from_ascii2netcdf(onshore_250_dem_name, use_cache=True, verbose=True)
    6161
    6262#creates pts file for onshore DEM
    63 dem2pts(on_offshore25_dem_name,
    64         easting_min=project_smf.eastingmin25,
    65         easting_max=project_smf.eastingmax25,
    66         northing_min=project_smf.northingmin25,
    67         northing_max= project_smf.northingmax25,
    68         use_cache=True, verbose=True)
     63dem2pts(onshore_250_dem_name, use_cache=True, verbose=True)
    6964
    7065print 'create offshore'
    71 G1 = Geospatial_data(file_name = project_smf.offshore_dem_name1 + '.xya')+\
    72      Geospatial_data(file_name = project_smf.offshore_dem_name2 + '.xya')+\
    73      Geospatial_data(file_name = project_smf.offshore_dem_name3 + '.xya')
     66G1 = Geospatial_data(file_name = project_smf.offshore_dem_name1 + '.txt')+\
     67     Geospatial_data(file_name = project_smf.offshore_dem_name2 + '.txt')+\
     68     Geospatial_data(file_name = project_smf.offshore_dem_name3 + '.txt')
    7469
    7570print 'create onshore'
    76 G2 = Geospatial_data(file_name = project_smf.on_offshore25_dem_name + '.pts')
     71G2 = Geospatial_data(file_name = project_smf.onshore_250_dem_name + '.pts')
     72
     73print 'create coastline'
     74G3 = Geospatial_data(file_name = project_smf.coast_line + '.txt')
    7775
    7876print 'add'
    79 G = G1 + G2
     77G = G1.clip(Geospatial_data(project_smf.polyAll)) + G2 + G3
    8078
    8179print 'export points'
Note: See TracChangeset for help on using the changeset viewer.