Description
This test case is provided for you to check your code before running the assigned problem and configuration for Computer Problem 5.Please note the settings below, and remember to change all settings back to the official problem(s) before running and handing in results.
- This is an advection test (only), for theta
- Grid layout: 201x41; ∆x=∆z=100 meters
- Quasi-compressible sound speed Cs = 100 meters/sec
- Potential temperature perturbation radius (x,z) = 4000, 1000
- Integration: time step ∆t = 0.25 sec
- You will disable
diffusion and pressure gradient/buoyancy forces for
these tests.
- These tests use Lax-Wendroff
advection.
- Split U test:
- Initial W, P=0. Initial U = +20m/s in bottom half of domain where k = [1:(nz-1)/2], -20 m/s on top half [(nz-1)/2+1:nz]
- Put +25K temperature perturbations at (x,z) = (10050, 50) and (10050, 4050) meters; plot.
- Run 1600 steps, plotting every 400. Theta
perturbation on top boundary moves left,
perturbation on bottom moves right.
- Split W test:
- Initial U, P=0. Initial W = +20m/s in left half of domain where i = [1:(nx-1)/2], -20 m/s on right half [(nx-1)/2+1:nx]
- Put +25K temperature perturbations at (x,z) = (50, 2050)
and (20050, 2050) meters; plot.
- Run 400 steps, plotting every 50. Theta
perturbation on left half should move up,
perturbation on right moves down.
Test case results
- Plots below use contour intervals of 2.0 (for U and W), 1.0 (for Theta'), and 25.0 (for P').
- Partial text output from these runs: Split
U, Split
W (links
fixed)
| Pgm
#5 test case B |
T=0 | T=100 | T=200 | T=300 | T=400 | |
|---|---|---|---|---|---|---|
| N=0 | N=400 | N=800 | N=1200 | N=1600 | ||
| Split U |
T |  |
 |
 |
 |
 |
| Pgm
#5 test case B |
T=0.0 | T=12.5 | T=25.0 | T=37.5 | T=50.0 | |
| N=0 | N=50 | N=100 | N=150 | N=200 | ||
| Split W |
T |  |
 |
 |
 |
 |