Last updated: 2022
The following table shows which developments have been or are implemented into cosmo and int2lm, per working group.
WG subtask |
Start date | Task status | Expected delivery |
Task subject |
Responsible person |
Target version |
|
---|---|---|---|---|---|---|---|
WG1 | |||||||
1.5 | 2016-01-01 | finish | 2016-01-01 | changes to SPPT scheme in order to improve stability | Walser (MeteoSwiss) | 5.05 | |
DocumentationProvide specific information on availability of the following documentation required by the COSMO source code management rules:
|
|||||||
1.2 | 2015-12-16 | finish | 2016-01-15 | removal of AOF reading | Schraff (DWD) | 5.04 | |
Descriptionremoval of option to read observations from AOF files, to be done when AOF reading is not used any more by any of the COSMO members / partners |
|||||||
1.1 | 2011-09-01 | test | 2015-12-01 | refined assimilation of T-2m observations | ARPA-Piemonte | N/A | |
DescriptionThe constant increase of computing power allows to run progressively higher resolution numerical weather prediction models in the operational chains. Italian reference operational chain for limited area numerical weather prediction is represented by the COSMO-I7 and COSMO-I2 congurations of the COSMO model, being the COSMO-I2 the finer resolution one, with a grid spacing of about 2.8 km. The initialisation stage of the COSMO chains is deficient in two points: on the one hand it is missing a proper initialisation of the land surface energy and water exchanges, on the other hand it does not make use of all the information coming from the high density networks of weather stations. These two problems are assessed with the proposal to implement an analysis cycle for the COSMO-I2. In order to cope with the problems at the land surface, an explicit coupling scheme between the existing COSMO data assimilation code, which directly influences only the atmospheric state, and the description of the soil state has been introduced. This coupling scheme, called FASDAS, acts as a new balance relation in the assimilation framework. |
|||||||
The code has been prepared for version COSMO v5.0. |
|||||||
Technical Issues (Coding Standards, 4-eyes Assurance)The modifications required by the inclusion of FASDAS have been made according to the COSMO coding standards. The code has ben given to the WG coordinator C. Schraff. The code has originally been developed by Marco Galli (now at U.S.A.M.). Matteo Giorcelli (ARPA Piemonte) reviewed the code in details. He also has rewritten some parts to reach full standards compliance and to ensure ”COSMO quality level”. The code referent person will be Massimo Milelli (ARPA Piemonte). |
|||||||
TestingTests performed on COSMO-I2 and COSMO-I7 configurations, over north-west Italy (firstly) and over the whole domain (secondly). Results have been published (links below). |
|||||||
DocumentationTwo articles in the COSMO Newsletter have been published: - model/documentation/newsLetters/newsLetter11/1_galli.pdf - model/documentation/newsLetters/newsLetter14/cnl14_02.pdf Also a few presentations have been shown during the GMs: - consortium/generalMeetings/general2011/wg1/wg1_3_galli_t2m_ts.pdf - consortium/generalMeetings/general2012/wg1-kenda/wg1_galli_t2m.pdf - consortium/generalMeetings/general2014/wg1-kenda/fasdas_giorcelli.pdf Another presentation given during the CUS 2012: |
|||||||
N/A | N/A | finish | 2012-12-13 | optional new QC check for surface pressure against lateral BC fields | Schraff (DWD) | 4.28 | |
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 |
|||||||
1.3 | N/A | finish | 2099-12-31 | assimilation of radar-derived rain rates by a 1DVAR + nudging approach | Poli (ARPA-EMR) | 0.00 | |
Description
|
|||||||
tried to resolve the problem of the bias between precipitation generated by COSMO model and the forward model of 1D-Var modifying some parameters in the 1D-Var algorithm
This task has been stopped (Sept. 2013). |
|||||||
1.4 | N/A | idea | N/A | optimisations from project OPCODE | Leuenberger (MeteoSwiss) | N/A | |
WG2 | |||||||
0.15 | 2015-01-01 | finish | 2015-07-01 | Increase of stability of the numerical diffusion scheme for 3D turbulence | Baldauf (DWD) | 5.03 | |
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. |
|||||||
0.14 | 2015-01-01 | work | 2016-03-01 | Kaas (2008)-conserving extension of the Semi-Lagrangian Advection | deMorsier (MeteoCH) | N/A | |
DescriptionThe current implemented Semi-Lagrangian tracer advection scheme shall be extended by the idea of the Kaas (2008) Tellus A paper. Aim: conserving advection scheme. |
|||||||
StatusThis work firstly has been started at the EMPA (ETH) institute. A first visit of Eigil Kaas has been taken place, but it was recognized that the implementation will take a longer time. Guy deMorsier has been willing to overtake this work. |
|||||||
0.11 | 2013-12-15 | finish | 2014-03-31 | Optional targeted diffusion to avoid cold pools in narrow valleys | Baldauf (DWD) | 5.01 | |
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 |
|||||||
0.13 | 2013-12-15 | finish | 2014-03-31 | Remove hacks in the tracer module | Fuhrer, Roches (MeteoCH) | 5.01 | |
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 |
|||||||
0.09 | 2013-01-01 | finish | 2014-06-30 | Extension of the Bott tracer advection scheme | Schneider (Univ. Bonn), Blahak | 5.06 | |
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 |
|||||||
0.10 | 2013-01-01 | finish | 2014-06-30 | Adaptation of the RK dynamical core for stochastic physics | Torrisi (CNMCA), Baldauf (DWD) | 5.01 | |
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 |
|||||||
N/A | 2012-12-04 | work | 2013-12-31 | High order horizontal schemes for the Runge-Kutta dynamical core | [CLM] Andreas Will | N/A | |
Description3rd to 6th order horizontal discretisation error convergence schemes are implemented in the Runge-Kutta dynamical core of COSMO. The steps planned are the following: 1. Theoretical comparison of the advection term formulations of Baldauf (standard COSMO adv-schemes), Wicker-Skamarock (as in WRF) and Morinishi (1998) (symmetric, conservative form). 2. Convergence properties of RK dynamical core with 2nd to 6th order advection (Baldauf) 2. 4th order scheme for the 2D mountain flow test case terms 3. 4th order scheme for propagating waves test case terms 4. 4th order scheme for real case dynamics terms 5. 3rd to 6th order schemes for (1) to (3) |
|||||||
Status1. CLM_CO informed 2. Releaed version used: COSMO_4.23 3. Implementation of 4th order scheme ongoing. |
|||||||
Technical IssuesCoding Standards, Technical Test Suite, 4-eyes Assurance (Michael Baldauf) |
|||||||
TestingSingle Test Cases: 2D mountain flow, 2D wave propagation (ongoing) Experiments: Standard evaluation run (planned) |
|||||||
Documentation1. Reporting: Results for 1,2 presened at COSMO/CLM Seminar 2012 2. Model Documentation 3. External Documentation |
|||||||
N/A | 2012-12-04 | idea | 2013-12-31 | 2-way nesting ECHAM-COSMO | [CLM] Andreas Will | N/A | |
DescriptionA two-way nesting between ECHAM6 and COSMO-CLM is being developed. The following working steps and COSMO model developments are planned: 1. OASIS-MCT interface for ECHAM6, NEMO and Community Land Model 2. vertical interpolation of ECHAM6 fields at every ECHAM time step in COSMO by implementation of the int2lm procedure, simplified to one step interpolation 3. exchange of the 3D fields with ECHAM5 |
|||||||
Status1. CLM_CO informed 2. Released version: cosmo_4.8_clm19 and int2lm_1.10_clm11 3. Simlified vertical interpolation implemented in int2lm, implementation in cosmo ongoing |
|||||||
Technical IssuesCoding Standards, Technical Test Suite, 4-eyes Assurance |
|||||||
Testing1. Single Test Cases: comparison int2lm original and simplified (successfully done) 2. Experiments: Central America to northern West-Atlantic region 10y period 1999-2009. |
|||||||
N/A | 2012-12-01 | idea | 2013-06-30 | Optimal relaxation function for lateral boundary conditions | [CLM] Andreas Will | N/A | |
DescriptionThe relaxation function for the lateral boundaries is optimized in order to reduce the reflection of waves at position rlwidth, which is the inner boundary of the relaxation zone. Hereto different exponents of the current exponential function and different functions of order (cos)**n are implemented and tested. Steps of development are: 1. Introduce and test different exponents of exp-function 2. Introduce and test different powers of cos-function |
|||||||
Technical Issues1. CLM_CO informed 2. Released version used: cosmo_4.18 |
|||||||
Testing1. Single Test Cases: 2 D mountain flow test case successfull with exp(-10 x) 3. Experiments: |
|||||||
N/A | 2012-09-30 | finish | 2013-04-30 | Consolidation work in the New fast waves solver | Baldauf (DWD) | 4.27 | |
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) |
|||||||
N/A | 2012-09-10 | work | 2100-12-01 | New implementation of 3D diffusion in terrain following coordinates | Oliver Fuhrer (MeteoCH) | N/A | |
DescriptionThe actual calculation of horizontal turbulent fluxes in terrain following coordinates suffers from stability restrictions in steep terrain. An alternative consits in the calculation on a z-coordinate. |
|||||||
0.07 | 2012-06-30 | finish | 2012-09-01 | Bug fix in advection operators (Bott-scheme) and cloud diag | Fuhrer (MeteoCH) | 4.23 | |
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 IssuesCoding 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. |
|||||||
TestingThe 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/. |
|||||||
DocumentationPresentation of results: N/A (bug) Provide extensions / modifications for the Documentation System: N/A (bug) Future Code responsible: Michael Baldauf |
|||||||
0.06 | 2012-03-31 | work | 2013-03-31 | Mahrer-discretization of horizontal pressure gradients in fast waves (current and new version) | deMorsier (MeteoCH), Baldauf | N/A | |
DescriptionMahrer (1984) MWR approach; possibly with some improvements: quadratic terms to improve extrapolation near the bottom boundary (Zaengl (2012) submitted to MWR) |
|||||||
Technical Issues4-eyes: in current fast waves solver (Baldauf in new fast waves solver (deMorsier) |
|||||||
0.08 | 2012-03-01 | finish | 2012-09-01 | calculation of the moisture divergence | Fuhrer (MeteoCH) | 4.23 | |
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 IssuesCoding 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 |
|||||||
Testing3) 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). |
|||||||
DocumentationPresentation of results: not applicable Changes to Model Documentation: not applicable Changes have been documented in the Release Notes Future Code responsible: Michael Baldauf |
|||||||
0.03 | 2012-02-15 | finish | 2012-09-01 | New tracer advection schemes: MPDATA and 'Bott2_Strang' | deMorsier (MeteoCH), Baldauf ( | 4.23 | |
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 IssuesCoding 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 |
|||||||
TestingAt DWD an experiment is ongoing, where the efficient Strang-splitting version is tested. |
|||||||
DocumentationCurrent state of MPDATA: see presentation at COSMO GM 2011 by G. deMorsier |
|||||||
0.04 | 2012-01-21 | finish | 2012-09-01 | Reorganisation of metrics stuff | Baldauf (DWD) | 4.23 | |
Descriptionconcentrate metric terms for the spherical, terrain following grid in an own module. This will facilitate the coupling of other dynamical cores (EULAG, ...). |
|||||||
Technical IssuesCoding 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 |
|||||||
TestingThis is only a collection of subroutines from a diversity of files into one common file. Therefore no change of results |
|||||||
0.01 | 2011-09-30 | finish | 2013-01-31 | New Fast Waves Solver for RK dynamical core | Baldauf(DWD) | 4.24 | |
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 IssuesCoding 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 |
|||||||
TestingAll 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. |
|||||||
DocumentationBaldauf (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. |
|||||||
0.02 | 2011-09-01 | finish | 2011-12-06 | Horizontal nonlinear Smagorinsky diffusion | Baldauf (DWD) | 4.21 | |
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 Issues4-eyes: Uli Schaettler Code responsibility: M. Baldauf |
|||||||
TestingExperiments with COSMO-DE during 01.-28. Feb. 2011 and 01.-31. Aug. 2011 were performed. Verification results are mainly neutral as expected. |
|||||||
0.05 | 2008-09-01 | work | 2013-09-30 | Conservative Dynamical Core | PP participants | N/A | |
Descriptionsee PP CDC project plan and PP CELO project plan |
|||||||
Technical Issues4-eyes: code responsibility: |
|||||||
Documentationpapers in Acta Geophysica (2011) 59/6: Michał Z. Ziemiański, Marcin J. Kurowski, Zbigniew P. Piotrowski, Bogdan Rosa und Oliver Fuhrer: Bogdan Rosa, Marcin J. Kurowski und Michał Z. Ziemiański: Marcin J. Kurowski, Bogdan Rosa und Michał Z. Ziemiański: |
|||||||
WG3a | |||||||
1.08 | 2017-10-31 | submit | 2020-08-31 | TERRA-related issues of ConSAT: implicit surface temperature, surface-proceses, canopy and snow | Raschendorfer | N/A | |
Description7. Implementation of an extra, massive roughness layer, which is imlicitly coupled to the soil and semi-implicitly coupled to the atmosphere via the temperature of rouhness elements. This includes the parameterization of transfer of heat, moisture and radiation through this layer as well as the heat storage of elevated roughness layer elements (canopy), which coveres the surface of the dense soil. For that purpose, the surface scheme of TERRA needs to be modified considerably, in particular with respect to the effect of evaporating or tranpirating surface fractions and the upper boundary condition for the implicit solution of the soil heat conduction equation. This should reduce the currently too excessive soil heat fluxes and possibly help to increase numerical stability related to near surface processes. That measure is also closely related to the treatment of snow in general and particularly with regard to snow-interception by roughness elements. 8. Formulation of surface processes as implicit functions of surface temperature. This includes the solution of coupled multi-layer heat budgets for the soil, a possibly fractional snow-pack on top of the compact soil and a roughness cover (canopy), where the latter may reduce to a massless skin layer. This implicit formulation is a prorequisite of for a consistent and numerically save description of near-surface processes, including various phase transitions of hydrometeors at the surface and the physical effect of elevated roughness elements. This formulation is also necessary in order to avoid excessive time-step oscillations of near-surface tempertures as large time-steps, as it is currently the case with the operational code. |
|||||||
Status notesIssue 7: A first working version of the elevated roughness layer (canopy) has been implemented into the previous (non-blocked) vesion 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. The implicit coupling of all heat-equations under the particular consideration of snow, however, requires a rather extensive refomulation of TERRA and associated modules; This implicit coupling is also closely related to a formulation of surface processes as implicit funtions of surface temperature, which seems to be necessary in order to avoid undesired temperature oscillations at large time steps. Therefore it has been decided to realize an implicit formulation of thermal surface processes as a first step based on the blocked and further tuned TERRA version running with ICON (s. Issue 8). Thereafter, the above mentioned development of a particular canopy layer is going to be adapted as an extension of the coupled skin-layer equation of the snow-free surface. Finally, a canopy layer (with possible snow-interception) is planned to be considered also above a snow-pack covering the compact soil. Issue 8: A first version of primary surface processes, which are formulated as an implicit function of surface temperature, is present in a private development branch for ICON. This includes major modifications in TERRA and its interface routines as well as related modifiactions in TURBTRAN including its interface routines. Since the official code had been written in terms of an explicit treatment of surface temperature, various related limitations and restrictions could now be removed or had to be adapted to an implicit treatment of surface processes. For that purpose large parts of the code had to be extensively altered or even to be rewritten completely, particularly with regard to phase transitions of precipitation reaching the ground, to the sequence of simulated processes, to the formulation of the single-layer snow model and to the treatment of dynamic snow tiles. Further, rather elaborate adaptations were necessary, in order to adapt various empirical parameterizations (that had been implemented into the ICON code in the meantime by G. Zängl) to the more rigorous implicit treatment. In this TERRA version, a single coupled system of linearized heat budgets is solved for the following three domains: a snow-free skin layer, a partial (single-layer) snow-cover and the rigid soil. With all these implementations we aim to keep as close as possible to the results of the operational formulation in order to have optimal control about the various formal, methododical and numerical modifications. Related physical extensions with expected positive impact to the model scores will be implemented as next steps. Various test-runs demonstrate that near-surface temperature oscillations are eliminated with only minor further impact to the results. Thes even holds, if the so called flux-limiter is deactivated, which is necesary in the operational version in order to avoid model crashes due to the mentioned oscillations, is now deactivated! Nevertheless, some further tests are pending. July 2019: The coupled system of linearized heat equations has been extended towards a multi-layer heat equation for snow with an individual, time-step dependent number of snow-layers for each grid point. This has been implemented as a ConSAT-contribution to PT SAINT, where the latter is aiming to introduce the sophisticated multi-layer snow hydrology from the snow-model SNOWPACK. After merging the private development branch (containing the implicit surface-temperature treatment) with the running ICON-development and after further testing, the result is expected to be committed to the official development branch of ICON soon (at least by the end of 2019). |
|||||||
TestingThe canopy-extension has been tested with COSMO for a cloud-free summer-case over Germany with very promising results, particularly with respect to 2m dew-point temperature. The various implementations for the implicit treatment of surface tempertature have been tested technically and through several cases with ICON, showing (so far) the expected behaviour, that means: Vanishing time-step oscillations of near surface temperature, even without the so far necessary flux-limiter in TERRA. The various empirical extensions and extra limitations (the latter in order to damp inherent numerical instability) that had been implemented in the operational TERRA before, could be adapted so as to obtain similar results with the implicit version compared to the operational version. |
|||||||
DocumentationProvide specific information on availability of the following documentation required by the COSMO source code management rules:
|
|||||||
3.07 | 2016-05-01 | work | 2016-09-30 | Implementation of ICON/COSMO two-moment microphysics scheme | Seifert/Blahak | 5.06 | |
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 |
|||||||
DocumentationProvide specific information on availability of the following documentation required by the COSMO source code management rules:
|
|||||||
1.07 | 2015-10-28 | submit | 2015-12-01 | Soil and surface fluxes: effective slope-angle correction factor (1/cos(slopeangle) | Blahak/De Morsier/Raschendorfe | 5.06 | |
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. |
|||||||
0.02 | 2015-10-28 | submit | 2015-12-01 | Stochastic boundary layer perturbations reflecting unresolved effects of variability of surface fluxes on convection initiation | Kober / Blahak | 5.06 | |
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) |
|||||||
2.01 | 2015-10-28 | finish | 2015-12-01 | Shallow convection: new closure based on conv. vertical velocity scale, modified entrainment, modified shallow conv. cloud cover | Boeing / Blahak | 5.05 | |
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) External Doc: shcu.pdf |
|||||||
1.06 | 2014-04-01 | finish | 2014-08-30 | Turbulence: Implementing and extending the restructured ICON-version of TURBDIFF | Raschendorfer | 5.05 | |
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:
|
|||||||
StatusDecember 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
|
|||||||
Testing
|
|||||||
Documentation
|
|||||||
3.06 | 2014-03-01 | finish | 2014-06-30 | Implementation of a unified COSMO-ICON microphysics | Kahler / Fuhrer / Lapillonn | 5.01 | |
DescriptionThe unified COSMO-ICON model microphysics includes different changes conducted by Felix Rieper and Gunther 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. |
|||||||
1.02 | 2014-02-01 | finish | 2014-12-31 | Turbulence: 3D-extensions for TURBDIFF within the concept of scale separation | Raschendorfer / Blahak | 5.01 | |
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 IssuesThe 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: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. |
|||||||
0.01 | 2013-12-01 | submit | 2014-07-31 | Providing COSMO-SC for general application | Raschendorfer / Schaettler | 5.05 | |
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 IssuesThe 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. |
|||||||
TestingThe 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). |
|||||||
DocumentationA documentation of the SC-framework is already available. However, it is not completely up to date. |
|||||||
1.05 | 2012-04-01 | finish | 2015-09-30 | Consolidation of surface-to-atmosphere transfer ConSAT | Raschendorfer | 5.05 | |
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 IssuesAll 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! |
|||||||
TestingAll 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. |
|||||||
DocumentationA full documentation of the TKE-based transfer scheme has started as well as the preparation of publications. |
|||||||
1.04 | 2012-03-31 | finish | 2014-06-30 | Turbulence: Implementing TKE-advection for TURBDIFF | Blahak | 5.01 | |
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 IssuesThe code modifications are along the coding standards and they have been inspected by M. Raschendorfer and J. Farstner. 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. |
|||||||
TestingThe 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. |
|||||||
DocumentationA 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. |
|||||||
3.05 | 2012-01-01 | finish | 2014-01-01 | Microphysics: Improved supercooled liquid water content | A. Seifert / U. Blahak | 5.01 | |
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. |
|||||||
StatusThe 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 Garsdorf (DWD). |
|||||||
1.03 | 2011-11-30 | idea | 2012-03-31 | Turbulence: Smagorinsky LES turbulence parameterization | N/A | N/A | |
DescriptionThis work package would allow to introduce this standard parameterization as a further optional turbulence scheme for very high resolution LES-type runs (dx < 1 km), for research and development purposes. Features: Full implementation of all metrical terms, which is an improvement over the already implemented Smagorinsky-type LES from Hans Herzog (itype_turb=7 or 8). Documentation: this parameterization will be described in the next COSMO News Letter; could serve as a basis for the COSMO Latex documentation. |
|||||||
Technical Issues4-eyes: will be done by Michael Baldauf |
|||||||
4.03 | 2011-10-31 | finish | 2012-01-31 | Radiation: Zenith Angle Corrections | Reinhardt (DWD/AGeoBw) | 4.22 | |
3.01 | 2011-08-31 | work | 2012-12-31 | Modified saturation adjustment | Blahak/Seifert | N/A | |
DescriptionNew saturation adjustment through isochoric processes instead of isobaric processes (with problems of mass conservation). Implementation in 4.16 as a test version. |
|||||||
Technical IssuesFor the pressure interpolation problem: use hydrostatically diagnosed pressure as the basis for interpolation. For the super-saturation problems: Change diagnosis of QVsat from the previous formulation to QVsat = E(T) / (rho*R_v*T) |
|||||||
TestingLonger term tests revealed subtle problems in conjunction with strong super-saturation in the upper troposphere and technical problems with the interpolation of model variables to p-surfaces, because it can now happen that the pressure is no more strictly decreasing with height, but has some small increasing "kinks" due to the pressure diagnosis after saturation adjustment. Solutions: known, but I had no time to implement them up to now. |
|||||||
3.02 | 2011-08-31 | finish | 2012-12-01 | Microphysics: 2-Moment Scheme | Seifert / Blahak | 4.25 | |
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 IssuesThe 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. |
|||||||
TestingThe scheme has been developed now since several years and has been tested extensively. It has lately been presented to a workshop at ETH Zurich. 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. |
|||||||
DocumentationThere are already 5 published papers (peer review and conferences), more are in preparation. |
|||||||
3.03 | 2011-08-31 | finish | 2013-08-31 | Microphysics: Melting Snow | Seifert / Frick (DWD) | N/A | |
DescriptionMelting of snow using explicit melt-water on snowflakes by Claudia Frick and Axel Seifert. Still under development. The initial development has been successfully completed and an implementation of the scheme which works in a 3D real case COSMO simulation is available. During the next months the scheme is going to be tested, evaluated and maybe some calibration of parameters will be done based on the case studies (aka tuning). Some journal papers on the new scheme is currently being written, and will hopefully be submitted within this year. |
|||||||
3.04 | 2011-08-31 | finish | 2013-03-31 | Microphysics: Cirrus clouds | Seifert / Kahler | N/A | |
DescriptionImproved cirrus cloud microphysics by Carmen Koehler and Axel Seifert. The initial development has been successfully completed and an implementation of the scheme which works in a 3D real case COSMO simulation is available. During the next months the scheme is going to be tested, evaluated and maybe some calibration of parameters will be done based on the case studies (aka tuning). Some journal papers on the new scheme is currently being written, and will hopefully be submitted within this year. |
|||||||
Is being tested in a unified COSMO-ICON version. |
|||||||
Technical Issues |
|||||||
1.02 | 2011-08-31 | idea | N/A | Turbulence: Clean up namelist switches | Raschendorfer | N/A | |
DescriptionClean up the different options: itype_turb, lturb, imode_turb Remove unused options, e.g. lexpcor, lprfcor, tmpcor Work started to re-structure the organization of turbulence because of unified COSMO-ICON physics package (cf. "Implementing and extending the restructured ICON-version of TURBDIFF"). |
|||||||
4.02 | 2011-05-31 | finish | 2012-09-01 | Changes in parameterization for ice cloud cover in radiation | Seifert | 4.23 | |
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. |
|||||||
4.01 | 2011-01-01 | work | 2013-12-01 | Improved cloud radiation coupling | Blahak / Ritter | N/A | |
Description- Optical properties of hydrometeors depending on their effective radius - Including snow, graupel and rain in the radiative calculations - Effective factor to consider subgrid scale variability within clouds - Radiative properties of SGS clouds - Improved consideration of aerosols for radiation and microphysics |
|||||||
StatusConsolidated test version implemented in COSMO 4.22 Continued in PP T2(RC)2. |
|||||||
TestingA number of new namelist parameters have to be investigated by extensive sensitivity studies based on 4 real cases (COSMO-DE). If possible, the number of new namelist parameters has to be reduced to the absolutely necessary minimum. This work is under way. Lots of tests are performed within PP T2(RC)2 by use of data at the observatories at Moskow State University and Lindenberg. |
|||||||
WG3b | |||||||
1.12 | 2014-09-01 | finish | N/A | Improved representation of the surface temperature including the effects of vegetation | J.P. Schulz (DWD) | N/A | |
Land surface processes have a significant impact on near-surface |
|||||||
The new scheme was implemented in TERRA offline and also in ICON-TERRA. See also here. |
|||||||
Not yet done |
|||||||
The new scheme was tested in TERRA in offline mode, using Lindenberg |
|||||||
Schulz, J.-P. and G. Vogel, 2017: An improved representation of the land surface temperature including the effects of vegetation in the COSMO model. Geophysical Research Abstracts, 19, EGU2017-7896. Viterbo, P., and A. C. M. Beljaars, 1995: An improved land surface parameterization scheme in the ECMWF model and its validation. J. Climate, 8, 2716–2748. |
|||||||
1.08 | 2013-01-01 | finish | 2015-12-31 | Revision of transpiration and root parameterization | G.Vogel (DWD) | N/A | |
Goal: For the sake of consistency, the parameterization of root characteristics (vertical root profile, rooting depth) should also be updated if satellite-based LAI and plant cover values are used in the NWP models. In particular, a more realistic model approach of the annual cycle of the water uptake processes by the roots is needed. The widely used parameterization introduced by Arora and Boer (2003) seems appropriate for this task. However, some adaptation work is needed for extending the current formulation of the transpiration calculation in the present TERRA model. Moreover, sensitivity studies will be made for testing its reliability under various soil and vegetation conditions. Description: In order to improve the vegetation impact on the energy and water transfer between surface and atmosphere, the annual cycles of the vegetation properties such as the LAI, the plant cover and the rooting depth should be carefully considered. According to offline evaporation studies with the TERRA module for the Falkenberg field site (DWD, MO Lindenberg) the rooting depth turned out to be the predominant quantity. Since many years simple relationships are used in the current weather prediction models which prescribe the vegetation yearly cycles in a fairly artificial but consistent manner. But, if the parameterised LAI and plant cover are replaced by climatologic SEAWIFS-based values (e.g. in the operational GME predictions), the treatment of the rooting depth should not remain unchanged; an inconsistent use of the SEAWIFS data ignores completely the close allometric relationships between leaf and root biomass. To also improve the root parameterisation it is proposed to adapt an approach by Arora and Boer (2003). It enables to describe the evolving root density profile, the cumulative root fraction and the rooting depth as a function of the varying root biomass. The needed root biomass is available from the literature for the relevant vegetation kinds. An appropriate scaling of the root biomass with the already used SEAWIFS-based NDVI-index provides then a maximum consistency between all vegetation properties. Tasks: > Conceptual studies for adapting the current transpiration scheme to the root parameterization proposed by Arora and Boer (2003) > Refinement of a coupling scheme of rooting profiles with satellite-based/parameterized LAI values > Tentative implementation of the modifications into the offline TERRA code (based on COSMO release 4.11) and evaluation with offline runs > Documentation of the modified parameterizations of rooting depth and transpiration Resources: Work done in collaboration with JP. Schulz (Uni Frankfurt). |
|||||||
Status presented at the COSMO GM in Sibiu, 09.2013: - The root parameterization should only be improved in the COSMO models together with the soil heat conduction and if the shading effect of the vegetation is considered. - A time-constant exponential root profile is sufficient in order to simulate the soil moisture development in the annual cycle. - At present, the root parameterisation by Arora and Boer (2003) fits the annual soil moisture cycle at best. Moreover, it is open for future develop-ments. |
|||||||
Preliminary study has been done and shown at the COSMO GM 2012 (see consortium/generalMeetings/general2012/wg3a-wg3b-utcs.htm). |
|||||||
1.11 | 2012-12-01 | finish | 2015-08-31 | Revision of TERRA to support HWSD data | J.Helmert (DWD) | 5.05 | |
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. |
|||||||
3.02 | 2012-10-01 | finish | N/A | Veg3D Coupled with OASIS | [CLM] M. Breil (KIT) | N/A | |
VEG3D coupled via OASIS3-MCT to COSMO_4.8_CLM19 and COSMO_4.21_clm2. |
|||||||
Implemented and in testing phase |
|||||||
1.05 | 2012-09-10 | work | 2015-08-31 | Revision of rainfall interception | J.Helmert (DWD) | N/A | |
Work done in collaboration with G.Vogel (DWD) The interception of rain water plays an important role in simulating soil moisture and evapo-transpiration, in particular in case of high vegetation. The following actions are planned:
|
|||||||
Offline tests have shown encouraging results. Some improvements from SURFEX have been implemented in ICON-TERRA, but snow interception is still missing and complicated. |
|||||||
Offline tests have shown encouraging results. |
|||||||
1.03 | 2012-03-01 | stop | N/A | Tile approach to support partial snow cover | E.Machulskaya (DWD) | N/A | |
See task 1.04. |
|||||||
The experiments with snow as a tile have shown that the overall scores of the model are not definitively better; in the regions with partial snow cover the 2m temperature is predicted more accurately, but in general the winter temperature became lower, which is not what we would like to have. Recently there was a discussion at the NWP Working Group at DWD concerning exactly the problem that the tile approach is designed to improve: this March 2013, in the regions with partial snow cover, there was a noticeable underestimation of the 2m temperature. This period will be re-computed with the tile approach, and if the results are positive, a decision could be taken to transfer these developments in the official code. This development is available in the ICON model only. It will be available to COSMO with the shared COSMO / ICON physics library; however, no time line exists yet for this particular feature (it requires substantial modifications of the interface in the COSMO code). |
|||||||
A partial implentation of the tile approach with two tiles for snow-covered and snow-free regions is implemented and being tested. |
|||||||
2.07 | 2012-01-01 | work | 2015-10-31 | Parameterization of urban effects | [CLM] H. Wouters (KU Leuven) | 6.00 | |
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
An action has started to incorporate this development in the official COSMO code (coordination U.Blahak) |
|||||||
DocumentationProvide 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. |
|||||||
4.02 | 2011-10-01 | finish | 2013-12-31 | Community Land Model coupled with OASIS | [CLM] E. Davin (ETHZ) | N/A | |
Community Land Model version 3.5 (CLM3.5) and version 4.0 (CLM4.0) coupled via OASIS3-MCT to COSMO4.8-CLM19. |
|||||||
Implemented and in testing phase. |
|||||||
1.02 | 2011-09-30 | stop | N/A | Multi-layers snow model | J.Helmert (DWD) | N/A | |
The multi-layer snow model differs mainly in two points from the current one-layer snow model. These are, 1) an arbitrary number of layers in snow instead of one bulk layer and 2) the possibility of water phase changes, existence of liquid water content, water percolation and refreezing within snowpack. The explicit vertical stratification (multi-layer structure) of various properties of snow (temperature, density etc.) allows a more correct representation of the temperature at the soil-snow and snow-atmosphere interface which is important for calculation of snow melting rate and surface turbulent fluxes. The accounting for liquid water and water phase changes within snowpack allows a more accurate calculation of the evolution of the snow properties, in particular, snow water-equivalent depth and snow density, which in turn determines snow heat conductivity. |
|||||||
A first reasonably stable implementation is available in COSMO 4.25 version; however, this implementation has still some problems (see below). Two small bugs were found during February 2012, which however did not affect the results (could be seen in the ICON model only so far); one more numerical stability issue has been corrected in September 2012. Additionally, one security check was added in January 2012 in order to eliminate the consequences of grib-degrib procedures that lead to the loosing of accuracy, internal inconsistency of values and finally to the crash of the model. One additional bug fix has been given to COSMO SCA for inclusion in a future release (in the case of partial snow cover, latent heat flux due to rain freezing was added not only to the snow surface heat balance, but to the heat balance of free of snow surface too. However, if soil temperature is above freezing point, this flux is zero over snow-free soil because rain is not freezingin this case – in some situations this can be important). The multi-layer snow model has been tested within the 3d COSMO model at DWD including the full assimilation cycle. As expected, the explicit representation of the vertical profile of temperature within snowpack leads to the decreasing of the snow surface temperature at least during nights which is detremental for current COSMO model scores. Furthermore, the presence of liquid water leads to slower melting of snow, and, as a result, to the lower surface temperatures in general. It increases the cold temperature bias during the winter, although the snow height itself is modeled better. The origin of the cold bias and the possibilities to reduce it have been investigated. [Juergen H., 23.01.2017] The scheme will be considered with high priority in the next weeks |
|||||||
A parallel experiment is running at DWD with all bugs corrected. The new snow model seems to maintain more snow as compared to the current snow model; further analysis of the results is needed. |
|||||||
1.01 | 2011-09-01 | finish | 2017-09-01 | Mire parametrization | I.Rozinkina (RHM) | 5.06 | |
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 |
|||||||
4.01 | 2009-01-01 | finish | N/A | Community Land Model coupled as subroutine | [CLM] E. Davin (ETHZ) | N/A | |
Community Land Model version 3.5 (CLM3.5) coupled to COSMO4.0 and COSMO4.8-CLM11 References: Davin, E. L., R. Stoeckli, E. B. Jaeger, S. Levis and S.I. Seneviratne (2011), COSMO-CLM2: A new version of the COSMO-CLM model coupled to the Community Land Model, Clim. Dyn., 37, 9, 1889-1907, doi: 10.1007/s00382-011-1019-z. |
|||||||
Implemented and evaluated. | |||||||
N/A | 2015-06-01 | finish | 2015-12-01 | DFI: Better treatment of clouds / precipitation | Lucio Torrisi | 5.04 | |
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 |
|||||||
N/A | 2015-03-01 | finish | 2015-07-31 | Consolidation of GRIB2 and usage at other centers | Darte Liermann, Ulrich | 5.03 | |
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. |
|||||||
N/A | 2013-12-05 | finish | 2014-02-05 | Online Trajectory Module | Annette Miltenberger (ETH) et | 5.01 | |
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
|
|||||||
N/A | 2013-11-11 | finish | 2013-12-20 | GRIB2 and general vertical coordinates | U. Schaettler | 5.01 | |
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. |
|||||||
N/A | 2013-02-10 | work | 2100-12-01 | Unified interface for OASIS3/OASIS3-MCT (land, ocean, GCM couplings) | [CLM] J. Brauch (DWD) | N/A | |
DescriptionA unified interface for the OASIS3/OASIS3-MCT coupler needs to be implemented to allow couplings with multiple components (land, ocean, GCM...). OASIS-MCT coupler is the new parallelised coupler supporting 3D field exchange. |
|||||||
0.09 | 2012-09-10 | finish | 2014-09-10 | GPU Version of COSMO | PP POMPA participants | 5.06 | |
0.11 | 2012-01-31 | finish | 2012-09-01 | Additional values in (long) meteograph output | Reinhardt (DWD) | 4.23 | |
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 IssuesCoding 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). |
|||||||
TestingOnly functionality tests are needed. |
|||||||
DocumentationThere was no need to update the User Guide, because the additional soil and surface variables are not explained in detail. |
|||||||
0.12 | 2012-01-31 | finish | 2012-09-01 | Smoothing of fields on p-/z-levels | Uli Blahak | 4.23 | |
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:
First modification:
The actions of l_z_filter / l_p_filter are the same as before |
|||||||
StatusThe modifications have been implemented in Version 4.23 and could be tested by everybody. No complaints up to now. |
|||||||
Technical IssuesCoding 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) |
|||||||
TestingOnly functionality tests needed. |
|||||||
DocumentationThe User Guide has been updated accordingly. |
|||||||
0.15 | 2012-01-15 | finish | 2012-09-01 | Time measurement for long simulations: Problems with memory | Uli Schaettler | 4.23 | |
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 IssuesThis is a pure technical change. The functionality has been tested with some standard runs. |
|||||||
TestingOnly check the functionality. |
|||||||
DocumentationNo need to update existing documentations. The change has been documented in the Release Notes. |
|||||||
N/A | 2011-10-31 | finish | 2012-12-01 | Asynchronous NetCDF I/O Strategy | Carlos Osuna (C2SM, ETH) | 4.25 | |
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 IssuesCoding 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. |
|||||||
0.13 | 2011-09-10 | finish | 2012-12-01 | I/O Optimization with MPI_ALL2ALL | Fuhrer (MCH) | 4.25 | |
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 IssuesCoding 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). |
|||||||
TestingFunctionality has been tested successfully. |
|||||||
DocumentationUser Guide has been updated accordingly. |
|||||||
0.14 | 2011-09-10 | finish | 2012-09-01 | Removal of (dead) code | Schaettler, et al | 4.23 | |
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 IssuesNot applicable; Just make sure that nobody is using that any more |
|||||||
TestingWorking configurations must not be touched. Up to now, nobody complaint. |
|||||||
DocumentationThe User Guide has been updated accordingly. |
|||||||
0.16 | 2011-09-10 | finish | 2012-09-01 | Use SST forecast from an ocean model as boundary condition | Lucio Torrisi | 4.23 | |
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 IssuesCoding Standards: fulfilled Technical Test Suite: not yet available, but some standard tests have been conducted without problems 4-eyes Assurance: done by Uli Schaettler |
|||||||
TestingNo changes of results. Only some functionality tests are necessary. |
|||||||
DocumentationUser Guide has been updated. |
|||||||
0.04 | 2011-09-10 | finish | 2012-12-01 | Implementation of a tracer module | Roches / Fuhrer (MeteoSwiss) | 4.25 | |
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 IssuesCoding 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. |
|||||||
TestingOnly functionality tests are necessary, as the code gives bit-identical results. |
|||||||
DocumentationA COSMO Technical Report on the new Tracer Module, which serves as technical documentation, is in preparation. It is expected in December. |
|||||||
0.07 | 2011-09-10 | finish | 2013-11-30 | Unified COSMO-ICON Physics | PP POMPA | 5.05 | |
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 IssuesMCH 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. |
|||||||
TestingMicrophysics: 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. |
|||||||
DocumentationExtensions for the documentation of the microphysics in the "Physical Parameterizations Documentation" is available and will be incorporated soon. |
|||||||
0.10 | 2011-09-10 | finish | 2012-09-01 | Physics-Dynamics Coupling | Oliver Fuhrer (MCH) | 4.23 | |
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 IssuesCoding 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 |
|||||||
TestingAs there are no changes to the results, only the technical Test Suite was necessary. |
|||||||
DocumentationChange is documented in the Release Notes. No other change of documentation necessary. |
|||||||
0.05 | 2011-06-30 | finish | 2013-05-31 | Implementation of grib-api | Schaettler (DWD) | 4.28 | |
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 IssuesCoding 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 |
|||||||
Testing27.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 |
|||||||
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. |
|||||||
0.06 | 2011-06-30 | finish | 2012-09-01 | Reconstruction of mpe_io.f90 | Prill (DWD) | 4.25 | |
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 IssuesCoding 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. |
|||||||
TestingThe 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. |
|||||||
DocumentationA 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. |
|||||||
0.02 | 2011-04-16 | finish | 2011-08-31 | Digital Filter Initialization | Schaettler | 4.20 | |
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 IssuesThe changes are still evaluated and tested by colleagues from the Bundeswehr. |
|||||||
TestingSee Technical Issues. |
|||||||
Technical IssuesThe changes have been implemented by colleagues from ETH Zurich and have been cross-checked and implemented by the Source Code Administrator. |
|||||||
TestingFew single test cases have been performed to ensure that the results have not been changed. |
|||||||
WG7 | |||||||
5.1 | 2016-09-01 | finish | 2019-08-31 | Development of COSMO-E | A. Walser | N/A | |
The COSMO-E ensemble is now operational at MCH: - 2.2 km hor. res., 60 vertical levels - 20 members - ICs from KENDA - BCs from ECMWF ENS - SPPT as model perturbation |
|||||||
Status:
|
|||||||
5.2 | 2016-09-01 | finish | 2019-08-31 | Development of COSMO-2I-EPS | C. Marsigli | N/A | |
COSMO-IT-EPS is in its development phase: - 2.8 km hor. res. (soon 2.2), 50 levels - 10 members - ICs now from downscaling (later KENDA) - BCs either from COSMO-ME-EPS or COSMO-LEPS - model perturbations: SPPT and parameters (see SPRED PP) - perturbation of soil mositure ICs (see SPRED PP) |
|||||||
Developments (see SPRED PP): - study the combination of SPPT and parameter perturbation - impact of soil mositure perturbation - use of ICs from KENDA |
|||||||
5.4 | 2016-09-01 | finish | 2019-08-31 | Development of TLE-MVE | A. Mazur | N/A | |
DescriptionTLE-MVE is in development phase: - 2.8 km hor. res. - 20 members - ICs and BCs are provided by 4 successive ICON deterministic run (intermediate step with COSMO 7 km) - no model perturbation - soil perturbations are applied to each member |
|||||||
Statussoil perturbation (see SPRED PP) |
|||||||
5.5 | 2016-09-01 | finish | 2019-08-31 | Development of COSMO-DE-EPS | C. Gebhardt | N/A | |
DescriptionCOSMO-DE-EPS is operational since 2012: - 2.8 km hor. res. 50 levels - 20 members - ICs and BCs from BCEPS - perturbed parameters - soil perturbation |
|||||||
StatusDevelopment (see SPRED PP): - random parameters - more parameter perturbed (energy) - EM scheme for model perturbation - test of ICs from KENDA (also in combination with BCEPS) |
|||||||
4.1 | 2014-03-01 | finish | 2015-08-31 | Test of the verification of ensemble forecasts with VERSUS, with WG5 | Tesini, Walser, Bundel | N/A | |
Exchanges of expertise about ensemble verification between WG7 and WG5 is on-going. The work is now concentrating, for WG7, in providing to WG5 and especially to VERSUS PL and developers, specifications about the further implementation and refinement of EPS verification in VERSUS. |
|||||||
On-going. |
|||||||
5.3 | 2013-05-01 | finish | 2013-12-10 | Evaluate ECMWF LAM-EPS BC test data-sets for the convection-permitting scale | A. Walser & C. Marsigli | N/A | |
The EPS-BCs test data-sets described in WP 7.1.5 are tested for the convection-permitting ensembles developed at MCH and ARPA-SIMC. A feedback is provided to ECMWF. |
|||||||
Completed the first test. Now is on-going a test of more frequent LBCs from ENS. |
|||||||
2.1 | N/A | finish | 2014-01-31 | Development of the 7km ensemble system for the Sochi olympic games | Andrea Montani | N/A | |
DescriptionA 7-km ensemble similar to COSMO-LEPS is implemented on the Sochi area (COSMO-RU-LEPS). This work is part of the CORSO Priority Project. A 10-member downscaling of the 12 UTC EPS is performed with the COSMO model, up to 72h forecast range. Products specifically designed for RHM have been prepared and disseminated. |
|||||||
Verification has been performed and presented at the last COSMO GM. |
|||||||
Technical IssuesFor some products, coding in GRIB2 require more work, since grib-api is not ready for this. |
|||||||
TestingTesting will be performed for Winter 2011-2012 and 2012-2013. |
|||||||
1.5 | N/A | finish | 2012-08-31 | Migration to grib-api of the COSMO-LEPS clustering code | Andrea Montani | N/A | |
DescriptionIn the clustering code EMOSLIB is still used. The code should be migrated to grib-api. Technical change. |
|||||||
Technical IssuesUse of grib-api, F90 code. |
|||||||
TestingCheck that results do not change after and before the migration. |
|||||||
KENDA | |||||||
N/A | 2018-10-31 | finish | 2019-02-27 | ensure reproducibility of redundancy check, adaptations for buoy and AMV reports | Schraff (DWD) | 5.06 | |
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. |
|||||||
DocumentationProvide specific information on availability of the following documentation required by the COSMO source code management rules:
|
|||||||
N/A | 2018-04-11 | work | 2019-08-31 | new data types: high-res (BUFR) radiosonde, radiosonde descent data, tower profile data; option for superobbing of vertical profiles; minor (technical) changes and bug fixes | Schraff (DWD) | 5.07 | |
Description
input file name , description of code type
The criterion for which profiles superobbing is applied depends on the new namelist variable 'av_reso'. The target superobbing layers depend on the new The redundancy check is adjusted to avoid supplementing a superobbed profile by data from the original report (with state 'merged').
Changes of the namelist variables:
Changes of the results:
|
|||||||
StatusAll modifications implemented, currently in private version of C. Schraff. |
|||||||
Technical IssuesCoding Standards: Yes. Technical Test Suite, 4-eyes Assurance |
|||||||
TestingAll changes tested in short experimental trials (technical tests). An experiment over several weeks is currently ongoing, with neutral results expected (without changing the selection of actively used observation types). |
|||||||
N/A | 2016-06-01 | finish | 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 | Schraff (DWD) | 5.05 | |
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. |
|||||||
DocumentationProvide specific information on availability of the following documentation required by the COSMO source code management rules:
|
|||||||
N/A | 2016-06-01 | finish | 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 | Schraff (DWD) | 5.05 | |
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. |
|||||||
N/A | 2013-01-30 | finish | 2013-11-08 | for OSSE capability: read obs from feedback files | Schraff (DWD) | 4.28 | |
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 |
|||||||
N/A | 2012-12-01 | finish | 2014-03-31 | read Meteosat radiances | Faulwetter (DWD) | 5.06 | |
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 . |
|||||||
N/A | 2012-12-01 | finish | 2014-06-30 | read / process GPS slant path delay data | Bender (DWD) | 5.05 | |
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. |
|||||||
DocumentationAvailability of the following documentation required by the COSMO source code management rules:
|
|||||||
7 | 2012-12-01 | work | 2014-03-31 | pattern generator for additive covariance inflation and stochastic physics | Tsyrulnikov (HMC) | N/A | |
Description
|
|||||||
4 | 2012-09-10 | finish | 2013-11-08 | modular observation operators for existing observation types | Schraff (DWD) | 4.28 | |
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 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 A new test is being started |
|||||||
DocumentationIModel Documentation : improved in-line documentation External Documentation : no changes required for scientific documentation, User's Guide |
|||||||
6 | 2012-05-31 | finish | 2012-12-13 | RTTOV obs operator + read / write satellite radiances | Messer (DWD) | 4.26 | |
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 IssuesCoding 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. |
|||||||
TestingThe 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. |
|||||||
DocumentationOffline documentation is available from DWD, FE12. The implementation is documented in the Release Notes of 4.26. The User Guide has been updated accordingly. |
|||||||
2 | 2012-03-01 | finish | 2012-12-01 | adaptions for sub-hourly analysis update with LETKF | Reich (DWD) | 4.24 | |
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 IssuesCoding 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
|
|||||||
DocumentationNew 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. |
|||||||
3 | 2011-09-10 | finish | 2014-04-16 | stochastic physics | Torrisi (CNMCA) | 5.01 | |
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 |
|||||||
StatusPrevious 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
|
|||||||
TestingTesting 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. |
|||||||
DocumentationInline 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. |
|||||||
N/A | 2011-03-31 | work | 2014-06-30 | read / process cloud analysis for KENDA | Schomburg (DWD) | N/A | |
Description
|
|||||||
5 | 2011-03-31 | finish | 2014-03-31 | radar obs operators + nudging of radial winds | Blahak / Stephan (DWD) | 5.06 | |
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 IssuesCoding 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 |
|||||||
TestingThe 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. |
|||||||
DocumentationScientific 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.
|
|||||||
N/A | 2011-03-31 | finish | 2014-04-16 | interface for radar obs operators | Blahak (DWD) | 5.01 | |
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. |
|||||||
1 | 2010-09-30 | finish | 2012-03-31 | new modular nudging code (with many new options) | Schraff (DWD) | 4.22 | |
Description
|
|||||||
Technical Issues
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
task completed. |
|||||||
Documentation
task completed |
|||||||
External people | |||||||
1 | 2012-02-01 | finish | 2012-09-01 | CLM Contributions to Unified Version 5.0 | Hans-Jurgen Panitz | 4.23 | |
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 IssuesCoding Standards: fulfilled Technical Test Suite: not yet available, but changes could be tested by all COSMO Partners 4-eyes Assurance: done by Uli Schaettler |
|||||||
TestingAll changes have been tested by the CLM Community and in addition by the COSMO partners. No complaints were raised. |
|||||||
DocumentationThe COSMO User Guide has been updated accordingly. The changes have been described in detail in the Release Notes of Version 4.23 |
|||||||
2 | 2010-03-31 | finish | 2013-04-30 | Modular Earth Submodel System (MESSy) | N/A | 4.27 | |
4 | N/A | idea | N/A | Fusion of COSMO-ART and M7 | N/A | N/A | |
3 | N/A | idea | 2012-02-28 | bufr2netcdf | Cesari and Patruno (ARPA-SIMC) | N/A | |
DescriptionLibrary and applications to use bufr/crex reports in COSMO data assimilation. The library provides featureful BUFR and CREX encoding and decoding. Bufr2netcdf convert bufr/crex messages to netcdf file format suitable for cosmo data input. Notice that the bufr2netcdf program is able to convert almost any BUFR message to Netcdf format, but the COSMO model is able to import only Netcdf files created from a BUFR which was coded according to the proper template for each kind of observation (e.g. WMO for SYNOP and TEMP). For template conversion a downlodable external software is needed (DB-all.e). Components are:
|
|||||||
Technical IssuesNetcdf observation data input format fault documentation so the software and test are based on reverse engenering of the existing software and files. External sofware should be provided for bufr template conversion: DB-all.e, high-level tools to work with point-based weather data, based on its physical interpretation. It can do conversion between report templates, data manipulation from C++, Fortran and Python using a database, which can also be explored graphically. The software is provided under the GNU-GPL license:
ARPA-SIMC holds the copyright and, if needed, may agree to relax the last requirement for COSMO partners (i.e. license it under LGPL). |
|||||||
TestingPackage include automatic tests. Bufr2netcdf has been tested on bufr data provided by DWD. It produces the same netcdf output as DWD Bufrx2netcdf for: Further testing is required to ensure that these differences are negligible. The proposed test to be done is: 1) run of COSMO with data assimilation (nudging) with operational-like bufr converted to netcdf by bufrX2netcdf |
Start date | Task status | Expected delivery |
Task subject |
Responsible person |
Target version |
---|---|---|---|---|---|
WG3a | |||||
2016-01-01 | finish | 2017-06-01 | Kinne (2015) climatology (AOT12,SSA12,ASY12 --> TAU_KINNE,SSA_KINNE,G_KINNE) | Marina Shatunova (RHM) | 2.05 |
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. |
|||||
WG3b | |||||
2013-12-01 | work | 2014-05-30 | Interpolation of soil moisture in the vicinity of glaciers | Guy DeMorsier | |
DescriptionWhen interpolating coarse grid COSMO to fine grid and in the case where the SOILTYP ice (1, therefore no soil moisture) of the coarse grid covers larger areas as the fine grid, the interpolation of the soil moisture variable W_SO is modified. In this case, the original coarse grid COSMO soil moisture comes from the nearest grid point which is not an ice point and the interpolated fine COSMO grid points have a meaningful moisture different from 0. This procedure could be used in a similar way for the case of rock (SOILTYP 2 which neither has any soil moisture). |
|||||
StatusWork on this issue has been given back to the developers. There are concerns, because the soil temperature is not modified accordingly, which in special cases could lead to cold pools in the COSMO-Model. |
|||||
WG6 | |||||
2017-01-16 | finish | 2017-01-31 | Remove clipping of T_S and T_SO to -0.8 deg C below (partially) snow covered grid points | Ulrich Blahak | 2.03 |
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. |
|||||
TestingDoes what it should do. Removes unrealistic cold soil temperature patches in partial snow covered areas. |
|||||
DocumentationDocumentation:
|
|||||
2016-12-05 | finish | 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 | Ulrich Blahak (DWD) | 2.03 |
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. |
|||||
DocumentationSee slides of U. Blahak (page 30 ff) |
|||||
2016-11-01 | finish | 2017-01-31 | New methods for profile adaption from coarse to fine orography in case of non-hydrostatic input models | Ulrich Blahak | 2.03 |
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. |
|||||
TestingExperiments at DWD showed that the initial noise in COSMO-DE is reduced over mountaineous regions. |
|||||
Documentation
|
|||||
2016-08-01 | finish | 2017-01-31 | Hydrostatic pressure computation instead of PP interpolation for non-hydrostatic input models | Ulrich Blahak | 2.03 |
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
|
|||||
2016-01-01 | finish | 2016-08-22 | Linear vertical interpolation as an option to splines | Ulrich Blahak (DWD) | 2.03 |
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. |
|||||
2016-01-01 | finish | 2016-08-22 | For GFS, possibility to use RELHUM instead of QV | Ulrich Blahak (DWD) | 2.03 |
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! |
|||||
TestingTested by Helmut Frank (DWD) and colleagues from RosHydromet. |
|||||
2016-01-01 | finish | 2016-08-22 | General 4D field to hold new 3D fields for easier integration into the code. | Ulrich Blahak (DWD) | 2.03 |
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. |
|||||
TestingIt was successfully used to implement the 11 CAMS aerosol mixing ratios from ECMWF in the framework of PP T2RC2. |
|||||
2016-01-01 | work | 2016-08-22 | MACC (CAMS) aerosol data from ECMWF | Harel Muskatel (IMS) | 2.03 |
DescriptionImplemented interpolation of MACC aerosol data to COSMO grid with the help of the new 4D field. 11 new fields. |
|||||
2016-01-01 | finish | 2016-08-22 | TERRA_URB: new EXTPAR fields AHF and FR_PAVED/ISA | Hendrik Wouters (KU Leuven) | 2.03 |
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. |
|||||
2016-01-01 | finish | 2016-08-22 | Field SSO_STDH needed for the new blocked COSMO turbulence code | Ulrich Blahak (DWD) | 2.03 |
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). |
|||||
TestingWorks correctly |
|||||
2016-01-01 | work | 2100-12-01 | For ICON input: Interpolation of convective precip rates for moisture enhancement in COSMO | Ulrich Blahak (DWD) | 2.03 |
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. |
|||||
2016-01-01 | finish | 2016-08-22 | ICON-nest, ICON-lam or ICON-sub input: domain checker | Ulrich Blahak (DWD) | 2.03 |
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 |
|||||
Technical IssuesNone. 4-eyes pending. |
|||||
TestingWorks as expected |
|||||
2014-01-15 | finish | 2014-05-30 | SLEVE2 coordinate | Guy DeMorsier | 2.01 |
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. Schär, 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 Schättler for reviewing. |
|||||
Technical IssuesThe code has been given to Uli Schättler (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-01-01 | finish | 2014-09-15 | Implementation of ICON to COSMO interpolation | Thorsten Reinhardt | 2.01 |
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" |
|||||
2013-12-01 | finish | 2014-05-30 | Consolidation of GRIB2 and usage of grib_api | Uli Schättler | 2.01 |
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. |
|||||
2013-06-28 | finish | 2013-08-19 | Implementation of grib_api for reading / writing | Uli Schättler | 1.22 |
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-06-28 | finish | 2013-08-19 | Multi-layer treatment of IFS soil | Davide Cesari | 1.22 |
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 |
|||||
finish | 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 | Ulrich Blahak (DWD) | 2.04 |