There have been several changes for this new version, but most of them are pure technical and do not change the results.
The interfaces to use COSMO-ART have been implemented now also in the Runge-Kutta dynamical core. The modules affected are
y_scalar_advect | Description | lsl_adv_qx | yef_adv_qx |
---|---|---|---|
'SL3_MF' | SL-Advection with tricubic interpolation, multiplicative filling-option | .TRUE. | --- |
'SL3_SFD' | SL-Advection with tricubic interpolation, selective filling diffusion-option | not available before | |
'Bott2' | Bott 2nd order finite-volume scheme | .FALSE. | 'Bott_2' |
'Bott2_Strang' | Bott 2nd order finite-volume scheme, with Strang splitting ('z-y-2x-y-z') | not available before | |
'Bott4' | Bott 4th order finite-volume scheme | .FALSE. | 'Bott_4' |
'Bott4_Strang' | Bott 4th order finite-volume scheme, with Strang splitting ('z-y-2x-y-z') | not available before | |
'vanLeer' | .FALSE. | 'vanLeer' | |
'vanLeer_Strang' | not available before | ||
'PPM' | .FALSE. | 'PPM' | |
'PPM_Strang' | not available before | ||
'MPDATA' | not yet available, work at MeteoCH | not available before |
In addition to the module src_sat_tbs.f90 (which is the interface to the RTTOV7 library) a new module src_sat_rttov.f90 has been introduced as interface to the RTTOV9 library. The synthetic satellite images can now also be computed using the RTTOV9 library.
In contrast to the RTTOV7 library, multiple-scattering effects are not computed explicitly, but are parametrised. Also, microphysical cloud properties are now taken into account better than before. For areas with a complex cloud coverage in the vertical there can be differences in the brightness temperatures up to some Kelvin. With low / high brightness temperatures (high / low clouds) the values will be lower / higher. The contrasts between lower and higher brightness temperatures will be intensified.
To choose between the libraries, a new Namelist parameter itype_rttov has been introduced in the group /SATCTL/:
Name | Meaning | Default | ||||
---|---|---|---|---|---|---|
itype_rttov | to specify the RTTOV version that should be used
|
7 |
Conditional compilation has been introduced for Nudging, the synthetic satellite images and the use of the NetCDF library. Not all users of the COSMO-Model do use these components, but have to compile and link them together with the additional libraries needed. Up to now we were offering "dummy"-files, but now we introduced conditional compilation (using the #ifdef pragma) for switching on / off these components during compilation time:
Note: | NetCDF can only be switched off, if also Nudging is switched off, because the treatment of observations in the Nudging is done with NetCDF and cannot be switched off then. |
If the source code is NOT compiled with these extensions, but the corresponding component is initiated during a model run, the model will stop with an appropriate error message.
Some technical adaptations were necessary to run the forward launching filter (ndfi=1), because the forecast then starts for nstart > 0, but some initializations (e.g. in the soil model) were done only for nstart = 0.
For the forward launching (ndfi=1), the soil processes are now taken into account. The initial values of all soil variables are saved and after the filtering step, all soil variables are initialized again as mean values of the initially saved values and the newly computed ones.
Note: | The reproducibility problem for the backward-forward initialization (ndfi=2) still is not solved. |
There have been changes to few diagnostic variables:
This tuning parameter will be tested for COSMO-DE-EPS:
Shallow convection parameterization becomes active only if cloud
thickness from cloud base to cloud top exceeds a threshold. To evaluate
this condition a parcel is launched. This threshold is typically set to
values between 20000.0 Pa and 30000.0 Pa with a COSMO_DE default of 25000.0 Pa.
(for comparison: IFS default is 20000.0 Pa set by Peter Bechtold).
Reasonable values are between 10000.0 and 45000.0 Pa.
In case of luvmasspoint=.TRUE. there were stripes from the processor grid in the output of AUMFL_S and AVMFL_S because of the averaging to the mass grid point. The boundaries of these variables must be exchanged before computing the averaging.
Call to subroutine potential_vorticity_rho: variable zvarlev is now passed within (1:ke) to this subroutine. (by Jean-Marie Bettems)
The computation of qc_rad and qi_rad was wrong, when using a coarser grid for the radiation. This has been fixed. (by Victor Venema, Uni Bonn)
There were the following changes for the Namelist variables:
Group | Name | Meaning | Default | |||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
/DYNCTL/ | lsl_adv_qx | eliminated | - | |||||||||||||||||||||
yef_adv_qx | eliminated | - | ||||||||||||||||||||||
y_scalar_advect | string to specify the horizontal advection scheme for scalar variables
in the Runge-Kutta dynamics. The string is case-insensitive, lower and upper
case letters are accepted.
|
'SL3_MF' | ||||||||||||||||||||||
/TUNING/ | thick_sc | limit for convective clouds to be "shallow" (in Pa) Range: 10000.0 ... 45000.0 |
25000.0 | |||||||||||||||||||||
/IOCTL/ | yform_write | eliminated | - | |||||||||||||||||||||
/GRIBOUT/ | yform_write | to specify output format for every output group separately. | 'grb1' | |||||||||||||||||||||
/SATCTL/ | itype_rttov | to specify the RTTOV version that should be used
|
7 | |||||||||||||||||||||
the following namelist variables are new and are only in effect for using RTTOV9 | ||||||||||||||||||||||||
sat_long_01 | position of first satellite (longitude) | -999.0 | ||||||||||||||||||||||
sat_long_02 | position of second satellite (longitude) | -999.0 | ||||||||||||||||||||||
extrp_type | Type of extrapolation above highest model level
|
0 | ||||||||||||||||||||||
iceshape | To specify, whether ice particles are
|
1 | ||||||||||||||||||||||
iwc2effdiam | Type of conversion of ice water content to
effective diameter of ice particles
|
4 |
Back to Contents
The modifications are purely technical, so that the results do not change.
Only few diagnostics are affected because of bug fixes / technical changes:
11. Changes of Results
VMAX_10M
because of computations on the boundaries;
more changes, if itype_diag_gusts=2 is used, because of a bug fix.
HZEROCL
changed search algorithm
SNOWLMT
changed search algorithm
ASWDIFU_S
introduced smoothing of this field after
computation of radiation on a coarse grid
SYNMSG
if RTTOV9 is used, the results will change
z-levels
fields on z-levels do change, because the lower boundary
condition for the vertical interpolation has changed.