Ignore:
Timestamp:
Oct 31, 2007, 4:56:02 PM (16 years ago)
Author:
duncan
Message:

Working on the anuga-2007 paper. Adding results from the UQ flume.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • anuga_validation/uq_flume_2007/run_dam.py

    r4750 r4773  
    11"""Script for running a dam break simulation of UQ's dam break tank.
     2
     3The simulation is based on the simulation Matt has presented and will
     4form part of the ANUGA validation paper.
     5
     6i want to compare the stage at .4m and the velocity at .45m.
    27
    38
    … …  
    3338def main():
    3439     
    35     slopes = [0]
    36     frictions = [0.01] #can enter multiple frictions eg [0.00, 0.01]
    37     inital_depths = [0.2] #can enter multiple initial depths eg [0.1, 0.2]
    38     gate_positions = [0.75]
     40    slope= 0
     41    friction = 0.01
     42    inital_depth = 0.2
     43    gate_position = 0.75
    3944   
    40     scenario(0, 0.01, 0.02, 0.75)
     45    return scenario(slope, friction, inital_depth, gate_position)
    4146           
    4247
    … …  
    4954    create_mesh.generate(project.mesh_filename,
    5055                         gate_position,
    51                          is_course=True) # this creates the mesh
    52                          #is_course=False) # this creates the mesh
     56                         #is_course=True) # this creates the mesh
     57                         is_course=False) # this creates the mesh
    5358
    5459    head,tail = path.split(project.mesh_filename)
    … …  
    101106        print 'finished'
    102107       
    103     points = [[0.4,0.2]]
     108    points = [[0.4,0.2], [0.45,0.2]]
    104109    #-------------------------------------------------------------------------
    105110    # Calculate gauge info
    … …  
    108113    interpolate_sww2csv( project.basename +".sww",
    109114                        points,
    110                         project.depth_filename + '.csv',
    111                         project.velocity_x_filename + '.csv',
    112                         project.velocity_y_filename + '.csv')
     115                        project.depth_filename ,
     116                        project.velocity_x_filename,
     117                        project.velocity_y_filename )
    113118
    114119
Note: See TracChangeset for help on using the changeset viewer.