Changeset 1838


Ignore:
Timestamp:
Sep 16, 2005, 11:49:45 AM (18 years ago)
Author:
ole
Message:

Introduced caching in extract_timeseries.py

File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/validation/Completed/LWRU2/extract_timeseries.py

    r1835 r1838  
    44import sys
    55from os import sep
     6from caching import cache
     7
    68sys.path.append('..'+sep+'..'+sep)
    79
     
    5961#Read model output
    6062#filename = 'output.sww'
    61 filename = 'lwru2.sww'
    62 f = file_function(filename,
    63                   quantities = 'stage',
    64                   interpolation_points = gauges,
    65                   verbose = True)
     63#filename = 'lwru2.sww'
     64filename = 'lwru2_.05s.sww'
     65
     66#f = file_function(filename,
     67#                  quantities = 'stage',
     68#                  interpolation_points = gauges,
     69#                  verbose = True)
     70
     71f = cache(file_function,filename,
     72          {'quantities': 'stage',
     73           'interpolation_points': gauges,
     74           'verbose': True},
     75          verbose = True)
    6676
    6777
     
    99109    ion()
    100110    hold(False)
    101     plot(reference_time, gauge_values[k], '-k',
    102          reference_time, model, '-b')
     111    plot(reference_time, gauge_values[k], 'r.-',
     112         reference_time, model, 'k-')
    103113    title('Gauge %s' %name)
     114    xlabel('time(s)')
     115    ylabel('stage (m)')   
     116    legend(('Observed', 'Modelled'), shadow=True, loc='upper left')
    104117    savefig(name, dpi = 300)
    105118
Note: See TracChangeset for help on using the changeset viewer.