Ignore:
Timestamp:
Apr 7, 2008, 3:55:51 PM (16 years ago)
Author:
duncan
Message:

near shore work, new limiters primary function, with build-up phase

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/near_shore_PMD/run_beach.py

    r5191 r5192  
    7777
    7878    # Structured mesh variables
    79     dx = 5.
     79    dx = 1.0
    8080    dy = dx
    8181    L = 2080.
     
    8787        outputdir_name += '_test'
    8888        yieldstep = 1
    89         finaltime = 100.
     89        finaltime = 200.
    9090        maximum_triangle_area=3000
    9191        thinner=True
    9292    else:
    9393        yieldstep = 1.
    94         finaltime = 1600
     94        finaltime = 1700 #1600
    9595        maximum_triangle_area = 100
    9696        thinner=True
     
    193193    Br = Reflective_boundary(domain)
    194194    Bd = Dirichlet_boundary([10.,0,0])
     195
     196   
     197    def wave_func(t):
     198        if t < 100:
     199            wave = [(1.2*t/100.0*sin(2*t*pi/10)), 0.0 ,0.0]
     200        else:
     201            wave = [(1.2*sin(2*t*pi/10)), 0.0 ,0.0]           
     202        return wave
     203   
    195204    Bwp = Transmissive_Momentum_Set_Stage_boundary(domain,
    196                    lambda t: [(1.2*sin(2*t*pi/10)), 0.0 ,0.0])
     205                   wave_func)
    197206    Bwp_velocity = Time_boundary(domain,
    198207                   lambda t: [(1.2*sin(2*t*pi/10)),
     
    206215                              'bottom': Br ,
    207216                              'left' : Br ,
    208                               'right' : Bwp } )
     217                              'right' : Bwp} )
    209218    else:
    210219        domain.set_boundary( {'wall': Br, 'wave': Bwp} )
     
    218227    for t in domain.evolve(yieldstep, finaltime):   
    219228        domain.write_time()
    220 
     229       
    221230    print 'That took %.2f seconds' %(time.time()-t0)
    222231    print 'finished'
     
    247256        friction=0.0,
    248257        is_structured=True,
    249         outputdir_name='structured_old_limiters')
     258        outputdir_name='structured_ramp_up2')
Note: See TracChangeset for help on using the changeset viewer.