Changeset 5411 for anuga_work/development/convergence_okushiri_2008/create_okushiri_truescale.py

Timestamp:

Jun 23, 2008, 1:01:17 PM (15 years ago)

Author:

Leharne

Message:

Updates to okushiri convergence study

File:

• anuga_work/development/convergence_okushiri_2008/create_okushiri_truescale.py (modified) (2 diffs)

anuga_work/development/convergence_okushiri_2008/create_okushiri_truescale.py

@@ -6 +6 @@
 
 from anuga.pmesh.mesh import *
-from anuga.pmesh.mesh_interface import create_mesh_from_regions
 from anuga.coordinate_transforms.geo_reference import Geo_reference
 from anuga.geospatial_data import Geospatial_data
@@ -112 +111 @@
 
 
-    #--------------------------------------------------------------------------
-    # Create the triangular mesh based on overall clipping polygon with a
-    # tagged
-    # boundary and interior regions defined in project_truescale.py along with
-    # resolutions (maximal area of per triangle) for each polygon
-    #--------------------------------------------------------------------------
-
-
-    base_resolution = 1 # Use this to coarsen or refine entire mesh. 
-
-    # Basic geometry and bounding polygon
-    xleft   = 0
-    xright  = 2179.2
-    ybottom = 0
-    ytop    = 1360.8
-
-    point_sw = [xleft, ybottom]
-    point_se = [xright, ybottom]
-    point_nw = [xleft, ytop]    
-    point_ne = [xright, ytop]
-
-    bounding_polygon = [point_se,
-                        point_ne,
-                        point_nw,
-                        point_sw]
-         
-    # Localised refined area for gulleys
-    xl = 1920
-    xr = 2120
-    yb = 640
-    yt = 920
-    p0 = [xl, yb]
-    p1 = [xr, yb]
-    p2 = [xr, yt]
-    p3 = [xl, yt]
     
-    gulleys = [p0, p1, p2, p3]
-    
-
-    # Island area and drawdown region (original)
-    #island_0 = [xleft + 2*(xright-xleft)/3+1.2, ytop-0.5]
-    #island_1 = [xleft + 2*(xright-xleft)/3+0.5, ybottom + 2*(ytop-ybottom)/3]
-    #island_2 = [xleft + (xright-xleft)/2+0.3, ybottom + 2*(ytop-ybottom)/3-0.3]
-    #island_3 = [xleft + (xright-xleft)/2+0.3, ybottom + (ytop-ybottom)/3+0.3]
-    #island_4 = [xleft + 2*(xright-xleft)/3+0.4, ybottom + (ytop-ybottom)/3-0.3]
-    #island_5 = [xleft + 2*(xright-xleft)/3+1.2, ybottom+0.2]
-    #island_6 = [xl-.01, yb]  #OK
-    #island_7 = [xl-.01, yt]  #OK      
-
-
-    # Island area and drawdown region
-    island_0 = [xleft + 2*(xright-xleft)/3+480, ytop-200]
-    island_1 = [xleft + 2*(xright-xleft)/3+200, ybottom + 2*(ytop-ybottom)/3]
-    island_2 = [xleft + (xright-xleft)/2+160, ybottom + 2*(ytop-ybottom)/3-120]
-    island_3 = [xleft + (xright-xleft)/2+160, ybottom + (ytop-ybottom)/3+120]
-    island_4 = [xleft + 2*(xright-xleft)/3+160, ybottom + (ytop-ybottom)/3-120]
-    island_5 = [xleft + 2*(xright-xleft)/3+480, ybottom+320]
-    island_6 = [xl-4, yb]  # Keep right edge just off the gulleys
-    island_7 = [xl-4, yt]
- 
-    island = [island_0, island_1, island_2,
-              island_3, island_4, island_5,
-              island_6, island_7]
-
-
-    # Region spanning half right hand side of domain just inside boundary (org)
-    #rhs_nw = [xleft + (xright-xleft)/3+1, ytop-0.02]
-    #rhs_sw = [xleft + (xright-xleft)/3+1, ybottom+0.02]
-    #rhs_se = [xright-0.02, ybottom+0.02]
-    #rhs_ne = [xright-0.02, ytop-0.02]
-
-    # Region spanning half right hand side of domain just inside boundary
-    rhs_nw = [xleft + (xright-xleft)/3+400, ytop-560]
-    rhs_sw = [xleft + (xright-xleft)/3+400, ybottom+200]
-    rhs_se = [xright-4, ybottom+80]
-    rhs_ne = [xright-4, ytop-80]        
-
-    rhs_region = [rhs_nw, rhs_ne, rhs_se, rhs_sw]
-
-     
-    # Interior regions and creation of mesh
-    interior_regions = [[rhs_region, 80],
-                        [island, 32*base_resolution],
-                        [gulleys, 3.2*base_resolution]]    
-
-    meshname = project_truescale.mesh_filename + '.msh'
-    m = create_mesh_from_regions(bounding_polygon,
-                                 boundary_tags={'wall': [0, 1, 3],
-                                                'wave': [2]},     
-                                 maximum_triangle_area=16000*base_resolution,
-                                 interior_regions=interior_regions,
-                                 filename=project_truescale.mesh_filename,
-                                 verbose=True)
-