1 | em_grav2d_x test case: |
---|
2 | |
---|
3 | The default version of this test case reproduces the |
---|
4 | 100 m resolution 2D gravity current. This test case |
---|
5 | was the focus of a workshop. The results of this workshop |
---|
6 | are published in |
---|
7 | |
---|
8 | Straka et al, NUMERICAL-SOLUTIONS OF A NONLINEAR DENSITY-CURRENT - A |
---|
9 | BENCHMARK SOLUTION AND COMPARISONS, |
---|
10 | INT J NUMER METH FL 17 (1): 1-22 JUL 15 1993 |
---|
11 | |
---|
12 | This test is 2d (x,z) and uses a fixed physical viscosity, hence the |
---|
13 | solution converges. It also uses a finite length domain with periodic |
---|
14 | conditions, so tests with non-zero horizontal translation are also |
---|
15 | appropriate. The reference shows tests are several different |
---|
16 | resolutions and these are quite illuminating |
---|
17 | |
---|
18 | There are a number of interesting features in this very non-linear |
---|
19 | solution (eddy structure, mins and maxes of temperature and velocity, |
---|
20 | propagation of the gravity current, etc.) |
---|
21 | |
---|
22 | Solution results from the WRF-ARW core can be found on the web at |
---|
23 | http://www.mmm.ucar.edu/individual/skamarock/test_cases/test_cases.html |
---|
24 | |
---|
25 | A number of versions test cases can be run using namelist.input files |
---|
26 | and input_sounding files provided in the test directory. |
---|
27 | |
---|
28 | Two input_sounding files are provided |
---|
29 | |
---|
30 | (1) input_sounding.um=0 |
---|
31 | (2) input_sounding.um=20 |
---|
32 | |
---|
33 | The first has no mean wind, the second has a constant horizontal |
---|
34 | environmental wind of 20 m/s. The no-mean wind case produces a |
---|
35 | symmetric gravity current, while the 20 m/s wind translates the |
---|
36 | current (see the web page for further details and examples). |
---|
37 | Copy the input_sounding.um=? into input_sounding to run the |
---|
38 | desired case. The input_sounding file is an ascii file, so |
---|
39 | it can be edited to change the mean wind speed to other values |
---|
40 | if desired. |
---|
41 | |
---|
42 | Three namelist.input files are in this directory and they provide |
---|
43 | for three different resolutions (dx, dz = 100, 200 and 400 m). |
---|
44 | Again, see the web page listed above for example solutions at |
---|
45 | the different resolutions. As with the sounding file, copy the |
---|
46 | namelist.input.??? into namelist.input to run the desired case. |
---|
47 | |
---|
48 | The default case is 100 m resolution and um=0. |
---|
49 | |
---|
50 | |
---|
51 | |
---|
52 | |
---|
53 | |
---|