...
Mass conservation
OpenIFS (with thanks to Michail Diamantakis)IFS) does not formally conserve mass. The change in mass is small but needs correcting if using OpenIFS for long integrations, or with chemistry/aerosol modelling.
The version of OpenIFS based on Cycle 38r1 has a limited set of "mass-fixers" available. More recent versions of IFS have improved options for mass conservation.
Types of mass fixers
In OpenIFS 38r1 there are only global mass fixers available. The first ensures conservation of the global model mass by using the spectral characteristics of the model. The second method (enabled by some code changes detailed below) is also global and works in gridpoint space only on tracer quantities.
Global mass fixer
For global mass conservation of air the following settings should be added in namelistsor enabled in the namelist file (fort.4):
Code Block |
---|
&NAMDYN LMASCOR=true, LMASDRY=true, &NAMCT0 NFRMASSCON=1, |
LMASCOR enables the global mass conservation by ensuring the global mass before and after a timestep remains the same.
LMASDRY ensures the mass conservation applies only to the dry air mass. This is done to ensure that (the above fixes the dry air mass only and leaves the fluctuation in the total mass due to moist parametrizations - this preserves the annual water cycle signal which would otherwise be removed by fixing the total mass. ERA-Interim also use these settings.)
The above turn on the spectral mass fixer.
Code Block |
---|
&NAMCT0
NFRMASSCON=1, |
...
is unaffected by the mass correction. If this option is set false
the mass correction is done to the total air mass (dry+water vapour). Note that the ERA-Interim reanalyses used LMASDRY=true.
NFRMASSCON is the frequency at which the mass fix is applied. The above example enables it at every timestep. It must be set as the default is set to a very large number (i.e. it's default is to effectively disable it).!
Mass correction for moisture and tracer variables
...
A further global mass fixer can also be applied to correct the moisture variables and any tracers (if used). Note that this fixer can be used together with the spectral mass fixer above or independently.
To enable it, ensure the following namelist variables are set:
Code Block |
---|
&NAMGFL YS_NL%LMASSFIX=true, YR_NL%LMASSFIX=true, YI_NL%LMASSFIX=true, YL_NL%LMASSFIX=true, YQ_NL%LMASSFIX=true, LTRCMFIX=true, |
the above activates and controls the tracer mass fixer. It applies it on spec The model applies the fix to the specific humidity, cloud water, cloud ice, cloud rain, cloud snow: YQ_NL, YL_NL, YI_NL, YR_NL, YS_NL IFS variables respectively.
...
Code Block |
---|
112 !LTRCMFP=.TRUE. ! proportional 113 !LTRCMFA_LAP=.FALSE. ! additive based on hor. Laplace 114 !LTRCMFA_VER=.FALSE. ! additive based on vert. Lapace 115 !LTRCMFA_DIF=.FALSE. ! additiva based on abs gradient |
These enable mass fixers that work in gridpoint space on the above variables.
...
Acknowledgement
Dr Michail Diamantakis, research department, ECMWF, provided much of the information on this page.
HTML |
---|
<script type="text/javascript" src="https://software.ecmwf.int/issues/s/en_UKet2vtj/787/12/1.2.5/_/download/batch/com.atlassian.jira.collector.plugin.jira-issue-collector-plugin:issuecollector/com.atlassian.jira.collector.plugin.jira-issue-collector-plugin:issuecollector.js?collectorId=5fd84ec6"></script> |
...