# source:anuga_core/source/anuga/config.py@3619

Last change on this file since 3619 was 3619, checked in by ole, 17 years ago

Config

File size: 3.1 KB
Line
1"""Module where global pyvolution model parameters are set
2"""
3
4
6from anuga_config import epsilon, default_boundary_tag
7
8
9
10#FIXME (Ole): More of these may need to be moved to anuga_config.py
11time_format = '%d/%m/%y %H:%M:%S'
12
13min_timestep = 1.0e-6 #Should be computed based on geometry
14max_timestep = 1.0e+3
15#This is how:
16#Define maximal possible speed in open water v_max, e.g. 500m/s (soundspeed?)
17#Then work out minimal internal distance in mesh r_min and set
18#min_timestep = r_min/v_max
19#
20#Max speeds are calculated in the flux function as
21#
22#lambda = v +/- sqrt(gh)
23#
24# so with 500 m/s, h ~ 500^2/g = 2500 m well out of the domain of the
25# shallow water wave equation
26#
27#The actual soundspeed can be as high as 1530m/s
28#(see http://staff.washington.edu/aganse/public.projects/clustering/clustering.html),
29#but that would only happen with h>225000m in this equation. Why ?
30#The maximal speed we specify is really related to the max speed
31#of surface pertubation
32#
33
34
35#v_max = 100 #For use in domain_ext.c
36sound_speed = 500
37
38
39max_smallsteps = 50  #Max number of degenerate steps allowed b4 trying first order
40
41manning = 0.3  #Manning's friction coefficient
42g = 9.81       #Gravity
43#g(phi) = 9780313 * (1 + 0.0053024 sin(phi)**2 - 0.000 0059 sin(2*phi)**2) micro m/s**2, where phi is the latitude
44#The 'official' average is 9.80665
45
46
47
48
49eta_w = 3.0e-3 #Wind stress coefficient
50rho_a = 1.2e-3 #Atmospheric density
51rho_w = 1023   #Fluid density [kg/m^3] (rho_w = 1023 for salt water)
52
53
54#Betas [0;1] control the allowed steepness of gradient for second order
55#extrapolations. Values of 1 allow the steepes gradients while
56#lower values are more conservative. Values of 0 correspond to
57#1'st order extrapolations.
58#
59# Large values of beta_h may cause simulations to require more timesteps
60# as surface will 'hug' closer to the bed.
61# Small values of beta_h will make code faster, but one may experience
62# artificial momenta caused by discontinuities in water depths in
63# the presence of steep slopes. One example of this would be
64# stationary water 'lapping' upwards to a higher point on the coast.
65#
66#
67#
68#There are separate betas for the w-limiter and the h-limiter
69#
70#
71#
72#
73#Good values are:
74#beta_w = 0.9
75#beta_h = 0.2
76
77
78
79beta_w = 0.9
80beta_h = 0.2
81CFL = 1.0  #FIXME (ole): Is this in use yet??
82           #(Steve) yes, change domain.CFL to
83           #make changes
84
85
86pmesh_filename = '.\\pmesh'
87
88
89import os, sys
90
91if sys.platform == 'win32':
93else:
95
96
97use_extensions = True    #Try to use C-extensions
98#use_extensions = False   #Do not use C-extensions
99
100use_psyco = True  #Use psyco optimisations
101#use_psyco = False  #Do not use psyco optimisations
102
103