Advanced Research WRF (ARW) Technical Note - MMM - University ...
Advanced Research WRF (ARW) Technical Note - MMM - University ...
Advanced Research WRF (ARW) Technical Note - MMM - University ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Real-Data Lateral Boundary Condition: Location of Specified and Relaxation Zones<br />
Specified<br />
Rows<br />
West<br />
Relaxation<br />
Rows<br />
Specified<br />
Columns<br />
North<br />
N=1<br />
N=2<br />
N=3<br />
N=4<br />
N=5<br />
1 2 3 4 5 5 4 3 2 1<br />
N=5<br />
N=4<br />
N=3<br />
N=2<br />
N=1<br />
South<br />
Relaxation<br />
Columns<br />
Figure 6.1: Specified and relaxation zones for a grid with a single specified row and column,<br />
and four rows and columns for the relaxation zone. These are typical values used for a specified<br />
lateral boundary condition for a real-data case.<br />
of the specified zone is run-time configurable, but is typically set to 1 (i.e., the last row and<br />
column along the outer edge of the coarse grid). The second region of the lateral boundary for<br />
the coarse grid is the relaxation zone. The relaxation zone is where the model is nudged or<br />
relaxed towards the large-scale forecast (e.g., rows and columns 2 through 5 in Fig. 6.1). The<br />
size of the relaxation zone is a run-time option.<br />
The specified lateral boundary condition for the coarse grid requires an external file, generated<br />
during the same pre-processing as the initial condition file. Let ψ be any prognostic value<br />
having a lateral boundary entry, after Davies and Turner (1977),<br />
East<br />
∂tψ n = F1(ψLS − ψ) − F2∆ 2 (ψLS − ψ), (6.1)<br />
where n is the number of grid points in from the outer row or column along the boundary<br />
(SpecZone + 1 ≤ n ≤ SpecZone + RelaxZone − 1; see Fig. 6.1) and ψLS is the large-scale value<br />
obtained by spatial and temporal interpolation from the analyses. ∆ 2 is a 5-point horizontal<br />
smoother applied along η-surfaces. The weighting function coefficients F1 and F2 in (6.1) are<br />
given by<br />
F1 = 1<br />
10∆t<br />
F2 = 1<br />
50∆t<br />
SpecZone + RelaxZone − n<br />
,<br />
RelaxZone − 1<br />
SpecZone + RelaxZone − n<br />
,<br />
RelaxZone − 1<br />
where n extends only through the relaxation zone (SpecZone+1 ≤ n ≤ SpecZone+RelaxZone−<br />
1). F1 and F2 are linear ramping functions with a maximum at the first relaxation row or column<br />
nearest the coarse grid boundary (just inside of the specified zone).<br />
41