Index: /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/dynamico.card
===================================================================
--- /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/dynamico.card	(revision 6895)
+++ /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/dynamico.card	(revision 6895)
@@ -0,0 +1,55 @@
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+RESOL_NBP=40
+RESOL_ATM_Z=79
+# create_etat0_limit=n is used for normal GCM run
+create_etat0_limit = n
+
+# The proposed experiment ensures consistency between nudging files and nudging_dynamico.xml.
+# If you change these nudging files, make sure you maintain this consistency. 
+guided_type=nudging
+
+# iflag_phys=1 is used for normal GCM run with continents (not aquaplanet)
+# iflag_phys is only used in combination with create_etat0_limit=y. See specific CREATE set up to know more.
+iflag_phys=1
+
+
+[InitialStateFiles]
+List=   (${ARCHIVE}/${config_UserChoices_CreateBasefolder}/${config_UserChoices_CreateName}/ICO/Output/Restart/${config_UserChoices_CreateName}_${year}_start.nc, start.nc)
+
+[BoundaryFiles]
+List=(${R_IN}/ATM/FORCING_LAM/${config_UserChoices_ForcingName}_${config_UserChoices_ForcingDate}.nc, nudging_forcing.nc)
+ListNonDel= ()
+
+[SmoothFiles]
+List=   ()
+
+[ParametersFiles]
+List=  	(${SUBMIT_DIR}/PARAM/DYNAMICO/earth_const.def, .),\
+   	(${SUBMIT_DIR}/PARAM/DYNAMICO/LAM/run_dynamico_lam_ERA5_France.def, run_dynamico.def),\
+	(${MODIPSL}/modeles/DYNAMICO/xml/DYNAMICO_XML/sponge_dynamico.xml, .),\
+	(${MODIPSL}/modeles/DYNAMICO/xml/DYNAMICO_XML/context_dynamico.xml, .),\
+   	(${MODIPSL}/modeles/DYNAMICO/xml/DYNAMICO_XML/context_input_dynamico.xml, .),\
+   	(${MODIPSL}/modeles/DYNAMICO/xml/DYNAMICO_XML/field_def_dynamico.xml, .),\
+	(${MODIPSL}/modeles/DYNAMICO/xml/DYNAMICO_XML/nudging_dynamico.xml, .),\
+   	(${MODIPSL}/modeles/DYNAMICO/xml/DYNAMICO_XML/file_def_dynamico.xml, .)
+
+
+[RestartFiles]
+List=   (restart.nc, restart.nc, start.nc)
+
+[OutputText]
+List=   (run_dynamico.def, earth_const.def, iodef.xml, context_dynamico.xml, context_input_dynamico.xml, field_def_dynamico.xml, file_def_dynamico.xml, out_icosa_lmdz.exe.out*, out_icosa_lmdz.exe.err*, xios_client_*)
+
+
+[OutputFiles]
+List=   (output_dcmip2016_native.nc,    ${R_OUT_ICO_O_H}/${PREFIX}_3H_output_dcmip2016_native.nc,     NONE), \
+        (output_dcmip2016_regular.nc,   ${R_OUT_ICO_O_H}/${PREFIX}_3H_output_dcmip2016_regular.nc,    NONE)
+
+
+#List=	( limit.nc,       ${R_OUT_ICO_O}/Boundary/${config_UserChoices_JobName}_clim_limit.nc,   NONE), \
+#	( start.nc,       ${R_OUT_ICO_O}/Restart/${config_UserChoices_JobName}_clim_start.nc,    NONE), \
+#	( startphy.nc,    ${R_OUT_ICO_O}/Restart/${config_UserChoices_JobName}_clim_startphy.nc, NONE)
+
Index: /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/lmdz.card
===================================================================
--- /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/lmdz.card	(revision 6895)
+++ /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/lmdz.card	(revision 6895)
@@ -0,0 +1,216 @@
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+# Physics package to use :
+# LMDZ_Physics=AP for standard/old physics(defalut), can be used with LMDZ4_AR5 or LMDZ5/trunk sources
+LMDZ_Physics=NPv6.2
+
+# Set ConfType to choose parameters for aerosols, solar and green house gazes.
+# The parameter file PARAM/config.def_$ConfType will be used.
+# ConfType=preind/actuel/annuel
+ConfType=annuel
+
+# Manage radiative aerosol effect and ozone
+flag_aerosol=0
+ok_ade=n
+ok_aie=n
+aerosol_couple=n
+read_climoz=0
+ok_cdnc=n
+flag_aerosol_strat=0
+
+# Set tau_gl depending on the forcing file for the sea-ice (use 2. for amip and 1. for hadisst)
+tau_gl=2.
+
+
+# Specify which output file to activate and at which output_level
+# Use an integer to specify the output level or NONE if to deactivate the file
+# output_level_filename=1     : activate the corresponding file with low output level
+# output_level_filename=10    : activate the corresponding file with high output level
+# output_level_filename=NONE  : the file will not be written
+
+output_level_histmth = 5
+output_level_histday = 3
+output_level_histhf  = NONE
+output_level_histmthNMC = NONE
+output_level_histdayNMC = NONE
+output_level_histhfNMC  = NONE
+output_level_histstrataer = NONE
+output_level_histstn      = NONE
+
+# Activate COSP in config.def and set corresponding output files
+LMDZ_COSP_OK=n
+output_level_histmthCOSP = NONE
+output_level_histdayCOSP = NONE
+output_level_histhfCOSP  = NONE
+
+
+[InitialStateFiles]
+# start.nc and startphy.nc are both set in dynamico.card to easier insure coherence of the files
+List=(${ARCHIVE}/${config_UserChoices_CreateBasefolder}/${config_UserChoices_CreateName}/ATM/Output/Restart/${config_UserChoices_CreateName}_${year}_startphy.nc, startphy.nc)
+
+[BoundaryFiles]
+#
+# With CMIP6 aerosols version 1
+# With 3D ozone
+List= (${ARCHIVE}/${config_UserChoices_CreateBasefolder}/${config_UserChoices_CreateName}/ATM/Output/Boundary/${config_UserChoices_CreateName}_${year}_limit.nc, limit.nc),\
+            (${R_IN}/ATM/AEROSOLS/CMIP6/v1/144x142/L${RESOL_ATM_Z}/aerosols${year}_from_inca.nc,   aerosols${year}.nc),\
+            (${R_IN}/ATM/SOLAR/CMIP6/v3.1corr/Historical/solarforcing_v3.1_daily_6bands_${year}.nc, solarforcing.nc),\
+            (${R_IN}/ATM/STRATAERO/CMIP6/v3/144x142/L${RESOL_ATM_Z}/taulwstrat.2D.${year}.nc, taulwstrat.2D.nc),\
+            (${R_IN}/ATM/STRATAERO/CMIP6/v3/144x142/L${RESOL_ATM_Z}/tauswstrat.2D.${year}.nc, tauswstrat.2D.nc)
+
+ListNonDel= (${R_IN}/ATM/GHG/CMIP6/historical/CO2_CMIP6_0000_2014.txt, CO2.txt),\
+            (${R_IN}/ATM/GHG/CMIP6/historical/CH4_CMIP6_0000_2014.txt, CH4.txt),\
+            (${R_IN}/ATM/GHG/CMIP6/historical/N2O_CMIP6_0000_2014.txt, N2O.txt),\
+            (${R_IN}/ATM/GHG/CMIP6/historical/CFC11eq_CMIP6_0000_2014.txt, CFC11.txt),\
+            (${R_IN}/ATM/GHG/CMIP6/historical/CFC12eq_CMIP6_0000_2014.txt, CFC12.txt),\
+            (${R_IN}/ATM/AEROSOLS/CMIP6/v1/144x142/L${RESOL_ATM_Z}/aerosols1850_from_inca.nc, aerosols.nat.nc)
+[SmoothFiles]
+List=   ()
+
+[ParametersFiles]
+List=   (${SUBMIT_DIR}/PARAM/LMDZ/physiq.def_${LMDZ_Physics}, physiq.def),       \
+        (${SUBMIT_DIR}/PARAM/LMDZ/vert.def_L${RESOL_ATM_Z},   vert.def),         \
+        (${MODIPSL}/modeles/LMDZ/DefLists/cosp_input_nl.txt, .),            \
+        (${MODIPSL}/modeles/LMDZ/DefLists/cosp_output_nl.txt, .),           \
+        (${SUBMIT_DIR}/PARAM/LMDZ/config.def_${ConfType},  config.def),          \
+        (${SUBMIT_DIR}/PARAM/LMDZ/guide.def, .),                                 \
+        (${SUBMIT_DIR}/PARAM/LMDZ/run.def, .),                                   \
+        (${SUBMIT_DIR}/PARAM/LMDZ/tracer.def, .),                               \
+        (${SUBMIT_DIR}/PARAM/LMDZ/LAM/context_lmdz.xml, context_lmdz.xml ),     \
+        (${MODIPSL}/modeles/LMDZ/DefLists/context_input_lmdz.xml, . ),      \
+        (${MODIPSL}/modeles/LMDZ/DefLists/field_def_lmdz.xml, . ),          \
+        (${MODIPSL}/modeles/LMDZ/DefLists/field_def_cosp1.xml, . ),          \
+	(${SUBMIT_DIR}/PARAM/LMDZ/file_def_histmth_lmdz.xml, . ),                \
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histmthCOSP_lmdz.xml, .),\
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histmthNMC_lmdz.xml, . ),\
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histday_lmdz.xml, . ),   \
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histdayCOSP_lmdz.xml, .),\
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histdayNMC_lmdz.xml, . ),\
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histhf_lmdz.xml, . ),    \
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histhfCOSP_lmdz.xml, .), \
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histhfNMC_lmdz.xml, . ), \
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histstn_lmdz.xml, . ),   \
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histins_lmdz.xml, . ),   \
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histstrataer_lmdz.xml, . ),\
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histdaystrataer_lmdz.xml, . ),\
+        (${MODIPSL}/modeles/LMDZ/DefLists/file_def_histLES_lmdz.xml, . )
+
+[RestartFiles]
+List=   (restartphy.nc, restartphy.nc,  startphy.nc)
+
+[OutputText]
+List=   (physiq.def, run.def, tracer.def, guide.def, config.def, vert.def, used_run.def, context_lmdz.xml, field_def_lmdz.xml, file_def_*lmdz.xml, debug.01, lmdz.prt)
+
+[OutputFiles]
+List=   (histmth.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,       Post_1M_histmth),       \
+        (histmthCOSP.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthCOSP.nc,   NONE),                  \
+        (histmthNMC.nc,    ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,    Post_1M_histmthNMC),    \
+        (histday.nc,       ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,       NONE),                  \
+        (histdayCOSP.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,   NONE),                  \
+        (histdayNMC.nc,    ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,    NONE),                  \
+        (histhf.nc,        ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,        NONE),                  \
+        (histhfCOSP.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,    NONE),                  \
+        (histhfNMC.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,     NONE),                  \
+        (histhf3h.nc,      ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,      NONE),                  \
+        (histhf3hm.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,     NONE)
+
+[Post_1M_histmth]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_centered, time_centered_bounds, aire)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFC, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, ave_t2m_daily_max, ave_t2m_daily_min, taux, taux_ter, taux_sic, taux_lic, taux_oce, tauy, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, nettop, pbase, ptop, topl0, topl, tops0, tops, topswad, topswai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, wind10m, od550aer, od550lt1aer, od550_STRAT, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl, swtoaas_ant, swtoacs_ant, swtoacf_ant, swsrfas_ant, swsrfcs_ant, swsrfcf_ant, swtoaas_nat, swtoacs_nat, swtoacf_nat, swsrfas_nat, swsrfcs_nat, swsrfcf_nat, solswai, solswad, swtoacf_zero, swsrfcf_zero, R_ecc, R_peri, R_incl, solaire, rsun1, co2_ppm, CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt, colO3_strat, colO3_trop)
+ChunckJob2D = 100Y 
+TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, tnt, dtdyn, dtphy, dtlschr, tntr, tntscpbl, tnhus, dqdyn, dqvdf, tnhusscpbl, dqphy, evu, mcd)
+ChunckJob3D = 100Y
+Seasonal=ON
+
+[Post_1M_histmthCOSP]
+Patches = ()
+GatherWithInternal = (lon, lat, pressure2, height, sza, time_counter, time_centered, time_centered_bounds, aire)
+TimeSeriesVars2D = (cllcalipso, clmcalipso, clhcalipso, cltcalipso, tclisccp, albisccp, ctpisccp)
+ChunckJob2D = OFF
+TimeSeriesVars3D = (clcalipso, clisccp2_c01, clisccp2_c02, clisccp2_c03, clisccp2_c04, clisccp2_c05, clisccp2_c06, clisccp2_c07, parasol_refl)
+ChunckJob3D = OFF
+Seasonal=OFF
+
+[Post_1M_histmthNMC]
+Patches = ()
+GatherWithInternal = (lon, lat, plev, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = (ta, tro3, zg, hus, hur, ua, va, wap)
+ChunckJob3D = NONE
+Seasonal=ON
+
+[Post_1D_histday]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, precip, pluc, plul, prw, psol, q2m, sens, sicf, slp, soll, sols, snow, SWdnSFC, t2m, t2m_max, t2m_min, topl, tops, tsol, u10m, v10m, wind10m, wind10max)
+ChunckJob2D = 50Y
+TimeSeriesVars3D = (lwcon, iwcon, mc, rhum, rneb, temp, theta, ovap, geop, vitu, vitv, vitw, pres, paprs)
+ChunckJob3D = OFF
+Seasonal=OFF
+
+[Post_1D_histdayCOSP]
+Patches = ()
+GatherWithInternal = (lon, lat, pressure2, height, sza, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = (cllcalipso, clmcalipso, clhcalipso, cltcalipso, tclisccp, albisccp, ctpisccp)
+ChunckJob2D = OFF
+TimeSeriesVars3D = (clcalipso, clisccp2_c01, clisccp2_c02, clisccp2_c03, clisccp2_c04, clisccp2_c05, clisccp2_c06, clisccp2_c07, parasol_refl)
+ChunckJob3D = OFF
+Seasonal=OFF
+
+[Post_1D_histdayNMC]
+Patches = ()
+GatherWithInternal = (lon, lat, plev, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = (ta, zg, hus, hur, ua, va, wap)
+ChunckJob3D = 10Y
+Seasonal=OFF
+
+[Post_HF_histhf]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = (cldt, psol, q2m, slp, precip, pluc, plul, t2m, tsol, u10m, v10m)
+ChunckJob2D = 50Y
+TimeSeriesVars3D = (temp, theta, ovap, vitu, vitv)
+ChunckJob3D = OFF
+Seasonal=OFF
+
+[Post_HF_histhfCOSP]
+Patches = ()
+GatherWithInternal = (lon, lat, height, sza, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = (cllcalipso, clmcalipso, clhcalipso, cltcalipso)
+ChunckJob2D = OFF
+TimeSeriesVars3D = (clcalipso, parasol_refl, cfad_lidarsr532_c01, cfad_lidarsr532_c02, cfad_lidarsr532_c03, cfad_lidarsr532_c04, cfad_lidarsr532_c05, cfad_lidarsr532_c06, cfad_lidarsr532_c07, cfad_lidarsr532_c08, cfad_lidarsr532_c09, cfad_lidarsr532_c10, cfad_lidarsr532_c11, cfad_lidarsr532_c12, cfad_lidarsr532_c13, cfad_lidarsr532_c14, cfad_lidarsr532_c15)
+ChunckJob3D = OFF
+Seasonal=OFF
+
+[Post_HF_histhfNMC]
+Patches = ()
+GatherWithInternal = (lon, lat, plev, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = (ta, zg, hus, hur, ua, va, wap)
+ChunckJob3D = 10Y
+Seasonal=OFF
+
+[Post_3H_histhf3hm]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = (flat, precip, pluc, snow, topl, topl0, SWupTOA, SWupTOAclr, SWdnTOA, SWupSFC, SWupSFCclr, SWdnSFC, SWdnSFCclr, LWupSFC, LWupSFCclr, LWdnSFC, LWdnSFCclr, bils, bils_tke, bils_diss, bils_ec, bils_kinetic, bils_enthalp, bils_latent, sens, cldt, zmax_th)
+ChunckJob2D = OFF
+TimeSeriesVars3D = ()
+ChunckJob3D = OFF
+Seasonal=OFF
+
+[Post_3H_histhf3h]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_centered, time_centered_bounds)
+TimeSeriesVars2D = (phis, tsol, t2m, q2m, ustar, u10m, v10m, psol)
+ChunckJob2D = OFF
+TimeSeriesVars3D = ()
+ChunckJob3D = OFF
+Seasonal=OFF
Index: /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/orchidee.card
===================================================================
--- /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/orchidee.card	(revision 6895)
+++ /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/orchidee.card	(revision 6895)
@@ -0,0 +1,102 @@
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+# Select the PFT and age class set-up. NOTE: the number of PFTs should match the 
+# number in the land cover and maps. The number of age classes does not affect the
+# maps that can be used.
+DefSuffix = 15pft.1ac
+
+# VEGET_UPDATE=0Y : no change in vegetation map. PFTmap.nc should be set only in InitialStateFiles/List. 
+# VEGET_UPDATE=1Y : change vegetation map every year. PFTmap.nc should be set only in BoundaryFiles/List.
+#
+# For VEGET_UPDATE=1Y, the PFTmap is read in January. In BoundaryFiles/List, copy the map for current year, use the variable ${year}.
+VEGET_UPDATE=1Y
+
+# Choice for routing scheme. Set ROUTING below to have following parameters set in orchidee.def : 
+# ROUTING=standard   => in orchidee.def RIVER_ROUTING=y and ROUTING_METHOD=standard
+# ROUTING=simple     => in orchidee.def RIVER_ROUTING=y and ROUTING_METHOD=simple
+# ROUTING=native     => in orchidee.def RIVER_ROUTING=y and ROUTING_METHOD=native (note: routing_simple.nc is used also for this case)
+# ROUTING=off        => in orchidee.def RIVER_ROUTING=n 
+ROUTING=simple
+
+# Specify output level for output files
+# Only the files sechiba_history.nc, sechiba_out_2.nc and sechiba_history_4dim.nc can be set here.
+# output_level_filename=0     : lowest level writing only variables needed for the monitoring
+# output_level_filename=12    : highest level which will output all variables
+# output_level_filename=NONE  : deactivate the file
+output_level_sechiba_history = 11
+output_level_sechiba_out_2 = NONE
+output_level_sechiba_history_4dim = 1
+
+# Specify output frequency for each file [1y, 1mo, 1d, 10800s, 1ts]
+# Settings using WriteFrequency in config.card are not longer used
+output_freq_sechiba_history = 1mo
+output_freq_sechiba_out_2 = 10800s
+output_freq_sechiba_history_4dim = 1mo
+
+
+[InitialStateFiles]
+List=       (${R_IN}/SRF/SOIL/soil_bulk_and_ph.nc, .), \
+            (${R_IN}/SRF/SOIL/soils_param_usda_with_clay_oxisol_coord.nc, soils_param.nc), \
+            (${R_IN}/SRF/SOIL/soils_param_zobler_coord.nc, soilcolor.nc), \
+            (${R_IN}/SRF/cartepente2d_15min4.nc, cartepente2d_15min.nc), \
+            (${R_IN}/SRF/ROUTING/floodplains4.nc, floodplains.nc), \
+            (${R_IN}/SRF/reftemp4.nc, reftemp.nc)
+	
+[BoundaryFiles]
+List= (${R_IN}/SRF/PFTMAPS/CMIP6/ESA-LUH2v2/historical/15PFT.v2023.1/PFTmap_${year}.nc, PFTmap.nc)
+ListNonDel= (${R_IN}/SRF/ROUTING/routing_simple.nc, .)
+
+[SmoothFiles]
+List=()
+
+[ParametersFiles]
+List=   (${SUBMIT_DIR}/PARAM/ORCHIDEE/orchidee.def_v4, orchidee.def)  ,\
+        (${SUBMIT_DIR}/PARAM/ORCHIDEE/orchidee_pft.def_${DefSuffix}_v4, orchidee_pft.def) ,\
+        (${MODIPSL}/modeles/ORCHIDEE_4/src_xml/file_def_orchidee.xml, .)  ,\
+        (${MODIPSL}/modeles/ORCHIDEE_4/src_xml/context_input_orchidee.xml, .) ,\
+        (${MODIPSL}/modeles/ORCHIDEE_4/src_xml/context_routing_orchidee.xml, .)   ,\
+        (${SUBMIT_DIR}/PARAM/ORCHIDEE/LAM/context_orchidee.xml_v4, context_orchidee.xml)   ,\
+        (${MODIPSL}/modeles/ORCHIDEE_4/src_xml/field_def_orchidee.xml, .)
+
+[RestartFiles]
+List=   (sechiba_rest_out.nc, sechiba_rest.nc, sechiba_rest_in.nc), \
+        (routing_restart.nc, routing_restart.nc, routing_start.nc, OPTIONAL)
+
+[OutputText]
+List=   (orchidee.def, orchidee_pft.def, used_orchidee.def, used_orchidee_pft.def, out_orchidee, context_orchidee.xml, context_input_orchidee.xml, field_def_orchidee.xml, file_def_orchidee.xml, context_routing_orchidee.xml, river_desc.nc)
+
+[OutputFiles]
+List=	(sechiba_history.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_history.nc, Post_1M_sechiba_history), \
+        (sechiba_history_4dim.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_history_4dim.nc, NONE), \
+        (sechiba_out_2.nc, ${R_OUT_SRF_O_H}/${PREFIX}_HF_sechiba_out_2.nc, NONE), \
+        (ORCHIDEE_interp_diag.nc, ${R_OUT_SRF_O_D}/${PREFIX}_ORCHIDEE_interp_diag.nc, NONE), \
+        (forcing_by_ORCHIDEE.nc, ${R_OUT_SRF_O_H}/${PREFIX}_forcing_by_ORCHIDEE.nc, NONE), \
+        (diag_routing.nc,   ${R_OUT_SRF_O_D}/${PREFIX}_1D_sechiba_routing.nc, NONE),\
+        (diag_routing_r.nc, ${R_OUT_SRF_O_D}/${PREFIX}_1D_sechiba_routing_r.nc, NONE)
+
+[Post_1M_sechiba_history]
+Patches = ()
+GatherWithInternal = (lon, lat, veget, time_counter, Areas, Contfrac, time_centered, time_centered_bounds)
+TimeSeriesVars2D = (nobiofrac, alb_nir, alb_vis, evap, fluxlat, fluxsens, netrad, qair, rain, runoff, snow, snownobio, snowf, frac_snow, subli, tair, temp_sol, tsol_max, tsol_min, drainage,mrsos, mrso, mrros, mrro, prveg, evspsblveg, evspsblsoi, tran, treeFrac, grassFrac, cropFrac, baresoilFrac, residualFrac, TWBR, riverflow, coastalflow, DelSoilMoist_daily, DelIntercept_daily, DelSWE_daily, delstock_routing_daily, snowmelt, LAImean, hydrographs, evapnu, Wind, lwdown, swdown)
+ChunckJob2D = 200Y
+TimeSeriesVars3D = (lai, maxvegetfrac, vegetfrac, ptn, nee, transpir, inter)
+ChunckJob3D = NONE
+Seasonal = ON
+
+# ALMAOUTPUT
+#TimeSeriesVars2D = (nobiofrac, SWnet, LWnet, Qh, Qle, Qg, Qf, DelSurfHeat, DelColdCont, Snowf, Rainf, Evap, Qs, Qsb, Qsm, DelSoilMoist, DelSWE, DelIntercept, AvgSurfT, RadT, Albedo, SWE, SoilMoist, SoilWet, SoilTemp, PotEvap, ESoil, ECanop, TVeg, RootMoist, SubSnow, ACond, SnowFrac, SAlbedo, SnowDepth, dis, GPP)
+#TimeSeriesVars3D = (ECanop, TVeg, GPP)
+
+
+[Post_HF_sechiba_out_2]
+Patches = ()
+GatherWithInternal = (lon, lat, time_counter, Areas, Contfrac, time_centered, time_centered_bounds)
+TimeSeriesVars2D = (mrsos, mrro)
+ChunckJob2D = 50Y
+TimeSeriesVars3D = (ptn)
+ChunckJob3D = 10Y
+Seasonal = OFF
+
+
Index: /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/stomate.card
===================================================================
--- /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/stomate.card	(revision 6895)
+++ /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/stomate.card	(revision 6895)
@@ -0,0 +1,82 @@
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+# SPINUP_ANALYTIC=y/n : Activate the spinup analytic option to solve the carbon in soil balance
+SPINUP_ANALYTIC=n
+
+# NINPUT_UPDATE=0Y : no change in nitrogen input maps. Nitrogen input maps should be set in InitialStateFiles/List
+# NINPUT_UPDATE=1Y : change in nitrogen input maps every year. Nitrogen input maps should be set only in BoundaryFiles/List.
+NINPUT_UPDATE=1Y
+
+# Specify output level for output files
+# Only the files stomate_history.nc and stomate_ipcc_history.nc can be set here.
+# output_level_filname=0     : lowest level writing only variables needed for the monitoring
+# output_level_filname=10    : highest level which will output all variables
+# output_level_filname=NONE  : deactivate the file
+output_level_stomate_history = 6
+output_level_stomate_ipcc_history = 1
+output_level_stomate_history_4dim = NONE
+output_level_stomate_fixed_dia = NONE
+
+# Specify output frequency for each file [1y, 1mo, 1d]
+# Settings using WriteFrequency in config.card are not longer used
+output_freq_stomate_history = 1mo
+output_freq_stomate_ipcc_history = 1mo
+output_freq_stomate_history_4dim = 1mo
+output_freq_stomate_fixed_dia = 1y
+
+[InitialStateFiles]
+List=   ()
+
+[BoundaryFiles]
+List=   (${R_IN}/SRF/NITROGEN/N_FERTILISATION/NMIP2/synthetic/historical/nfer_crop_nh4_${year}.nc, Nfert_ammo_cropland.nc) ,\
+        (${R_IN}/SRF/NITROGEN/N_FERTILISATION/NMIP2/synthetic/historical/nfer_crop_no3_${year}.nc, Nfert_nitr_cropland.nc) ,\
+        (${R_IN}/SRF/NITROGEN/N_FERTILISATION/NMIP2/synthetic/historical/nfer_pas_nh4_${year}.nc,  Nfert_ammo_pasture.nc) ,\
+        (${R_IN}/SRF/NITROGEN/N_FERTILISATION/NMIP2/synthetic/historical/nfer_pas_no3_${year}.nc,  Nfert_nitr_pasture.nc) ,\
+        (${R_IN}/SRF/NITROGEN/N_DEPOSITION/CMIP6.2/historical/nhx__${year}.nc, Nammonium.nc)\
+        (${R_IN}/SRF/NITROGEN/N_DEPOSITION/CMIP6.2/historical/noy__${year}.nc, Nnitrate.nc)
+
+ListNonDel= (${R_IN}/SRF/NITROGEN/BNF/bnf_1850.nc, bnf.nc) ,\
+	    (${R_IN}/SRF/CO2/CO2_1700_2022_TRENDYv2023.txt, CO2.txt)
+
+[SmoothFiles]
+# global_forest_management exists for years 1700-2015
+List= (${R_IN}/SRF/FOREST_MANAGEMENT/TRUNK_FM/historical/v1/global_forest_management_${year}.nc, FMmap.nc, 1:1:114) ,\
+      (${R_IN}/SRF/FOREST_MANAGEMENT/TRUNK_FM/historical/v1/global_forest_management_2015.nc, FMmap.nc, 115:1:)
+
+[ParametersFiles]
+List=   ()
+
+
+[RestartFiles]
+List=	(stomate_rest_out.nc, stomate_rest.nc, stomate_rest_in.nc)
+
+
+[OutputText]
+List=   ()
+
+[OutputFiles]
+List=	(stomate_history.nc, ${R_OUT_SBG_O_M}/${PREFIX}_1M_stomate_history.nc, Post_1M_stomate_history), \
+        (stomate_ipcc_history.nc, ${R_OUT_SBG_O_M}/${PREFIX}_1M_stomate_ipcc_history.nc, Post_1M_stomate_ipcc_history), \
+        (stomate_history_4dim.nc, ${R_OUT_SBG_O_M}/${PREFIX}_1M_stomate_history_4dim.nc, NONE), \
+        (stomate_fixed_dia.nc, ${R_OUT_SBG_O_Y}/${PREFIX}_1Y_stomate_fixed_dia.nc, NONE)
+
+[Post_1M_stomate_history]
+Patches= ()
+GatherWithInternal= (lon, lat, time_counter, Areas, CONTFRAC)
+TimeSeriesVars2D= (T2M_MONTH,RESOLUTION_X,RESOLUTION_Y,FLUX_PROD_S_LCC_c,FLUX_PROD_M_LCC_c,FLUX_PROD_L_LCC_c,FLUX_PROD_S_HARVEST_c,FLUX_PROD_M_HARVEST_c,FLUX_PROD_L_HARVEST_c)
+ChunckJob2D = 100Y
+TimeSeriesVars3D= (LAI_MEAN,VEGET_MAX,NPP,GPP,HET_RESP,MAINT_RESP,GROWTH_RESP,TOTAL_M_n,TOTAL_M_c,TOTAL_BM_LITTER_c,TOTAL_BM_LITTER_n,ADAPTATION,REGENERATION,SOIL_ACTIVE_c,SOIL_ACTIVE_n,SOIL_SLOW_c,SOIL_SLOW_n,SOIL_PASSIVE_c,SOIL_PASSIVE_n,SOIL_SURF_c,SOIL_SURF_n,LITTER_STR_AB_c,LITTER_STR_AB_n,LITTER_STR_BE_c,LITTER_STR_BE_n,LITTER_MET_AB_c,LITTER_MET_AB_n,LITTER_MET_BE_c,LITTER_MET_BE_n,LITTER_WOD_AB_c,LITTER_WOD_AB_n,LITTER_WOD_BE_c,LITTER_WOD_BE_n,TOTAL_TURN_c,TOTAL_TURN_n,LEAF_M_n,LEAF_M_c,N_UPTAKE_NO3,N_UPTAKE_NH4,NH3_EMISSION,N2_EMISSION,N2O_EMISSION,NOX_EMISSION,NBP_pool_c,NBP_flux_c,MBC_NBP1_c,MBC_NBP1_n,MBC_NBP2_c,MBC_NBP2_n)
+ChunckJob3D = 100Y
+Seasonal=ON
+
+[Post_1M_stomate_ipcc_history]
+Patches= ()
+GatherWithInternal= (lon, lat, time_counter, Areas, CONTFRAC)
+TimeSeriesVars2D=(cVeg, cLitter, cSoil, cProduct, cMassVariation, lai, gpp, ra, npp, rh, fFire, fHarvest, fLuc, nbp_c, nbp_n, fVegLitter, fLitterSoil, cLeaf, cWood, cRoot, cMisc, cLitterAbove, cLitterBelow, cSoilFast, cSoilMedium, cSoilSlow, landCoverFrac, treeFracPrimDec, treeFracPrimEver, c3PftFrac, c4PftFrac, rGrowth, rMaint, nppLeaf, nppWood, nppRoot, nep_c, nep_n, fWoodharvest)
+ChunckJob2D = 100Y
+TimeSeriesVars3D=()
+ChunckJob3D = NONE
+Seasonal=ON
+
Index: /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/xios.card
===================================================================
--- /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/xios.card	(revision 6895)
+++ /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/COMP/xios.card	(revision 6895)
@@ -0,0 +1,29 @@
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List=	()
+
+[BoundaryFiles]
+List=	()
+ListNonDel= ()
+
+[SmoothFiles]
+List=	()
+
+[ParametersFiles]
+List=   (${SUBMIT_DIR}/PARAM/XIOS/iodef.xml, .)
+
+[RestartFiles]
+List=	(NONE)
+
+[OutputText]
+List=   (xios_client_*, xios_server_*, out_xios.x.err.*, out_xios.x.out.*)
+
+[OutputFiles]
+List=   ()
+
+
+
Index: /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/config.card
===================================================================
--- /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/config.card	(revision 6895)
+++ /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_ERA5_LAM.orch4/config.card	(revision 6895)
@@ -0,0 +1,151 @@
+# $Id$
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+JobName=TEST-amip-ERA5-LAM.01.orch4 
+LongName=""
+TagName=ICOLMDZOR
+#D- Choice of experiment in EXPERIEMENTS directory
+ExpType=ICOLMDZOR/LAM/amip_ERA5_LAM.orch4
+ExperimentName=amip_LAM
+SpaceName=TEST
+
+# Use CreatBasefolder and CreateName to define the folder for the input files previously created by an expermient CREATE
+# The final folder will be ${ARCHIVE}/CreateBasefolder/CreateName
+# $ARCHIVE is added in dynamico.card and lmdz.card before the CreateBasefolder because it is not possible to add it at this level in config.card 
+CreateBasefolder=IGCM_OUT/ICOLMDZOR/DEVT/CE0L
+CreateName=CREATE-amip-ERA5-LAM.01
+
+ForcingName=ERA5
+ForcingDate=062000
+
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Begin and end of job 
+#D-- "YYYY-MM-DD"
+#D-- The year for this experiement has no significance as the SST/SIC are mean values over a whole period. 
+#D-- Change the DateEnd if you change CalendarType, set DateEnd=2000-12-31 for leap or noleap.
+DateBegin=2000-06-01
+DateEnd=2000-06-05
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=5D
+#============================
+#-- Compression level for netcdf output files
+CompressionLevel=2
+#============================
+# Source following file with module settings
+EnvFile=${SUBMIT_DIR}/../ARCH/arch.env
+#============================
+# OptMode indicates the optimization mode choosen during compilation
+# This variable is used in the executable name
+OptMode=prod
+
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ICO= (dynamico, DYNAMICO)
+ATM= (lmdz, LMDZ)
+SRF= (orchidee, ORCHIDEE)
+SBG= (stomate, ORCHIDEE)
+IOS= (xios, XIOS)
+
+#========================================================================
+#D-- Executable -
+[Executable]
+ICO= (icosa_lmdz_orch4_${OptMode}.exe, icosa_lmdz.exe, 48MPI)
+ATM= ("","")
+SRF= ("", "")
+SBG= ("", "")
+IOS= (xios_server_${OptMode}.exe, xios.x, 1MPI)
+
+#========================================================================
+#D-- Restarts -
+[Restarts]
+#D- If you want a GENERAL RULE FOR RESTARTS, put this flag to 'y'
+OverRule=n
+#D- Last day of the experience used as restart
+RestartDate=
+#D- Define restart simulation name
+RestartJobName=
+#D- Path Server Group Login
+RestartPath=
+# Attention login depend de la machine
+
+
+#========================================================================
+#D-- NO Post -
+[Post]
+PackFrequency=1Y
+TimeSeriesFrequency=NONE
+SeasonalFrequency=NONE
+MetricsPCMDI=FALSE
+
+
+#========================================================================
+#D-- ICO -
+[ICO]
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=
+# Define restart simulation name
+RestartJobName=
+RestartPath=
+
+#========================================================================
+#D-- ATM -
+[ATM]
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=
+# Define restart simulation name
+RestartJobName=
+RestartPath=
+
+#========================================================================
+#D-- SRF -
+[SRF]
+# If config_Restarts_OverRule == 'n' next 4 params are read
+Restart= n
+# Last day of the experience used as restart for this component if Restart=y
+RestartDate=2319-12-31
+#D- Define restart simulation name for this component
+RestartJobName=CM607-LR-pdCtrl-01
+#D- Path Server Group Login
+RestartPath=${R_IN}/RESTART/IPSLCM6/DEVT/pdControl
+
+#========================================================================
+#D-- SBG - STOMATE
+[SBG]
+# If config_Restarts_OverRule == 'n' next 4 params are read
+Restart= n
+# Last day of the experience used as restart for this component if Restart=y
+RestartDate=2319-12-31
+#D- Define restart simulation name for this component
+RestartJobName=CM607-LR-pdCtrl-01
+#D- Path Server Group Login
+RestartPath=${R_IN}/RESTART/IPSLCM6/DEVT/pdControl
+
+#========================================================================
+#D-- IOS -
+[IOS]
+# If config_Restarts_OverRule == 'n' next 4 params are read
+Restart= n
+#-- Last day of the experience used as restart for this component if Restart=y
+RestartDate=
+# Define restart simulation name for this component
+RestartJobName=
+RestartPath=
+# Old component name for restart (if empty, use new name)
+
+
Index: /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_LAM/COMP/lmdz.card
===================================================================
--- /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_LAM/COMP/lmdz.card	(revision 6894)
+++ /CONFIG/UNIFORM/v7/IPSLCM7/EXPERIMENTS/ICOLMDZOR/LAM/amip_LAM/COMP/lmdz.card	(revision 6895)
@@ -55,5 +55,4 @@
 # With 3D ozone
 List= (${R_OUT}/${config_UserChoices_TagName}/${config_UserChoices_SpaceName}/CE0L/${config_UserChoices_CreateName}/ATM/Output/Boundary/${config_UserChoices_CreateName}_${year}_limit.nc, limit.nc),\
-            (${R_IN}/ATM/AEROSOLS/CMIP6/v1/144x142/L${RESOL_ATM_Z}/aerosols1850_from_inca.nc,   aerosols.nat.nc),\
             (${R_IN}/ATM/AEROSOLS/CMIP6/v1/144x142/L${RESOL_ATM_Z}/aerosols${year}_from_inca.nc,   aerosols${year}.nc),\
             (${R_IN}/ATM/SOLAR/CMIP6/v3.1corr/Historical/solarforcing_v3.1_daily_6bands_${year}.nc, solarforcing.nc),\
Index: /CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/PARAM/ORCHIDEE/LAM/context_orchidee.xml_v4
===================================================================
--- /CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/PARAM/ORCHIDEE/LAM/context_orchidee.xml_v4	(revision 6895)
+++ /CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/PARAM/ORCHIDEE/LAM/context_orchidee.xml_v4	(revision 6895)
@@ -0,0 +1,472 @@
+<!-- =========================================================================================================== -->
+<!-- ORCHIDEE context                                                                                            -->
+<!-- context_orchidee.xml : Configuration file for ORCHIDEE for production of output files using XIOS            -->
+<!--                                                                                                             -->
+<!-- $HeadURL: svn://forge.ipsl.fr/orchidee/trunk/ORCHIDEE/src_xml/context_orchidee.xml $                -->
+<!-- $Date: 2023-08-29 11:56:15 +0200 (Tue, 29 Aug 2023) $                                                     -->
+<!-- $Revision: 8149 $                                                                                           -->
+<!-- =========================================================================================================== -->
+<context id="orchidee">
+  
+  
+  <!-- =========================================================================================================== -->
+  <!-- Definition of all existing variables                                                                        -->
+  <!-- DO NOT CHANGE THIS FILE                                                                                     -->
+  <!-- =========================================================================================================== -->
+  <field_definition src="./field_def_orchidee.xml"/>
+  
+  
+  <!-- =========================================================================================================== -->
+  <!-- Definition of output files                                                                                  -->
+  <!-- Definition of variables included in the different files                                                     -->
+  <!-- CHANGE THIS FILE BY ADDING THE FILES AND/OR VARIABLES YOU WANT TO PRODUCE                                   -->
+  <!-- Only variables and groups existing in field_def_orchidee.xml can be used                                    -->
+  <!-- =========================================================================================================== -->
+  <file_definition src="./file_def_orchidee.xml"/>
+  
+  
+  <!-- =========================================================================================================== -->
+  <!-- Definition of horizontal domain                                                                             -->
+  <!-- =========================================================================================================== -->
+  <domain_definition>
+
+    <!-- Model horizontal domain -->
+    <domain id="domain_landpoints"/>
+
+    <!-- Standard output horizontal domain:                                                                         -->
+    <!-- The same as the model grid is used by default for regular longitude-latitude grid. If the model is running -->
+    <!-- on an unstructred grid, this domain_landpoints will be changed to domain_landoints_regular, see below.     -->   
+    <domain id="domain_landpoints_out" domain_ref="domain_landpoints"/>
+
+    <!-- Interpolated output horizontal domain:                                                                   -->
+    <!-- This domain is used for output when the modele is running on an unstructured grid. For that case,        -->
+    <!-- the ORCHIDEE fortran code will switch so that domain_landpoints_out make reference to this domain.       -->
+    <domain id="domain_landpoints_regular" ni_glo="300" nj_glo="232" type="rectilinear"  >
+      <generate_rectilinear_domain lat_start="38" lat_end="55" lon_start="-9" lon_end="13"/>
+      <interpolate_domain order="1" renormalize="true" detect_missing_value="true"/>
+    </domain>
+
+    <!-- Output horizontal domain used for CMIP6 simulations:                                                    -->
+    <!-- This domain is only used when running with the workflow CMIP6 and dr2xml_orchidee.xml files.            -->
+    <domain id="greordered"  domain_ref="domain_landpoints">
+      <reorder_domain invert_lat="true" shift_lon_fraction="0.5" min_lon="0" max_lon="360" />
+    </domain>
+
+  </domain_definition>
+
+  <!-- =========================================================================================================== -->
+  <!-- Definition of vertical axis and extra dimensions                                                            -->
+  <!-- =========================================================================================================== -->
+  <axis_definition>
+    <!-- Vertical axis and extra dimensions -->
+    <axis id="nvm" name="veget" standard_name="model_level_number" long_name="Vegetation types" unit="1"/>
+    <axis id="nvmap" name="vegetmap" standard_name="model_level_number" long_name="Vegetation types in input file" unit="1"/>
+    <axis id="nlut" name="nlut" standard_name="landusetype4" long_name="Number of land use types specified by LUMIP" unit="1"/>
+    <axis id="ncarb" name="ncarb" standard_name="ncarb" long_name="Number of soil carbon pools" unit="1"/>
+    <axis id="nlaip1" name="laiax" standard_name="model_level_number" long_name="Interface of LAI layers" unit="1"/>
+    <axis id="ngrnd" name="solth" standard_name="model_level_number" long_name="Thermal soil levels" unit="m"/>
+    <axis id="nstm" name="soiltyp" standard_name="model_level_number" long_name="Soil tiles" unit="1"/>
+    <axis id="nnobio" name="nobio" standard_name="model_level_number" long_name="Other surface types" unit="1"/>
+    <axis id="albtyp" standard_name="model_level_number" long_name="Albedo types" unit="1"/>
+    <axis id="nlevels_tot" standard_name="model_level_number" long_name="Total number of vertical canopy levels" unit="1"/>
+    <axis id="nslm" name="solay" standard_name="model_level_number" long_name="Hydrological soil levels" unit="1"/>
+    <axis id="nscm" name="nscm" standard_name="model_level_number" long_name="Soil type" unit="1"/>
+    <axis id="nsnow" name="snowlev" standard_name="model_level_number" long_name="Snow levels" unit="1"/>
+    <axis id="nice" name="icelev" standard_name="model_level_number" long_name="Ice levels" unit="1"/>
+    <axis id="nbhtu" name="nbhtu" standard_name="model_level_number" long_name="HTU dimension for routing" unit="1"/>
+    <axis id="nbasmon" name="nbasmon" standard_name="model_level_number" long_name="HTU monitoring dimension" unit="1"/>
+
+    <!-- Vertical axis and extra dimensions in stomate -->
+    <axis id="P10" standard_name="model_level_number" long_name="Pool 10 years" unit="1"/>
+    <axis id="P100" standard_name="model_level_number" long_name="Pool 100 years" unit="1"/>
+    <axis id="P11" standard_name="model_level_number" long_name="Pool 10 years + 1" unit="1"/>
+    <axis id="P101" standard_name="model_level_number" long_name="Pool 100 years + 1" unit="1"/>
+    <axis id="nshort" standard_name="model_level_number" long_name="short-lived product pool" unit="yr"/>
+    <axis id="nmedium" standard_name="model_level_number" long_name="medium-lived product pool" unit="yr"/>
+    <axis id="nlong" standard_name="model_level_number" long_name="long-lived product pool" unit="yr"/>
+    <axis id="nshort1" standard_name="model_level_number" long_name="short-lived product pool" unit="yr"/>
+    <axis id="nmedium1" standard_name="model_level_number" long_name="medium-lived product pool" unit="yr"/>
+    <axis id="nlong1" standard_name="model_level_number" long_name="long-lived product pool" unit="yr"/>
+    <axis id="ncirc" standard_name="model_level_number" long_name="circumference classes" unit="1"/>
+    <axis id="ncut" standard_name="model_level_number" long_name="types of cuttings" unit="1"/>
+    <axis id="nparts" standard_name="model_level_number" long_name="biomass components" unit="1"/>
+    <axis id="nlctypes" standard_name="model_level_number" long_name="land cover types" unit="1"/>
+    <axis id="noutdiaclass" standard_name="model_level_number" long_name="number of output classes" unit="1"/>	
+    <axis id="lai_time" />
+    <axis id="time12" />
+    <axis id="time1" />
+  </axis_definition>
+
+
+  <!-- =========================================================================================================== -->
+  <!-- Definition of model and output grids                                                                        -->
+  <!-- A grid is a combination of the horizontal domain and optionally one or several axes.                        -->
+  <!-- All grids are specified twice: the grid used by the model "grid_axis", and the grid used for output         -->
+  <!-- "grid_axis_out", where "axis" represents one or several axes used in the grid, see below.                   -->
+  <!-- When the model runs on the default regular longitude-latitude grid, "grid_axis_out" is the same as          -->
+  <!-- "grid_axis", no interpolation will be done. When the model runs on an unstructured grid, grid_axis_out is   --> 
+  <!-- different from grid_axis and XIOS will do the interpolation.                                                -->
+  <!-- =========================================================================================================== -->
+  <grid_definition>
+
+    <!-- Definition of horizontal grid, 2D grid -->
+    <grid id="grid_landpoints">
+      <domain domain_ref="domain_landpoints" />
+    </grid>
+    <grid id="grid_landpoints_out">
+      <domain domain_ref="domain_landpoints_out" />
+    </grid>
+
+    <!-- Definition of grids with fields having one extra dimension (one vertical axis or another dimension), 3D grids -->
+    <grid id="grid_nvm">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nvm" />
+    </grid>
+    <grid id="grid_nvm_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nvm" />
+    </grid>
+
+    <grid id="grid_nvmap">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nvmap" />
+    </grid>
+    <grid id="grid_nvmap_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nvmap" />
+    </grid>
+
+    <grid id="grid_nlut">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nlut" />
+    </grid>
+    <grid id="grid_nlut_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nlut" />
+    </grid>
+    
+    <grid id="grid_ncarb">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="ncarb" />
+    </grid>
+    <grid id="grid_ncarb_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="ncarb" />
+    </grid>
+    
+    <grid id="grid_nstm">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nstm" />
+    </grid>
+    <grid id="grid_nstm_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nstm" />
+    </grid>
+
+    <grid id="grid_nslm">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nslm" />
+    </grid>
+    <grid id="grid_nslm_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nslm" />
+    </grid>
+
+    <grid id="grid_nscm">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nscm" />
+    </grid>
+    <grid id="grid_nscm_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nscm" />
+    </grid>
+
+    <grid id="grid_ngrnd">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="ngrnd" />
+    </grid>
+    <grid id="grid_ngrnd_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="ngrnd" />
+    </grid>
+
+    <grid id="grid_ngrnd_nvm">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="ngrnd" />
+      <axis axis_ref="nvm" />
+    </grid>
+    <grid id="grid_ngrnd_nvm_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="ngrnd" />
+      <axis axis_ref="nvm" />
+    </grid>
+
+    <grid id="grid_albtyp">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="albtyp" />
+    </grid>
+    <grid id="grid_albtyp_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="albtyp" />
+    </grid>
+
+    <grid id="grid_nnobio">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nnobio" />
+    </grid>
+    <grid id="grid_nnobio_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nnobio" />
+    </grid>
+
+    <grid id="grid_nlaip1">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nlaip1" />
+    </grid>
+    <grid id="grid_nlaip1_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nlaip1" />
+    </grid>
+
+    <grid id="grid_nsnow">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nsnow" />
+    </grid>
+    <grid id="grid_nsnow_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nsnow" />
+    </grid>
+
+    <grid id="grid_P10">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="P10" />
+    </grid>
+    <grid id="grid_P10_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="P10" />
+    </grid>
+    
+    <grid id="grid_P11">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="P11" />
+    </grid>
+    <grid id="grid_P11_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="P11" />
+    </grid>
+
+    <grid id="grid_P100">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="P100" />
+    </grid>
+    <grid id="grid_P100_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="P100" />
+    </grid>
+
+    <grid id="grid_P101">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="P101" />
+    </grid>
+    <grid id="grid_P101_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="P101" />
+    </grid>
+
+    <grid id="grid_ncut">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="ncut" />
+    </grid>
+    <grid id="grid_ncut_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="ncut" />
+    </grid>
+
+    <grid id="grid_nshort">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nshort" />
+    </grid>
+    <grid id="grid_nshort_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nshort" />
+    </grid>
+
+    <grid id="grid_nmedium">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nmedium" />
+    </grid>
+    <grid id="grid_nmedium_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nmedium" />
+    </grid>
+
+    <grid id="grid_nlong">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nlong" />
+    </grid>
+    <grid id="grid_nlong_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nlong" />
+    </grid>
+
+    <grid id="grid_nshort1">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nshort1" />
+    </grid>
+    <grid id="grid_nshort1_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nshort1" />
+    </grid>
+
+    <grid id="grid_nmedium1">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nmedium1" />
+    </grid>
+    <grid id="grid_nmedium1_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nmedium1" />
+    </grid>
+
+    <grid id="grid_nlong1">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nlong1" />
+    </grid>
+    <grid id="grid_nlong1_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nlong1" />
+    </grid>
+
+    <grid id="grid_time12">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="time12" />
+    </grid>
+    <grid id="grid_time12_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="time12" />
+    </grid>
+
+    <grid id="grid_nlctypes">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nlctypes" />
+    </grid>
+    <grid id="grid_nlctypes_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nlctypes" />
+    </grid>
+
+    <grid id="grid_nice">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nice" />
+    </grid>
+    <grid id="grid_nice_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nice" />
+    </grid>
+
+    <grid id="grid_nbhtu">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nbhtu" />
+    </grid>
+    <grid id="grid_nbhtu_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nbhtu" />
+    </grid>
+
+    <grid id="grid_nbasmon">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nbasmon" />
+    </grid>
+    <grid id="grid_nbasmon_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nbasmon" />
+    </grid>
+
+    <!-- Definition of grids with fields having two extra dimensions, 4D grids -->
+    <grid id="grid_nslm_nstm_read">
+      <domain domain_ref="domain_landpoints" >
+        <interpolate_domain order="1"/>
+      </domain>
+      <axis axis_ref="nslm" />
+      <axis axis_ref="nstm" />
+    </grid>
+    <grid id="grid_nslm_nstm"> 
+       <domain domain_ref="domain_landpoints" /> 
+       <axis axis_ref="nslm" /> 
+       <axis axis_ref="nstm" />
+    </grid> 
+    <grid id="grid_nslm_nstm_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nslm" />
+      <axis axis_ref="nstm" />
+    </grid>
+
+    <grid id="grid_nvm_nstm">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nvm" />
+      <axis axis_ref="nstm" />
+    </grid>
+    <grid id="grid_nvm_nstm_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nvm" />
+      <axis axis_ref="nstm" />
+    </grid>
+
+    <grid id="grid_nvm_nslm">
+      <domain domain_ref="domain_landpoints" />
+      <axis axis_ref="nvm" />
+      <axis axis_ref="nslm" />
+    </grid>
+    <grid id="grid_nvm_nslm_out">
+      <domain domain_ref="domain_landpoints_out" />
+      <axis axis_ref="nvm" />
+      <axis axis_ref="nslm" />
+    </grid>
+
+    <grid id="grid_nvm_ncirc">
+      <domain id="domain_landpoints" />
+      <axis id="nvm" />
+      <axis id="ncirc" />
+    </grid>
+    <grid id="grid_nvm_ncirc_out">
+      <domain id="domain_landpoints_out" />
+      <axis id="nvm" />
+      <axis id="ncirc" />
+    </grid>
+
+    <grid id="grid_nvm_nlevelstot">
+      <domain id="domain_landpoints" />
+      <axis id="nvm" />
+      <axis id="nlevels_tot" />
+    </grid>
+    <grid id="grid_nvm_nlevelstot_out">
+      <domain id="domain_landpoints_out" />
+      <axis id="nvm" />
+      <axis id="nlevels_tot" />
+    </grid>
+
+    <grid id="grid_nvm_nparts">
+      <domain id="domain_landpoints" />
+      <axis id="nvm" />
+      <axis id="nparts" />
+    </grid>
+    <grid id="grid_nvm_nparts_out">
+      <domain id="domain_landpoints_out" />
+      <axis id="nvm" />
+      <axis id="nparts" />
+    </grid>
+
+    <grid id="grid_nvm_noutdiaclass">
+      <domain id="domain_landpoints" />
+      <axis id="nvm" />
+      <axis id="noutdiaclass" />
+    </grid>
+    <grid id="grid_nvm_noutdiaclass_out">
+      <domain id="domain_landpoints_out" />
+      <axis id="nvm" />
+      <axis id="noutdiaclass" />
+    </grid>
+
+  </grid_definition>
+  
+  
+</context>
