source: production/onslow_2006/make_gauges.py @ 3158

""" Make a line of gauges from outside boundary to near centre polygon
for Onslow study to look at MOST wave and compare to ANUGA output
"""
import project
from pylab import plot, xlabel, ylabel, title, ion, axis, savefig, close
from utilities.polygon import poly_xy
from Numeric import arange

x_bound, y_bound = poly_xy(project.polyAll)

# nominated point in centre region for Onslow
x2 = 305000.0
y2 = 7605000.0

d0 = project.d0
d1 = project.d1
d2 = project.d2
d3 = project.d3

d = [d0, d1, d2, d3]
ion()
fid = open(project.gauge_comparison,'w')
s = 'Easting,Northing,Name \n'
fid.write(s)
count = 0
for i, point in enumerate(d):
    x = point[0]
    y = point[1]
    m = (y-y2)/(x-x2)
    c = y2-m*x2
    if m > 0:
        plotx = arange(240000,340000,1000.0)
        for j, xpoint in enumerate(plotx):
            name = 'Point%d' %count
            count += 1
            s = '%.2f,%.2f,%s \n' %(xpoint,m*xpoint+c,name)
            fid.write(s)
    else:
        plotx = arange(240000,340000,5000.0)
        for j, xpoint in enumerate(plotx):
            name = 'Point%d' %count
            count += 1
            s = '%.2f,%.2f,%s \n' %(xpoint,m*xpoint+c,name)
            fid.write(s)
    plot(x_bound, y_bound, plotx,m*plotx+c,'+')

axis([240000,340000, 7580000, 7800000])
savefig('boundarycomparison')
close('all')