Changeset 6271


Ignore:
Timestamp:
Feb 4, 2009, 1:51:49 PM (10 years ago)
Author:
ole
Message:

Attempt to get Patong running with wider extent to minimise boundary effects.
Mesh is still as original though. This needs to fixed.

Location:
anuga_work/production/patong
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/patong/build_patong.py

    r6269 r6271  
    8282
    8383print 'clipping outside extents'
    84 G1_clip = G1.clip(project.extent_elev_dir4, verbose=True)
    85 G3_clip1 = G3.clip_outside(project.extent_elev_dir4, verbose=True)
     84#G1_clip = G1.clip(project.extent_elev_dir4, verbose=True)
     85G3_clip1 = G3.clip_outside(project.extent_elev_dir1, verbose=True)
    8686G3_clip2 = G3_clip1.clip_outside(project.extent_elev_dir2, verbose=True)
    8787G4_clip = G4.clip_outside(project.extent_elev_dir3, verbose=True)
    8888
    8989#print 'hello', dir(G3_clip1), dir(G3_clip1), dir(G4_clip)
    90 print'add all geospatial objects'
    91 G = G1_clip + G2 + G3_clip1 + G3_clip2 + G4_clip
     90print 'Add all geospatial objects'
     91G = G1 + G2 + G3_clip2 + G4_clip
    9292
    93 print'clip combined geospatial object by bounding polygon'
     93print 'Clip combined geospatial object by bounding polygon'
    9494G_clipped = G.clip(project.poly_all)
    9595
    96 print'export combined DEM file'
    97 if access(project.topographies_dir,F_OK) == 0:
    98     mkdir (project.topographies_dir)
     96print 'Oxport combined DEM file'
    9997G_clipped.export_points_file(project.combined_dir_name + '.pts')
    10098G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC
  • anuga_work/production/patong/project.py

    r6264 r6271  
    4848finaltime=15000         # final time for simulation 15000
    4949
    50 setup='final'  # Final can be replaced with trial or basic.
     50setup='trial'  # Final can be replaced with trial or basic.
    5151               # Either will result in a coarser mesh that will allow a
    5252               # faster, but less accurate, simulation.
     
    8989elevation2 = 'patong_10m_small_grid_for_anuga_sub_Project.txt'
    9090elevation3 = 'patong_bay_1s_grid_Project.txt'
    91 elevation4 = 'andaman_3s_grid_Clip_Project.txt'
     91elevation4 = 'andaman_3s_grid_Clip2_Project.txt'
    9292
    9393# gauges - used in get_timeseries.py
Note: See TracChangeset for help on using the changeset viewer.