ModelTask1Settingupthebasestate.docx

Model Task 1: Setting up the base state ATM 562 Fall 2015 Fovell (see course notes, Chapter 9) 1 Overview Construct the base state (function of z alone) for five prognostic variables (u, w, q, qv, and p) and also r. The Weisman and Klemp (1982) sounding will be adopted. q and qv functions of z will be provided, and p and r will be computed. The grid will be staggered, using Arakawa’s “C” grid arrangement. Fake points above and below the model will facilitate handling of the boundary conditions. “C” grid arrangement (s = scalar) k+1 k+1/2 k k-1/2 k-1 NOTE: u(i,k), w(i,k) and s(i,k) not same point! ∆z ∆x Vertical grid Fortran The surface resides at the k = 2 level for w. k = 2 is also first real scalar level, so height of this level above is zT = (k-1.5)∆z, or 0.5∆z above ground C++ and other zero-based index languages The surface resides at the k = 1 level for w. k = 1 is also first real scalar level, so height of this level above is zT = (k-0.5)∆z, or (still) 0.5∆z above ground For this example problem, we take NZ = 40 and ∆z = 700 m W-K sounding Base state potential temperature (zT = scalar height [temperature] above ground; zTR = tropopause height above ground [12 km]; qTR = tropopause pot. temp. [343 K]; TTR = tropopause temp. [213 K]; g = 9.81 m/s2 ; cpd = 1004 J/kg/K). Note this is not qv. Base state water vapor mixing ratio can be specified as: Real and fake points For Fortran, the real points in the vertical for a scalar are k = 2, nz-1, with k=2 scalar level 0.5∆z above surface. Once we define mean potential temperature and mixing ratio (which I will call tb and qb) for the real points, we need to also fill in the fake points. Note the k=1 fake point is below the ground! We will presume the values 0.5∆z below the ground = those 0.5∆z above ground. That is, we assume zero gradient. With tb and qb, we can compute tbv, or mean virtual potential temperature, for all real and fake points. Derived quantities Given mean q, qv, we will compute the base sta