Ignore:
Timestamp:
Oct 27, 2006, 5:34:53 PM (18 years ago)
Author:
ole
Message:

Okushiri mesh using new mesh interface

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/okushiri_2005/extract_timeseries.py

    r3861 r3878  
    4545
    4646# old limiters
    47 expected_covariances = {'Boundary': 5.288601162783020386e-05,
    48                         'ch5': 1.167001054284431472e-04,
    49                         'ch7': 1.121474766904651861e-04,
    50                         'ch9': 1.249244820847215335e-04}
     47#expected_covariances = {'Boundary': 5.288601162783020386e-05,
     48#                        'ch5': 1.167001054284431472e-04,
     49#                        'ch7': 1.121474766904651861e-04,
     50#                        'ch9': 1.249244820847215335e-04}
     51#
     52#expected_differences = {'Boundary': 8.361144081847830638e-04,
     53#                        'ch5': 3.423673831653336816e-03,
     54#                        'ch7': 2.799962153549145211e-03,
     55#                        'ch9': 3.198560464876740433e-03}
    5156
    52 expected_differences = {'Boundary': 8.361144081847830638e-04,
    53                         'ch5': 3.423673831653336816e-03,
    54                         'ch7': 2.799962153549145211e-03,
    55                         'ch9': 3.198560464876740433e-03}
     57
     58expected_covariances = {'Boundary': 5.288392008865989237e-05,
     59                        'ch5': 1.166748190444680592e-04,
     60                        'ch7': 1.121816242516757850e-04,
     61                        'ch9': 1.249543278366777640e-04}
     62
     63expected_differences = {'Boundary':  8.373150808730501615e-04,
     64                        'ch5': 3.425914311580337875e-03,
     65                        'ch7': 2.802327594773105189e-03,
     66                        'ch9': 3.198733498646373370e-03}
    5667
    5768
     
    94105# Read and interpolate model output
    95106#--------------------------------------------------
     107
     108import sys
     109if len(sys.argv) > 1:
     110    sww_filename = sys.argv[1]
     111else:   
     112    sww_filename = project.output_filename
     113   
    96114#f = file_function('okushiri_new_limiters.sww',   #The best so far
    97115#f = file_function('okushiri_as2005_with_mxspd=0.1.sww',
    98 f = file_function(project.output_filename,
     116f = file_function(sww_filename,
    99117                  quantities='stage',
    100118                  interpolation_points=gauge_locations,
Note: See TracChangeset for help on using the changeset viewer.