|
Last change
on this file since 8427 was
8427,
checked in by davies, 13 years ago
|
|
Adding the trapezoidal channel validation test, and editing the ANUGA manual
|
|
File size:
592 bytes
|
| Line | |
|---|
| 1 | from numpy import sqrt |
|---|
| 2 | from Numeric import zeros, Float |
|---|
| 3 | from config import g |
|---|
| 4 | from scipy import sin, cos, tan, arcsin, arccos, arctan, pi |
|---|
| 5 | |
|---|
| 6 | L = 200.0 |
|---|
| 7 | N = 400 |
|---|
| 8 | |
|---|
| 9 | cell_len = L/N |
|---|
| 10 | points = zeros(N+1, Float) |
|---|
| 11 | for i in range(N+1): |
|---|
| 12 | #points[i] = i*cell_len |
|---|
| 13 | points[i] = i*cell_len - 0.5*L |
|---|
| 14 | |
|---|
| 15 | |
|---|
| 16 | |
|---|
| 17 | friction_slope = 0.05 #tan(delta) #0.05 # NOTE THAT friction_slope must less than bed_slope |
|---|
| 18 | bed_slope = 0.1#0.002 #tan(theta) #0.1 # bottom slope, positive if it is increasing bottom. |
|---|
| 19 | thet = arctan(bed_slope) |
|---|
| 20 | |
|---|
| 21 | F2 = g*cos(thet)*cos(thet)*friction_slope |
|---|
| 22 | m = -1.0*g*bed_slope + F2 # auxiliary variable |
|---|
Note: See
TracBrowser
for help on using the repository browser.
