Version du 27 avril 2023 à 14:05

How to : create mass flux files for inversion models

The lmdz model allows to output mass flux files that can be used as inversion model forcing files. This option is drive by the parameter offline

Note : Actually, this option is not working with XIOS in the dynamic

option offline in LMDZ

To use this option you need :

• to put parameter offline to yes
• add the reference to <file_definition src="./file_def_phystoke_lmdz.xml"/> in context_lmdz.xml
• choose the output storage frequency with parameter offline_time. For example, for a storage all the 3hours (1day/8) we will choose offline_time=8
• if you are running with modipsl, modify lmdz.card to save mass flux files :

 [OutputText]
 List=   (.... , wg.nc, tetag.nc, phig.nc, pbarvg.nc, pbarug.nc, masseg.nc)
 [OutputFiles]
 List=   (.... ),                  \
         (phystoke.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_phystoke.nc,      NONE)

• if you are not running with modipsl, save the output files : wg.nc, tetag.nc, phig.nc, pbarvg.nc, pbarug.nc, masseg.nc and phystoke.nc

Variables descriptions

       phiss.standard_name='surface geopotentiel'
       levs.standard_name='Model Level Number'
       upwds.standard_name='KE sat. up. mass flux'
       dnwds.standard_name='KE sat. down. mass flux'
       fmths.standard_name='Mass flux thermals'
       enths.standard_name='Ent therm.'
       fmths.units='kg m2 s-1'
       enths.units='kg m2 s-1'
       upwds.units='kg m2 s-1'
       dnwds.units='kg m2 s-1'
       ts.standard_name='t'
       das.standard_name='adiabatic ascent'

       phis.standard_name='geopotentiel instantane'
       phiss.standard_name='Surface Geopotential Height'
       aires.standard_name='Grid area'
       masses.standard_name='Mass'
       levs.standard_name='Model Level Number'
       phis.standard_name='geopotentiel instantane'
       pbarus.standard_name='flx de masse zonal'
       ws.standard_name='flx de masse vert'
       tetas.standard_name='temperature potentielle'
       pbarus.units = 'kg m2/s'
       ws.units = 'kg m2/s'
       tetas.units = '-'
       phis.units = '-'
       masses.units='kg'