Changeset 3159


Ignore:
Timestamp:
Jun 15, 2006, 4:53:13 PM (18 years ago)
Author:
sexton
Message:

adding MOST output to onslow report

Location:
production/onslow_2006
Files:
31 added
3 edited

Legend:

Unmodified
Added
Removed

  • production/onslow_2006/make_gauges.py

    r3158 r3159  
    66from utilities.polygon import poly_xy
    77from Numeric import arange
     8from pmesh.create_mesh import convert_points_from_latlon_to_utm
    89
    910x_bound, y_bound = poly_xy(project.polyAll)
    1011
     12MOST_south = project.south
     13MOST_north = project.north
     14MOST_east = project.east
     15MOST_west = project.west
     16p0 = [MOST_south, MOST_west]
     17p1 = [MOST_south, MOST_east]
     18p2 = [MOST_north, MOST_east]
     19p3 = [MOST_north, MOST_west]
     20MOSTpoly, zone = convert_points_from_latlon_to_utm([p0, p1, p2, p3])
     21
     22MOSTymax = 0
     23MOSTymin = 7800000.0
     24MOSTxmax = 0
     25MOSTxmin = 7800000.0
     26for i, point in enumerate(MOSTpoly):
     27    x = point[0]
     28    y = point[1]
     29    if y > MOSTymax: MOSTymax = y
     30    if y < MOSTymin: MOSTymin = y
     31    if x > MOSTxmax: MOSTxmax = x
     32    if x < MOSTxmin: MOSTxmin = x
     33
     34MOSTymax = MOSTymax-3000.0
     35print 'MOST clipping: ', MOSTxmin, MOSTxmax, MOSTymin, MOSTymax
    1136# nominated point in centre region for Onslow
    1237x2 = 305000.0
    1338y2 = 7605000.0
     39testy = 7625000.0
    1440
    1541d0 = project.d0
    … …  
    2147ion()
    2248fid = open(project.gauge_comparison,'w')
    23 s = 'Easting,Northing,Name \n'
     49s = 'Easting,Northing,Name,Elevation \n'
    2450fid.write(s)
    2551count = 0
     52miny = 7580000.0
     53maxy = 7800000.0
     54xplot = []
     55yplot = []
    2656for i, point in enumerate(d):
    2757    x = point[0]
    … …  
    3060    c = y2-m*x2
    3161    if m > 0:
    32         plotx = arange(240000,340000,1000.0)
     62        plotx = arange(MOSTxmin,MOSTxmax,1000.0)
    3363        for j, xpoint in enumerate(plotx):
    34             name = 'Point%d' %count
    35             count += 1
    36             s = '%.2f,%.2f,%s \n' %(xpoint,m*xpoint+c,name)
    37             fid.write(s)
     64            ypoint = m*xpoint+c
     65            if ypoint > testy and ypoint <= d0[1]:
     66            #if ypoint >= testy and ypoint <= MOSTymax:
     67                if xpoint > MOSTxmin and xpoint < MOSTxmax:
     68                    xplot.append(xpoint)
     69                    yplot.append(ypoint)
     70                    name = 'Point%d' %count
     71                    count += 1
     72                    s = '%.2f,%.2f,%s,%.2f \n' %(xpoint,ypoint,name,0.0)
     73                    fid.write(s)
    3874    else:
    39         plotx = arange(240000,340000,5000.0)
     75        plotx = arange(MOSTxmin,MOSTxmax,5000.0)
    4076        for j, xpoint in enumerate(plotx):
    41             name = 'Point%d' %count
    42             count += 1
    43             s = '%.2f,%.2f,%s \n' %(xpoint,m*xpoint+c,name)
    44             fid.write(s)
    45     plot(x_bound, y_bound, plotx,m*plotx+c,'+')
     77            ypoint = m*xpoint+c
     78            if ypoint > testy and ypoint < MOSTymax: #testy
     79                if xpoint > MOSTxmin and xpoint < MOSTxmax:
     80                    xplot.append(xpoint)
     81                    yplot.append(ypoint)
     82                    name = 'Point%d' %count
     83                    count += 1
     84                    s = '%.2f,%.2f,%s,%.2f \n' %(xpoint,ypoint,name,0.0)
     85                    fid.write(s)
    4686
    47 axis([240000,340000, 7580000, 7800000])
     87    plot(x_bound, y_bound, xplot, yplot,'b+')
     88
     89axis([240000,340000, miny, maxy])
    4890savefig('boundarycomparison')
    4991close('all')

  • production/onslow_2006/report/computational_setup.tex

    r3158 r3159  
    7878wave amplitude and speed.
    7979
    80 \begin{figure}[hbt]
     80\input{MOST_input_onslow}
    8181
    82   \centerline{ \includegraphics[width=100mm, height=75mm]
    83               {../report_figures/boundarycomparison.png}}
    8482
    85   \caption{Point locations used to illustrate form of tsunami wave.}
    86   \label{fig:boundarypoints}
    87 \end{figure}
    8883
    89 %\input{latexoutput_MOST}
    90 
    91 {\bf run compare\_timeseries.py for some boundary gauges to show
    92 what the source provides to the boundary for this study.}
    93 
Note: See TracChangeset for help on using the changeset viewer.