Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Section


Column

In this article we explain how to prepare and configure OpenIFS 43r3 for a nudged simulation. Therein the model needs to read meteorological parameters at the grid scale from pre-computed external forcing files. These forcing files have to be created prior to the nudged OpenIFS model experiment and this process is also described here. 

Info

Please note that nudging in OpenIFS is an experimental research tool and therefore may change between model versions. The description on this page refers to nudging in gridpoint space which also allows regional nudging. This differs from nudging in spectral space which constrains the model in the longer wave numbers on a global scale. 

For further assistance on nudging and configuring OpenIFS please post your question in the OpenIFS User Forums or alternatively email openifs-support@ecmwf.int

Newtonian Relaxation

OpenIFS uses initial and boundary conditions to calculate its own model dynamics, i.e. meteorological variables that are resolved on the grid scale. It is however possible to constrain the model dynamics with external data. Newtonian relaxation, sometimes referred to as "nudging", is a simple form of data assimilation which allows the user to constrain or "force" the model's meteorological fields with reanalysis data. This is sometimes referred to as running the model in "offline" mode. In nudged configuration the model's dynamics is continually nudged towards the meteorological reanalysis independent of the run length of the experiment. 

This method relaxes the model state towards gridded re-analysis data (or towards output from another atmospheric model, or gridded observational data) by adding a non-physical relaxation term to the model equations (Jeuken et al., 1996). In the equation below X represents a prognostic model variable and Fmodel the model forcing which determines the evolution of X. The relaxation term G (Xobs - X) includes the relaxation coefficient G (in s-1) which determines the "tightness" of the nudging.

Mathdisplay
\frac{\partial X}{\partial t} = F_{model}(X) + G(X_{obs} - X)

This method can be useful, for instance, in sensitivity studies which aim to isolate the model physics or chemistry while preventing feedbacks to the model dynamics. Another example for its use is to align a climate model simulation closer to historic meteorology for comparison with measurements. 


Column
width20%


Panel
bgColorwhite
titleBGColor#f0f0ff
titleOn this page...

Table of Contents
indent15px



...

Model Code Changes

Section

The source code changes to run OpenIFS 43r3 in nudged configuration are under version control in the OpenIFS Git repository in branch feature/nudging.In OpenIFS 43r3v1 the source code required for nudging is included by default. 

Column
width215px

In order to control the model's nudging configuration the following source code files are relevant:

  • namrlx.nam.h  –  the namelist block which appears in the fort.4 file. These variables allow to customise the model's nudging behaviour without recompiling the model source code.
  • yomrlx.F90  – in this module the  variables in the namelist block are defined.
  • surlx.F90  – this is the setup routine for the namelist variables, which includes their initialisation with default values and consistency checks.
  • updrlxref.F90  – this subroutine updates the reference fields for the relaxation.

The code in the nudging branch of OpenIFS 43r3 43r3v1 differs from that in IFS insofar that the variables CLPSHRLX and CLPGGRLX have been added to the namelist. These two variables allow to specify the location of the directory with the forcing files on the file system.  If these are not set then the model will expect to find the forcing files in the experiment directory.

The relaxation to the analyses involves code in further routines (e.g. yomsrlx.F90, relaxgp.F90, stepo.F90, stepo_oops.F90, su0yomb.F90, and updtim.F90) however this part does not need to be changed if the existing nudging capability of the model is to be used.

...