1 | """Module where global euler model parameters are set |
---|
2 | """ |
---|
3 | |
---|
4 | epsilon = 1.0e-12 |
---|
5 | |
---|
6 | default_boundary_tag = 'exterior' |
---|
7 | |
---|
8 | |
---|
9 | time_format = '%d/%m/%y %H:%M:%S' |
---|
10 | |
---|
11 | min_timestep = 1.0e-6 #Should be computed based on geometry |
---|
12 | max_timestep = 1000 |
---|
13 | |
---|
14 | |
---|
15 | v_max = 100 #For use in domain_ext.c |
---|
16 | sound_speed = 500 |
---|
17 | |
---|
18 | max_smallsteps = 50 #Max number of degenerate steps allowed b4 trying first order |
---|
19 | |
---|
20 | |
---|
21 | |
---|
22 | #Betas [0;1] control the allowed steepness of gradient for second order |
---|
23 | #extrapolations. Values of 1 allow the steepes gradients while |
---|
24 | #lower values are more conservative. Values of 0 correspond to |
---|
25 | #1'st order extrapolations. |
---|
26 | # |
---|
27 | # Large values of beta_h may cause simulations to require more timesteps |
---|
28 | # as surface will 'hug' closer to the bed. |
---|
29 | # Small values of beta_h will make code faster, but one may experience |
---|
30 | # artificial momenta caused by discontinuities in water depths in |
---|
31 | # the presence of steep slopes. One example of this would be |
---|
32 | # stationary water 'lapping' upwards to a higher point on the coast. |
---|
33 | # |
---|
34 | # |
---|
35 | # |
---|
36 | #There are separate betas for the w-limiter and the h-limiter |
---|
37 | # |
---|
38 | # |
---|
39 | # |
---|
40 | # |
---|
41 | #Good values are: |
---|
42 | #beta_w = 0.9 |
---|
43 | #beta_h = 0.2 |
---|
44 | |
---|
45 | |
---|
46 | |
---|
47 | beta = 0.9 |
---|
48 | |
---|
49 | CFL = 1.0 #FIXME (ole): Is this in use yet?? |
---|
50 | #(Steve) yes, change domain.CFL to |
---|
51 | #make changes |
---|
52 | |
---|
53 | |
---|
54 | pmesh_filename = '.\\pmesh' |
---|
55 | |
---|
56 | |
---|
57 | import os, sys |
---|
58 | |
---|
59 | if sys.platform == 'win32': |
---|
60 | #default_datadir = 'C:\grohm_output' |
---|
61 | default_datadir = '.' |
---|
62 | else: |
---|
63 | #default_datadir = os.path.expanduser('~'+os.sep+'grohm_output') |
---|
64 | default_datadir = '.' |
---|
65 | |
---|
66 | |
---|
67 | use_extensions = True #Try to use C-extensions |
---|
68 | #use_extensions = False #Do not use C-extensions |
---|
69 | |
---|
70 | use_psyco = True #Use psyco optimisations |
---|
71 | #use_psyco = False #Do not use psyco optimisations |
---|
72 | |
---|
73 | |
---|
74 | optimised_gradient_limiter = True #Use hardwired gradient limiter |
---|
75 | |
---|
76 | #Specific to Euler |
---|
77 | minimum_allowed_density = 1.0e-3 #Water depth below which it is considered to be 0 |
---|