# Changeset 5212 for anuga_work/production/tonga/run_fangauta.py

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

update_fanguata

File:
1 edited

Unmodified
Removed
• ## anuga_work/production/tonga/run_fangauta.py

 r5194 import time import sys from math import radians # Related major packages from math import sin, pi Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/60), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/60)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/180), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/180)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=32400,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/300), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/300)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=43200,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/420), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/420)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=54001,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/540), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/540)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=64800,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/660), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/660)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=75600,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/780), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/780)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=86400,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/900), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/900)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=97200,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/1020), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/1020)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=108000,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/1140), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/1140)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time() for t in domain.evolve(10,finaltime=11800,skip_initial_step=True): Bw = Time_boundary(domain = domain, f=lambda t: [0.5*sin(t*2*pi/1260), 0, 0]) f=lambda t: [0.5*sin(t*radians(2*pi/1260)), 0, 0]) domain.set_boundary({'land':Bt,'mouth':Bw}) domain.write_time()
Note: See TracChangeset for help on using the changeset viewer.