Changeset 8433


Ignore:
Timestamp:
May 18, 2012, 4:20:29 PM (12 years ago)
Author:
davies
Message:

Cleaning validation

Location:
trunk/anuga_core/validation_tests/Tests/Analytical
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/validation_tests/Tests/Analytical/deep_wave/plotme.py

    r8431 r8433  
    1818pyplot.ylim((-1.5,2.0))
    1919pyplot.legend()
     20pyplot.xlabel('Time (s)')
     21pyplot.ylabel('m')
    2022pyplot.savefig('wave_atten.png')
  • trunk/anuga_core/validation_tests/Tests/Analytical/deep_wave/results.tex

    r8431 r8433  
    99\begin{center}
    1010\includegraphics[width=0.9\textwidth]{wave_atten.png}
    11 \caption{Stage over time in the centre of the parabola}
     11\caption{Stage over time at 3 points in space}
    1212\label{stagewave}
    1313\end{center}
  • trunk/anuga_core/validation_tests/Tests/Analytical/parabolic_basin_1D/parabolaplot.py

    r8429 r8433  
    7373pyplot.plot(p2.time,p2.stage[:,ppp], color='green', label='numerical')
    7474pyplot.legend()
     75pyplot.xlabel('time (s)')
     76pyplot.ylabel('Stage (m)')
    7577pyplot.savefig('Stage_centre_v2.png')
    7678#       pyplot.savefig('runup_x_velocities.png')
     
    7981pyplot.plot(p2.time,p2.stage[:,ppp2], color='green', label='numerical')
    8082pyplot.legend(loc=10)
     83pyplot.xlabel('time (s)')
     84pyplot.ylabel('Stage (m)')
    8185pyplot.savefig('Stage_centre_v3.png')
    8286
     
    9094pyplot.plot(p2.x[v], -omega*A*numpy.sin(omega*p2.time[pltind])*(w>p2.elev[v]), label='Analytical')
    9195pyplot.legend()
     96pyplot.xlabel('x (m)')
     97pyplot.ylabel('x velocity (m/s)')
    9298pyplot.title('Velocity near to time=3T')
    9399pyplot.savefig('Vel_3T_v2.png')
     
    101107pyplot.plot(p2.x[v], -omega*A*numpy.sin(omega*p2.time[pltind])*(w>p2.elev[v]), label='Analytical')
    102108pyplot.legend()
     109pyplot.xlabel('x (m)')
     110pyplot.ylabel('x velocity (m/s)')
    103111pyplot.title('Velocity near to time=3.25T')
    104112pyplot.savefig('Vel_3_5T_v2.png')
  • trunk/anuga_core/validation_tests/Tests/Analytical/runup_sinusoid/runup_sinusoidplot.py

    r8412 r8433  
    4747pyplot.quiver(p1.x,p1.y,p1.xvel[t1,:],p1.yvel[t1,:])
    4848pyplot.title('The maximum VERTEX speed is '+ str(p1.vel[t1,:].max()) + ' m/s at time '+ str(p1.time[t1])+'s')
     49pyplot.xlabel('x (m)')
     50pyplot.ylabel('y (m)')
    4951pyplot.savefig('vel_sinu3s_v2.png')
    5052
     
    5557pyplot.quiver(p2.x,p2.y,p2.xvel[t1,:],p2.yvel[t1,:])
    5658pyplot.title('The maximum CENTROID speed is '+ str(p2.vel[t1,:].max()) + ' m/s at time ' + str(p1.time[t1]) + ' s')
     59pyplot.xlabel('x (m)')
     60pyplot.ylabel('y (m)')
    5761pyplot.savefig('vel_sinu3s_v2_cent.png')
    5862
     
    6367pyplot.quiver(p1.x,p1.y,p1.xvel[t2,:],p1.yvel[t2,:])
    6468pyplot.title('The maximum VERTEX speed is '+ str(p1.vel[t2,:].max()) + ' m/s at time ' + str(p1.time[t2]) +' s')
     69pyplot.xlabel('x (m)')
     70pyplot.ylabel('y (m)')
    6571pyplot.savefig('vel_sinu30s_v2.png')
    6672
     
    7177pyplot.quiver(p2.x,p2.y,p2.xvel[t2,:],p2.yvel[t2,:])
    7278pyplot.title('The maximum CENTROID speed is '+ str(p2.vel[t2,:].max()) + ' m/s at time' + str(p1.time[t2])+ ' s')
     79pyplot.xlabel('x (m)')
     80pyplot.ylabel('y (m)')
    7381pyplot.savefig('vel_sinu30s_v2_cent.png')
    7482
  • trunk/anuga_core/validation_tests/Tests/Analytical/trapezoidal_channel/plotme.py

    r8427 r8433  
    8383# But in my tests, they are not equal
    8484pyplot.legend( ('computed velocity', 'Analytical velocity', 'computed depth', 'analytical depth'), loc=4)
     85pyplot.xlabel('Down channel distance (m)')
     86pyplot.ylabel(' Generic scale (m or m/s)')
    8587pyplot.savefig('trapz_velocity_downstream_l0_eq_1_EL.png')
    8688
     
    101103pyplot.legend( ('computed velocity (m/s)', 'analytical velocity (m/s)', 'bed elevation (m)', 'computed stage (m)', 'analytical_stage (m)') ,loc=10)
    102104pyplot.title('Velocity (analytical and numerical) and Stage:' + '\n' +'Central channel regions (470 to 500m)' +'\n')
     105pyplot.xlabel('Cross-channel distance (m)')
     106pyplot.ylabel(' Generic scale (m or m/s)')
    103107pyplot.savefig('trapz_velocity_cross_channel_l0_eq_1_EL.png')
    104108
     
    120124pyplot.legend( ('computed velocity (m/s)', 'analytical velocity (m/s)', 'bed elevation (m)', 'computed stage (m)', 'analytical_stage (m)') ,loc=10)
    121125pyplot.title('Velocity (analytical and numerical) and Stage:' + '\n' +'Central channel regions (470 to 500m)' +'\n')
     126pyplot.xlabel('Cross-channel distance (m)')
     127pyplot.ylabel(' Generic scale (m or m/s)')
    122128pyplot.savefig('trapz_velocity_cross_channel_l0_eq_1b_EL.png')
    123129
     
    138144pyplot.legend( ('computed velocity (m/s)', 'analytical velocity (m/s)', 'bed elevation (m)', 'computed stage (m)', 'analytical_stage (m)') , loc=10)
    139145pyplot.title('Velocity (analytical and numerical) and Stage:' + '\n' +'Central channel regions (870 to 900m)' +'\n')
     146pyplot.xlabel('Cross-channel distance (m)')
     147pyplot.ylabel(' Generic scale (m or m/s)')
    140148pyplot.savefig('trapz_velocity_cross_channel_l0_eq_1c_EL.png')
Note: See TracChangeset for help on using the changeset viewer.