Last updated: 20 May 2014
The following table shows which developments have been or are implemented into cosmo and int2lm, per version.
Note: Changes from older to newer versions are documented in each "releaseNotes" field, as links to "Release Notes" documents. All these documents can also be collectively found.
Target version |
Version status | Development start |
Planned delivery |
Delivery date |
---|---|---|---|---|
6.00 | work | 2019-07-01 | 2019-12-31 | N/A |
COSMO-Model Version 6.0 will be the last and ultimate unified release of the COSMO-Model. It will contain all the loose ends, i.e. the developments which at the moment are ongoing and will be finished by the end of 2019. NOTE: Deadline for submitting contributions to COSMO-Model 6.0 is 31st of December 2019. |
||||
2015-10-31 | Parameterization of urban effects | WG3b | work | [CLM] H. Wouters (KU Leuven) |
Based on in-depth urban climate modeling research
Urban upgrade of TERRA-ML to TERRA-URB
Integration of TERRA-URB in COSMO
|
||||
'Offline' evaluations were performed for urban sites of Marseille, Toulouse and Basel COSMO-CLM/TERRA 'Online' evaluation for Flanders, Belgium Status
An action has started to incorporate this development in the official COSMO code (coordination U.Blahak) |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
The TERRA-URB scheme and its effects are documented in the following paper currently under review for GMD: Wouters, H., Demuzere, M., Blahak, U., Fortuniak, K., Maiheu, B., Camps, J., Tielemans, D., and van Lipzig, N. P. M.: Efficient urban canopy parametrization for atmospheric modelling: description and application with the COSMO-CLM model (version 5.0_clm6) for a Belgian Summer, Geosci. Model Dev. Discuss., doi:10.5194/gmd-2016-58, minor revisions, 2016. |
||||
2017-03-01 | Higher order discretization in the Runge-Kutta dynamical core | WG2 | work | Will, Ogaja (Univ. Cottbus) |
DescriptionThe current Runge-Kutta dynamical core will be extended in respect to the discretization order of the horizontal operations. The 5th order advection shall be replaced by a 4th order symmetric advection scheme (Morinishi et al ,1998, JCP). Similarly in the fast waves solver the horizontal currently 2nd order operations are replaced by 4th order discretizations. First tests with this new version show an improved behaviour of the kinetic energy spectrum. Furthermore climate runs over 20 years have been done in a stable manner without any artificial horizontal diffusion. One problem which must still be solved is an excessive underestimation of precipitation during summer with this new scheme. |
||||
StatusA first version based on COSMO 5.0 has been delivered in April 2016. However, it could not deal with the COSMO-DE setup. A fixed version has been delivered in Nov. 2016. |
||||
5.06 | refVersion | 2018-05-06 | 2019-03-29 | 2019-10-18 |
Two subversions have been implemented on the way to Version 5.06. There are release notes for every sub-version. This version 5.06 has been tested with the COSMO NWP Test Suite before the official release.
|
||||
2014-03-31 | read Meteosat radiances | KENDA | finish | Faulwetter (DWD) |
Description
Modification for V5.06:
|
||||
StatusRobin Faulwetter produces a new official version of the external library for satellite data. Based on this, Axel Hutt upgrades the adaptions of the COSMO code from V5.02 to V5.03, and then Uli Schättler continues the upgrade to V5.05 . |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2014-03-31 | radar obs operators + nudging of radial winds | KENDA | finish | Blahak / Stephan (DWD) |
Description
|
||||
StatusSome (further) code clean-up will be done by Uli Blahak before the code will enter V5.06. Later on (after the code is introduced in the official version), some (even) further code cleanup is planned, e.g. splitting of the large src_obs_radar.f90 into several modules. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding standards fulfilled (except for some in-line comments in German, which will be partly replaced by comments in English, but probably not all by the time of bringing it into V5.06 - note that the code consists of 20'000 code lines). 4-eyes assurance for much of the code (commonly written by Yuefeil Zeng, Dorit Jerger, Uli Blahak). The code has been applied by various users, and the latest bug report was two year ago. This shows that the code is very stable, indeed |
||||
Testing (Single Test Cases, Experiments)The operator has been thoroughly tested and already applied by various users and in various studies, incl. the peer-reviewed paper of Bick et al, 2016, QJRMS. The latest bug report was two year ago which shows that the code is very stable. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Scientific documentation:
Presentation of results: E.g.:
User's Guide: There are 3 parts (to be) written:
These descriptions has been found sufficient by the several users of the code to get started, apply the operator and do experiments.
|
||||
2014-06-30 | Extension of the Bott tracer advection scheme | WG2 | finish | Schneider (Univ. Bonn), Blahak |
DescriptionAn improvement (socalled deformational correction) of the Bott-schemes is currently done at the University Bonn (W. Schneider, A. Bott) to circumvent the necessity of the "true" Strang splitting (A. Bott (2010) AtmRes). The traditional strang splitting stretched over two timesteps is still necessary. A further improvement is the use of newer discretization coefficients for the gradient operators. The current COSMO implementation uses coefficients from a very old paper of Bott (1989). |
||||
StatusPreliminary version is available for 4th oder advection with the options "no flux limiter", "positive definite flux limiter" and "shape preserving" (Schneider, Bott), which runs efficiently on cache-based machines. Optimization of this new version for NEC SX9 has been done for the 4th order positive definite option (Vectorization) (Blahak). This new Version is currently under testing at DWD. Upgrade of this method to COSMO 5.4b1 has been done. A COSMO-DE experiment has been started. If successful, this new scheme might be a condidate for COSMO 5.07. |
||||
Technical IssuesThe code is the original code of Bott and Schneider taken from their own model and coupled separately to COSMO, totally different from the already implemented Bott schemes. Unfortunately, this implementation is not very efficient, because there the model fields are divided into an array for the interior domain and separate fields for the 4 lateral boundary zones, which necessitates the copying/mapping of the COSMO fields to these fields and back in every timestep for every tracer. Also, all grid lengths and velocities are normalized somehow. In case of successfull meteorological testing, this should be changed in that the method is re-coded in the COSMO "standard" way of coding the Bott-Schemes (if possible). This requires a new global array for summing up and storing the deformational correction term, which is subtracted again at the end of the advection algorithm. To be tested/checked further: Coding Standards, Technical Test Suite, 4-eyes Assurance |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2014-09-10 | GPU Version of COSMO | WG6 | finish | PP POMPA participants |
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2015-12-01 | Soil and surface fluxes: effective slope-angle correction factor (1/cos(slopeangle) | WG3a | submit | Blahak/De Morsier/Raschendorfe |
DescriptionTo account for the effect of an enhanced area for fluxes and their divergences at the surface and in the soil in case of sloping terrain (enlarged area but same grid box volume to heat up), an effective factor 1/cos(slope) is multiplied to the exchange coefficients tcm and tch. This factor has been mentioned in the PhD-Thesis of M. Buzzi (2008), based on Müller and Scherrer (1995), but it cancels out in the surface energy balance. The only terms where it does not cancel out (and this is not mentioned by Buzzi), are the fluxes of momentum, heat and moisture on the atmospheric side at the ground (transfer scheme). Within the atmosphere, where turbulent fluxes are assumed to be parallel to the gradient of the grid scale quantity (flux-gradient-relations), this factor is implicitly taken into account when taking into account the full 3D turbulent fluxes. In the soil, heat gradients are mainly perpendicular to the surface, so also here the factor has to be taken into account, whereas for soil water fluxes where gravity is an important factor, this is not so clear. Therefore, we include the factor in the heat conductivity, but exclude it for hydraulic conductivity at the moment. |
||||
StatusCode has been implemented into a test version based on COSMO 5.3 and is currently being tested. This measure belongs to the generalized boundary-layer-approximation already implemented in the surface-to-atmosphere transfer scheme (considering already the surface enlargement by sub-grid scale slopes) and should now be revied in order to implement the effect of grid scale slopes consistent to the planned implementation of "the vertically resoled roughness layer" and its vertical profiles of a surface area index. Alredy with the parameterization of a mixed roughness layer, which is being impemented in order to be coupled impicitly with the air above and the soil below, the general effect of the surface enlargement will be included. Thus a consistent representation of the mean surface slope effect will be implemented in connection with this deveopment. |
||||
Technical IssuesCoding Standards are fullfilled. Technical Test Suite has still to be run. 4-eyes principle has yet to be insured |
||||
TestingSingle Test Cases, Experiments: First tests with the present version and 3D-diffusion show only a rather small impact for the Alps only, mainly in terms of upward long-wave radiation. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2015-12-01 | Stochastic boundary layer perturbations reflecting unresolved effects of variability of surface fluxes on convection initiation | WG3a | submit | Kober / Blahak |
Description
In Short:
ttens = ttens + const * random_number * sqrt(tetl_tetl)
- random number seeding and -generation |
||||
StatusOriginal code provided by K. Kober, LMU Munich. Code technically re-written and implemented into Cosmo 5.3 test version. Modified random number seeding. Implemented random number field into restart. Has the potential to increase convective activity, but some more tuning seems to be necessary. Although first vefication reslults are not yet convincing, the code is going to be implemented, as it is already of value for ensemble methods and can be controlled by namelist switches. In a test-version, the neeed output of turbulence statistics has also been integrated into the new blocked SUB 'turbdiff'. Howeve, this additional modification needs to be introduced into the official code as well. |
||||
Technical IssuesCoding Standards fulfilled Technical Test Suite not yet 4-eyes Assurance not yet |
||||
TestingSingle Test Cases, Experiments First test simulations for convective periods do not yet provide a significant increase in precipitation. The results of a 1-month verification based on COSMO 5.2 (August 2015) are now yet convincing. |
||||
DocumentationPresentation of Results, Model Documentation, External Documentation Publication by Kirsin Kober et al: Kober, K., and C. Craig, 2016: Physically Based Stochastic Perturbations (PSP) in the Boundary Layer to Represent Uncertainty in Convective Initiation, J. Atmos. Sci., 73, 2893-2911. PPTx by Ulrich Blahak (Description, new namelist parameters, experiments) |
||||
2016-09-30 | Implementation of ICON/COSMO two-moment microphysics scheme | WG3a | work | Seifert/Blahak |
DescriptionThe two-moment microphysics scheme has been rewritten and strner modulaized for the ICON code and thus is already formulated in block-data stucture. This code being physcally identical with a previous implementation will now be portet to COSMO. |
||||
StatusCode is being committed to the code administrator and needs to be implmented into the next COSMO-version |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2017-03-15 | New communication interfaces | POMPA | work | Oliver Fuhrer (MeteoSwiss) |
DescriptionThe existing interface for halo-updates is not flexible enough in order to make it work on GPUs and allows for dynamically allocated fields to be exchanged. A new interface has been introduced in the code which is based on the GCL library (C++) and can handle halo-updates both on CPU and GPU. Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesThe work is finished and has been implemented in COSMO 5.0. |
||||
Technical IssuesCoding standard have been observed. Technical testsuite has been run using 3 compilers, 5 different systems and both in single and double precision, as well as in CPU and GPU mode. 4-eyes assurance has been done by Xavier Lapillonnne and Andrea Arteaga. |
||||
TestingResults are bit-reproducible with respect to reference version using the exchg_boundaries() interface on CPU. The new interfaces are being used regularly in COSMO-1 and COSMO-E and for the nightly builds on Jenkins. |
||||
DocumentationProcess documentation is done. No other documentation required, since the code modifications are within a #ifdef and the standard implementation is not being modified. |
||||
2017-05-15 | Data assimilation in SP | POMPA | work | Xavier Lapillonne (MeteoSwiss) |
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2017-06-01 | C++ Dycore for 5.06 | POMPA | work | Pascal Spoerri (C2SM) |
Description
|
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2017-06-15 | Physical parametrizations to GPU | POMPA | work | Xavier Lapillonne (MeteoSwiss) |
DescriptionUpdate SON 2016 : work on Turblence and sso (not completed) - Radiation (5.03) - Microphysics (5.06) - Turbulence (5.06) - Soil (5.06) - Shallow Convection (5.04b) - Tiedke Convection (5.06) - sso (5.06) - flake (5.06) Other not planned - sea ice |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2017-06-15 | LHN to GPU | POMPA | work | Xavier Lapillonne (MeteoSwiss) |
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2017-06-15 | Nudging to GPU | POMPA | work | Xavier Lapillonne (MeteoSwiss) |
DescriptionSome parts which can not be run in parallel have been modified for the GPU:
Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2017-06-15 | Output to GPU | POMPA | work | Xavier Lapillonne (MeteoSwiss |
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2017-06-15 | Static memory allocations | POMPA | work | Xavier Lapillonnne (MeteoSwiss |
DescriptionCost of memory allocation is very high on GPU, all arrays used on the GPU need to be allocated once before the time loop. In order to address this issue automatic arrays used in routine running on GPU will be made allocatable. An option will be provided to allocate these arrays on the fly at every time step so as to reduce memory usage when running on CPU. When running on GPU the allocation will be done once before the time loop. Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesThe developments are implemented in COSMO 5.0 |
||||
Technical IssuesThe coding standards have been observed, although static memory allocations may not be considered good coding style on a CPU-based system. Nevertheless, allocating/deallocating memory on every timestep can incurr an overhead also on traditional hardware. Technical testsuite has been passed using 3 compilers, 5 systems and in single/double precision. 4-eyes ensurance is done by Stefan Ruedisuehli and Christoph Schraff/Kaus Stephan (in the assimilation). |
||||
TestingThese modifications are being run regularly for COSMO-1. Results are bit-identical with respect to a version without modifications. |
||||
DocumentationNo documentation is required, purely technical change. Process documentation will be delivered with code changes to SCA. |
||||
2017-06-15 | New lateral BC module | POMPA | work | Oliver Fuhrer (MeteoSwiss) |
DescriptionThe lateral boundary conditions in the current model version are inconsistent, hard to understand and implemented differently throughout the code. See presentations at COSMO GM: ../../consortium/generalMeetings/general2014/wg6-pompa/wg6_lbc.pdf ../../consortium/generalMeetings/general2012/wg6-pompa/boundary_conditions.pdf Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesThe work is done and has been implemented in COSMO v5.0 |
||||
Technical IssuesThe coding standards have been observed. Technical testsuite has been passed using 3 different compilers, 5 different systems and both single/double precision as well as CPU/GPU execution. 4-eyes assurance has been done by Carlos Osuna and Anne Roches. |
||||
TestingVerification results between the pompa code and the offfical COSMO 5.0 have shown neutral effects of these changes. |
||||
DocumentationThe new interface has been presented at the COSMO GM's and has a very lightweight user-interface which is very easy to understand. The in code documentation explains what is done where directly in the code. The process documentation is available. |
||||
2017-06-15 | Port of lmorg and other "glue code" to GPU | POMPA | work | Xavier Lapillonne (MeteoSwiss) |
DescriptionPort using OpenACC |
||||
Technical IssuesWork done in COSMO v5.0 |
||||
Technical IssuesCoding standards have been observed. Technical testsuite has been run using 3 compilers, 5 different systems, in single/double precision and on CPU and GPU. 4-eyes assurance has been done by Stefan Ruedisuehli and Christoph Schraff |
||||
TestingResults are bit-reproducible. Version with modifications is being run with COSMO-1 regular runs and nightly with Jenkins. |
||||
DocumentationOn technical refactoring. No further documentation except process documentation required. |
||||
2017-09-01 | Mire parametrization | WG3b | finish | I.Rozinkina (RHM) |
Mires (organic soil) are extremely moist landscapes, which occupy relatively large areas in Scandinavia, North-West of European Russia and Siberia. Due to their unique physical and hydrological properties mires as well as lakes (Mironov et al., Boreal Environment Research, 2009) deserve a specific description in land surface schemes. The goal of this task is to incorporate a mire parameterization into the TERRA land surface scheme. It is planned to investigate the influence of mire parameterization on the components of the heat and water balance simulated by TERRA and compare them with available observations.
|
||||
Features of the mire hydrological and thermal regime have been introduced in the TERRA code and are being tested in production at RHM over Wester Siberia. There is now a final ICON branch for the Mire scheme that could in principle used also for COSMO. Final test will be made at DWD. |
||||
Juergen Helmert has seen the code. |
||||
Those modifications were tested with the single column model and the overall performance is reasonable and matches the observation. A verification of the 3D COSMO version on the basis of SYNOP stations in Siberia and North European Russia where mires are numerous has been done. These modifications are now tested in the production model at RHM. |
||||
A paper has been submitted to a special issue of 'Water Resources Research', this will be used for the scientific documentation of the model Input for the COSMO User Guide will be prepared by Inna |
||||
2019-02-27 | ensure reproducibility of redundancy check, adaptations for buoy and AMV reports | KENDA | finish | Schraff (DWD) |
Description
Changes of the namelist variables:
Changes of the results: None for the use of conventional data, unless 'lredn_repro' is switched on. |
||||
StatusAll changes implemented. All changes tested in short experimental trials, handed to CNMCA for successful testing of the adaptations of the buoy and AMV processing. |
||||
Technical IssuesCoding Standards: YES Technical Test Suite: By Uli. |
||||
TestingAll changes tested in short technical experimental trials. New option ensuring reproducibility avoids a crash of the LETKF if the ensemble members run with different domain decompositions. Adaptations for buoy and AMV data further tested at CNMCA where the new code allows successful processing of these data. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
5.05 | release | 2016-04-05 | 2017-04-30 | 2018-02-23 |
Version 5.05 is implemented in several substeps and includes the following sub-versions (with contents listed). There are release notes for every sub-version. In addition a summary of Namelist changes and changes to the results will be added shortly. This version 5.05 is now tested with the COSMO NWP Test Suite before an official release to the public can be done (expected in May 2018).
|
||||
2013-11-30 | Unified COSMO-ICON Physics | WG6 | finish | PP POMPA |
DescriptionIn the future, DWD wants to unify the physical packages between the COSMO-Model and the new global model ICON in a way, that if the same package is used from both, the COSMO-Model and the ICON, then also the same source code has to be used (which is at the moment not the case between COSMO-Model and GME, for example). For such a unification, several issues have to be addressed:
There are some other advantages of such an approach: When passing 2D structures to the parameterizations, the user can make the choice to pass rather few grid points (or short vectors), which is good for cache based machines, or to pass long vectors (good for vector processors) to the subroutines. All grid points have to be grouped to blocks then, which allows for an easy parallelization over these blocks (e.g. for OpenMP, but could also be good for implementation on GPUs). Work on the physical parameterizations has started to convert them to the correct memory layout. During this work, several issues occured at MeteoSwiss:
|
||||
Work done since 2011: MeteoSwiss ported most of the physical parameterizations (which are necessary to run COSMO-1) to GPUs using the block structure version. A framework was developed to copy the necessary COSMO fields from the (i,j,k) data structure to the blocked structure and back. This framework has been modified early 2014 to make its usage more comfortable and straightforward. But developments were also going on in the ICON physics, so that most packages have now diverged. The new ICON versions therefore have to be brought back to the COSMO-Model and tested there. Microphysics (Work done in 2014):
Radiation (Work done in 2015, May - July)
|
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)MCH developed the "copy-to-and-from block data structure" environment and most of the interfaces for the different parameterizations, when porting the code to the GPUs. The work has been followed and monitored by DWD. When updating the parameterizations to the new ICON versions, the still missing interfaces will be developed. |
||||
Testing (Single Test Cases, Experiments)Microphysics: The modifications to the microphysics have been tested in several experiments at DWD. A first experiment for COSMO-EU has already been done in Summer 2013 based on model version 4.26. In May / June 2014 more experiments have been conducted based on model version 5.0.2. Some basic verification results are summarized here. Radiation: Since results are bit-identical to old version, only the technical test suite has been run, which confirms the bit-identity. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Extensions for the documentation of the microphysics in the "Physical Parameterizations Documentation" is available and will be incorporated soon. |
||||
2014-06-30 | read / process GPS slant path delay data | KENDA | finish | Bender (DWD) |
DescriptionObservation operator for the use of GNSS zenith (ZTD) and slant total delay (STD) data in KENDA. This includes the writing of NetCDF feedback files by the COSMO model which can then be read by the LETKF. The current version can read two types of data:
In addition, a file is read containing the geoid undulation which is required to convert the heights above the geoid (used for the observations) into ellipsoidal heights (used by the model). The inclusion of the current version of this observation operator into V5.05 is important for DWD, because DWD requires to run this in passive mode as soon as possible in the pre-operational suite to prepare a meeting with GFZ delegates in December. |
||||
StatusThe code is ready and sent to Uli. |
||||
Technical IssuesCoding Standards are fulfilled, some technical tests have been performed (code was running stably in BACY for short period (few days) with STD and/or STD or without), the interface has been cross-checked briefly. |
||||
TestingSingle observation sensitivity experiments showed expected behaviour, and BACY trials using the full data set showed some positive impact of precipitation forecasts compared to a control with use of only conventional obs (both in comparisons with and without LHN). Some plots should be available here before end of August. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Availability of the following documentation required by the COSMO source code management rules:
|
||||
2014-07-31 | Providing COSMO-SC for general application | WG3a | submit | Raschendorfer / Schättler |
DescriptionThe single column test bed (COSMO-SC) had been developed years ago by M. Raschendorfer and consists of a couple of additional modules facilitating SC-runs with the parameterizations schemes used in the COSMO model. The framework provides a flexible configuration of a SC model run that can be forced by various SC measurement data (e.g. from meteorological observatories) in order to perform a kind of component testing of singular parameterization schemes. In order to perform SC runs, these modules need to be compiled with a sub-set of COSMO modules that branches from module organize_physics and includes also the related NAMELIST settings. As we want to produce a common library of physical parameterizations for COSMO and ICON, the SC framework automatically would be a test-bed for future ICON development. Nevertheless, the framework modules can also be linked to a subset of ICON modules in principal, what however would require some major adaptions. The last running SC configuration was based on COSMO version 4.24. In particular the new tracer module is not yet compatible to the SC framework, which involves some major modifications in these modules and some minor ones in the original COSMO code. We aim to prepare a running SC configuration based on version 5.0 and to include the modules of the SC framework to the official COSMO code, what means that the SC test-bed is available for all COSMO users, and that modifications of the original COSMO code need to guarantee the functionality of the SC runs. |
||||
StatusThe SC modules have been adapted to COSMO version 5.03_beta and a few related modifications have been implemented into the test version including the modifications related to "Turbulence: Implementing and extending the restructured ICON-version of TURBDIFF", "Implementing TKE-advection for TURBDIFF" and "Implementing TKE-advection for TURBDIFF". The standard workbench and makefile of Uli Schättler has been extended, so that the additional SC-code is always included now. Optionally either the paralles binary for the full 3D COSMO or a sequential binary for the COSMO-SC can be compiled. In particular all common physic schemes already implemented in COSMO can be tested by SC-runs, which includes component testing by means of an accordingly designed measurement forcing of SC-runs. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)The code of the additional SC-framework itself is rather comprehensive and its complete inspection is impossible. Moreover, this code is purely sequential and is not optimized with respect to performance, as it has been designed as a tool for testing and developing physical parameterizations only. However it written in rather modular and the code-lines are very inensively commented. It has been used intensively for various component testing experiments. |
||||
Testing (Single Test Cases, Experiments)The functionality of the SC-framework in combination with COSMO-5.0 has been tested by some users at DWD and I. Cerenzia at ARPA-SIMC. A previous version has been tested also by M. Buzzi (MeteoSwiss). |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)A documentation of the SC-framework is already available. However, it is not completely up to date. |
||||
2014-08-30 | Turbulence: Implementing and extending the restructured ICON-version of TURBDIFF | WG3a | finish | Raschendorfer |
DescriptionThe standard turbulence scheme of the COSMO (TURBDIFF) is the default scheme for ICON as well and a lot of adaptations have been introduced to the ICON-version including a major reorganization of the code.
This task aims to merge the different development into a common TURBDIFF module valid for COSMO and ICON:
|
||||
Status December 2016:
June 2016:
December 2015:
August 2014: After a first attempt has been introduced into the single column (SC) test-bed (COSMO-SC) based on version 4.29, which didn't contain all the different development completely and was still based on 2 horizontal dimensions of variable arrays,
Based on COSMO version 5.1a
Comments:
|
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)
|
||||
Testing (Single Test Cases, Experiments)
|
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)
|
||||
2015-08-31 | Revision of TERRA to support HWSD data | WG3b | finish | J.Helmert (DWD) |
The COSMO model uses operationally soil type data from the FAO Digital Soil Map of the World (DSMW) with a rather coarse resolution of 5 arc minutes compared to other external parameters (e.g., land use, orography). With the advent of the Harmonized World Soil Database (HWSD), a new global data set is available that provides a resolution of about 1 km (30 arc seconds by 30 arc seconds) and allows characterization of selected soil parameters (e.g., organic carbon, water storage capacity, soil depth, textural class and granulometry) [1]. Using recent database extensions in COSMO-CLM [2], the aim of the task is to use the HWSD in the EXTPAR system and adapt the soil model TERRA for the use of the high-resolution soil types. The benefit for the COSMO model is the advanced consideration of surface heterogeneities that can be used in future versions of TERRA. The implementation in the COSMO model uses pedotransfer functions to translate the fractions of sand, silt, clay and bulk density in physical properties used by TERRA. In the current implementation, only the shallow soil layer (0-30cm) is considered. This behaviour is controlled by the new namelist switch itype_soil (itype_soil=2). [1] http://www.iiasa.ac.at/Research/LUC/External-World-soil-database/HWSD_Documentation.pdf |
||||
The HWSD data are ready to be used as TERRA soil types. They are already used in the ART scheme for dust emission. |
||||
N/A |
||||
2015-09-30 | Consolidation of surface-to-atmosphere transfer ConSAT | WG3a | finish | Raschendorfer |
ConSAT is a coordinated and planned, sequential as well as partly interactive and dynamicaly adapted list of measures aiming to generally improve the description of Surface-to-Atmosphere Transfer (SAT), through which particularly the simulated diurnal cycle of near-surface variables can be strongly influenced. From 2013 till 2016 ConSAT was organized as sequence of PTs. More details about the aims, measures and results during this phase can be found in the related documents of the PTs. Currently ConSAT is continued at the level of WG tasks. The focus of ConSAT is the coupling between the rigid surface and the atmosphere and thus refers to the numerical simulation of atmospheric and non-atmospheric processes near the surface, which mainly are represented by the modules TURBDIFF (boundary layer turbulence), TURBTRAN (calculation of bulk atmospheric resistances between surface and the lowermost atmospheric model level) and TERRA (physical processes at the rigid surface). Since a consistent numerical description of this complex requires the consideration of each of these domains, ConSAT tasks necessarily are also dealing with surface processes being represented in TERRA. The work started in the framework of the COSMO-model, where the Single-Column framework (COSMO-SC) was an important working tool. Within 2017 the ConSAT-development moved to the ICON model, since COSMO could not realy be upgraded so as to contain the same treatment of surface processes (at least as far as TERRA is concerned). Furthermore, some related methodological development in ICON (such as surface tiles and dynamical snow tiles) have not been transferred to COSMO. Nevertheless, COSMO-SC will be used again, as soon as modules containing the related physical parameterizations can really be interchanged without further adaptations between COSMO and ICON. So far, the ConSAT-work takes place completely in ICON. ConSAT is related to the following implementation taks:
|
||||
StatusSome of the mentioned work (issues 1 to 6) have already been picked up in the development of the ICON-version of TURBDIFF (cf. "Implementing and extending the restructured ICON-version of TURBDIFF"). The full ICON-version of TURBDIFF is used now in COSMO exclusively (since version 5.04e). However, at DWD, TURBDIFF is still configured in order to provide similar results as compared to the previous TURBDIFF-version. This is realized by a couple of parameter settings. Most of them are activated by the NAMELIST-setting "loldturb=T", which can be used since COSMO-version 5.04f. July. 2019: In a particular development branch for ICON more recent development is being implemented and tested (see issue 3 and 8) Issue 1: An extended interpolation of the vertical profile function of turbulent velocity scale through the transfer layer (with a new hyperbolic branch for stable stratification) has been implemented in the common TURBDIFF version as well as an option to run the surface scheme without the input of turbulent diffusion coefficients from the lowermost model half level. The first implementation has been intensively inspected by Ines Cerencia based on SC-simulations and tower measurements and evidenced that the scheme is now also able to reproduce the strong decoupling of the homogeneous surface under strongly stable conditions, if minimal diffusion coefficients can be reduced accordingly. The second implementation needs some more inspection. Issue 2: Already included in the common TURBDIFF version. Not all measures are switched on operationally, even in ICON. A modified calculation of stability function is active, but the matrix conditioning seems to be dispensable. The new positive definete semi-implicit solution of the TKE-equation is not yet active. It promises to reduce some further security measures for stable stratification and will be tested in near future. A Ri-number dependent correction-factor for our minimal diffusion coefficients reduces their impact in the boundary layer and is already active in ICON. Issue 3: Already included in the common TURBDIFF as far as the separated horizontal shear interaction is concerned. In contrast to COSMO this interaction term is already active in ICON and produces more turbulence along frontal zones throughout the whole troposphere. The near surface interaction with non-turbulent thermal circulations is now dependent on the standard deviation of SSO and reduces nocturnal SAT for flat terrain. July. 2019: A new formulation of the circulation term is ready in a private development branch for ICON, which needs to be further tested and validated. Issue 4: Implemented in the common TURBDIFF and active in ICON. Issue 5: Already included in the common TURBDIFF version and runnin in ICON. In order to ease surface tiling, all input of the turbulence scheme from the transfer scheme schould be formulated in terms of fluxes in future. An according reformulation is planned for the next COSMO years. Issue 6: Is implemented in the common TURBDIFF version and running in ICON (as far as the application in the transfer layer is concerned). It remains to substitue the grid scale saturation adjustment by the turbulent adjustment, providing a possibly important additional heat source that is not yet considered. Issue 7: A generalized boundary layer approximation has been formulated and a related extension of the turbulence scheme has been developed. An implementation is planned for a medium time range. Issue 8: A first working version of the elevated roughness layer (canopy) has been implemented into the previous (non-blocked) version of TERRA and the new (blocked) version of TURBDIFF. In this version, however, a full implicit thermal coupling between roughness-elements, soil and (particularly) a fractional snow-cover is not yet considered. With this prototype the presence of snow is still excluded. The canopy-extension has ability to simulate the mean diurnal cycle of T_2m and Td_2m for a clear-sky summer day almost perfectly. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)All implementations of official code are part of the common TURBDIFF. Thus the related notes in task "Implementing and extending the restructured ICON-version of TURBDIFF" apply here as well. The prototypic canopy-extension is still present only in a test-version of COSMO! |
||||
Testing (Single Test Cases, Experiments)All implementations are being tested with the new common TURBDIFF version. All those being operational active in ICON have already been verified there. Mainly the extended formulation of vertical profile functions (Issue 1) the new positive definite solution of the TKE equation (Issue 2) and the prototypic canopy extension (Issue 8) have been tested by SC simulations and 3D test cases only. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)A full documentation of the TKE-based transfer scheme has started as well as the preparation of publications. |
||||
2015-12-01 | Shallow convection: new closure based on conv. vertical velocity scale, modified entrainment, modified shallow conv. cloud cover | WG3a | finish | Boeing / Blahak |
DescriptionModifications of the shallow convection parameterization based on the work of S. Boeing at MeteoSwiss. Cloud base muss flux now based on convective velocity scale instead of moisture convergence. Modified entrainment parameterization and modified shallow cloud cover. |
||||
StatusCode provided by S. Boeing and transferred from the GPU-implementation of MeteoSwiss to a Cosmo 5.3 test version by U. Blahak. New namelist switch y_conv_closure with the 2 possibilities "standard" and "Boeing" to activate the changes. |
||||
Technical IssuesCoding Standards are fullfilled, but additional nested IF branches are introduced within several outer loops over gridpoints in src_conv_shallow.f90 and make the code not very readable and perhaps also inefficient due to vectorization problems. Has to be re-checked. In cosmo-prerelease the SELECT CASE is now used and with OPENACC the performance of standard is the same. Technical Test Suite: Ok. 4-eyes Assurance: By Ulrich Blahak. |
||||
TestingExperiments made with COSMO-1 for 2 periods (Winter 2014-20159 and May/June 2015) show little impact. Results in Verification_shallow_convection.pdf |
||||
DocumentationPresentation of Results, Model Documentation, External Documentation A short description of by S. Boing in LaTeX is available. Talk at COSMO GM 2014: Evaluation of kilometer scale COSMO simulations using LES and observations. Latest results from Turb-i-sim. (p.19-33 of GMtkesv.pdf) External Doc: shcu.pdf |
||||
2016-01-01 | changes to SPPT scheme in order to improve stability | WG1 | finish | Walser (MeteoSwiss) |
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-06-15 | Code refactoring and bugfixes in Fortran for 5.5 | POMPA | finish | Pascal Spoerri |
DescriptionWhile porting the Fortran dynamical core to C++, some bugfixes or refactorings to improve clarity have been introduced in the Fortran dynamical core. Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-08-23 | technical options for LETKF/MEC: UV10M station selection; reader of DWD national stations; simple IAU; analysis as initial LBC; extrapolation of lateral boundary fields | KENDA | finish | Schraff (DWD) |
DescriptionAll changes are new options, or relate to options that are currently not used operationally:
|
||||
Status
|
||||
Technical IssuesCoding Standards: YES Technical Test Suite: By Uli. 4-eyes Assurance: By Uli for IAU / LBC |
||||
Testing(Single Test Cases, Experiments:) Short technical tests have been performed to make sure that the code run properly. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-08-23 | technical options for LETKF/MEC: UV10M station selection; reader Schraff (DWD) of DWD national stations; simple IAU; analysis as initial LBC; extrapolation of lateral boundary fields | KENDA | finish | Schraff (DWD) |
DescriptionAll changes are new options, or relate to options that are currently not used operationally:
|
||||
Status
|
||||
Technical IssuesCoding Standards: YES Technical Test Suite: By Uli. 4-eyes Assurance: By Uli for IAU / LBC |
||||
Testing(Single Test Cases, Experiments:) Short technical tests have been performed to make sure that the code run properly. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-09-15 | Block module | POMPA | finish | Xavier Lapillonne (MeteoSwiss) |
DescriptionBlock module need to be extended to handle tracers and fields of different types (int/bool) as well. |
||||
Technical IssuesWork is ongoing |
||||
Technical IssuesCoding standards will be respected. Technical testsuite will be applied using 3 compilers, 5 systems and in single/double precision. 4-eyes assurance will be done by Oliver Fuhrer |
||||
TestingResults are bit-identical with respect to a reference version. Testing in regular COSMO-1 runs of MeteoSwiss. |
||||
DocumentationOnly process documentation required. Block module is documented by a PowerPoint presentation. |
||||
2016-12-15 | C++ Dycore | POMPA | finish | Pascal Spoerri (C2SM) |
DescriptionUpdate for 12/2016: * New error handlig * Completed missing ftues * Dycore documentaton -- C++ dynamical core is integrated into the code via a library (libDycore.a). Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesThis development is ready for integration. |
||||
Technical IssuesCoding standard have been observed. Technical testsuite is being passed on 3 compilers, 5 systems and in both single/double precision every day. 4-eyes principle is assured by Xavier Lapillonnne and Andrea Arteaga. |
||||
TestingThis development does not change results of standard model setup (bit-reproducible). The results when lcpp_dycore=T are within rounding error as compared to the Fortran dycore (ensured via testsuite). The code runs on CPU and GPU since over one month and a standard verification will be done by MeteoSwiss beginning of January 2015. |
||||
DocumentationThe C++ dynamical core implements _exactly_ the same numerics as the Fortran dynamical core and thus the documentation already exists. |
||||
2016-12-15 | Change of application domain in relaxation | POMPA | finish | Oliver Fuhrer (MeteoSwiss) |
DescriptionRelaxation currently also sets lateral boundary conditions (BCs) on prognostic variables (time interpolated values from _bd fields). This is a bugfix since the dynamical core destroys the lateral BCs. If correct BCs are set in the dynamical core, this is no longer required. Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesThe work is done. Implemented in Version 5.0 of COSMO. |
||||
Technical IssuesCoding standards have been respected. If correct BCs are applied after the dynamical core (see dynamics refactoring task further up in the list), the results are bit-reproducible. The technical testsuite has been run using 3 different compilers, 5 different systems, in single/double precision, on CPU/GPU. The 4-eyes principle has been assured by Carlos Osuna. |
||||
TestingTests show an improvement in the fields at the boundary of the model domain. This is likely due to the changes in BCs in the dynamical core. Since results are bit-reproducible once correct BCs are applied after the dynamics, no further testing is required. The modified version of the code is used for the regular COSMO-1 runs which have been validated and proven to give neutral verification scores with respect to the reference version. |
||||
DocumentationProcess documentation is available. |
||||
5.04 | develop | 2015-09-30 | 2016-01-31 | 2016-03-10 |
The version 5.04 is a development version, which means that it is not distributed officially. There are no changing of results for the operational applications compared to version 5.03, but there are some changes in the Namelists. Take a look to the |
||||
2014-12-31 | Tracking and copying of boundary fields | POMPA | finish | Oliver Fuhrer (MeteoSwiss) |
DescriptionFields which are read from the BC files (lbff*) have to be tracked and copied to GPU Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesWork done and implemented in COSMO 5.0 |
||||
Technical IssuesCoding standard are observed. Technical testsuite passed using 3 different compilers, 5 different systems, in single/double precision, on CPU/GPU. 4-eyes assurance done by Xavier Lapillonne. |
||||
TestingResults are bit-reproducible. Version of code with modifications runs for COSMO-1 regular runs and nightly in the Jenkins suite of tests. |
||||
DocumentationSmall technical change. No further documentation except process documentation required. |
||||
2015-12-01 | DFI: Better treatment of clouds / precipitation | WG6 | finish | Lucio Torrisi |
DescriptionThe subroutine dfi_initializations.f90 has been modified to have a better treatment of clouds / precipitation during the diabatic DFI step. Filtering the qx-variables tends to smooth the structure of the initial state clouds. To mitigate this problem, two new namelist parameter have been introduced (working with the option ndfi=1):
For consistency with clouds / precipitation treatment also another namelist variable has been introduced:
If the default values are used, the results will not be changed. The recommended settings that were tested (results shown at COSMO-GM 2015) are:
|
||||
Technical IssuesCode has been modified and tested in Version 5.01 |
||||
Technical IssuesCoding Standards: fulfilled Technical Test Suite: 4-eyes Assurance: will be checked by the SCA |
||||
Testing |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-01-15 | removal of AOF reading | WG1 | finish | Schraff (DWD) |
Description
to be done when AOF reading is not used any more by any of the COSMO members / partners |
||||
Testing (Single Test Cases, Experiments)Fill information about testing |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-01-30 | Code refactoring and bugfixes in Fortran dycore 5.4 | POMPA | finish | Oliver Fuhrer (MeteoSwiss) |
DescriptionWhile porting the Fortran dynamical core to C++, some bugfixes or refactorings to improve clarity have been introduced in the Fortran dynamical core. Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesDevelopments are ready and have been implemented into 5.0. |
||||
Technical IssuesCoding standards have been observed. Technical testsuite passes using 3 compilers, 5 systems and in single/double precision. 4-eyes assurance has been done by Andrea Arteaga. |
||||
TestingThe modifications are used in regular COSMO-1 runs and visual feedback is that the boundaries of the model look much beter as compared to the standard version. Standard verificiation has been done by MeteoSwiss and scores are neutral. |
||||
DocumentationModifications have been discussed with Michael Baldauf and various occasions. Model documentation does not require modification. |
||||
2016-01-31 | Serialization | POMPA | finish | Oliver Fuhrer (MeteoSwiss) |
DescriptionIn order to validate the C++ dynamical core against the Fortran version and in order to ensure that the two versions are in-synch (and to rapidly identify code regions which are no longer in-synch) we need to introduce a small set of comments that help serialize data for unit-test validation. Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesFinished and ready for delivery to SCA. |
||||
Technical IssuesThese are only comments, so coding standard are fulfilled. No technical testsuite necessary, since only comments. 4-eyes assurance trivial and done by SCA when integrating the comments. |
||||
TestingNot required. |
||||
DocumentationChanges have been discussed and accepted by Ulrich Schaettler. |
||||
5.03 | release | 2015-06-01 | 2015-09-30 | 2015-12-02 |
The modifications in this version are grouped in the following categories and are documented in more detail in the Release Notes.
More information to special tasks can be found below. |
||||
2015-01-31 | Re-ordering of assimilation / relaxation | POMPA | finish | Oliver Fuhrer (MeteoSwiss) |
DescriptionRelaxation is part of C++ dynamical core. In the Fortran code the assimilation is done in between the dynamics and the relaxation. The reason for this is historical, since time-filtering – which is only used for LF dycore – is done in the same call as relaxation. Impact on code:
Impact on build process:
Impact on results:
Involved code owners (with effort):
|
||||
Technical IssuesWork done. Implemented in version 5.0. |
||||
Technical IssuesCoding standards have been observed. Technical testsuite has been run using 3 different compilers on 5 different systems. It has been run in single/double precision and in CPU as well as GPU mode. 4-eyes assurance has been done by Xavier Lapillonne. |
||||
TestingIt is not expected that switching the assimilation and relaxation leads to huge differences in the results, since both are simply relaxing the fields against a reference on each timestep and switching the order of operations is only a second order effect. The version with switched assimilation and relaxation is being used for the regular COSMO-1 runs at MeteoSwiss. Verification results between the pompa code and the offfical COSMO 5.0 have shown neutral effects of these changes. |
||||
DocumentationThe switch of order is only a very small change in lmorg.f90. Since the order of operations is not documented anywhere in the already existing documentation, there is no documentation to update. The process documentation is ready for delivery. |
||||
2015-07-01 | Increase of stability of the numerical diffusion scheme for 3D turbulence | WG2 | finish | Baldauf (DWD) |
DescriptionThe 1D vertical diffusion (i.e. our standard diffusion scheme for T, p, qx and the velocity components) is unconditionally stable by the vertically implicit treatment. However, the 3D extension for terrain-following coordinates suffers from numerical instability of the metric terms in tilted terrain as was pointed out by W. Langhans and O. Fuhrer. A recent stability analysis by M. Baldauf indicates that a significant increase in stability may be achieved by using more terms in the vertically implicit scheme and by use of some off-centering for them. 3D diffusion (or 3D turbulence) is not relevant for the currently operational model applications in the COSMO consortium with resolutions in the range 2 ... 14 km. It is probably even not relevant for dx ~1 km (e.g. COSMO-1 at MeteoCH). But it is expected that 3D turbulence becomes important for sub-km model runs. |
||||
StatusThe stability analysis for both the scalar and the vector diffusion is ready. A reviewed publication is under progress. The coding work has been finished and mainly consists in The validity of this new discretization has been tested by comparisons with exact analytic solutions over steep terrain (see 'Documentation'). The new code is based on the current version COSMO 5.2 and is ready for implementation in COSMO 5.3. |
||||
Technical Issuescode written according to standards: fulfilled 4-eyes assurance: by ??? responsible person: Michael Baldauf |
||||
TestingThe main testing consists in comparisons with analytically available solutions (see 'Documentation'). These tests have been successfully performed. Furthermore, single standalone runs with the operational COSMO-DE setup (2.8. km resolution) have been made. As expected, effects of 3D diffusion are quite small on this scale and will only be more prominent for much finer model resolutions. Nevertheless, the runs didn't show unusual behaviour of the new implementation. |
||||
DocumentationThe theory of 3D diffusion is documented in The new development is described in the presentation The newest test results, which show a very good agreement for the above mentioned code version, can be found in the document tests_of_3D_diffusion_Baldauf_v2.pdf. |
||||
2015-07-31 | Consolidation of GRIB2 and usage at other centers | WG6 | finish | Dörte Liermann, Ulrich |
DescriptionWith the earlier implementation of grib_api (use of local shortname concept, local use sections) it was not possible for centers /= DWD (WMO Code 78) to use GRIB2. Modifications have been implemented therefore, which now allow all centers to run with grib_api (in io_metadata.f90):
Also there have been changes to the grib2 definition files that have to be used together with grib_api. These changes are first implemented for the DWD definition files (definitions.edzw) of grib_api-1.13.1. local.250.def -> local.78.def |
||||
Technical IssuesThe modifications have been implemented in a pre-release of version 5.03, which has been given to MeteoSwiss for testing. This pre-release version has to work with grib_api-1.13.1 and the corresponding definition- and sample-files. |
||||
Technical IssuesCoding Standards: are fulfilled Technical Test Suite: runs every day now by Mr. Jenkins at MCH for a GRIB2 MCH case. 4-eyes Assurance: implemented by D. Liermann, checked by U. Schättler |
||||
TestingMeteoSwiss made tests with INT2LM and the COSMO-Model. Other centers can now also test with this version. |
||||
DocumentationA special COSMO Web Page for documenting GRIB2 and its implementation in INT2LM and the COSMO-Model is (still) under development. |
||||
2016-01-15 | Radiation to GPU | POMPA | finish | Xavier Lapillonne (MeteoSwiss) |
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-01-15 | Re-ordering of microphysics | POMPA | finish | Xavier Lapillonne (MeteoSwiss) |
DescriptionIn order to move all physical parametrizations into a single block (in order to increase memory locality) the microphysics should be shifted to the beginning of the timeloop. In order to avoid multiple copies of field used in the physics an optional compile switch (#idef GSP_FIRST) is introduced to call the microphysics at the beginning of the physics. If at all possible, an additional namelist switch (instead of the macro) would be preferable for testing purposes. |
||||
Technical IssuesDone. Implemented in COSMO v5.0 |
||||
Technical IssuesCoding standard have been respected. Technical testsuite has been run using 3 different compilers and 5 different systems. It has been run in single/double precision and on CPU/GPU. 4-eyes principle assured by Oliver Fuhrer |
||||
TestingThe switched order version is used for the regular runs of COSMO-1. The standard verification has been done by MeteoSwiss and has shown neutral results. |
||||
DocumentationNo further documentation than process documentation is required. The move of the microphysics is simply a cyclic permutation of the operations that are being done on a timestep of the model. |
||||
5.02 | develop | 2015-02-01 | 2015-05-01 | 2015-05-21 |
This version contains mainly bug fixes and an update of the code to process radiances. More details can be found in the Release Notes. This is only a DWD-internal development version, which will not be released. |
||||
5.01 | release | 2014-06-01 | 2014-09-30 | 2014-11-30 |
The following contributions have been implemented in COSMO-Model 5.1. More details can be found in the Release Notes
|
||||
2013-12-20 | GRIB2 and general vertical coordinates | WG6 | finish | U. Schaettler |
DescriptionINT2LM 2.0 and COSMO-Model 5.0 are now able to read and write model data with the GRIB2 format using the ECMWF grib_api. Another change we plan when going to GRIB2, is the introduction of the new general vertical coordinate. This brings some changes in the way we run our operational models:
There are now some ideas how these problems can be handled. These ideas were discussed within WG 6:
This is something we will implement and test here at DWD during November / December |
||||
Status of the work08. November:
15. November:
06. December:
14. April 2014:
12. May 2014
|
||||
Technical IssuesCoding Standards: Coding standards are fulfilled Technical Test Suite: the new versions passed technical tests. 4-eyes Assurance: The changes have been implemented by Uli Schaettler and have been cross-checked by DWD colleagues also involved in grib_api and GRIB2 implementation (Helmut Frank, et al.) |
||||
TestingThe new (DWD special) versions INT2LM 2.0.2 and COSMO-Model 5.0.2 are running in the parallel production suites. Verification results will be available end of May / beginning of June. There were 2 different parallel runs, where DWD special versions (COSMO-Model 5.0.1.1, 5.0.2.1 and INT2LM 2.0.1.1, 2.0.2.1) were tested. These versions contain the changes necessary for working with grib_api using GRIB1 and / or GRIB2.
The following pictures show the COSI Index for forecast days 1, 2 and 3: It can be seen that the two different runs are nearly identical. The differences that occur are due to different start times (and different cut-off times) for the operational and the parallel suite at DWD.
Again, the COSI index for forecast days 1, 2 and 3 is shown and it can be seen that there is no significant difference. |
||||
DocumentationA special web page for the COSMO Web is under construction. It will be extended in the next weeks. |
||||
2014-01-01 | Microphysics: Improved supercooled liquid water content | WG3a | finish | A. Seifert / U. Blahak |
DescriptionTo improve the forecast of aircraft icing, it is helpful that the microphysics schemes have some predictive skills concerning the super-cooled liquid water content in mixed-phase clouds. Currently, this is drastically underestimated in the operational COSMO 1-moment schemes (at the same time overestimating ice and snow content), and model based super-cooled liquid cannot be used for predictive purposes. The aim of this work is therefore to identify the cause(s) of the underestimation and, if possible, improve the situation, without degrading the precipitation forecasts. A better prediction of supercooled liquid can also influence the cloud radiative feedback, because water droplets have very different optical properties compared to ice particles. |
||||
Status The cause for underestimation of super-cooled liquid within the framework of the COSMO 1-moment schemes has been identified and it is a too efficient Bergeron-Findeisen-process. The reasons behind this are: - Too aggressive ice initiation at comparatively warm temperatures T > -10 deg C: IN are a simple function of temperature, and too much IN are assumed at these temperatures - Quite aggressive drop freezing, initiating even more ice. - Too coarse vertical resolution in operational setups within cloud layers, therefore no ability to simulate the often observed super-cooled liquid water layers at the top of stratus clouds. These problems could be mitigated by reducing the IN(T)-function at higher temperatures, by reducing the freezing probability function for droplets at higher temperatures, and by introducing the new parameterization of Richard Forbes (ECMWF) of super-cooled liquid layers at the top of stratus clouds. |
||||
Technical IssuesMethodology and code has been developed by F. Rieper (DWD) and is available in a test version. The code has been cross-checked by Carmen Köhler and has been implemented in Version 5.01 |
||||
TestingThe effects of the changes have been tested by approx. 20 case studies using COSMO-EU and observational data (retrievals of cloud phase and SLWC) from the Lindenberg cloud-net facility. All cases showed more or less improvement of SLWC in the model, some of them very clear. Overall, the changes were beneficial. Month-long experiments (Jan. 2013, June 2013) have been conducted and showed somewhat positive impacts (reduction of RMS) in the winter period for pressure and temperature, and slight improvements of precipitation scores. The summer period (convective precipitation) behaved more or less neutral in COSMO-EU. Note that there have been no model changes in the convection parameterization! All in all, the long-term tests demonstrated that there are no negative side-effects on the operational COSMO-EU NWP at DWD, on the contrary, there is an improvement especially in winter time. |
||||
DocumentationThe methodology and the results of the case studies and long-term tests have been presented at the COSMO General Meeting 2013. The presentation is available from Ulrich Blahak (DWD). Model Documentation: has to be adapted. External Documentation: the parameterization of Richard Forbes is documented in a journal paper. The implementation and testing of the supercooled liquid water representation in COSMO-EU has been submitted by Felix Rieper and Uli Görsdorf (DWD). |
||||
2014-02-05 | Single precision version of COSMO (WG6) | POMPA | finish | Oliver Fuhrer |
Description
|
||||
StatusThe single precision version is on the desk of the SCA and is pending for introduction into model version 5.1
|
||||
Technical IssuesThe code modifications are documented in the COSMO Newsletter 14 (Chapter 5.1).
|
||||
TestingThe results of all testing is documented in the COSMO Newsletter 14 (Chapter 5.1). Difference in double precision mode are only allowed within rounding errors of the calculations. But they are expected due to reordering/rewriting of certain operations to be safer in single precision mode. Results in single precision mode have been be verified extensively on a module-per-module basis and for real-case applications (3 months of COSMO-7 and COSMO-E). It is not guaranteed that single precision will work for all model setups and users and each user will have to carefully consider the precision requirements of his/her application. |
||||
DocumentationPresentation at the General Meeting in Sibiu (2013) and at previous POMPA parallel sessions and COSMO User Seminars. The Article in the COSMO Newsletter n°14 documents how to use the single precision version for users who are interested. For users which will remain with double precision, there are no implications. |
||||
2014-02-05 | Online Trajectory Module | WG6 | finish | Annette Miltenberger (ETH) et |
DescriptionA module to calculate online trajectories has been implemented into the nonhydrostatic limited-area weather prediction and climate model COSMO. Whereas offline trajectories are calculated with wind fields from model output, which is typically available every one to six hours, online trajectories use the simulated resolved wind field at every model time step (typically less than a minute) to solve the trajectory equation. As a consequence, online trajectories much better capture the short-term temporal fluctuations of the wind field, which is particularly important for mesoscale flows near topography and convective clouds, and they do not suffer from temporal interpolation errors between model output times. |
||||
Status21. November: the WG6 coordinator has been informed about the ongoing activity 28. November: the proposed activity has been accepted by TAG 11. December: presentation to SMC |
||||
Technical Issues
|
||||
TestingAs a first application of the new COSMO-model module, an Alpine north foehn event in summer 1987 has been simulated with horizontal resolutions of 2.2, 7 and 14 km. |
||||
Documentation
|
||||
2014-03-31 | Optional targeted diffusion to avoid cold pools in narrow valleys | WG2 | finish | Baldauf (DWD) |
DescriptionAround 06. Dez. 2013, a strong and unrealistic cool down in a narrow Alpine valley occured in the COSMO-DE model setup at DWD. (this caused even a model crash in one of the COSMO-DE ensemble members in the parallel routine) Similar problems have been reported by MeteoSwiss, too. The upwind advection operator of 5th order causes such cold pools. There exist several solutions for this problem: 1.) use the 3rd order advection operator (but this probably not accurate enough for convection-permitting runs (?)) 2.) It seems to occur less frequently with the dynamical bottom boundary condition (ldyn_bbc=.TRUE.) (but the reason is not clear) 3.) avoid this undershooting by a targeted diffusion like in ICON (G. Zaengl) (which exactly cures the artificial behaviour of the advection operator) 4.) the current available limiter in the horizontal advection operator could be used, too. However, in the current form it takes more computation time and (as reported by MeteoCH) it cannot cure every cold pool problem. The 3rd item in fact cures the problem and changes only the few grid cells where a cold pool occurs. It is highly a efficient method (acts at most in the lowest 3 layers) This targeted diffusion may have a very small impact to results. However, in most simulations, no critical grid cell will be found and consequently results will not be changed. |
||||
StatusThe targeted diffusion to avoid cold pools has been implemented and is ready for COSMO 5.1. |
||||
Technical Issuescode written according to standards: fulfilled 4-eyes assurance: by Oliver Fuhrer responsible person: Michael Baldauf |
||||
TestingThe 'targeted diffusion' cured the above described cold pools. Extensive numerical experiments at DWD for both COSMO-DE (2.8 km) and COSMO-EU (7 km) setups have been done. In the attached presentation, verification results of this action (together with the 'reformulation in the calculation of the divergence damping coefficient' and the 'adaptation of RK for stochastic physics') are shown: the results are neutral (as expected). |
||||
DocumentationA new section 8.4.3 is ready for the COSMO Sci. Doc part I. (it will be available with the next release of this document). User's guide: no adaptation necessary process documentation (history/changes log file): is updated |
||||
2014-03-31 | Remove hacks in the tracer module | WG2 | finish | Fuhrer, Roches (MeteoCH) |
DescriptionWhen implementing the new Tracer Module (Roches, Fuhrer (2013) COSMO Tech. Rep.) several decisions had to be made concerning specialized treatment of tracer species. In the first implementation of the Tracer Module, those had been treated with the Metadata functionality of the Tracer Module. However the question arises, if some (or the most) of these special treatments are still necessary. In particular, the following special treatments (‘hacks’) should be removed (from Roches, Fuhrer (2013) COSMO Tech. Rep.): Fixes concerning the Leapfrog dynamical core: CLP_10E-12: clip qi to zero, if qi < 10-12 Fixes concerning the Runge-Kutta dynamical core: ADD_CLP_ADV: add clipping for sedimenting moisture species (qr, qs, qg) at the end of the advection routine. BD_0GRAD_FORCED: for the species qi, qr, qs, qg either boundary values are read from a file or the boundary condition (BC) grad=0 is used. In the original code version one of these is done in any case despite the fact, that one can prescribe also other boundary conditions. DAMP_FORCED: For precipitating species, Rayleigh damping is done in any case, even when a grad=0 BC is prescribed. It is quite difficult to decide, if these measurements are still necessary (they have often been implemented during the development phase and sometimes accidentally remained in the code, even when they are not longer necessary). This is in particular the case if the original developer of the code is not longer available. Therefore, one needs real case test runs to decide, if anything strange happens or if the forecast quality suffers. |
||||
Statushacks which should be removed are identified. Code is ready for COSMO 5.1 |
||||
Technical Issues4-eyes Assurance: Oliver Fuhrer, Anne Roches |
||||
TestingTests during a longer period have been done at MeteoCH (Roches, Fuhrer). The documentation of these tests is available. Results are neutral. |
||||
DocumentationThe test results are documented in a report The whole tracer module is described in A. Roches, O. Fuhrer (2012) COSMO technical report no. 20. Sci. Doc. part I: no changes necessary User's guide: no adaptation necessary Process documentation (history/changes log file): not yet available |
||||
2014-04-16 | stochastic physics | KENDA | finish | Torrisi (CNMCA) |
DescriptionIntroduction of stochastic perturbation of the physics tendencies (SPPT). Its application will be optional, selectable by a namelist switch. For an outline of the basic method, see also: ../../consortium/generalMeetings/general2012/wg1-kenda/torrisi_stochastic_physics.pdf |
||||
Status Previous questions on how to perturb which tendencies (e.g. microphysics) have been sorted out. Namelist switches have been added to also perturb qc, qi resp. qs, qr tendencies Revised SPPT version (with i.e. removal of erroneous perturbations of Coriolis term) has been tested over 1 month; SPPT gives positive results In order to allow for applying SPPT in a (rapid) DA cycle without destroying the temporal correlations of the random fields: an option has been implemented recently to make all random fields in subsequent model runs reproducible. Originally, a minor technical modification to avoid using vertical coordinate parameters (which are not available any more with Grib-2 input) has been planned. However, recent communication with Uli Schaettler revealed that a modification is going to be implemented in propriety DWD Version V5.0.2. in mid / end of May which makes the vertical coordinate parameters available even for Grib-2 input. A more uniform approach will be implemented thereafter nevertheless to avoid the use of these parameters everywhere in the (physics and data assimilation parts of the) code. The bottom line is that the current code of SPPT is the final candidate to enter in V5.1 without any further changes. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)
|
||||
Testing (Single Test Cases, Experiments)Testing has shown small, but positive impact: - ../../consortium/generalMeetings/general2013/wg7-kenda.htm : talk by Andre Walser, for COSMO-E forecasting component, showing small but positive impact (larger positive impact than from using perturbed parameters, as done in COSMO-DE-EPS) - ../../consortium/generalMeetings/general2013/wg1.htm : talk by Daniel Leuenberger, in LETKF assimilation with COSMO-2, showing small but positive impact - ../../consortium/generalMeetings/general2012/wg1-kenda/torrisi_stochastic_physics.pdf : talk by Lucio Torrisi, in COSMO-ME forecasting component, showing small but positive impact Recent Testing has been performed at MeteoSwiss for COSMO-E and has been reported at the COSMO-User Seminar. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Inline version history is written. Changes log-file is availableThis includes a description of all the new namelist parameters. Update of the COSMO User's Guide: Shall be postponed, since Uli Schaettler is going to change the format of the namelist tables in this document. When the new format is being implemented, the new namelist parameters will be added, based on the changes log-file. This is planned to be done before V5.1 will be available. Scientific documentation: A short section has been written (Lucio Torrisi) The scheme has been outlined in a talk by Lucio Torrisi held at a 'Stochastic Physics' workshop at DWD on 26 Nov. 2013. |
||||
2014-04-16 | interface for radar obs operators | KENDA | finish | Blahak (DWD) |
DescriptionOnly technical changes: Add interface to radar observation operators; this in particular includes the calls for the cpu time measurement. Use of ifdef e.g. for the call to the radar obs operator and even the cpu time measurement. |
||||
StatusAll is ready, code has been handed to Uli Schaettler |
||||
Technical Issues
|
||||
TestingWorks technically. |
||||
DocumentationModification comments included in code; no external documentation required. |
||||
2014-06-30 | Adaptation of the RK dynamical core for stochastic physics | WG2 | finish | Torrisi (CNMCA), Baldauf (DWD) |
DescriptionIn the current development of the stochastic physics packages, it his desired to perturb only physical tendencies and not dynamic tendencies (i.e. no perturbations of advection, fast waves or Coriolis forces). Furthermore, one targets at a proper perturbation of terms, which are a the border between physics and dynamics (in particular latent heating). This development probably must be done in several steps. The first step is a simple shift of the call of the Coriolis term in src-runge_kutta.f90. This will be go into COSMO 5.1. Other steps concerning the latent heating will follow. |
||||
StatusProgramming work for COSMO 5.1 is done. |
||||
Technical Issuescode written according to standards: fulfilled 4-eye principle: Lucio Torrisi responsible person: Michael Baldauf |
||||
TestingFrom a theoretical viewpoint there is no objection in shifting the call of the Coriolis term from the current position behind the call of the vertical implicit diffusion (neither in changing stability nor accuracy properties). Several stand alone runs have shown that the position of the 'call coriolis' has a negligible influence on the results, as expected (Baldauf). Tests in two longer standing tests with COSMO-DE and COSMO-EU didn't show any significant influence to the scores (tested together with the 'targeted diffusion to avoid cold pools' and 'Reformulation in the calculation of the divergence damping coefficient'). |
||||
DocumentationChange of position in 'call coriolis': this very small change needs no extra documentation. (Possibly there will be a larger common documentation on the whole stochastic physics model changes) Sci. Doc part I: update of section 8.1.1 done (available with the next release) User's Guide: no change necessary process documentation (history/changes log file): is updated |
||||
2014-06-30 | Turbulence: Implementing TKE-advection for TURBDIFF | WG3a | finish | Blahak |
DescriptionTKE-advection can be activated by choosing lprog_tke=.TRUE.. This is being working for years in the framework of the Runge-Kutta Core (l2tls=.TRUE.), if one of the alternative turbulence schemes are selected, which had been introduced mainly for LES purposes (itype_turb=5...8). In the course of a stepwise generalization of the standard Raschendorfer scheme (TURBDIFF, active for itype_turb=3) TKE-advection should be possible for this scheme as well. Thus this task can be regarded as a part of "Implementing 3D-extionsions for TURBDIFF". The implementation should be along the lines of tracer advection and thus be applicable for semi-Lagrangian advection, flux-form density based advection (Bott et al.) or the traditional formulation with divergence correction. Since the standard turbulence scheme employs a prognostic equation of the turbulent velocity scale q (square root of the turbulent stress tensor trace) and not of TKE (being half of the squared q), q is going to be advected, where it is TKE for the LES schemes. Both of them are contained in the array 'tke', which is defined on half levels. In case of the flux form advection, we do it without a separate density advection performed at half levels for the time being. In order not to apply the optional exponential time filter for q (controlled by the parameter 'tkesmot') not to the advection tendencies, an additional array for the advection tendencies of q is used. We expect benefits mainly for higher horizontal resolution. Regarding to the forecast of eddy dissipation rate (EDR) for aviation some further benefit is possible, in particular if the additional scale interaction terms related to SSO wake production, separated horizontal shear eddies or convective currents are considered in the prognostic TKE equation. |
||||
StatusFirst implementation ready in a test code including also horizontal shear production and horizontal diffusion (applicable for the standard scheme). The test code is based on COSMO version 5.0 and has been tested to work formally correct. The modifications are going to get part of version 5.1 to be available for more rigorous testing for high resolution model runs within the whole WG3a community. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)The code modifications are along the coding standards and they have been inspected by M. Raschendorfer and J. Förstner. The implemented extensions are optional and allow for a more stringent selection of 3D-aspects for the standard turbulence scheme. Their operational benefit will be tested afterwards. |
||||
Testing (Single Test Cases, Experiments)The extensions have been tested by U. Blahak. Technical correctness was confirmed by a simple idealized test. Meteorological tests have been performed for a single real test case (COSMO-DE) and for an idealized LES-like setup with a much smaller grid spacing. Longer experiments are needed for investigating possible positive impacts already with the currently operational model resolution, but not undertaken so far. All tests are documented in tkeadv.pdf. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)A draft documentation of the implementation and preliminary results are described in a LaTeX document: tkeadv.pdf Longer experiments are needed for investigating possible positive impacts already with the currently operational model resolution, but not undertaken so far. |
||||
2014-06-30 | Implementation of a unified COSMO-ICON microphysics | WG3a | finish | Köhler / Fuhrer / Lapillonn |
DescriptionThe unified COSMO-ICON model microphysics includes different changes conducted by Felix Rieper and Günther Zaengl: Cloud ice sedimentation, an alteration in the sticking efficiency of cloud ice, the treatment of supercooled liquid water vapor, the reduction in the freeezing rate of rain and some bug fixes related to the values of physical parameters. These changes have been implemented and can be controlled by switches in the code. Default is false. The code will be adopted from ICON, where it is already written in the block data sturcture. |
||||
StatusThe code is ready in a test version and is being tested for real cases. |
||||
Technical IssuesCoding Standards, Technical Test Suite, 4-eyes Assurance The code has been implemented along the coding standards and all modifications have been checked by at least one additional persson. |
||||
TestingSingle Test Cases, Experiments The new common microphysics module has been tested within the ICON framework with some remaining problems in the ICON stratosphere, when running the supercooled liquid water version. A comparison using the operational setup at MeteoSwiss along a 15 days period showed only minor differences. |
||||
DocumentationPresentation of Results, Model Documentation, External Documentation A short description of the main modifications and the related NAMELIST-switches is in cluded in a document of C. Köhler. A 1st-verification of the new code has been performed by MeteoSwiss for a 3-week period starting at April 12 2014. In this parallel experiment all the new development (except some bug fixes) has been switched off. Thus mainly the technical adaptations (e.g. the code in block-data structure) has been inspected. As expected, the verification results are neutral. Now also verification results from DWD's experiments with COSMO-EU and COSMO-DE are available. |
||||
2014-12-31 | Turbulence: 3D-extensions for TURBDIFF within the concept of scale separation | WG3a | finish | Raschendorfer / Blahak |
DescriptionSo far, the standard Raschendorfer-scheme (TURBDIFF, active for itype_turb=3) is based on the closed 2-nd order equations for turbulent statistical moments, which uses a prognostic equation for TKE (and in future also optional prognostic equations for scalar variances - cf. "extension to the TKESV configuration") and a system of diagnostic linear equations for the remaining 2-nd order moments. This all together is solved in horizontal boundary approximation HBLA similar to the classical Mellor/Yamada scheme, what makes it to be a pure single column (SC) scheme. We are now aiming to extend this scheme towards a 3D-scheme that can be applied for the targeted convection-permitting model resolution without losing its applicability for coarser model resolution. For a first step we want to treat the prognostic TKE-equation as a 3D-equation, considering full 3D shear-production, advection and diffusion of TKE, where the remaining linear system is still solved in HBLA except the fact that the vertical shear forcing term is substituted by the 3D-term. The arising horizontal momentum flux densities are expressed by extending the flux gradient representation of the SC solution to the horizontal directions employing the resulting diffusion coefficients as isotropic properties. The same representation is also used for horizontal diffusion of prognostic variables (including TKE). In a second step we want to include the separated horizontal shear mode that has been introduced before as an optional scale-interaction production term in our TKE-equation. For that purpose, this production term is treated as an sink term for the separated horizontal shear mode being balanced by horizontal shear related to horizontal grid scales. A derived additional horizontal diffusion coefficient can be added to the isotropic turbulent one for calculation of the overall horizontal diffusion. We further want to modify the parameterization of the turbulent master length scale to be valid also for grid cells with a vertical dimension larger than the horizontal one. For controlling the options, we use the NAMELIST switch l3dturb that has already been introduced for the alternative turbulence schemes for LES purposes (itype_turb=5...8) and the selector itype_sher: l3dturb=.TRUE.: horizontal shear and horizontal diffusion active (also for itype_turb=3) itype_sher=1: using isotropic turbulent diffusion coefficients and related shear production only itype_sher=2: using the separated horizontal shear as an additional TKE source itype_sher=3: considering also the additional horizontal diffusion by the separated horizontal shear mode The activation of TKE-advection is controlled by the switch lprog_tke, similar as for the LES turbulence shemes (cf. "Implementing TKE-advection for TRUBDIFF"). Note that both, TKE-advection and horizontal diffusion of all prognostic variables (including TKE) is only provided, if the Runge-Kutta Core (l2tls=.TRUE.) is active. We expect benefits along frontal zones and along steep orography with increasing significance for increased horizontal model resolution. |
||||
StatusFirst implementation ready in a test code including also TKE-advection (applicable for the standard scheme). The test code is based on COSMO version 5.0 and has been tested to work formally correct. The modifications are going to get part of version 5.1 to be available for more rigorous testing for high resolution model runs within the whole WG3a community. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)The code modifications are along the coding standards and they has been inspected by U. Blahak. The implemented extensions are optional and allow for a more stringent selection of 3D-aspects for the standard turbulence scheme. Their possible operational benefit (even for the currently operational model resolution) will be tested afterwards. |
||||
Testing (Single Test Cases, Experiments):Technical correctness has been checked by U. Blahak (cf. tkeadv.pdf). Further meteorological testing is required with special attention to higher horizontal resolution. An additional uncertain scaling parameter for the description of the separated horizontal shear mode can perhaps be estimated by fitting kinetic energy spectra for high resolution runs. |
||||
DocumentationPresentation of Results, Model Documentation, External Documentation The whole extended scheme will be documented as soon as the common turbulence module for COSMO and ICON is ready. The 3D-extension will not yet be used operationally and has been implemented in order to facilitate test runs with this 3D-extionsion of the operational scheme. Explainations of the needed NAMELIST-switches can be taken form the "descrciption" section. |
||||
5.00 | refVersion | 2012-07-01 | 2012-12-21 | 2013-11-08 |
The contents of COSMO-Model 5.0 has been accepted and approved by the COSMO Steering Committee at its meeting during the COSMO General Meeting in Sibiu in September 2013. Since then there have been only some technical changes and updates. COSMO-Model 5.0 now is the Reference Version of the COSMO Consortium. COSMO-Model 5.0 will be distributed to all COSMO-partners during November 2013. |
||||
4.30 | test | 2013-10-11 | 2013-11-01 | 2013-11-08 |
This version contains the very last technical changes and fixes for the new release COSMO-Model 5.0.: Release Notes This version has been distributed to all COSMO-partners, the CLM Community and the COSMO-ART group for further testing on November 08th. It will be distributed to all Licence takers and other partners in the next weeks. |
||||
4.29 | test | 2013-08-12 | 2013-09-30 | 2013-10-02 |
This version contains a unification of the COSMO and the MESSy tracer structures. Many of the "ifdef MESSY" constructs could be avoided now. In addition, there are a few other technical changes and bug fixes (see Release Notes) This version has been distributed to all COSMO-partners, the CLM Community and the COSMO-ART group for further testing on September 26th |
||||
4.28 | test | 2013-03-01 | 2013-06-14 | N/A |
The major developments of this version are documented in extra entries below. Besides there have been the following adaptations (see the Release Notes for further documentation):
This new version has been distributed to all COSMO-Partners, the CLM Community and the COSMO-ART group for further testing on July 15th. |
||||
2012-12-13 | optional new QC check for surface pressure against lateral BC fields | WG1 | finish | Schraff (DWD) |
DescriptionNew option to add a quality control check for surface pressure obs against the fields that provide the lateral boundary fields. At DWD, this means checking against interpolated GME fields. New namelist parameter 'qcflbcp': scales QC theshold for check against COSMO fields to obtain threshold for new check. |
||||
Technical IssuesIt is coded and tested in V4_22, and ported to V4_27. It is part of the task 'modular observation operators for existing obs types' (see KENDA list), and is being submitted together with that task. |
||||
Technical IssuesCoding Standards : fulfilled. Technical Test Suite : - 4-eyes Assurance : together with task 'modular observation operators for existing obs types' (see KENDA list), partly Andreas Rhodin and Uli Schaettler |
||||
TestingSingle Test Cases : - target met: very positive results in 3 cases with series of erroneous pressure obs from a single buoy - target met: no negative impact in first Christmas storm 1999 (no further rejection of good data) Experiments: - target met: 1-month test period (April 2012) with neutral results. |
||||
DocumentationPresentation of Results : see talk at COSMO GM 2013 in Lugano, WG1 Parallel Session Additional results because of recent problems in operational COSMO-EU. Model Documentation : in-line documentation External Documentation : being written |
||||
2013-05-31 | Implementation of grib-api | WG6 | finish | Schaettler (DWD) |
DescriptionSince some years, the new GRIB standard (GRIB2), defined by WMO, is ready to be used. ECMWF developed an application programmers interface (grib_api), which can be used within programs to read, create and manipulate grib messages for GRIB2, but still also for GRIB1. At DWD it has been decided to use grib_api as the official GRIB library in the future. Also the COSMO partners agreed to use grib_api. Therefore, grib_api will be implemented in the COSMO-Model and the INT2LM in the next few months. At DWD the implementation is based on using the concept of the shortnames and the tables and definitions provided by grib_api. These tables and definitions have to be coordinated with the COSMO-partners. |
||||
StatusWork has started with delay in March 2013 to implement grib_api. It is implemented in addition to the DWD Grib-library, so the user can choose, whether Grib1 shall be written using the DWDLIB or whether the grib_api (for Grib1 or for Grib2) shall be used. 27.05.13: First implementation finished, Functionality tests started; General vertical coordinate still missing. 28.06.13: Added general vertical coordinate (new level type 150). But there are still issues with grib_api. Need Version grib_api 1.11.0 to process the new keys for the general vertical coordinate. 12.07.13: Several tests have been done and finished. There are still problems with writing the general vertical coordinate in ensemble mode; but new model version released now anyhow (for testing) |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: have been met Technical Test Suite: a DWD-pre-implementation of the test suite has been run.
4-eyes Assurance: code is cross-checked by DWD colleagues with grib_api experience |
||||
Testing (Single Test Cases, Experiments)27.05.13: Functionality tests started, but will go on during June 2013. 28.06.13: Functionality tests using grib_api are positive. 16.07.13: The new releases of INT2LM and the COSMO-Model are put to DWD's parallel suite. Experiments will be started soon using grib_api and GRIB2 |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)A COSMO-web page has been set up to document the implementation and GRIB2 usage. This page is not yet visible (still work in progress). It will be enhanced in the next few weeks. User Guide has been updated. |
||||
2013-11-08 | modular observation operators for existing observation types | KENDA | finish | Schraff (DWD) |
Description
|
||||
StatusIt is coded and tested in V4_22, and ported to V4_27. It is part of the task 'optional new QC check for surface pressure against lateral BC fields' (see WG1 list), and is being submitted together with that task. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards : fulfilled. Technical Test Suite : - 4-eyes Assurance : together with task 'modular observation operators for existing obs types' (see KENDA list), partly Andreas Rhodin and Uli Schaettler |
||||
Testing (Single Test Cases, Experiments)1-month test with V4_22 and only some of these technical modifications A new test is being started |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)IModel Documentation : improved in-line documentation External Documentation : no changes required for scientific documentation, User's Guide |
||||
2013-11-08 | for OSSE capability: read obs from feedback files | KENDA | finish | Schraff (DWD) |
DescriptionDevelop option to read observations from feedback ('fof') files instead of (or in addition to) 'cdfin' files, to use the simulated observation values from the feedback files as observations (from the nature run of an OSSE), and to perturb them randomly. |
||||
StatusIt is coded and tested in V4_22, and ported to V4_27. It is part of the task 'optional new QC check for surface pressure against lateral BC fields' (see WG1 list), and is being submitted together with that task |
||||
Technical IssuesCoding Standards : fulfilled. Technical Test Suite : - 4-eyes Assurance : together with task 'modular observation operators for existing obs types' (see KENDA list), partly Andreas Rhodin and Uli Schaettler |
||||
Testing1-month test with V4_22 and only some of these technical modifications Single test runs with all modifications based on V4_27 done, a longer test period is about to start |
||||
DocumentationModel Documentation : improved in-line documentation External Documentation : no changes required for scientific documentation, User's Guide |
||||
4.27 | test | 2012-12-01 | 2013-02-28 | 2013-03-21 |
Besides the consolidation work in the new fast-waves solver and the MESSy interfaces, version 4.27 now contains:
Also there are few technical modifications. More information on all the changes can be found in the Release Notes The new version has been distributed to all COSMO-Partners, the CLM Community and the COSMO-ART group for further testing on January 21st. |
||||
2013-04-30 | Consolidation work in the New fast waves solver | WG2 | finish | Baldauf (DWD) |
DescriptionAdd features in the New fast waves solver (see item 0.1 'Revised Runge-Kutta Core') which are contained in the old fast_waves_rk: 1.) Upper sponge Layer of Klemp et al (2008) 2.) Radiation Condition (with an additional switch for the selctive use of x- and/or y-direction) 3.) lw_freeslip functionality These are switchable functionalities (the NAMELIST switches exist already) |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2013-04-30 | Modular Earth Submodel System (MESSy) | Ext | finish | N/A |
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
4.26 | test | 2012-10-29 | 2012-12-01 | 2012-12-12 |
This version is mainly a consolidation version. Several problems have been detected in the developments of the former versions, which have been fixed now. But also new issues have been added: The new issues:
And several things have been modified / adapted:
More information can be found in the Release Notes. The updated User Guide is available on the COSMO Web Page. The new version has been distributed to all COSMO-Partners, the CLM Community and the COSMO-ART group for further testing on December 12th. |
||||
2012-12-13 | RTTOV obs operator + read / write satellite radiances | KENDA | finish | Messer (DWD) |
Description
|
||||
StatusThe modifications have been implemented in Version 4.26 (using ifdef RTTOV10). With setting -DRTTOV10 it is now possible to use the RTTOV10-library for producing the synthetic satellite images and to perform satellite observation processing. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled Technical Test Suite: Not yet available, but technical functionality has been tested with DWD COSMO-EU and COSMO-DE setup. 4-eyes Assurance: by Messer and Schaettler. |
||||
Testing (Single Test Cases, Experiments)The use of RTTOV10 is optional and could be tested now. Because of performance problems with the vectorization of the RTTOV10 library on the SX-9, it is not yet used operationally at DWD. See some pictures from runs with the different versions. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Offline documentation is available from DWD, FE12. The implementation is documented in the Release Notes of 4.26. The User Guide has been updated accordingly. |
||||
4.25 | test | 2012-07-01 | 2012-08-31 | 2012-09-28 |
The Version 4.25 again contains many different changes. Despite of the planning, there have again been some bug fixes, which do change the results of all model applications. Besides the items that are listed in the development table, the following issues have been implemented:
More information can be found in the Release Notes. This version has been distributed to all COSMO-Partners and the CLM Community at the end of September. In October / November there was a longer debugging and consolidation phase, because problems have been reported for the restart functionality and also for both asynchronous I/O modules. |
||||
2012-09-01 | Reconstruction of mpe_io.f90 | WG6 | finish | Prill (DWD) |
DescriptionThe module mpe_io.f90 is a stand-alone module of the NWP suite, containing an asynchronous parallel implementation for reading and writing GRIB1 data from / to disk. It is used in the COSMO-Model and also in INT2LM. At DWD it is also used in the GME. Lately, a bug has been detected for writing the ready-files (see below) in asynchronous mode. At the same time, some other optimization issues treated.
Current implementation work:
It turned out, that some interfaces of mpe_io had to be extended. Side Note: Within HP2C and the PP POMPA, work on a scaling, asynchronous I/O not only for GRIB, but also for NetCDF is ongoing. This work on mpe_io.f90 is not an extra parallel work, but just to fix the bug with ready files and to improve I/O of GRIB1 files immediately. The work is coordinated within PP POMPA. |
||||
StatusAll above mentioned items have been implemented into a new module mpe_io2.f90. In this way we can adapt the single programs (INT2LM, GME) one after the other to the changes in the asynchronous GRIB I/O. The modifications have been implemented in Version 4.25 and functionality has been tested successfully. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled Technical Test Suite: technical tests have been performed to prove functionality. 4-eyes Assurance: All above mentioned items have been implemented into mpe_io2.f90 and will be cross-checked by Uli Schaettler. |
||||
Testing (Single Test Cases, Experiments)The reconstructed module has been tested with the COSMO-Model in the experimentation system of DWD. A problem when using many compute processors has been reported by DWD and by MeteoSwiss and could be fixed in Version 4.26. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)A LaTeX documentation of mpe_io2.f90 is available. Together with a (still to be developed) similar documentation for the asynchronous NetCDF I/O this should be published as a COSMO Technical Report to serve as technical documentation of (asynchronous) I/O. |
||||
2012-12-01 | I/O Optimization with MPI_ALL2ALL | WG6 | finish | Fuhrer (MCH) |
DescriptionAn additional method has been implemented for gathering the subdomains of a field from all other processors using MPI_ALLTOALL. A Namelist switch itype_gather has been implemented to choose the special method. Which method behaves might depend on the computer used. |
||||
StatusThe modifications have been implemented in Version 4.25 and functionality has been tested successfully. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled. Technical Test Suite: The MeteoSwiss prototype of the technical test suite has been run successfully. 4-eyesAssurance: done by the Source Code Administrator (Uli Schaettler). |
||||
Testing (Single Test Cases, Experiments)Functionality has been tested successfully. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)User Guide has been updated accordingly. |
||||
2012-12-01 | Asynchronous NetCDF I/O Strategy | WG6 | finish | Carlos Osuna (C2SM, ETH Z� |
DescriptionAn asynchronous solution for output of Netcdf files will be implemented similar to the asynchronous Grib I/O. The focus will first be on Output. Depending on the available time and resources, also prefetching of input will be considered. |
||||
StatusA new module netcdf_io.f90 has been implemented in the source code of the COSMO-Model and existing modules (src_input.f90, src_output.f90, organize_data.f90) have been adapted to properly use that module. All changes have been implemented in Version 4.25. When I/O is done without asynchronous processors, nothing has been changed in the behaviour of the model. In that way, the new strategy can now be tested by the CLM community (or NetCDF users). |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: have been met Technical Test Suite: the prototype suite at MeteoSuisse has been used to test the developments 4-eyes Assurance: The code has been reviewed and checked by the CLM Source Code Administrators. |
||||
TestingFunctionality tests have been performed. Now the users of the model should test with their configurations. |
||||
DocumentationA presentation on the implementation and results has been given at a COSMO User Workshop in January 2012. |
||||
2012-12-01 | Microphysics: 2-Moment Scheme | WG3a | finish | Seifert / Blahak |
DescriptionTwo-moment microphysical parameterization for mixed-phase clouds to improve the explicit representation of clouds and precipitation. The scheme predicts the evolution of mass as well as number densities of the five hydrometeor types cloud droplets, raindrops, cloud ice, snow and graupel. The working implementation is based on a COSMO 4.20 test version. |
||||
StatusOnly the interfaces for the 2-moment microphysics have been implemented in COSMO Version 4.25. This facilitates the implementation of the code itself into a COSMO-Model version considerably. As this is only a technical change, only moderate tests have been necessary. The interfaces have been implemented with ifdef TWOMOM_SB |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)The code of the 2-moment microphysics is currently not according to the COSMO standards (more than 1 module per source code file; comments partly in German) and needs polishing. This is the reason why it has not been taken over to the official version at the moment. For the interfaces, all requirements are fulfilled: Coding Standards: are fulfilled Technical Test Suite: not yet available, but functionality of code without 2-moment microphysics has not been touched. 4-eyes Assurance (only for interfaces): done by SCA Uli Schaettler. |
||||
Testing (Single Test Cases, Experiments)The scheme has been developed now since several years and has been tested extensively. It has lately been presented to a workshop at ETH Zürich. Further optimization to run on the SX9 machines has been done. On this machine the scheme runs faster by a factor of 4 - 5 compared to the previous code, and it is slower by "only" a factor of 4 compared to the standard graupel scheme. It needs about half the time of the dynamical core. For non-vector machines, we would expect similar behaviour but this needs to be tested. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)There are already 5 published papers (peer review and conferences), more are in preparation. |
||||
2012-12-01 | Implementation of a tracer module | WG6 | finish | Roches / Fuhrer (MeteoSwiss) |
DescriptionCurrently, every module that uses passive tracers (e.g. COSMO-ART, microphysics) handles this individually. This detoriates the maintainability of the code (there is a lot of redundant code) and makes the introduction of a new tracer field tedious. Therefore, a general implementation of tracers is developed. There will be 2 new source files (data_tracers.f90 with constants and type definitions, src_tracers.f90 with the methods to implement the framework for handling tracers). The tracers will be stored in a contiguous memory block of the form trcr(i,j,k,ntracer,nt). The memory structure of current tracers in the model is not consistent. Therefore a choice has to be made about the data layout. The tracer module will keep as much information private as possible, in order to allow future changes to the internal data structure. Additional metadata may be added to the above strucutre as needed. The functionality of the tracer module will contain
|
||||
StatusThe new tracer module is available in COSMO-Model Version 4.25. All humidity variables are now treated as tracers, which led to a significant reduction in source code of the dynamics (instead of re-writing code for every qx-Variable, there now is a loop over the tracers). All modifications have been done in a way that the code gives bit-identical results. The new version has already been given to the COSMO-ART group and they started to adapt their codes to that new structure. Feedback and some necessary adaptations of the new tracer module have already been implemented in Version 4.26. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled Technical Test Suite: The new version has been tested with the MeteoSwiss Technical Test Suite 4-eyes Assurance: the code has been checked by Uli Blahak and Uli Schaettler. |
||||
Testing (Single Test Cases, Experiments)Only functionality tests are necessary, as the code gives bit-identical results. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)A COSMO Technical Report on the new Tracer Module, which serves as technical documentation, is in preparation. It is expected in December. |
||||
4.24 | test | 2012-05-01 | 2012-06-22 | 2012-06-22 |
Besides the changes listed in the Model development list, the following technical issues have been implemented. More information can be found in the Release Notes.
This version contains all (foreseen) modifications for the new unified release COSMO-Model 5.0, which do change the results of the simulations. All outstanding modifications for Version 5.0 are purely technical (or should at least be). It has been distributed to all COSMO-Partners and the CLM-Community for testing at the end of June. Feedback from Roshydromet, MeteoSwiss and the CLM Community has already been incorporated to Version 4.25. |
||||
2012-12-01 | adaptions for sub-hourly analysis update with LETKF | KENDA | finish | Reich (DWD) |
Description
|
||||
StatusThe modifications have been implemented in Version 4.24. The extension of file names has been implemented in a way that it is backward compatible:
Introducing the new file names can therefore be done smoothly by every partner. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled Technical Test Suite: not yet available, but technical functionality has been tested by running LETKF experiments and for the old file names by running the DWD COSMO-EU and COSMO-DE setup. 4-eyes assurance: done by Christoph Schraff and Uli Schaettler |
||||
Testing (Single Test Cases, Experiments)
|
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)New way of specifying ydate_ini has been documented in the COSMO User Guide The implementation of the modifications has been documented in the Release Notes. |
||||
2013-01-31 | New Fast Waves Solver for RK dynamical core | WG2 | finish | Baldauf(DWD) |
DescriptionMost of the operational setups of the COSMO model now use the so-called Runge-Kutta time integration scheme (Wicker, Skamarock, 2002). As in the original proposal of Wicker, Skamarock (2002, 1998) the fast waves are treated horizontally with a backward-forward scheme and vertically implicit to allow larger values for the small time step. An additional filter process must stabilize this whole time-splitting procedure;usually a divergence damping is used (e.g. Skamarock, Klemp, 1992) in the fast waves solver. 2. Use of the divergence operator in strong conservation form 3. Isotropic treatment of the artificial divergence damping References: |
||||
StatusThe modifications have been implemented in Version 4.24. The new solver now runs in the parallel suite of DWD for COSMO-EU and COSMO-DE. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled (checked by SCA Ulrich Schaettler) Technical Test Suite: not yet available, but functionality has been checked: fast waves part of the code runs 30% slower on NEC (increase in total runtime is 5%), no numbers yet for other platforms; numbers to be confirmed when code is finalised. Also on scalar x86-platforms the code runs slower. Only on Itanium platform (SGI Altix) it is nearly twice as fast as the old fast waves solver (reported from tests by RosHydromet). 4-eyes Assurance: done by Andreas Will Code Responsibility: M. Baldauf |
||||
Testing (Single Test Cases, Experiments)All Standard idealised test cases look fine (linear gravity wave in channel (Skamarock, Klemp (1992)), several Mountain flows (Bonaventura (2000)), cold bubble (Straka et al (1993)) test, moist warm bubble (Weisman, Klemp (1982)) test, ...) Real-case stand-alone test runs for several COSMO-EU (7 km), COSMO-DE (2.8km), and high-resolution COSMO-DE (2.2 km, 2km) are performed. Synoptic and upper air verification for winter and summer periods are available for COSMO-DE (2.8 km) runs. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Baldauf (2013): COSMO Technical Report No. 21 is available Results for several idealised test cases were presented at the COSMO user seminar 2011, Langen and the SRNWP-workshop, 2011, Bad Orb. Results for realistic simulations and time periods were presented at the COSMO User Seminar 2012 The User Guide (documentation of new switch itype_fast_waves) has been updated. |
||||
4.23 | test | 2012-03-17 | 2012-04-30 | 2012-04-30 |
In this version many developments have been incorporated that are planned for Version 5.0. Besides the items from the model development list, the following technical issues have been implemented:
From the CLM Community, the following contributions have been added:
More information can be found in the Release Notes. Version 4.23 has been implemented and distributed to the CLM-Community beginning of May. There was some technical feedback, which has been incorporated in Version 4.24 |
||||
2012-09-01 | Additional values in (long) meteograph output | WG6 | finish | Reinhardt (DWD) |
DescriptionAdditional values, which are calculated in the radiation, should be written to the (long) meteograph output to enhance the diagnostic output of the model. This will change the format of the long meteograph output, therefore all users have to be informed before. The new values are:
|
||||
StatusThe modifications have been implemented in Version 4.23 and could be tested by everybody. No complaints up to now. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: fulfilled. This is a pure technical change. It is only necessary to check the functionality. Technical Test Suite: not yet available; Functionality tested with DWD COSMO-EU and COSMO-DE setup. 4-eyes assurance: done by the Source Code Administrator (Uli Schaettler). |
||||
Testing (Single Test Cases, Experiments)Only functionality tests are needed. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)There was no need to update the User Guide, because the additional soil and surface variables are not explained in detail. |
||||
2012-09-01 | Smoothing of fields on p-/z-levels | WG6 | finish | Uli Blahak |
DescriptionAt the moment it is only possible to smooth all fields in the output list for p- or z-levels, or none of them. Because some of the output fields should not be smoothed (e.g. W), the smoothing should only be done for the fields, that need it (e.g. PMSL) Up to now, the following means were available to smooth special fields:
|
||||
StatusThe modifications have been implemented in Version 4.23 and could be tested by everybody. No complaints up to now. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: fulfilled. This is a pure technical change, but needs additional Namelist switches. Technical Test Suite: not yet available; Functionality tested with DWD COSMO-EU and COSMO-DE setup 4-eyes assurance: done by the Source Code Administrator (Uli Schaettler) |
||||
Testing (Single Test Cases, Experiments)Only functionality tests needed. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)The User Guide has been updated accordingly. |
||||
2012-09-01 | Bug fix in advection operators (Bott-scheme) and cloud diag | WG2 | finish | Fuhrer (MeteoCH) |
DescriptionThe use of the same field for more than one formal parameters in subroutine arguments can cause compilation errors with some compilers. (not clear if this is a compiler bug or a program bug) This issue must be combined with item "New tracer advection ..." |
||||
StatusThe modified interfaces have been implemented into the official version (4.23). This version has been provided for testing to the COSMO testing partners and to the CLM Community. No problems have been reported. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards are fulfilled. 4-eyes assurance has been done by Michael Baldauf and also by Uli Schaettler. Technical Test Suite: The MeteoSwiss Version of the Technical Test Suite has been applied successfully. |
||||
Testing (Single Test Cases, Experiments)The modified code gives bit-identical results (only the interfaces of certain subroutines have been changed). This has been tested using the operational COSMO-2 and COSMO-7 model runs. For certain compilers (NAG, Lahey Fujitsu) differences may appear when strongly optimizing the code, the new version should be correct. For the Cray XT4 the model performance was not significantly influenced by the bugfix and was within the measurement uncertainty. Code responsible: Michael Baldauf https://cosmo.cscs.ch/cosmo/branch/mch/olifu/bugs/bug01_arguments/results/. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Presentation of results: N/A (bug) Provide extensions / modifications for the Documentation System: N/A (bug) Future Code responsible: Michael Baldauf |
||||
2012-09-01 | calculation of the moisture divergence | WG2 | finish | Fuhrer (MeteoCH) |
DescriptionThe diagnosis of dqvdt was intermixed with the computation of vertical diffusion of qv and qc in the subroutine complete_tendencies_qvqcqi_tke (and the subroutine slow_tendencies for the Leapfrog core). This should be separated into two different subroutines (for Runge-Kutta) or different code blocks (for Leapfrog) in order to make the code more modular and extensible in the future. This requires changes to the source files slow_tendencies.f90, src_runge_kutta.f90 and src_slow_tendencies_rk.f90. The changes are of pure technical nature and do not influence the model results. |
||||
StatusThe changes have been implemented in a private version at MeteoSwiss and tested with the operational MeteoSwiss setup. The modified code has been submitted to the Source Code Administrator and has been implemented in Version 4.23 |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled Technical Test Suite: not yet available; operational MeteoSwiss setup has been tested 4-eyes Assurance: code has been reviewed by Michael Baldauf and Uli Schaettler |
||||
Testing (Single Test Cases, Experiments)3) Testing of the changes according to the Quality Control System of the COSMO partner or of the external partner that implements them. Testing: The modified code has been checked for the MeteoSwiss operational setup (COSMO-2 and COSMO-7) and proven to give bit-reproducible results. If an optimization is made in the LF core (removal of "+qvtens-qvtens") rounding errors can lead to differences in the range of the numerical precision. These changes only influence the computation of dqvdt and thus are only visible if the convection scheme is active (COSMO-7) or in the DQVDT diagnostic output. The impact on model performance is negligible (timings of LF and RK simulations attached) and cannot be discerned from the jitter on the machine. In the results, the modified version of the code consistenly ran slightly faster than the original version, which is not understood and probably not statistically significant. The 4-eyes principle has been applied (Oli and Anne) but a review of the code responsible and WG2 chair (Michael) would nevertheless be optimal. In summary, no significant changes are expected, bit-reproducibility can be achieved (always for RK, but only with non-optimal code version in the LF case). |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Presentation of results: not applicable Changes to Model Documentation: not applicable Changes have been documented in the Release Notes Future Code responsible: Michael Baldauf |
||||
2012-09-01 | Removal of (dead) code | WG6 | finish | Schaettler, et al |
DescriptionThere are several parts in the model. which are not used and not maintained any more. These parts should be removed from the code. The following items could be removed soon:
|
||||
StatusThe items listed above have been removed in Version 4.23. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Not applicable; Just make sure that nobody is using that any more |
||||
Testing (Single Test Cases, Experiments)Working configurations must not be touched. Up to now, nobody complaint. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)The User Guide has been updated accordingly. |
||||
2012-09-01 | Changes in parameterization for ice cloud cover in radiation | WG3a | finish | Seifert |
DescriptionThe code in src_radiation to reduce the cloud cover of ice clouds in the upper troposphere has been modified:
|
||||
StatusCode has been tested in 2011 by Axel Seifert in private experiments with neutral / good results. The modifications have been implemented in Version 4.23. |
||||
Technical IssuesCoding Standards: are fulfilled Technical Test Suite: not yet available 4-eyes Assurance: has been reviewed by Uli Blahak |
||||
TestingLonger experiments have been conducted at DWD with neutral verification results. |
||||
DocumentationThe problems and the code changes have been discussed in a (draft) paper. The implementation has been documented in the Release Notes. No changes of documentation necessary. |
||||
2012-09-01 | CLM Contributions to Unified Version 5.0 | Ext | finish | Hans-Jürgen Panitz |
DescriptionFor the next unified version of the COSMO-Model there are some contributions from the CLM Community. All these contributions are either of a pure technical nature or are optional (for the NWP mode) and can be switched on/off by Namelist parameters. Therefore the results of NWP applications are not affected by these contributions. This is a short summary of the contributions:
The changes are documented and described in more detail in the Release Notes of Version 4.23 |
||||
StatusThe change has been implemented into Version 4.23 and has been given to the CLM Community for further testing. Additional feedback has been implemented in Versions 4.24 / 4.25. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: fulfilled Technical Test Suite: not yet available, but changes could be tested by all COSMO Partners 4-eyes Assurance: done by Uli Schaettler |
||||
Testing (Single Test Cases, Experiments)All changes have been tested by the CLM Community and in addition by the COSMO partners. No complaints were raised. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)The COSMO User Guide has been updated accordingly. The changes have been described in detail in the Release Notes of Version 4.23 |
||||
2012-09-01 | Time measurement for long simulations: Problems with memory | WG6 | finish | Uli Schaettler |
DescriptionProblems have been reported with the memory usage, when running long simulations. The model crashes at the end of the forecast, when the timing informations from the processors are gathered for output. This is due to the fact, that all hours are gathered separately, so the amount of memory can be too big then. Solution: For itype_timing = 1/3: A loop over all hours has been implemented to gather values only from one hour and process them in PE 0 For itype_timing = 2/4: The sum over all hours is computed in the PEs and only the sums are gathered from the different PEs. |
||||
StatusThe change has been implemented into Version 4.23 and could be tested by all COSMO Partners and the CLM Community. No complaints were raised. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)This is a pure technical change. The functionality has been tested with some standard runs. |
||||
Testing (Single Test Cases, Experiments)Only check the functionality. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)No need to update existing documentations. The change has been documented in the Release Notes. |
||||
2012-09-01 | Use SST forecast from an ocean model as boundary condition | WG6 | finish | Lucio Torrisi |
DescriptionThe usual procedure for numerical weather prediction for few days is Usually in numerical weather prediction for few days, the Sea Surface Temperature (SST) is held constant over the ocean. Only in climate runs, it is updated as are all the other slowly varying external parameters (as leaf area index, root depth, etc.). But sometimes it is useful to have an updated SST without updating the other external parameters. Therefore a new logical switch is implemented in the group /IOCTL/, to update only the SST over the ocean during a NWP simulation. |
||||
Status
The change has been implemented into Version 4.23 and could be tested by all COSMO Partners and the CLM Community. No complaints were raised. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: fulfilled Technical Test Suite: not yet available, but some standard tests have been conducted without problems 4-eyes Assurance: done by Uli Schaettler |
||||
Testing (Single Test Cases, Experiments)No changes of results. Only some functionality tests are necessary. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)User Guide has been updated. |
||||
2012-09-01 | Physics-Dynamics Coupling | WG6 | finish | Oliver Fuhrer (MCH) |
DescriptionMove computation of total physical tendencies to the physics. Up to now this has been done seperatly in the different dynamical cores. But the proper place is in the module organize_physics.f90, after all tendencies have been computed. This is a pure technical change without changing of results. |
||||
StatusThe change has been implemented into Version 4.23 and could be tested by all COSMO Partners and the CLM Community. No complaints were raised. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: fulfilled Technical Test Suite: The prototype of the Test Suite by MeteoSwiss has been run without problems. 4-eyes Assurance: done by Uli Schaettler |
||||
Testing (Single Test Cases, Experiments)As there are no changes to the results, only the technical Test Suite was necessary. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Change is documented in the Release Notes. No other change of documentation necessary. |
||||
2012-09-01 | New tracer advection schemes: MPDATA and 'Bott2_Strang' | WG2 | finish | deMorsier (MeteoCH), Baldauf ( |
DescriptionMPDATA is currently developed in the PP 'Conservative dynamical core'. It has the advantage of combinig tracer mass conservation and to be a fully 3-dimensional advection scheme; however the order of accuracy is less than that of the existing schemes 'Bott' (direction splitted scheme) and Semi-Lagrange (non-conservative). The Strang-splitting version of Bott is more stable than the 'standard' Bott version, but for the price of about 60% more computing time. This could be avoided by applying Strang-splitting only in the lowest levels. To this purpose some reprogramming of the calling routines of the Bott-scheme seems to be advisable. Currently almost the same (large!) piece of code is contained 4 times in src_advection_rk.f90. This almost identical code can be reduced to one subroutine. |
||||
StatusThis technical rewriting and a preliminary version of the more efficient Strang-splitting version is implemented now in COSMO 4.23. MPDATA must be extended into a fully 3D-version before it can go into the official version. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled Technical Test Suite: not yet available; technical tests have been performed to prove functionality. 4-eyes: for MPDATA by M. Baldauf (DWD) for revised Bott2_Strang: by G. deMorsier (MeteoCH) code responsibility: M. Baldauf |
||||
Testing (Single Test Cases, Experiments)At DWD an experiment is ongoing, where the efficient Strang-splitting version is tested. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Current state of MPDATA: see presentation at COSMO GM 2011 by G. deMorsier |
||||
2012-09-01 | Reorganisation of metrics stuff | WG2 | finish | Baldauf (DWD) |
Descriptionconcentrate metric terms for the spherical, terrain following grid in an own module. This will facilitate the coupling of other dynamical cores (EULAG, ...). |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)Coding Standards: are fulfilled. The developers of FieldExtra have been consulted to discuss the interfaces. Technical Test Suite: not yet available; technical tests have been performed to guarantee functionality. 4-eyes: U. Schaettler (DWD) code responsible: M. Baldauf |
||||
Testing (Single Test Cases, Experiments)This is only a collection of subroutines from a diversity of files into one common file. Therefore no change of results |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
4.22 | release | N/A | N/A | 2012-01-31 |
2012-01-31 | Radiation: Zenith Angle Corrections | WG3a | finish | Reinhardt (DWD/AGeoBw) |
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2012-03-31 | new modular nudging code (with many new options) | KENDA | finish | Schraff (DWD) |
Description
|
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)
task completed, except (as at end April):: completion of 4-eyes assurance (Daniel Leuenberger, MeteoSwiss: check interfaces first, and based on the results, decide on what other parts of the code should rigorously be check by means of the 4-eyes principle) |
||||
Testing (Single Test Cases, Experiments)
task completed. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)
task completed |
||||
4.21 | release | N/A | N/A | 2011-12-06 |
2011-12-06 | Horizontal nonlinear Smagorinsky diffusion | WG2 | finish | Baldauf (DWD) |
DescriptionSometimes the COSMO model aborts due to horizontal shear instabilities. Some of these crashes can be avoided by the 4th order 'artificial horizontal diffusion' with a prescribed constant diffusion coefficient. But in rare events this diffusion is not strong enough and a more physically based diffusion mechanism is necessary. The nonlinear Smagorinsky diffusion (Smagorinsky (1963) MWR) determines the diffusion coefficient by the horizontal shear (and tension) strain and therefore acts in particular to reduce too strong horizontal shear. It is switched on by the DYNCTL namelist parameter 'l_diff_Smag=.TRUE.' (otherwise it has no impact to the results). An internal parameter, the Smagorinsky constant, is currently set to c_smag=0.03. This value is chosen to prevent shear instabilities in COSMO-DE applications, but otherwise to influence as less as possible the verification scores. |
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)4-eyes: Uli Schaettler Code responsibility: M. Baldauf |
||||
Testing (Single Test Cases, Experiments)Experiments with COSMO-DE during 01.-28. Feb. 2011 and 01.-31. Aug. 2011 were performed. Verification results are mainly neutral as expected. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Powerpoint-presentation with verification results held at 08. Nov. 2011 at 'Routine meeting' at DWD. New section 8.4.2 in COSMO Sci. Doc. part I Article by M. Baldauf and G. Zaengl in COSMO-Newsletter no. 12 (2012) |
||||
4.20 | release | N/A | N/A | 2011-08-31 |
Release Notes | ||||
2011-08-31 | Digital Filter Initialization | WG6 | finish | Schaettler |
DescriptionIn tests for historical storm events, driven by ERA Reanalyses, it turned out that the digital filter initialization used together with the Runge-Kutta scheme did not deliver reproducible results. But for the Runge-Kutta scheme an adiabatic backward integration is not possible at all. So the option ndfi=2 (backward-forward initialization) has been disabled, if the Runge-Kutta scheme (l2tls=.TRUE.) is used (in Version 4.20). To run the forward initialization (ndfi=1) with the Runge-Kutta scheme, some technical adaptations were necessary. These have been implemented in Version 4.18. Status of work: This issue has been fully implemented in Version 4.20:
|
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)The changes are still evaluated and tested by colleagues from the Bundeswehr. |
||||
Testing (Single Test Cases, Experiments)See Technical Issues. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Still missing: update User Guide to explain that no backward filtering can be done with Runge-Kutta dynamical core. |
||||
2011-08-31 | OASIS3 coupling | WG6 | finish | N/A |
DescriptionAn interface to the OASIS3 coupler is implemented to couple the Community Land Model (CLM) and an ocean model. Status of work: Necessary changes in the cosmo-model:
|
||||
Technical Issues (Coding Standards, Technical Test Suite, 4-eyes Assurance)The changes have been implemented by colleagues from ETH Zürich and have been cross-checked and implemented by the Source Code Administrator. |
||||
Testing (Single Test Cases, Experiments)Few single test cases have been performed to ensure that the results have not been changed. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
4.19 | release | N/A | N/A | 2011-08-01 |
Release Notes | ||||
2011-07-29 | Conditional compilation for external libraries | WG6 | finish | Schaettler (DWD) |
DescriptionThe cosmo-model and the int2lm use several external libraries: DWD Grib1 library, NetCDF library, a miscellaneous library, RTTOV libraries. Up to now, "dummy"-routines have been provided to satisfy external references, if these libraries are not available. The following conditional compilations have been provided, to get rid of these dummy-routines:
Status of Work:
Still missing: -DNOMPI (for creating a sequential binary without MPI) |
||||
Testing (Single Test Cases, Experiments)Few single test cases have been performed to ensure that the results have not been changed. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)The User Guide has been updated accordingly. Proper description of -DGRIBDWD still outstanding. |
||||
4.18 | release | N/A | N/A | 2011-05-25 |
Release Notes | ||||
4.17 | release | N/A | N/A | 2011-02-24 |
Release Notes | ||||
4.16 | release | N/A | N/A | 2010-12-06 |
Release Notes | ||||
4.15 | release | N/A | N/A | 2010-11-19 |
4.14 | release | N/A | N/A | 2010-06-14 |
Release Notes | ||||
4.13 | release | N/A | N/A | 2010-05-11 |
4.12 | release | N/A | N/A | 2010-05-11 |
Release Notes | ||||
4.11 | release | N/A | N/A | 2009-11-30 |
Release Notes | ||||
4.10 | release | N/A | N/A | 2009-09-11 |
4.09 | release | N/A | N/A | 2009-07-16 |
4.08 | release | N/A | N/A | 2009-02-16 |
4.07 | release | N/A | N/A | 2008-12-12 |
4.06 | release | N/A | N/A | 2008-09-24 |
Release Notes | ||||
4.05 | release | N/A | N/A | 2008-09-10 |
4.04 | release | N/A | N/A | 2008-07-16 |
Release Notes | ||||
4.03 | release | N/A | N/A | 2008-02-25 |
Release Notes | ||||
4.02 | release | N/A | N/A | 2007-12-05 |
Release Notes | ||||
4.01 | release | N/A | N/A | 2007-12-04 |
4.00 | release | N/A | N/A | 2007-07-26 |
Target version |
Version status | Development start |
Planned delivery |
Delivery date |
---|---|---|---|---|
2.05 | release | 2017-02-24 | N/A | 2017-03-31 |
There were some technical changes and bug fixes without changing the results, apart from preventing crashes. This release also contains the major changes that have been introduced in the earlier test versions 2.02a, 2.02a1, 2.03. Please consult the release notes of these versions: Version 2.04 contains no changes compared to 2.03 and has been just re-named to be consistent with DWD internal version control, without changing the content. There are some further slides documenting the changes and their effects of the vertical profile adaption (nonhydrostatic input models), the pressure computation (nonhydrostatic input models) and the soil temperature height correction bugfix (ICON-input): |
||||
2017-06-01 | Kinne (2015) climatology (AOT12,SSA12,ASY12 --> TAU_KINNE,SSA_KINNE,G_KINNE) | WG3a | finish | Marina Shatunova (RHM) |
DescriptionImplement new fields for the Kinne (2015) aerosol climatology, which are read from the external parameters and written to the analysis. |
||||
StatusOngoing work. Various new external parameters have been obtained from Stefan Kinne by Natalia Chubarova (Moscow State Univ.), have been prepared to be used in COSMO by Marina Shatunova (RHMS) and have been submitted to EXTPAR as netcdf files. Once these are available from EXTPAR, final implementation in int2lm is possible. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2.04 | test | 2017-02-14 | 2017-02-14 | 2017-02-14 |
This version is just a re-naming of version numbers to be consistent with DWD internal version control, without changing the content. |
||||
2017-02-14 | This version is just a re-naming of version numbers to be consistent with DWD internal version control, without changing the content | WG6 | finish | Ulrich Blahak (DWD) |
2.03 | test | 2017-01-31 | N/A | 2017-01-31 |
The version 2.03 contains the changes of the two DWD-internal versions 2.02a and 2.02a1. Thus, the release notes of 2.03 consist of these three parts: There are some further slides documenting the changes and their effects of the vertical profile adaption (nonhydrostatic input models), the pressure computation (nonhydrostatic input models) and the soil temperature height correction bugfix (ICON-input): |
||||
2016-08-22 | Linear vertical interpolation as an option to splines | WG6 | finish | Ulrich Blahak (DWD) |
DescriptionLinear vertical interpolation for some variables can now be chosen more easily in the code (src_vert_inter_lm.f90) by new internal switch in subroutine call: vert_interp (qr_lm, 'qr', ... --> vert_interp (qr_lm, 'qr', 'linear‘‚ … (or ‚ctspline‘) but no actual change yet (I‘m thinking about changing it for moisture variables and tracers in the future). No change of the results yet. |
||||
Technical IssuesInternally implmentet, but not activated so far. Would require an additional namelist parameter or a hardwired choice in src_gribtabs.f90. |
||||
TestingTechnically tested for a single case, where it worked correctly. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-08-22 | For GFS, possibility to use RELHUM instead of QV | WG6 | finish | Ulrich Blahak (DWD) |
DescriptionIn case of GFS input data possibility for input of RELHUM data. Previously only QV was read, but since 01/2015 GFS data has switched to RELHUM. Huge change of results, because before, qv_lm was set to the value for 1 % relative humidity, if no QV present in input data! |
||||
StatusFinished |
||||
Technical IssuesNone |
||||
TestingTested by Helmut Frank (DWD) and colleagues from RosHydromet. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-08-22 | General 4D field to hold new 3D fields for easier integration into the code. | WG6 | finish | Ulrich Blahak (DWD) |
DescriptionGeneral 4D field to hold new 3D fields for easier integration into the code. Works a bit like the tracer structure in COSMO, but much more rudimentary. It is meant only for the code developers to easily test and integrate new fields. Purely technical, does not change results. |
||||
StatusFinished |
||||
Technical IssuesNone |
||||
TestingIt was successfully used to implement the 11 CAMS aerosol mixing ratios from ECMWF in the framework of PP T2RC2. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-08-22 | MACC (CAMS) aerosol data from ECMWF | WG6 | work | Harel Muskatel (IMS) |
DescriptionImplemented interpolation of MACC aerosol data to COSMO grid with the help of the new 4D field. 11 new fields. |
||||
StatusFinished |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-08-22 | TERRA_URB: new EXTPAR fields AHF and FR_PAVED/ISA | WG6 | finish | Hendrik Wouters (KU Leuven) |
DescriptionTERRA_URB needs the two new external parameter fields AHF and FR_PAVED/ISA from the COSMO external parameters. These have been implemented in int2lm. New namelist switch lterra_urb. New grib-codings have been implemented in the DWD grib_api in the local use section. They have also been submitted to WMO, but this will take time. |
||||
StatusWorks technically , but grib codings are still local and preliminary. Int2lm uses shortnames 'aermrXX' (where XX stands for the number of the species, 01 to 11) for grib2 and local DWD-table 206, ee=101...111 for grib1. In the COSMO, these fields are only usable in the PP T2RC2 test version of the code. |
||||
Technical Issuesgrib coding of the CAMS aerosol fields 'AERMRXX' |
||||
TestingWorks correctly using our local preliminary grib coding. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-08-22 | Field SSO_STDH needed for the new blocked COSMO turbulence code | WG6 | finish | Ulrich Blahak (DWD) |
DescriptionFor the new blocked turbulence code, the field SSO_STDH is needed. This field is now required in any case and must be present in the external parameter files. Therefore, new external parameter files should be generated (from ASTER orography). |
||||
StatusFinished |
||||
Technical IssuesNone |
||||
TestingWorks correctly |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-08-22 | ICON-nest, ICON-lam or ICON-sub input: domain checker | WG6 | finish | Ulrich Blahak (DWD) |
DescriptionFor ICON input which is not a global dataset (ICON-LAM, ICON-NEST or ICON-SUB) the int2lm simply crashes with a segfault if some COSMO domain points are outside the ICON domain. Two different methods for domain checks have now be implemented: 1) A graphical ASCII-ART map showing the ICON and COSMO domains and their overlap. This is very handy to detect domain problems. 2) In case of ICON-SUB, a checker for the rotated lat/lon domain specifications of ICON-SUB and COSMO |
||||
StatusFinished |
||||
Technical IssuesNone. 4-eyes pending. |
||||
TestingWorks as expected |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2016-12-10 | Option to blend the soil temperature profile from the "normal" adaption method to new orography to a climatological-only height correction for deeper levels | WG6 | finish | Ulrich Blahak (DWD) |
DescriptionNew option for climatological height correction of deep soil temperatures for the multi-layer soil model New switch "lmultlay_deepsoil_clim_hcorr" (logical, default: .TRUE.) to specify if blending to a climatological height correction of soil temperatures for deeper soil layers is desired. The normal method is to preserve the temperature differences to the lowest atmospheric temperature across the interpolation, and this method is then only applied in the upper soil levels with a decreasing weight with increasing depth. This should reduce the danger that, e.g., local extremes in near-surface air temperatures, which may have been artificially created by the PBL-profile adaption from coarse to fine orography, create also artificial extremes in deep-soil temperatures where they have a long memory. The climatological height correction is based on a constant temperature gradient of 0.007 K/m with respect to orography difference coarse minus fine levels. In this respect it is similar to the already existing option "lt_cl_corr" for the climatological deep-soil temperature t_cl in the old 2-layer soil model. |
||||
StatusImplemented |
||||
Technical IssuesCoding Standards fulfilled, 4-eyes assurance pending |
||||
TestingNot really possible at DWD. Should be investigated at a center where COSMO is started from interpolated analyses without DA on a daily basis. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2017-01-31 | Hydrostatic pressure computation instead of PP interpolation for non-hydrostatic input models | WG6 | finish | Ulrich Blahak |
DescriptionFor non-hydrostatic input models (ICON, COSMO, UM, CM), up to now the pressure deviation was vertically interpolated rather computed from the basic hydrostatic equation. Especially over mountaineous terrain this lead to very noisy pressure fields in the horizontal, which generated a lot of sound- and gravity waves after COSMO model start. To mitigate this, now the pressure is hydrostatically integrated based on the interpolated T and QV profiles in each grid column, so that this noise is considerably reduced. |
||||
StatusFor the hydrostatic pressure computation, existing subroutines from COSMO src_artifdata.f90 have been implemented and used. No technical problems and issues were encountered. Work and implementation are finished. |
||||
Technical IssuesCoding standards are fulfilled. Technical testing has been done at DWD and other COSMO centers and found no problems. |
||||
TestingTests at DWD showed that, when started from an interpolated analysis (ICONEU in this case) the initial noise in COSMO simulations (DPSDT, WA500) decreased by 15 % on average in the entire COSMO-DE domain. Most of this reduction comes from the Alpine region, which in turn means that in this region the noise reduction is even larger. Small positive effects are also found through the boundary data alone, but COSMO-DE does not significantly "slice" through mountain ranges at its domain boundaries. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Documentation:
|
||||
2017-01-31 | New methods for profile adaption from coarse to fine orography in case of non-hydrostatic input models | WG6 | finish | Ulrich Blahak |
DescriptionFor non-hydrostatic input models, the methods for adapting the profiles from the coarse input orography to the fine orography are ok in general, but still there is room for improvement. First, the PP profile is no longer needed (because of the new hydrostatic pressure computation instead of PP interpolation), second new methods to stretch/compress the other profiles into deeper valleys or over higher mountains have been developed. These produce less initial noise when COSMO model is started from an interpolated coarse model analysis over mountaineous regions. |
||||
StatusFinished |
||||
Technical IssuesNone |
||||
TestingExperiments at DWD showed that the initial noise in COSMO-DE is reduced over mountaineous regions. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Documentation:
|
||||
2017-01-31 | Remove clipping of T_S and T_SO to -0.8 deg C below (partially) snow covered grid points | WG6 | finish | Ulrich Blahak |
DescriptionPreviously the surface temperature T_S has been clipped to a maximum of -0.8 deg C for (partially) snow covered soil (W_SNOW > 0 and T_SNOW <= 2 deg C). However, this clipping propagates also down into the soil via conservation of vertical differences when computing the grid-mean soil temperatures and may cause artificial cold patches of soil in situations with very little snow cover. Clearly this is not appropriate. Also, if there is an insulating effect of the snow on soil temperature, it should already be contained in the input soil temperatures. Therefore, this artificial clipping has been removed. |
||||
StatusFinished |
||||
Technical IssuesNone |
||||
TestingDoes what it should do. Removes unrealistic cold soil temperature patches in partial snow covered areas. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Documentation:
|
||||
2100-12-01 | For ICON input: Interpolation of convective precip rates for moisture enhancement in COSMO | WG6 | work | Ulrich Blahak (DWD) |
DescriptionICON has a very active convection scheme, which means that initial and boundary conditions derived from ICON have generally less moisture than for example COSMO-EU. This leads to a precipitation gap near the inflow boundary and a considerable spinup zone for precipitation. To reduce this effect in a practical way, the 3D convective rain and snow rates from ICON are interpolated and written to the analysis and boundary condition files so that the COSMO model will be able to convert these rates to mass densities (division by fall speed) and add them either to QR and QS or to QV. |
||||
StatusImplemented in int2lm. PRR_CON and PRS_CON can be interpolated either as 2D or 3D fields. But the corresponding COSMO development is still under testing. |
||||
Technical IssuesRequires ICON to write out two more 3D fields, which is a heavy burden on the disk space! Therefore, the COSMO part should rely on the 2D fields and use a climatological relative profile instead of the full 3D fields. |
||||
Documentation (Presentation of Results, Model Documentation, External Documentation)Provide specific information on availability of the following documentation required by the COSMO source code management rules:
|
||||
2.02 | release | 2014-10-01 | 2015-05-01 | 2015-06-30 |
2.01 | release | 2014-04-01 | 2014-09-01 | 2014-11-30 |
2014-05-30 | Consolidation of GRIB2 and usage of grib_api | WG6 | finish | Uli Schaettler |
DescriptionAfter implementing the new ECMWF grib_api to read and write GRIB1 and GRIB2 data into INT2LM, tests have started at DWD to port all model I/O to GRIB2. These tests showed that still not all GRIB2 meta data are set properly to run a forecast suite, where INT2LM, data assimilation cycle, forecast model and ensemble mode have to be adjusted. Therefore the implementation has to be adjusted. Also, changes in the source code of INT2LM and the COSMO-Model are necessary to deal with the non-existent vertical coordinates. It has been agreed that the height field HHL and the full pressure P are transfered from INT2LM to the COSMO-Model and also within the COSMO data assimilation cycle with higher grib packing rate (24 bits per value). Several modules have to be adapted to eliminate the usage of the vertical coordinate parameters. |
||||
StatusNov 2013: Tests (and work) started at DWD April 2014: Consolidated code has been implemented in DWD versions of INT2LM. July 2014: DWD switched on operational use of GRIB2 in all models. | ||||
Technical IssuesThe code has been implemented according to the COSMO-Standards. Several colleagues at DWD, who were also involved in coding GRIB2 in other models, reviewed the code. Since July, the code runs operational at DWD in a special DWD version. |
||||
TestingThe code has been technically tested in experiments and in the DWD parallel suite, before it was used operationally since July 2014. |
||||
DocumentationA special COSMO page for GRIB2 is in preparation. |
||||
2014-05-30 | SLEVE2 coordinate | WG6 | finish | Guy DeMorsier |
DescriptionA generalized SLEVE (Smooth LEvel VErtical) coordinate with a modified vertical decay of the topographic signature with height is implemented as a new vertical coordinate type in INT2LM and in the COSMO-Model. It can be chosen with the parameter ivctype=4. For a detailed description see Leuenberger, D., M. Koller and C. Schaer, 2010: An improved formulation of the SLEVE coordinate. Mon. Wea. Rev., 138 (9), 3683-3689, DOI: 10.1175/2010MWR3307.1 Basic changes to the INT2LM are in the module vgrid_refatm_utils.f90, in the Subroutines reference_atmosphere_2 and reference_atmosphere_BVconst, where the vertical coordinate values zbk and zbk2 are computed according to the modified formulation: zbk (k) = SINH( (vc_type%vcflat/svc1)**zn - (vc_type%vert_coord(k)/svc1)**zn ) / SINH( (vc_type%vcflat/svc1)**zn ) With zn=1, the original SLEVE coordinate is computed. For the modified SLEVE a value of zn=1.35 is chosen. | ||||
StatusUp to February 2014: The modified SLEVE coordinate has been implemented and tested in local MeteoSwiss INT2LM implementations. March 2014: The code has been given to the Source Code Administrator and to Uli Schaettler for reviewing. |
||||
Technical IssuesThe code has been given to Uli Schaettler (DWD) to review the implementation (4-eyes principle). All coding standards are fulfilled. |
||||
TestingThe modified SLEVE coordinate is only activitated if the Namelist variable ivctype=4 is set. For other values the results are not changed. Tests have been performed at MeteoSwiss. |
||||
DocumentationAn update of the INT2LM User Guide has been provided. |
||||
2014-09-15 | Implementation of ICON to COSMO interpolation | WG6 | finish | Thorsten Reinhardt |
DescriptionInterpolation of the data from DWD's new global model ICON to the COSMO-grid is implemented. ICON uses a similar icosahedral grid as GME, but the implementation is rather difficult. GME still uses regular 2D Fortran arrays to store the variables, but ICON chose a full irregular approach using indirect addressing. Therefore, the computation of the ICON grid is rather complex, and it is not computed during model run, but only pre-computed grids are read when the program starts. INT2LM takes over the necessary routines from ICON to deal with this irregular grid and do the horizontal interpolations to the COSMO grid. A new module src_icon_interpol.f90 is implemented for that issue. |
||||
StatusICON to COSMO interpolation has been implemented in the first half of 2014. Since then it is tested. A preliminary version of INT2LM has been distributed to all COSMO partners mid of October, who whose GME as driving model, to test ICON to COSMO interpolation for their setups. |
||||
Technical IssuesThe COSMO Coding standards are fulfilled. The implementation has been cross-checked by U. Blahak and U. Schättler. |
||||
TestingThe interpolation is tested in DWD's parallel suite. |
||||
DocumentationExtension of INT2LM documentation: How to use ICON data for the COSMO-Model has been documented in a new chapter of the INT2LM User Guide: "Driving Models for the COSMO-Model" |
||||
2.00 | refVersion | 2013-10-12 | 2013-11-01 | 2013-11-08 |
INT2LM 1.24 has been renamed to INT2LM 2.1. For the Release Notes just see Release Notes of 1.24 |
||||
1.24 | test | 2013-10-12 | 2013-11-01 | 2013-11-01 |
This version will be renamed to INT2LM 2.0 |
||||
1.23 | test | 2013-08-12 | 2013-09-30 | 2013-10-02 |
1.22 | test | 2013-04-01 | 2013-06-28 | 2013-07-11 |
2013-08-19 | Implementation of grib_api for reading / writing | WG6 | finish | Uli Schaettler |
DescriptionSince some years, the new GRIB standard (GRIB2), defined by WMO, is ready to be used. ECMWF developed an application programmers interface (grib_api), which can be used within programs to read, create and manipulate grib messages for GRIB2, but still also for GRIB1. At DWD it has been decided to use grib_api as the official GRIB library in the future. Also the COSMO partners agreed to use grib_api. Therefore, grib_api will be implemented in the COSMO-Model and the INT2LM in the next few months. At DWD the implementation is based on using the concept of the shortnames and the tables and definitions provided by grib_api. These tables and definitions have to be coordinated with the COSMO-partners. In INT2LM, grib_api has been implemented in several steps:
|
||||
Status27.05.13: the new general vertical coordinate has been implemented, but there are still issues with grib_api. Version grib_api 1.11.0 is necessary to compile and link INT2LM with grib_api. 28.06.13: implementation finished; functionality tests ongoing 17.06.13: while restructuring the treatment of the vertical coordinate parameters, a bug in the module src_vert_inter_lm has been found:
11.07.13: new version released as a test version |
||||
Technical IssuesCoding Standards: have been met Functionality tests have been performed 4-eyes Assurance: will be done by DWD colleagues with grib_api experience |
||||
TestingSingle Test Cases: have been performed. For GME => COSMO-EU there are no changes of results. The results for COSMO-EU => COSMO-DE, the results do change because of the bug fix in src_vert_inter_lm.f90 |
||||
DocumentationA COSMO-web page has been set up to document the implementation and GRIB2 usage. This page is not yet visible (still work in progress). It will be enhanced in the next few weeks. User Guide has been updated. | ||||
2013-08-19 | Multi-layer treatment of IFS soil | WG6 | finish | Davide Cesari |
DescriptionWith the current interpolation scheme for the IFS soil model, if you use l_multi_layer_lm=.TRUE. and yinput_model='IFS', the soil is first interpolated to the 2-layer soil model levels, and successively interpolated (actually extrapolated) to the multi-layer levels, so input temperature below the COSMO T_CL level gets lost and deep temperature is unrealistic. With a new direct interpolation all the input layers should influence the output, so higher precision and less extrapolation is obtained. |
||||
StatusA new direct interpolation scheme has been implemented. |
||||
Technical IssuesCoding Standards: are met Technical Tests: have been performed The new interpolation is activated when
which was a forbidden combination before. When l_multi_layer_in=.FALSE. the old behaviour is retained. So now the new version can be tested. |
||||
TestingSingle Test Cases, Experiments: can now be done by all interested partners. |
||||
DocumentationSee the release notes of Version 1.22 |
||||
1.21 | test | 2012-12-01 | 2013-03-01 | 2013-03-22 |
1.20 | test | 2012-06-25 | 2012-08-31 | 2012-09-28 |
1.19 | test | N/A | 2012-06-06 | N/A |
1.18 | release | N/A | N/A | 2011-03-11 |
Release Notes | ||||
1.17 | release | N/A | N/A | 2011-03-11 |
Release Notes | ||||
1.16 | release | N/A | N/A | 2010-12-20 |
Release Notes | ||||
1.15 | release | N/A | N/A | 2010-12-10 |
Release Notes | ||||
1.14 | release | N/A | N/A | 2010-11-19 |
Release Notes | ||||
1.13 | release | N/A | N/A | 2010-10-13 |
Release Notes | ||||
1.12 | release | N/A | N/A | 2010-06-14 |
Release Notes | ||||
1.11 | release | N/A | N/A | 2010-04-23 |
Release Notes | ||||
1.10 | release | N/A | N/A | 2009-12-17 |
Release Notes | ||||
1.09 | release | N/A | N/A | 2009-09-03 |
1.08 | release | N/A | N/A | 2009-05-29 |
Release Notes | ||||
1.07 | release | N/A | N/A | 2007-11-26 |
Release Notes | ||||
1.06 | release | N/A | N/A | 2007-09-07 |
1.05 | release | N/A | N/A | 2007-07-09 |
1.04 | release | N/A | N/A | 2006-01-03 |
1.03 | release | N/A | N/A | 2005-12-12 |
1.02 | release | N/A | N/A | 2005-07-22 |
1.01 | release | N/A | N/A | 2005-04-11 |