Ignore:
Timestamp:
Apr 14, 2008, 3:20:21 PM (17 years ago)
Author:
herve
Message:

update_fanguata

File:
1 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/tonga/run_fangauta.py

    r5194 r5212  
    2222import time
    2323import sys
     24from math import radians
    2425
    2526# Related major packages
    … …  
    175176    from math import sin, pi
    176177    Bw = Time_boundary(domain = domain,
    177                                 f=lambda t: [0.5*sin(t*2*pi/60), 0, 0])
     178                                f=lambda t: [0.5*sin(t*radians(2*pi/60)), 0, 0])
    178179       
    179180    domain.set_boundary({'land':Bt,'mouth':Bw})
    … …  
    193194                           
    194195        Bw = Time_boundary(domain = domain,
    195                                         f=lambda t: [0.5*sin(t*2*pi/180), 0, 0])
     196                                        f=lambda t: [0.5*sin(t*radians(2*pi/180)), 0, 0])
    196197        domain.set_boundary({'land':Bt,'mouth':Bw})
    197198        domain.write_time()
    … …  
    201202    for t in domain.evolve(10,finaltime=32400,skip_initial_step=True):
    202203        Bw = Time_boundary(domain = domain,
    203                                         f=lambda t: [0.5*sin(t*2*pi/300), 0, 0])
     204                                        f=lambda t: [0.5*sin(t*radians(2*pi/300)), 0, 0])
    204205        domain.set_boundary({'land':Bt,'mouth':Bw})
    205206        domain.write_time()
    … …  
    209210    for t in domain.evolve(10,finaltime=43200,skip_initial_step=True):
    210211        Bw = Time_boundary(domain = domain,
    211                                         f=lambda t: [0.5*sin(t*2*pi/420), 0, 0])
     212                                        f=lambda t: [0.5*sin(t*radians(2*pi/420)), 0, 0])
    212213        domain.set_boundary({'land':Bt,'mouth':Bw})
    213214        domain.write_time()
    … …  
    217218    for t in domain.evolve(10,finaltime=54001,skip_initial_step=True):
    218219        Bw = Time_boundary(domain = domain,
    219                                         f=lambda t: [0.5*sin(t*2*pi/540), 0, 0])
     220                                        f=lambda t: [0.5*sin(t*radians(2*pi/540)), 0, 0])
    220221        domain.set_boundary({'land':Bt,'mouth':Bw})
    221222        domain.write_time()
    … …  
    225226    for t in domain.evolve(10,finaltime=64800,skip_initial_step=True):
    226227        Bw = Time_boundary(domain = domain,
    227                                         f=lambda t: [0.5*sin(t*2*pi/660), 0, 0])
     228                                        f=lambda t: [0.5*sin(t*radians(2*pi/660)), 0, 0])
    228229        domain.set_boundary({'land':Bt,'mouth':Bw})
    229230        domain.write_time()
    … …  
    233234    for t in domain.evolve(10,finaltime=75600,skip_initial_step=True):
    234235        Bw = Time_boundary(domain = domain,
    235                                         f=lambda t: [0.5*sin(t*2*pi/780), 0, 0])
     236                                        f=lambda t: [0.5*sin(t*radians(2*pi/780)), 0, 0])
    236237        domain.set_boundary({'land':Bt,'mouth':Bw})
    237238        domain.write_time()
    … …  
    241242    for t in domain.evolve(10,finaltime=86400,skip_initial_step=True):
    242243        Bw = Time_boundary(domain = domain,
    243                                         f=lambda t: [0.5*sin(t*2*pi/900), 0, 0])
     244                                        f=lambda t: [0.5*sin(t*radians(2*pi/900)), 0, 0])
    244245        domain.set_boundary({'land':Bt,'mouth':Bw})
    245246        domain.write_time()
    … …  
    249250    for t in domain.evolve(10,finaltime=97200,skip_initial_step=True):
    250251        Bw = Time_boundary(domain = domain,
    251                                         f=lambda t: [0.5*sin(t*2*pi/1020), 0, 0])
     252                                        f=lambda t: [0.5*sin(t*radians(2*pi/1020)), 0, 0])
    252253        domain.set_boundary({'land':Bt,'mouth':Bw})
    253254        domain.write_time()
    … …  
    257258    for t in domain.evolve(10,finaltime=108000,skip_initial_step=True):
    258259        Bw = Time_boundary(domain = domain,
    259                                         f=lambda t: [0.5*sin(t*2*pi/1140), 0, 0])
     260                                        f=lambda t: [0.5*sin(t*radians(2*pi/1140)), 0, 0])
    260261        domain.set_boundary({'land':Bt,'mouth':Bw})
    261262        domain.write_time()
    … …  
    264265    for t in domain.evolve(10,finaltime=11800,skip_initial_step=True):
    265266        Bw = Time_boundary(domain = domain,
    266                                         f=lambda t: [0.5*sin(t*2*pi/1260), 0, 0])
     267                                        f=lambda t: [0.5*sin(t*radians(2*pi/1260)), 0, 0])
    267268        domain.set_boundary({'land':Bt,'mouth':Bw})
    268269        domain.write_time()
Note: See TracChangeset for help on using the changeset viewer.