!-----------------------------------!-----------------------------------
 Two simulations of this ensemble were run with double radiation calls : 
!-----------------------------------!-----------------------------------
   CM619-LR-piClim-histall-2 : for CO2 only 
       (see PARAM/config.def_annuel : co2_ppm_per = 2.8432e+02 =co2_preind)

    CM619-LR-piClim-histall-3 : for well-mixed GHG and aerosols 
       (see PARAM/config.def_annuel : all GHG "_per" values = preind)

!-----------------------------------!-----------------------------------
  A t t e n t i o n   i s   d r a w n  to the following files, 
 modified to allow easy (de)activation of IRF(-related) variable output :
!-----------------------------------!-----------------------------------

 A/ Customized *.xml files (stored in PARAM) 
!-------------------------

 ===> context_lmdz.xml_2RadCall : 
  Includes new grid "presnivs0srf" to extract only the surface (level 1) of 3D variables.
  Includes new output file "histmthForIRFtropo" with 3D radiation vars for GHG, and *_tropopause variables

 ===> field_def_lmdz.xml_2RadCall : 
  Groups of variables are created for :

     aerosols (TOA and SRF), with standard or CMIP-style variable names : 
        group_ref="fields_irf_aer_2D_lmdz"  
                or "fields_irf_aer_2D_CMIPstyle"

     GHG-2D (TOA and SRF), for 4xCO2, or another GHG combination (typically preindustrial) :
        group_ref="fields_rad_4co2_TOA"  and "fields_rad_4co2_SRF",
                or "fields_rad_ghg_TOA"  and "fields_rad_ghg_SRF"

     GHG-3D (for IRFs at tropopause level);
        group_ref="fields_rad_4co2_3D"
                or "fields_rad_ghg_3D"

Output of these variables can be activated/disabled in file*xml as a group ("enabled= TRUE/FALSE" for the group), instead of changing the "level" for each variable individually. 

 ===> file_histmth_lmdz.xml_2RadCall : 
   Allows (de)activating output for 2D groups "TOA" and "SRF"

 ===> file_def_histmth_lmdz.xml_2RadCall_ForIRFtropo : 
  Contains GHG-3D groups (4co2 and ghg) and the *_tropopause variables.
  The corresponding output file "histmthForIRFtropo" can be erased after producing TS_MO for all these variables in the post-processing.


 B/ COMP/lmdz.card 
!-----------------

1/ Points to the modified *xml files :
[ParametersFiles]
...
        (${SUBMIT_DIR}/PARAM/context_lmdz.xml_2RadCall, context_lmdz.xml ),            \
        (${SUBMIT_DIR}/PARAM/field_def_lmdz.xml_2RadCall, field_def_lmdz.xml ),          \
        (${SUBMIT_DIR}/PARAM/file_def_histmth_lmdz.xml_2RadCall, file_def_histmth_lmdz.xml ),   \
        (${SUBMIT_DIR}/PARAM/file_def_histmth_lmdz.xml_2RadCall_ForIRFtropo, file_def_histmthForIRFtropo_lmdz.xml ),   \
...

2/ Includes "histmthForIRFtropo" in the Output file list :
[OutputFiles]
List=….
   (histmthForIRFtropo.nc,  ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthForIRFtropo.nc,  Post_1M_histmthForIRFtropo),       \
   ….

3/ 2D and 3D timeseria are requested in postprocessing for all these variables needed to compute IRF at TOA, SRF and tropopause :
see sections [Post_1M_histmth] and the newly introduced [Post_1M_histmthForIRFtropo]