Ignore:
Timestamp:
Jul 11, 2008, 4:15:59 PM (14 years ago)
Author:
duncan
Message:

Current Hinwood scenario - added plotting of froude number

File:
1 edited

Legend:

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

    r5461 r5494  
    5959         metadata_dic,
    6060         boundary_path=None,
    61          friction=0.01,
     61         friction=0.012,  # planed wood. http://www.lmnoeng.com/manningn.htm
    6262         outputdir_name=None,
    6363         run_type=0):
     
    6565   
    6666    basename = 'zz_' + metadata_dic['scenario_id']
     67    end_tag = '_limiterD'
    6768    if run_type == 1:
    68         outputdir_name += '_test_D'
     69        outputdir_name += '_test' + end_tag
    6970        yieldstep = 1.0
    7071        finaltime = 15.
     
    7273       
    7374    elif run_type == 2:
    74         outputdir_name += '_test_long_time'
     75        outputdir_name += '_test_long_time' + end_tag
    7576        yieldstep = 0.5
    7677        finaltime = None
     
    7879       
    7980    elif run_type == 3:
    80         outputdir_name += '_yieldstep_0.1_tri_area_0.01_D'
     81        outputdir_name += '_yieldstep_0.1_tri_area_0.01' + end_tag
    8182        yieldstep = 0.1
    8283        finaltime = None       
    8384        maximum_triangle_area=0.01
    8485    elif run_type == 4:
    85         outputdir_name += '_good_tri_area_0.01_D'
     86        outputdir_name += '_good_tri_area_0.01' + end_tag
    8687        # this is not a test
    8788        # Output will go to a file
     
    9192        maximum_triangle_area=0.01
    9293    elif run_type == 5:
    93         outputdir_name += '_good_tri_area_0.001_D'
     94        outputdir_name += '_good_tri_area_0.001' + end_tag
    9495        # this is not a test
    9596        # Output will go to a file
     
    117118    #return pro_instance
    118119    if finaltime is None:
    119         finaltime = boundary_final_time
     120        finaltime = boundary_final_time - 0.1 # Edge boundary problems
    120121    # Boundary file manipulation
    121122    if boundary_path is None:
     
    168169    domain.set_datadir(pro_instance.outputdir)
    169170    domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
    170     domain.set_minimum_storable_height(0.001)
    171     #domain.set_store_vertices_uniquely(True)  # for writting to sww
     171    domain.set_minimum_storable_height(0.0001)
     172
     173    domain.set_default_order(2) # Use second order spatial scheme
     174    domain.set_timestepping_method('rk2')
     175    domain.use_edge_limiter = True
     176    domain.tight_slope_limiters = True
     177   
     178    domain.beta_w      = 0.6
     179    domain.beta_uh     = 0.6
     180    domain.beta_vh     = 0.6
     181   
    172182
    173183    #-------------------------------------------------------------------------
     
    206216
    207217    # It seems that ANUGA can't handle a starttime that is >0.
    208     domain.starttime = 1.0
     218    #domain.starttime = 1.0 #!!! what was this doing?
    209219    for t in domain.evolve(yieldstep, finaltime):   
    210220        domain.write_time()
     
    230240                            pro_instance.outputdir + "velocity_x_" + id,
    231241                            pro_instance.outputdir + "velocity_y_" + id,
    232                             pro_instance.outputdir + "stage_" + id)
     242                            pro_instance.outputdir + "stage_" + id,
     243                            pro_instance.outputdir + "froude_" + id)
    233244 
    234245    return pro_instance
     
    259270                             run_type = run_type,
    260271                             outputdir_name=run_data['scenario_id'])
    261         gauges_for_slope(pro_instance.outputdir,[run_data])
     272        #gauges_for_slope(pro_instance.outputdir,[run_data])
Note: See TracChangeset for help on using the changeset viewer.