Changeset 4386


Ignore:
Timestamp:
Apr 17, 2007, 2:23:21 PM (17 years ago)
Author:
duncan
Message:

new UA friction scenario - elevation calculated by a function, initial stage of 0.06

Location:
anuga_work/development/friction_UA_flume_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/friction_UA_flume_2006/create_mesh.py

    r4351 r4386  
    77
    88
    9 def generate(mesh_filename, is_coarse = False):
     9def generate(mesh_filename, maximum_triangle_area=0.01):
    1010    """
    1111    Generate mesh for University of Aberbeen dam break flume.
     
    6868    # this is the location of the slope region.
    6969    slope.setTag("slope")
     70    #slope.setTag("dam")
    7071   
    71     if is_coarse:
    72         m.generate_mesh(maximum_triangle_area=0.01)
    73     else:
    74         m.generate_mesh(maximum_triangle_area=0.0001)
     72    m.generate_mesh(maximum_triangle_area=maximum_triangle_area)
    7573
    7674    m.export_mesh_file(mesh_filename)
  • anuga_work/development/friction_UA_flume_2006/parallel_run_dam.py

    r4351 r4386  
    6565        yieldstep = 1
    6666        finaltime = 15.1
     67        maximum_triangle_area=0.01
    6768    else:
    6869        yieldstep = 0.02
    6970        finaltime = 15.1
     71        maximum_triangle_area=0.0001
     72       
     73        maximum_triangle_area=0.001
    7074       
    7175
     
    107111    # Create the triangular mesh
    108112    #-------------------------------------------------------------------------
    109    
     113
     114    # this creates the mesh
    110115    #gate_position = 12.0
    111116    create_mesh.generate(mesh_filename,
    112                          is_coarse=is_trial_run) # this creates the mesh
     117                         maximum_triangle_area=maximum_triangle_area)
    113118
    114119    head,tail = path.split(mesh_filename)
     
    157162
    158163    for t in domain.evolve(yieldstep, finaltime):
     164   
    159165        domain.write_time()
    160    
    161166        print 'That took %.2f seconds' %(time.time()-t0)
    162167        print 'finished'
     
    199204                 0.019, 0.02,
    200205                 0.0225, 0.025, 0.03, 0.035, 0.04] # 17 processes
    201     #frictions = [0.0]
     206    frictions = [0.0, 0.014]
    202207    for i, friction in enumerate(frictions):
    203208        if i%size == rank:
    204209            print "I am processor %d of %d on node %s" %(rank, size, node)
    205             main(friction, is_trial_run = True, outputdir_name='friction_UA_new_meshII')
     210            main(friction, is_trial_run = True,
     211                 outputdir_name='friction_UA_new_meshII_0.06_stageIII')
Note: See TracChangeset for help on using the changeset viewer.