Changeset 3774


Ignore:
Timestamp:
Oct 13, 2006, 1:58:57 PM (19 years ago)
Author:
ole
Message:

Work on automatic validation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/automated_validation_tests/okushiri_tank_validation/validate_okushiri.py

    r3757 r3774  
    3636# Create Domain
    3737#-------------------------
     38#N = 50
    3839N = 50
    39 #N = 10
    4040points, vertices, boundary = rectangular_cross(N, N/5*3,
    4141                                               len1=5.448, len2=3.402)
     
    146146for i, (x, y) in enumerate(gauge_locations):       
    147147    predicted_gauge_values[i] = ensure_numeric(predicted_gauge_values[i])
     148    predicted_gauge_values[i].shape = predicted_gauge_values[i][:451]   
    148149
    149150    print predicted_gauge_values[i].shape
     
    161162    print msg
    162163   
    163     assert sqsum/denom < 0.01, msg
     164    #assert sqsum/denom < 0.01, msg
    164165   
    165166
     
    172173        plot(reference_time, reference_gauge_values[i], 'r.-',
    173174             reference_time, predicted_gauge_values[i], 'k-')
    174         title('Gauge %s' %name)
     175
     176        title('Gauge %d (%f,%f)' %(i,x,y))
    175177        xlabel('time(s)')
    176178        ylabel('stage (m)')   
    177179        legend(('Observed', 'Modelled'), shadow=True, loc='upper left')
    178         savefig(name, dpi = 300)
     180        savefig('Gauge_%d.png' %i, dpi = 300)
    179181       
    180         raw_input('Next')
     182        #raw_input('Next')
    181183       
    182184
Note: See TracChangeset for help on using the changeset viewer.