Changeset 3426


Ignore:
Timestamp:
Jul 31, 2006, 3:35:38 PM (19 years ago)
Author:
duncan
Message:

check in code that can compare jagged catchment boundary with Dirichlet boundary.

Location:
development/dam_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • development/dam_2006/create_mesh.py

    r3421 r3426  
    4848                      point_edge_bottom,
    4949                      point_dam_bottom]
     50    #if False:
     51    if False:
     52        dict['segments'] = [[0,1], [1,2], [2,3], [3,4],
     53                            [4,5], [5,6], [6,7 ],[7,0],  #The outer border
     54                            [2,7],         #dam Separator
     55                            [3,6]]         #edge Separator
     56   
     57        dict['segment_tags'] = ['wall',
     58                                'wall',
     59                                'wall',
     60                                'wall', #back of the dam
     61                                'wall',
     62                                'wall',
     63                                'wall',
     64                                'wall',
     65                                '',
     66                                '']           #Interior
    5067
    51     dict['segments'] = [[0,1], [1,2], [2,3], [3,4],
    52                         [4,5], [5,6], [6,7 ],[7,0],  #The outer border
    53                         [2,7],         #dam Separator
    54                         [3,6]]         #edge Separator
     68    else:
     69        dict['segments'] = [[0,1], [1,2], [2,3],
     70                            [6,7 ],[7,0],  #The outer border
     71                            [2,7],         #dam Separator
     72                            [3,6]]         #edge Separator
    5573   
    56     dict['segment_tags'] = ['wall',
    57                             'wall',
    58                             'wall',
    59                             'wall', #back of the dam
    60                             'wall',
    61                             'wall',
    62                             'wall',
    63                             'wall',
    64                             '',
    65                             '']           #Interior
     74        dict['segment_tags'] = ['wall',
     75                                'wall',
     76                                'wall',
     77                                'wall',
     78                                'wall',
     79                                '',           #Interior
     80                                'edge']           #Interior
    6681
    6782       

  • development/dam_2006/run_dam.py

    r3423 r3426  
    9999
    100100    Br = Reflective_boundary(domain)
    101     Bd = Dirichlet_boundary([tide,0,0])
    102 
    103     domain.set_boundary( {'wall': Br, 'wave': Br} )
     101    Bd = Dirichlet_boundary([0,0,0])
     102    #Bd = None
     103    domain.set_boundary( {'wall': Br, 'edge': Bd} )
    104104
    105105    #-------------------------------------------------------------------------
    … …  
    108108    t0 = time.time()
    109109
    110     for t in domain.evolve(yieldstep = 0.1, finaltime = 25):
     110    for t in domain.evolve(yieldstep = 0.01, finaltime = 1):
    111111        domain.write_time()
    112112   
Note: See TracChangeset for help on using the changeset viewer.