Changeset 4044


Ignore:
Timestamp:
Nov 29, 2006, 2:26:23 PM (17 years ago)
Author:
duncan
Message:

make friction a parameter

File:
1 edited

Legend:

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

    r4033 r4044  
    3131
    3232
    33 def main():
     33def main(friction):
    3434
    3535    #--------------------------------------------------------------------------
     
    4747    myid = 0
    4848    numprocs = 1
    49     start_screen_catcher(project.outputtimedir, myid, numprocs)
     49    #start_screen_catcher(project.outputtimedir, myid, numprocs)
    5050
    5151    print 'USER:    ', project.user
     
    5454    #-------------------------------------------------------------------------
    5555
    56     gate_position = 2.0
     56    gate_position = 0.85
    5757    create_mesh.generate(project.mesh_filename,
    5858                         gate_position,
     
    9090       
    9191    domain.set_quantity('stage', elevation_tilt)
    92     domain.set_quantity('friction', 0.03)
     92    domain.set_quantity('friction', friction)
    9393    domain.set_quantity('elevation',elevation_tilt)
    9494
    9595    print 'Available boundary tags', domain.get_boundary_tags()
    96     domain.set_region('dam','stage',0.15,
     96    domain.set_region('dam','stage',0.20,
    9797                                 location = 'unique vertices')
    9898    #domain.set_region(Set_region('dam','stage',0.2,
     
    108108    t0 = time.time()
    109109
    110     for t in domain.evolve(yieldstep = 0.05, finaltime = 10):
     110    for t in domain.evolve(yieldstep = 0.01, finaltime = 30):
    111111        domain.write_time()
    112112   
     
    118118              [gate_position - 0.45,0.2],
    119119              [gate_position - 0.35,0.2],
    120               [gate_position - 0.25,0.2],
    121               [0.55,0.2]
     120              [gate_position - 0.25,0.2]
    122121              ]
    123122
     
    133132#-------------------------------------------------------------
    134133if __name__ == "__main__":
    135     main()
     134    main(0.000)
    136135   
Note: See TracChangeset for help on using the changeset viewer.