Changeset 4911


Ignore:
Timestamp:
Jan 8, 2008, 10:05:27 AM (16 years ago)
Author:
sexton
Message:

updates to plotting of bearing

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

    r4910 r4911  
    10211021    speed = zeros((n0,m,p), Float)
    10221022    bearings = zeros((n0,m,p), Float)
     1023    due_east = 90.0*ones(n0, Float)
     1024    due_west = 270.0*ones(n0, Float)
    10231025    depths = zeros((n0,m,p), Float)
    10241026    eastings = zeros((n0,m,p), Float)
     
    10621064                       +gaugeloc+'.csv'
    10631065            fid_out = open(thisfile, 'w')
    1064             s = 'Time, Stage, Momentum, Speed, Elevation, xmom, ymom \n'
     1066            s = 'Time, Stage, Momentum, Speed, Elevation, xmom, ymom, Bearing \n'
    10651067            fid_out.write(s)
    10661068            #### generate quantities #######
     
    10771079                    else:
    10781080                        vel = m / (depth + 1.e-6/depth)
    1079                     #bearing = calc_bearing(uh, vh)                   
     1081                    bearing = calc_bearing(uh, vh)                   
    10801082                    model_time[i,k,j] = (t + starttime)/scale #t/60.0
    10811083                    stages[i,k,j] = w
     
    10851087                    momenta[i,k,j] = m
    10861088                    speed[i,k,j] = vel
    1087                     #bearings[i,k,j] = bearing
     1089                    bearings[i,k,j] = bearing
    10881090                    depths[i,k,j] = depth
    10891091                    thisgauge = gauges[k]
    10901092                    eastings[i,k,j] = thisgauge[0]
    1091                     s = '%.2f, %.2f, %.2f, %.2f, %.2f, %.2f, %.2f\n' %(t, w, m, vel, z, uh, vh)
     1093                    s = '%.2f, %.2f, %.2f, %.2f, %.2f, %.2f, %.2f, %.2f,\n' %(t, w, m, vel, z, uh, vh, bearing)
    10921094                    fid_out.write(s)
    10931095                    if t == 0:
     
    12251227                        units = 'm / sec'
    12261228                    if which_quantity == 'bearing':
    1227                         due_east = 90.0*ones(shape(model_time[0:n[j]-1,k,j],Float))
    1228                         due_west = 270.0*ones(shape(model_time[0:n[j]-1,k,j],Float))
    1229                         plot(model_time[0:n[j]-1,k,j], bearings, '-',
    1230                              model_time[0:n[j]-1,k,j], due_west, '-.',
    1231                              model_time[0:n[j]-1,k,j], due_east, '-.')
     1229                        #due_east = 90.0*ones(shape(model_time[0:n[j]-1,k,j],Float))
     1230                        #due_west = 270.0*ones(shape(model_time[0:n[j]-1,k,j],Float))
     1231                        plot(model_time[0:n[j]-1,k,j], bearings[0:n[j]-1,k,j], '-',
     1232                             model_time[0:n[j]-1,k,j], due_west[0:n[j]-1], '-.',
     1233                             model_time[0:n[j]-1,k,j], due_east[0:n[j]-1], '-.')
    12321234                        units = 'degrees from North'
    12331235                        ax = axis([time_min, time_max, 0.0, 360.0])
Note: See TracChangeset for help on using the changeset viewer.