Handout (725.2 kB) Handout (2.3 MB)
However, the observations (reflectivity, radial velocity, polarisation parameters) are not directly comparable to the prognostic variables of the model. In order to, on one hand, enable radar data assimilation in the framework of the above-mentioned LETKF-assimilation system and, on the other hand, to facilitate comparisons of numerical simulations with radar observations in the context of cloud microphysics verification, a comprehensive modular radar forward operator is under development. This operator calculates the radar observables reflectivity, radial wind and (some) polarisation parameters from the prognostic model output. Given the planned range of applications, it has to be applicable on supercomputer architectures in an operational environment and efficiency is a major design criterion, which requires for good parallelization and vectorization properties of the code.
The operator consists of several modules, each of which handles a special physical process (scattering, extinction, microwave propagation, etc.). Each of these modules offers different formulations to choose from, which enables variable possibilities to configure the forward operator with respect to the often conflicting constraints of physical accuracy and computational costs. Here, a main goal is to find in some sense “optimal” configurations for both applications (data assimilation and model verification). For example: 1) the radar beam can be considered to propagate as a simple ray or with the actual volume averaging characteristics (beam function). 2) The beam bending can be either derived from a 4/3 earth radius concept or from the actual simulated vertical gradient of the refractive index of air. 3) Radar reflectivity may be calculated from the full Mie-theory or from various (more efficient) approximations, and attenuation effects may be taken into account or not. The strategy is to start from the "full detailed" operator and, by subsequent neglect, simplification and/or approximation of single physical effects, the "best" balance between computational efficiency and physical accuracy for each of the two applications is investigated. Later on, impact studies on radar data assimilation and studies on model microphysics verification will round out the project.
Currently the project is in the stage of finishing the development of the "full" operator. The presentation will show first results.
![](images/flag.png)
![](images/blue_ribbon.png)