1 / 12

NDdiag – diagnostics on pressure levels from the Unified Model

NDdiag – diagnostics on pressure levels from the Unified Model. Peter Panagi Dept. Meteorology Univ. Reading 0118 378 8794 p.panagi@reading.ac.uk. What’s NDdiag for?. Re-maps the data onto constant pressure surfaces (user-specified)

wynn
Download Presentation

NDdiag – diagnostics on pressure levels from the Unified Model

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. NDdiag – diagnostics on pressure levels from the Unified Model Peter Panagi Dept. Meteorology Univ. Reading 0118 378 8794 p.panagi@reading.ac.uk UWERN Feb 13th 2004

  2. What’s NDdiag for? • Re-maps the data onto constant pressure surfaces (user-specified) • Computes a set of user-selected diagnostics (from a pre-defined set) on those pressure levels • Takes Unified Model data, on (height-based) model-levels, from the New Dynamics (V5.0 or higher) • Works with UK Mesoscale, other limited area and Global model domains UWERN Feb 13th 2004

  3. Whats’required • We MUST have the following fields UWERN Feb 13th 2004

  4. I/O Formats • Input and output data are in pp-format • Fields are stored sequentially as a function of model-level, eg u(1),u(2)….u(38) • Binary, read in as (Fortran-90) READ(file)IPP(45),RPP(19) NX = IPP(19) NY = IPP(18) READ(file)DATA(NX,NY) UWERN Feb 13th 2004

  5. Where to get data • Output from your own run of the UM • UWERN support programmer Chang-Gui Wang c.g.wang@reading.ac.uk • BADC http://badc.nerc.ac.uk UWERN Feb 13th 2004

  6. What does NDdiag calculate? • A large number of pre-defined (approx. 80) pressure-level diagnostics, which has evolved over time to meet requirements. These have been driven primarily by synoptic and mesoscale research, e.g. Geostrophic diagnostics, pot. temps., divergence, PV, Q-vectors, vorticity, frontal diagnostics, frontogenesis diagnostics, static stability, Richardson, q-g omega eq. inversion, shear, vertical gradients UWERN Feb 13th 2004

  7. How does NDdiag work? • Requires 2 files as input • NDdiag.CONTROL • NDdiag.INPUT • Generates as output • NDdiag.OUTPUT • LATLON.ASCII (optionally) UWERN Feb 13th 2004

  8. More details • NDdiag.CONTROL #NDdiag.CONTROL # Pressure levels 1000,50,50 #lat/lon of output grid – code as 0,0,0,0 for full domain 0,0,0,0 # linear recursive filter 0 #Required diagnostics 90001 geop ht m 90016 Dry PV 90041 theta-w K 90083 frontogenesis tilt K^2 m-2 s-1 UWERN Feb 13th 2004

  9. NDdiag structure • Fortran 90/95 BUT • ALLOCATABLE ARRAYS in SUBROUTINE argument list. This is an extension that will become part of the Fortran 2000 standard. • Uses MODULES, INTERFACES • Fixed name I/O files (to make my life easier!) UWERN Feb 13th 2004

  10. Compilation & Running • Compiles with • SUN Forte, NAG, HP • make NDdiag • But NOT PGI, or other Linux-based (such as on Newton) • Can be run from a script • NDdiag.scr <your pp-file> • Typical executable times • Mesoscale :1-2 mins./time step, for all diagnostics • Global : about 8 times longer UWERN Feb 13th 2004

  11. Getting hold of NDdiag • http://www.met.rdg.ac.uk/~panagi/NDdiag • This web page contains • NDdiag.tar.gz • NDdiag.ps • Instructions for unpacking, compiling and running • Continuous updates on modifications and bug fixes UWERN Feb 13th 2004

  12. History • MDIAG, written in f77, is still used by many people. It was written to work on output from the hydrostatic UM (V4.5 and previous) on hybrid model-levels. • Some diagnostics have not yet been re-written for NDdiag e.g. CAPE, some frontal • Available from http://www.met.rdg.ac.uk/~panagi/mdiag UWERN Feb 13th 2004

More Related