Version du 5 juillet 2023 à 13:10

How to : create mass flow files for inversion models

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

Note:

• we are working with LMDZ after rev 4608
• Actually, this option is not working with XIOS in the dynamic (we are using writefield to write dynamic output files)

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_H}/${PREFIX}_3H_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'