source: inundation/pmesh/create_gong.py @ 3031

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

from mesh import *
from coordinate_transforms.geo_reference import Geo_reference

#-------------------------------------------------------------
if __name__ == "__main__":


    bellambi_northing = 6195000.0
    north_beach_northing = 6183000.0
    north_beach_easting = 308500
    border = 12000.0
    west = 301000.0 # From the map 
    inner_east = 310000.0 # From the map 
    inner_west_south = 303000.0 # From the map 
    inner_west_north = 307000.0  # From the map 
    east = 328421.0
    sea_boundary = (east+inner_east)/2
    inner_center_easting =  ( inner_east +inner_west_north)/2
    inner_center_northing = (bellambi_northing + north_beach_northing)/2

    # from 100m dem
    geo = Geo_reference(xllcorner = inner_center_easting,
                        yllcorner = inner_center_northing,
                        zone = 56)
    print "***********************"
    print "geo ref", geo
    print "***********************"
    
    m = Mesh(geo_reference=geo)
    
    factor = 1000 #low res 10000, high res 1000
    
    poly = [[west, north_beach_northing - border],  #sw
                         [west, bellambi_northing + border],  #nw
                         [sea_boundary,
                          bellambi_northing + border],  #n_sea_boundary
                         [east, bellambi_northing + border], #ne
                         [east, north_beach_northing - border],    #se
                         [sea_boundary,
                          north_beach_northing - border]    #s_sea_boundary
                         ]
    tags = {'wall':[0,1,2,3,4,5]}

    m.add_region_from_polygon(poly, tags, max_triangle_area=factor*1000)
    poly =  [
        [inner_west_south,north_beach_northing], #sw
        [inner_west_north,bellambi_northing],   #nw
        [inner_east,bellambi_northing], #ne
        [inner_east,north_beach_northing] #se
        ]
    m.add_region_from_polygon(poly, max_triangle_area=factor*100)
    

    #Even more inner points
    flagstaff_x = 308000
    flagstaff_y = 6189000
    west_x = 305000
    west_y = 6185000
    
    north_beach_northing = 6183000.0
    
    poly = [
                         [flagstaff_x,flagstaff_y], #n
                         [west_x,west_y],   #w
                         [north_beach_easting,north_beach_northing], #s
                         ]
    m.add_region_from_polygon(poly, max_triangle_area=factor*10)
    
    #m.generateMesh()
    m.generateMesh("pzq28.0za1000000a") 
    
    m.export_mesh_file("wollongong_4regions.msh")
    #m.export_mesh_file("wollongong_highres_a28.tsh")