Version du 21 juillet 2023 à 10:47

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 in phystoke.nc

• Name="t_stok" long_name="Temperature" unit="K"
• "coefh_stok" long_name="coeff de melange dans la couche limite" unit="-"
• "fm_th_stok" long_name="Mass flux thermals" unit="kg m-2.s-1"
• "en_th_stok" long_name="Thermal plume entrainment" unit="kg m-1.s-1"
• "da_stok" long_name="adiabatic ascen" unit="kg m-2.s-1"
• "mp_stok" long_name="precipitating downdraft mass flux under hydrostatic approximation" unit="kg m-2.s-1"
• "wght_stok" long_name="weights of the layers feeding convection" unit="-"
• "upwd_stok" long_name="saturated updraft mass flux" unit="kg m-2.s-1"
• "dnwd_stok" long_name="saturated downdraft mass flux" unit="kg m-2.s-1"
• "phi_lev"[0-79] "mass flow environmental mixing" unit="kg m-2 s-1"
• "yu1_stok" long_name="zonal surface wind" unit="m.s-1"
• "yv1_stok" long_name="wind meridian surface" unit="m.s-1"

All of these variables are an average of offline_time frequency, on the horizontal grid with vertical level klev (and not presnivs)

Variables descriptions in dynamic

• "phig" long_name="Surface geopotential height" unit="m2 s-2"
• "massg" long_name="Air mass (averaged over 3 hours)" unit="kg"
• "tetag" long_name="Potential temperature (averaged over 3 hours)" unit="K"
• "pbarug" long_name="Zonal advected mass fluxes (averaged over 3 hours)" unit="kg m-2 s-1"
• "pbarvg" long_name="Meridional advected mass fluxes (averaged over 3 hours)" unit="kg m-2 s-1"
• "wg" long_name="Vertical advected mass flux (averaged over 3 hours)" unit="kg m-2 s-1"