source: production/gippsland_2005/create_mesh.py @ 2378

"""Create a mesh of the Gippsland Lakes coast based on 100m
   bathymetric data.
"""

#import os, sys
#sys.path.append('..')



# add storm_surge dir to your pythonpath
from pmesh.mesh import *
from pyvolution.coordinate_transforms.geo_reference import Geo_reference

#Lower right corner
xlrcorner = 601018
ylrcorner = 5797629
    
def create_mesh(max_area, mesh_file=None, triangles_in_name = False):
    global xlrcorner
    global ylrcorner
    
    #Lower left corner
    xllcorner = 570066
    yllcorner = 5773910

    #Upper left corner
    xulcorner = 553571
    yulcorner = 5796228

    #Upper right corner
    xurcorner = 583675
    yurcorner = 5820018

    #Lower right corner
    #xlrcorner = 601018
    #ylrcorner = 5797629


    #create geo reference object for use with mesh
    geo = Geo_reference(xllcorner = xllcorner,
                        yllcorner = yllcorner,
                        zone = 56)
    print "***********************"
    print "geo ref", geo
    print "***********************"
    
    m = Mesh(geo_reference=geo)

    #Boundary
    dict = {}
    dict['points'] = [[xllcorner, yllcorner],   #sw
                      [xulcorner, yulcorner],
                      [xurcorner, yurcorner],
                      [xlrcorner, ylrcorner]]
    
    dict['segments'] = [[0,1], [1,2], [2,3], [3,0]] # the outer border
    dict['segment_tags'] = ['side',
                            'back',                            
                            'side',
                            'ocean']
    m.addVertsSegs(dict)

    m.generateMesh('q28.0z', maxArea = max_area)
    triangle_count = len(m.getTriangulation())
    if mesh_file is None:    
        return m, triangle_count
    else:
        if triangles_in_name is True:
            mesh_file = mesh_file[:-4] + '_' + str(triangle_count) \
                        + mesh_file[-4:]
        m.export_mesh_file(mesh_file)
        return mesh_file, triangle_count


#-------------------------------------------------------------
if __name__ == "__main__":
    _, triangle_count = create_mesh(1000000)
    print "triangle_count",triangle_count