03.04.2013 Views

Advanced Research WRF (ARW) Technical Note - MMM - University ...

Advanced Research WRF (ARW) Technical Note - MMM - University ...

Advanced Research WRF (ARW) Technical Note - MMM - University ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!