Deutscher Wetterdienst
Status of the Support Activitiesand WG 6 Activities
COSMO Working Group 6: Reference Version and Implementation
Working Group Coordinator: Ulrich Schättler
06.09.2010 COSMO General Meeting 2010 2
Contents
Support Activities
Status of COSMO-Software
Status of the Web
Model Documentation
Working Group Activities
Grib2 Coding Issues
Source Code Management
SRNWP Interoperability Programme
Deutscher Wetterdienst
Status of COSMO Software
06.09.2010 COSMO General Meeting 2010 4
Status of COSMO Software
COSMO-Model versions
Version Implementation Date Status
4.9 15.07.09 released
4.10 11.09.09 tested
4.11 27.11.09 released
4.12 13.04.10 tested
4.13 31.05.10 released
4.14 14.06.10 in VCS (at DWD)
06.09.2010 COSMO General Meeting 2010 5
COSMO-Model 4.10 / 4.11 (released)
Dynamics:
Corrected use of cloud ice tendency from Tiedtke scheme in RK
Option for generalized Asselin-Filter
Option for preventing cold pools in steep valleys (ltadv_limiter) in combination with theta advection (itheta_adv=2) only
Option for selecting bottom boundary condition for vertical velocity
Option for reference atmosphere with constant Brunt-Väisälä frequency for idealized cases only
06.09.2010 COSMO General Meeting 2010 6
COSMO-Model 4.10 / 4.11 (released)
Physics
New sea-ice scheme
New multi-layer snow model
Changes to TKE-scheme
implicit vertical diffusion; 3D and horizontal correction of windshear production; separate horizontal shear mode; wake turbulence due to SSO; stability corrections for turbulent length scale
Radiation
New optional aerosol climatology
New external parameters for aerosol and surface emissivity
Soil-Model
New optional external parameter for minimum stomata resistance of plants
06.09.2010 COSMO General Meeting 2010 7
COSMO-Model 4.12 (tested)
Dynamics
Modifications to bottom boundary condition
Modifications to horizontal diffusion:
Removed the 2D mask field; diffusion factors for interior and boundary zone
Treat diffusion of pressure now separate from temperature
New routines to compute potential vorticity for output
Physics
Modifications to sea-ice scheme and multi-layer snow model
Latent Heat Nudging
Reintroduction of spatial quality function
Generalisation of the use of radar observation input frequency
More robust check in input of radar observations
06.09.2010 COSMO General Meeting 2010 8
COSMO-Model 4.13 (released) / 4.14 (in VCS)
4.13 (31.05.10)
Technical modifications to prepare for a new source code control system (SVN)
Introduction of an info-module which gives information about the binary, compiler options, etc, for every model run
Reconstruction of the Makefile: put Fopts, ObjFiles, ObjDependencies to separate files and keep extra files for every compiler in new directory LOCAL
4.14 (14.06.10)
Microphysics: Take density correction of terminal fall velocity into account. This leads to an increase of sedimentation velocity in the atmosphere
06.09.2010 COSMO General Meeting 2010 9
Status of COSMO Software
INT2LM versions
Version Implementation Date Status
1.9 03.03.09 released
1.10 16.12.09 released
1.11 23.04.10 in VCS (at DWD)
1.12 14.06.10 in VCS (at DWD)
06.09.2010 COSMO General Meeting 2010 10
INT2LM 1.10 (released)
Option lum2lm to process data from the Unified Model (UKMO)
Option to interpolate t_ice and h_ice from GME grid (only)
Option to compute FLake cold start variables
Decomposition independent computation of interpolation weights for reproducibility
Note: There still is a reproducibility problem for the „match“ interpolation of surface / soil variables for isolated land- or sea-points
06.09.2010 COSMO General Meeting 2010 11
INT2LM 1.11 / 1.12 (in VCS at DWD)
1.11 (23.04.10)
Technical modifications to prepare for a new source code control system (SVN)
Introduction of an info-module which gives information about the binary, compiler options, etc, for every model run
Reconstruction of the Makefile: put Fopts, ObjFiles, ObjDependencies to separate files and keep extra files for every compiler in new directory LOCAL
1.12 (14.06.10)
Possibility to switch the reference atmosphere in the option llm2lm
Vectorization of subroutine moist-split
06.09.2010 COSMO General Meeting 2010 12
Plans for the Near Future
Technical
Work Status
Reconstruction of mpe_io Prototype available
Implementation of ifdefs for external libraries not started
Choice of output format for every Namelist group not started
Implementation of grib_api started for INT2LM
New Physics Interface Prototype for Microphysics
Sem
i-Technical
COSMO-ART for Runge-Kutta under development
2-moment Microphysics under development
Generic treatment of idealized cases under development
Implementation of generic tracer variables needs coordination
Tile / Mosaic approach prototype available; definition of fluxes?
Full 3D strang splitting implemented at MCH; should be discussed in WG2
OASIS prototype for OASIS 3.0; new work for 4.0
Urban Module technical difficulties
Deutscher Wetterdienst
Status of COSMO Web
Deutscher Wetterdienst
The site now provides a
news feed for its updates
When the traditional “updates log” is
changed, subscribers get a notice of the new
items
06.09.2010 COSMO General Meeting 2010 15
New pages have been created for exchanging observation archives
both the browsing
and the downloading procedure is
web based
www.cosmo-model.org/srnwp/content
06.09.2010 COSMO General Meeting 2010 16
People interested for changes in PP taskLists can get informed fast;
When a PP page gets edited, a mail is send to the new cosmo-news mailist
06.09.2010 COSMO General Meeting 2010 17
The observation-archive browsing/downloading/uploading mechanism was extended to main site.
You can browse the site’s directories looking for non-web-page documents, going up and down the folders tree
www.cosmo-model.org/docView
06.09.2010 COSMO General Meeting 2010 18
www.cosmo-model.org/forum
New forum software installed
Possibilities for personal pages, discussions and distributed documentation
06.09.2010 COSMO General Meeting 2010 19
Conclusions
Although some contents are still too old, the site has seen some major changes
New features have to be communicated to the Users
We have to see how we can use the browsing/uploading/downloading and the Forum-pages
And (after all): The old contents have to be updated or removed
Deutscher Wetterdienst
Model Documentation
06.09.2010 COSMO General Meeting 2010 21
Documentation
Model Documentation
Dynamics and Numerics, Physics, Assimilation:
no updates have been done (at least none have been reported)
but: work in progress for Runge-Kutta Documentation
User Guides for INT2LM and the COSMO-Model have been updated to the latest versions in paper form. The web page tool has not been updated.
06.09.2010 COSMO General Meeting 2010 22
Documentation
COSMO Newsletter No. 10 has been published in January 2010 with (only) 10 contributions
But we had 6 new Technical Reports: TR 14 - 19
COSMO Priority Project "Further Developments of the Runge-Kutta Time Integration Scheme" (RK); Final Report
COSMO Priority Project "Tackle Deficiencies in Quantitative Precipitation Forecast" (QPF); Final Report
COSMO Priority Project "INTERP"; Final Report
Description of some convective indices, implemented in the COSMO model
Statistical Analysis of high-resolution COSMO Ensemble forecasts, in view of Data Assimilation
Seven years of activity in the field of mesoscale ensemble forecasting by the COSMO-LEPS system: main achievements and open challenges
Deutscher Wetterdienst
Other WG 6 Activities
Grib 2 Coding Issues
Source Code Management
06.09.2010 COSMO General Meeting 2010 24
Grib2 Coding Issues
Some COSMO centres have now decided to use grib_api from ECMWF as standard Grib library (DWD, MeteoSwiss, ARPA-SIM, …)
grib_api is already implemented into FieldExtra and work is in progress to implement it into INT2LM. Implementation into the COSMO-Model will follow.
Some issues have to be discussed (still ongoing):
Coding the vertical coordinate / grid of the COSMO-Model
What about the Reference Atmosphere
Coordination of usage of Grib2
06.09.2010 COSMO General Meeting 2010 25
Coding the Vertical Coordinate We code all these parameters as „vertical coordinate parameters“
ivctype Type of vertical coordinate (100: new reference atmosphere):
1: Pressure based hybrid coordinate
2: Height based hybrid coordinate
3: Height based hybrid SLEVE coordinate
ke Number of vertical levels
p0sl, t0sl, dt0lp, vcflat,
Parameters for the (old) reference atmosphere
delta_t, h_scal,
bv_ref
Additional parameters for the (new) reference atmosphere
svc1, svc2, nfltvc
Parameters for the SLEVE coordinates
vcoord(1:ke) Vertical coordinate parameters
06.09.2010 COSMO General Meeting 2010 26
Coding the Vertical Coordinate But only the list of vertical coordinate parameters should be there
Moreover, Grib only knows about „pairs of vertical coordinate parameters“ in combination with the surface pressure:
pressure = ak + bk * surface pressure (typeOfLevel = 119) or
height = ak + bk * surface height (typeOfLevel = 118; new in Grib2)
After the discussion we had in the last months, we proposed a new vertical coordinate type: typeOfLevel = 150
Do not specify vertical coordinate parameters together with a 2D field
but directly specify a 3D field (as a Grib message) with the height of all grid points
The list of vertical coordinate parameters can be used to specify additional information about the vertical grid:
NLAY: number of vertical layers (Note: # levels = # layers + 1)
NGRID: version number of vertical grid
06.09.2010 COSMO General Meeting 2010 27
Coding the Vertical Coordinate But what to do with the parameters for the reference atmospheres?
These parameters do not belong to the vertical grid and they are only necessary to interpret the field „pressure deviation from reference pressure“ (because they are needed to compute the reference pressure)
Using typeOfLevel = 150, we could still code them in the list of vertical coordinate parameters, because this is more or less „local use“ then.
But it is not nice
Are there different possibilities? Section 2, for example?
06.09.2010 COSMO General Meeting 2010 28
Source Code Management
COSMO-Standard for Source Code Development
Try to define clear rules for the development (design, coding rules, documentation, etc.)
And also give detailed rules for updating the Reference Version
This has been demanded several times
Now that the community is growing and other communities are also working with the COSMO-Model, this is inevitable
Discussion is going on and we will deliver a „Standard“ until the next COSMO User Seminar
Deutscher Wetterdienst
The SRNWP Interoperability Programme
06.09.2010 COSMO General Meeting 2010 30
Idea of Interoperability
To support future collaborative effort in European Limited Area Modelling by:
Defining a standard output format
Providing software tools (adaptors) to convert LAM output to the standard format
Enabling European LAMs to run from appropriate alternative model fields
To display or use in NMHS post-processing system
To provide operational backup
06.09.2010 COSMO General Meeting 2010 31
GME IFSUM ARPEGE
COSMO UM_local ALADIN HIRLAM
Fieldextra
Read (Write) Grib 2and transform to COSMO-Model grid:only 2D fields
For Parameter List 1:ps, psml, t2m, td2m, …
Plans for the Adaptor
06.09.2010 COSMO General Meeting 2010 32
GME ARPEGE
COSMO ALADIN HIRLAM
INT2LM
For Parameter List 2:u, v, w, t, p, qx, …
Read (Write) Grib 2and transform to COSMO-Model grid:2D fields, full atmosphere,soil parameters
to be donestatus unclear
GFS?GSM?
Plans for the Adaptor
06.09.2010 COSMO General Meeting 2010 33
lum2lm (Grib1): Temperature (lowest COSMO Level) from
GME IFS
UMAemet
UMDWD
06.09.2010 COSMO General Meeting 2010 34
lum2lm (Grib1): PP from reference pressure (lowest level) from
GME IFS
UMAemet
UMDWD
06.09.2010 COSMO General Meeting 2010 35
lum2lm (Grib1): Wind speed (lowest level) from
GME IFS
UMAemet
UMDWD
06.09.2010 COSMO General Meeting 2010 36
In December 2009, the “Interoperability Team” met in Reading for discussing last years work and do some planing for 2010.
The first 2 deliverables were finally agreed: The Standard Format shall be the native grid and native parameters Documentation supplied with the standard format should be the
minimum required Some Grib 2 Issues were discussed: There is a problem with coding the
vertical coordinate / vertical grid of COSMO in Grib 2 (see tomorrows session on Postprocessing)
Workshop in Reading
06.09.2010 COSMO General Meeting 2010 37
Work in 2010 Vertical coordinate / vertical grid of the COSMO-Model
A COSMO-internal discussion about coding the vertical grid in Grib2 has been initiated early 2010
For the Gal-Chen height based hybrid coordinate it is possible to code the grid with typeOfLevel=118, the corresponding orography (as a 2D field) and vertical coordinate parameters
This is NOT possible for e.g. the SLEVE coordinate
The COSMO discussion led to a proposal to the WMO Commission for Basic Systems (IPET-DRC: Inter-Programme Expert Team on Data Representation and Codes). A new vertical coordinate type is proposed, where no vertical coordinate parameters together with a 2D field are given, but a full 3D field with the heights of every grid point is specified.
The proposal has been discussed also with Enrico Fucile (ECMWF)
06.09.2010 COSMO General Meeting 2010 38
Work in 2010 grib_api
grib_api has been established at some COSMO centres to be the official Grib(2) library (DWD, MeteoSwiss, ARPA-SIM Bologna)
Work is in progress to improve performance on NEC machines (together with NEC, MeteoFrance and ECMWF)
It has been implemented in FieldExtra for reading and writing. FieldExtra therefore can now produce the parameters of the standard output format list directly
Work is in progress, to implement it into INT2LM
06.09.2010 COSMO General Meeting 2010 39
Work in 2010 New Grib 2 test data set and documentation
A new Grib2 test data set, which has been produced by FieldExtra, has been provided. This still has the old (Grib1-style) coding of the vertical coordinate
No Grib2 test data has been provided with typeOfLevel=118, because this might not be the default case any more (if proposal is accepted)
No surface test data (according to the guidelines from the ET on Surface Aspects) has been provided up to now
The minimum required level documentation has not been updated according to the plans from last workshop
describe native grids
describe fields in the standard output list: what field is it, what does it represent in the generating model
06.09.2010 COSMO General Meeting 2010 40
06.09.2010 COSMO General Meeting 2010 41
06.09.2010 COSMO General Meeting 2010 42
Thank you
very much
for your
attention