Voir la source de HowTo: create mass flux files

=  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 <code>offline</code>


Note: 
* we are working with LMDZ after [https://trac.lmd.jussieu.fr/LMDZ/changeset/4608/ 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 <code>offline</code> to yes 
*  add the reference to  <code><file_definition src="./file_def_phystoke_lmdz.xml"/></code> in <code>context_lmdz.xml</code> 
*  choose the output storage frequency with parameter <code>offline_time</code>. For example, for a storage all the 3hours (1day/8) we will choose <code>offline_time=8</code>
*  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''' 

== Physics mass flow Variables descriptions (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) 

== Dynamics mass flow Variables descriptions (wg.nc, tetag.nc, phig.nc, pbarvg.nc, pbarug.nc, masseg.nc)  == 
* "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"


[[Category:HowTo]]