TM 3 updates TM 3 new preprocessing chemical

TM 3 updates & TM 3 + new preprocessing (chemical forecast) & Diffusion module (TM 3/TM 5) & Photolysis module (TM 3/TM 5) TM meeting, IMAU, 10 June 2003 Michiel van Weele

TM 3 with new preprocessing (1) - HDF meteo fields TM 5 pre-processing - OD/ERA 40 data on 2 deg X 3 deg X 60 layers - Free to select (reduce) #levels within the code - Possibility to reduce horizontal resolution globally, e. g. , to 4 deg X 6 deg - Fluxes summed, other fields mass-averaged - Vertical diffusion coefficients calculated online - Surface fields flexible on 1 X 1 or 2 X 3 resolution - New Edgar emissions - Application: assimilation(, forecast, toy-model)

TM 3 with new preprocessing (2) 1. First experiences: 2. Online diffusion works fine (and is bit different) 3. cloud levels from entrain-/detrainments 4. Regular courant violations with 30 min timestep both in u and w directions => divide by 2 or 3 5. 3 day forecast 31 layers takes ~3 h (SUN) 6. severe problems with vertical transport

Resolution: VG 31 Old pre-processing New pre-procesing

Module Diffusion (1) subroutine bldiff subroutine calc_Kv subroutine pblhght subroutine difcoef Input: p, T, z, qv H, LE, ustar (not use surface stress) uwind, vwind (can use mass fluxes? ) Output: kvh, pblh Kvh = Vertical diffusion coefficient for heat Pblh = Planetary boundary layer height

Module Diffusion (2) subroutine calc_Kv (Holtslag and Boville, 1993) only shear production, mixing length stability functions in flux-Ri number full domain subroutine pblhght (Vogelezang and Holtslag, 1996) from bulk Ri number (surface heatflux, shear) subroutine difcoef (Beljaars and Viterbo, 1998) simila surface layer = 0. 1 * pblh Kvh=Kvm/Pr (Pr = Prandt number)

Kvh (on-line) vs Kzz (preprocessing) At the surface, 25 May 2003

Kvh (on-line) vs Kzz (preprocessing) 25 May 2003

Module Photolysis (work in progress) For tropospheric and stratospheric chemistry in TM 3 and TM 5 Use: psurf, temp. , clouds, ozone, surface albedo, (later: also aerosols, other trace gases) Absorption and scattering are separated - Absorption spectrally resolved (= exp. decay) - Scattering in limited #bands (with Landgraf scheme improved for low Sun) - Radiative transfer codes: LIDORT (van Oss) with 2, 4 or 6 streams Parameterization Krol and van Weele (fast)

Still to do with TM 3 Correct preprocessing forecasts Testing/improving diffusion module (Dirk? ) Writing of photolysis module Change TM 3 standard from VG/FG to 2 X 3/4 X 6 Use new surface fields, new EDGAR emissions Flexible chemistry module TM 3/TM 5