Changeset 6271

Timestamp:

Feb 4, 2009, 1:51:49 PM (15 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:

• build_patong.py (modified) (1 diff)
• project.py (modified) (2 diffs)

anuga_work/production/patong/build_patong.py

@@ -82 +82 @@
 
 print 'clipping outside extents'
-G1_clip = G1.clip(project.extent_elev_dir4, verbose=True) 
-G3_clip1 = G3.clip_outside(project.extent_elev_dir4, verbose=True) 
+#G1_clip = G1.clip(project.extent_elev_dir4, verbose=True) 
+G3_clip1 = G3.clip_outside(project.extent_elev_dir1, verbose=True) 
 G3_clip2 = G3_clip1.clip_outside(project.extent_elev_dir2, verbose=True)
 G4_clip = G4.clip_outside(project.extent_elev_dir3, verbose=True)
 
 #print 'hello', dir(G3_clip1), dir(G3_clip1), dir(G4_clip)
-print'add all geospatial objects'
-G = G1_clip + G2 + G3_clip1 + G3_clip2 + G4_clip
+print 'Add all geospatial objects'
+G = G1 + G2 + G3_clip2 + G4_clip
 
-print'clip combined geospatial object by bounding polygon'
+print 'Clip combined geospatial object by bounding polygon'
 G_clipped = G.clip(project.poly_all)
 
-print'export combined DEM file'
-if access(project.topographies_dir,F_OK) == 0:
-    mkdir (project.topographies_dir)
+print 'Oxport combined DEM file'
 G_clipped.export_points_file(project.combined_dir_name + '.pts')
 G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC

anuga_work/production/patong/project.py

@@ -48 +48 @@
 finaltime=15000         # final time for simulation 15000
 
-setup='final'  # Final can be replaced with trial or basic.
+setup='trial'  # Final can be replaced with trial or basic.
                # Either will result in a coarser mesh that will allow a
                # faster, but less accurate, simulation.
@@ -89 +89 @@
 elevation2 = 'patong_10m_small_grid_for_anuga_sub_Project.txt'
 elevation3 = 'patong_bay_1s_grid_Project.txt'
-elevation4 = 'andaman_3s_grid_Clip_Project.txt'
+elevation4 = 'andaman_3s_grid_Clip2_Project.txt'
 
 # gauges - used in get_timeseries.py