Ignore:
Timestamp:
Jul 25, 2008, 4:36:05 PM (16 years ago)
Author:
duncan
Message:

Current Hinwood scenario calculating RMSD's

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/Hinwood_2008/run_dam.py

    r5503 r5577  
    6262         outputdir_name=None,
    6363         run_type=0,
     64         width=1.0,
     65         use_limits=True,
    6466         end_tag = '_limiterD'):
    6567
     
    6870   
    6971    if run_type == 1:
    70         outputdir_name += '_test' + end_tag
    7172        yieldstep = 1.0
    7273        finaltime = 15.
    7374        maximum_triangle_area=0.1
     75        outputdir_name += '_test'
    7476       
    7577    elif run_type == 2:
    76         outputdir_name += '_test_long_time' + end_tag
    7778        yieldstep = 0.5
    7879        finaltime = None
    7980        maximum_triangle_area=0.01
     81        outputdir_name += '_test_long_time'
    8082       
    8183    elif run_type == 3:
    82         outputdir_name += '_yieldstep_0.1_tri_area_0.01' + end_tag
    8384        yieldstep = 0.1
    8485        finaltime = None       
    8586        maximum_triangle_area=0.01
     87        #outputdir_name += '_yieldstep_0.1'
     88       
    8689    elif run_type == 4:
    87         outputdir_name += '_good_tri_area_0.01' + end_tag
    8890        # this is not a test
    8991        # Output will go to a file
     
    9294        finaltime = None       
    9395        maximum_triangle_area=0.01
     96        #outputdir_name += '_good'
     97       
    9498    elif run_type == 5:
    95         outputdir_name += '_good_tri_area_0.001' + end_tag
    9699        # this is not a test
    97100        # Output will go to a file
     
    100103        finaltime = None       
    101104        maximum_triangle_area=0.001
    102      
     105        #outputdir_name += '_good'
     106       
     107    elif run_type == 6:
     108        # this is not a test
     109        # Output will go to a file
     110        # The sww file will be interpolated
     111        yieldstep = 0.01
     112        finaltime = None       
     113        maximum_triangle_area=0.0001
     114        #outputdir_name += '_good'
     115       
     116    if use_limits is True:
     117        outputdir_name += '_lmts'
     118    else:
     119        outputdir_name += '_nolmts'
     120    outputdir_name += '_wdth_' + str(width)
     121    outputdir_name += '_z_' + str(friction)
     122    outputdir_name += '_ys_' + str(yieldstep)
     123    outputdir_name += '_mta_' + str(maximum_triangle_area)
     124    outputdir_name += end_tag
     125   
    103126    metadata_dic = set_z_origin_to_water_depth(metadata_dic)   
    104127       
     
    150173    # this creates the mesh
    151174    #gate_position = 12.0
    152     create_mesh.generate(mesh_filename, metadata_dic,
     175    create_mesh.generate(mesh_filename, metadata_dic, width=width,
    153176                         maximum_triangle_area=maximum_triangle_area)
    154177
     
    172195    domain.set_minimum_storable_height(0.0001)
    173196
    174     domain.set_default_order(2) # Use second order spatial scheme
    175     domain.set_timestepping_method('rk2')
    176     domain.use_edge_limiter = True
    177     domain.tight_slope_limiters = True
    178    
    179     domain.beta_w      = 0.6
    180     domain.beta_uh     = 0.6
    181     domain.beta_vh     = 0.6
     197    if use_limits is True:
     198        domain.set_default_order(2) # Use second order spatial scheme
     199        domain.set_timestepping_method('rk2')
     200        domain.use_edge_limiter = True
     201        domain.tight_slope_limiters = True
     202       
     203        domain.beta_w      = 0.6
     204        domain.beta_uh     = 0.6
     205        domain.beta_vh     = 0.6
    182206   
    183207
     
    223247    print 'finished'
    224248
    225     flume_y_middle = 0.5
     249    flume_y_middle = 0.0
    226250    points = []
    227251    for gauge_x in metadata_dic['gauge_x']:
     
    259283    #from plot import plot
    260284
    261 
     285    # 1 is fast and dirty
    262286    # 4 is 0.01
    263287    # 5 is 0.001
     288    # 6 is 0.0001
     289   
     290    #run_type = 1
    264291    run_type = 5
    265292    #for run_data in [scenarios[5]]:
    266293    #scenarios = scenarios[2:]
    267294    #scenarios = [scenarios[0]]
     295    width = 1.0
     296    width = 0.1
    268297    for run_data in scenarios:
    269298        pro_instance = main( run_data['scenario_id'] + '_boundary.tsm'  ,
    270299                             run_data,
    271                              run_type = run_type,
     300                             width=width,
     301                             run_type=run_type,
    272302                             outputdir_name=run_data['scenario_id'],
    273                              end_tag='_limiterD')
     303                             use_limits=False,
     304                             friction=0.012,
     305                             end_tag='_G')
    274306        #gauges_for_slope(pro_instance.outputdir,[run_data])
Note: See TracChangeset for help on using the changeset viewer.