Changeset 5344


Ignore:
Timestamp:
May 20, 2008, 11:35:54 AM (16 years ago)
Author:
Leharne
Message:

Scaling the validation data to truescale

Location:
anuga_work/development/convergence_okushiri_2008
Files:
1 added
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/convergence_okushiri_2008/compare_timeseries.py

    r5343 r5344  
    2828#-------------------------
    2929
    30 finaltime = 22.5
    31 timestep = 0.05
     30finaltime = 450
     31timestep = 1
    3232
    33 gauge_locations = [[0.000, 1.696]] # Boundary gauge
    34 gauge_locations += [[4.521, 1.196],  [4.521, 1.696],  [4.521, 2.196]] #Ch 5-7-9
     33gauge_locations = [[0.000, 678.4]] # Boundary gauge
     34gauge_locations += [[1808.4, 478.4],  [1808.4, 678.4],  [1808.4, 878.4]] #Ch 5-7-9
    3535gauge_names = ['Boundary', 'ch5', 'ch7', 'ch9']
    3636
     
    4444#-------------------------
    4545
    46 print 'Reading', project.boundary_filename
    47 fid = NetCDFFile(project.boundary_filename, 'r')
     46print 'Reading', project_truescale.boundary_filename
     47fid = NetCDFFile(project_truescale.boundary_filename, 'r')
    4848input_time = fid.variables['time'][:]
    4949validation_data['Boundary'] = fid.variables['stage'][:]
    5050
    5151reference_time = []
    52 fid = open(project.validation_filename)
     52fid = open(project_truescale.validation_filename)
    5353lines = fid.readlines()
    5454fid.close()
     
    6262    for j, key in enumerate(gauge_names[1:]):  # Omit boundary gauge
    6363        value = float(fields[1:][j])           # Omit time
    64         validation_data[key].append(value/100) # Convert cm2m
     64        validation_data[key].append(value)
    6565
    6666
     
    8282    sww_filename = sys.argv[1]
    8383else:   
    84     sww_filename = project.output_filename
     84    sww_filename = project_truescale.output_filename
    8585   
    8686f = file_function(sww_filename,
     
    102102
    103103print 'Max runup elevation: ', q
    104 print 'Max runup elevation (scaled by 400): ', q*400
    105104print 'Max runup location:  ', loc
    106 
    107 #from create_okushiri import gulleys
    108 #assert is_inside_polygon(loc, gulleys)
    109105
    110106
Note: See TracChangeset for help on using the changeset viewer.