Index: CONFIG/UNIFORM/v5/LMDZ_v5/AA_make
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/AA_make	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/AA_make	(revision 1934)
@@ -0,0 +1,72 @@
+#- $Id$
+#-
+#- AA_make for LMDZ_v5 configuration
+#- 
+
+all : 
+	if [ -s ./.resol ] ; then $(M_K) `head -1 .resol |cut -c 8-` ; else $(M_K) LMD9695-L39 ; fi
+
+LMD4443 : libioipsl lmdz44x43x19 verif	
+	echo "noORCAxLMD4443" >.resol
+	echo "RESOL_ATM_3D=44x43x19" >>.resol
+
+LMD5655 : libioipsl lmdz56x55x19 verif	
+	echo "noORCAxLMD5655" >.resol
+	echo "RESOL_ATM_3D=56x55x19" >>.resol
+
+LMD9671 : libioipsl lmdz96x71x19 verif	
+	echo "noORCAxLMD9671" >.resol
+	echo "RESOL_ATM_3D=96x71x19" >>.resol
+
+LMD9695 : libioipsl lmdz96x95x19 verif
+	echo "noORCAxLMD9695" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+
+LMD9695-L39 : libioipsl lmdz96x95x39 verif
+	echo "noORCAxLMD9695-L39" >.resol
+	echo "RESOL_ATM_3D=96x95x39" >>.resol
+
+LMD144142 : libioipsl lmdz144x142x19 verif
+	echo "noORCAxLMD144142" >.resol
+	echo "RESOL_ATM_3D=144x142x19" >>.resol
+
+LMD144142-L39 : libioipsl lmdz144x142x39 verif
+	echo "noORCAxLMD144142-L39" >.resol
+	echo "RESOL_ATM_3D=144x142x39" >>.resol
+
+
+libioipsl : ../../modeles/IOIPSL/src
+	(cd ../../modeles/IOIPSL/src ; $(M_K) -f Makefile)
+
+lmdz44x43x19:
+	$(M_K) lmdz RESOL_LMDZ=44x43x19
+
+lmdz56x55x19:
+	$(M_K) lmdz RESOL_LMDZ=56x55x19
+
+lmdz96x71x19:
+	$(M_K) lmdz RESOL_LMDZ=96x71x19
+
+lmdz96x95x19:
+	$(M_K) lmdz RESOL_LMDZ=96x95x19
+
+lmdz96x95x39:
+	$(M_K) lmdz RESOL_LMDZ=96x95x39
+
+lmdz144x142x19:
+	$(M_K) lmdz RESOL_LMDZ=144x142x19
+
+lmdz144x142x39:
+	$(M_K) lmdz RESOL_LMDZ=144x142x39
+
+lmdz: ../../modeles/LMDZ
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -d $(RESOL_LMDZ) -cosp true -v false -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_$(RESOL_LMDZ)_phylmd_para.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -d $(RESOL_LMDZ) -cosp true -v false -parallel mpi -arch $(FCM_ARCH) gcm ; cp bin/gcm_$(RESOL_LMDZ)_phylmd_para.e ../../bin/gcm.e ; )
+
+
+verif: ../../bin
+	ls -l ../../bin
+
+clean :
+	(rm -rf $(LIBDIR)/* ; rm -f $(BINDIR)/* ; rm -rf ../../modeles/LMDZ/libo/* ../../modeles/LMDZ/.lock ; )
+
Index: CONFIG/UNIFORM/v5/LMDZ_v5/AA_make.ldef
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/AA_make.ldef	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/AA_make.ldef	(revision 1934)
@@ -0,0 +1,11 @@
+#- $Id$
+#---------------------------------------------------------------------
+#- Initialisations for IPSLCM5 configuration
+#---------------------------------------------------------------------
+SHELL = /bin/sh
+#---------------------------------------------------------------------
+LIBDIR = ../../lib
+BINDIR = ../../bin
+MODDIR = $(LIBDIR)
+#---------------------------------------------------------------------
+#-
Index: CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_amip/COMP/lmdz.card
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_amip/COMP/lmdz.card	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_amip/COMP/lmdz.card	(revision 1934)
@@ -0,0 +1,46 @@
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+
+[BoundaryFiles]
+List=   (${R_INIT}/ATM/${config_UserChoices_TagName}/AMIP/amipbc_sst_360x180_${year}.nc, amipbc_sst_1x1.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/AMIP/amipbc_sic_360x180_${year}.nc, amipbc_sic_1x1.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Ozone/HYBRIDE/v2.clim/tro3_${year}.new.nc, climoz.nc)
+
+ListNonDel= (${R_INIT}/ATM/${config_UserChoices_TagName}/Albedo.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/ECDYN.nc.20020101, ECDYN.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/ECPHY.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Relief.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Rugos.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/landiceref.nc, .)
+### Please, add the o2a.nc file to fit with the coupled configuration associated
+#	(${R_INIT}/ATM/IPSLCM5A/ORCA2.3x${RESOL_ATM}/o2a.nc, o2a.nc)
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}_${LMDZ_Physics}, physiq.def), \
+	(${SUBMIT_DIR}/PARAM/gcm.def_${RESOL_ATM_3D},    gcm.def),    \
+        (${SUBMIT_DIR}/PARAM/output.def_low        ,  output.def),    \
+	(${SUBMIT_DIR}/PARAM/config.def_actuel     ,  config.def),    \
+	(${SUBMIT_DIR}/PARAM/run.def, .), \
+	(${SUBMIT_DIR}/PARAM/traceur.def, .)
+
+
+[RestartFiles]
+List= (NONE)
+
+[OutputText]
+List=   (physiq.def, gcm.def, run.def, config.def, traceur.def, output.def)
+
+[OutputFiles]
+List=	( limit.nc,       ${R_OUT_ATM_O}/Boundary/${config_UserChoices_JobName}_${year}_limit.nc,       NONE), \
+	( climoz_LMDZ.nc, ${R_OUT_ATM_O}/Boundary/${config_UserChoices_JobName}_${year}_climoz_LMDZ.nc, NONE), \
+	( start.nc,       ${R_OUT_ATM_O}/Restart/${config_UserChoices_JobName}_${year}_start.nc,        NONE), \
+	( startphy.nc,    ${R_OUT_ATM_O}/Restart/${config_UserChoices_JobName}_${year}_startphy.nc,     NONE), \
+	( grilles_gcm.nc, ${R_OUT_ATM_O}/Grid/${config_UserChoices_JobName}_${year}_grilles_gcm.nc,     NONE)
+
Index: CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_amip/config.card
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_amip/config.card	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_amip/config.card	(revision 1934)
@@ -0,0 +1,77 @@
+# $Id$
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=ELI-${RESOL_ATM_3D}
+LongName="create_etat0_limit for interannuel amip run"
+TagName=LMDZ
+#D- Choice of experiment in EXPERIEMENTS directory
+ExpType=LMDZ/CREATE_amip
+#============================
+#D-- leap, noleap, 360d
+CalendarType=noleap
+#D-- Begin and end of job 
+#D-- "YYYY-MM-DD"
+DateBegin=1979-01-01
+DateEnd=2005-12-31
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=1Y
+#============================
+PeriodNb=50
+
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ)
+
+#========================================================================
+#D-- Executable -
+[Executable]
+#D- For each component, Real name of executable, Name of executable for oasis
+ATM= (create_etat0_limit.e, create_etat0_limit.e)
+
+#========================================================================
+#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-- ATM -
+[ATM]
+#
+WriteFrequency="1Y"
+# 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=
+# Old component name for restart (if empty, use new name)
+OldName=
+
+
+#========================================================================
+#D-- NO Post -
+[Post]
+RebuildFrequency=NONE
+TimeSeriesFrequency=NONE
+SeasonalFrequency=NONE
+#========================================================================
Index: CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_clim/COMP/lmdz.card
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_clim/COMP/lmdz.card	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_clim/COMP/lmdz.card	(revision 1934)
@@ -0,0 +1,46 @@
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+
+[BoundaryFiles]
+List= ()
+
+ListNonDel= (${R_INIT}/ATM/${config_UserChoices_TagName}/Albedo.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/ECDYN.nc.20020101, ECDYN.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/ECPHY.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Relief.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Rugos.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/landiceref.nc, .) \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Ozone/HYBRIDE/v2.clim/tro3_1995.new.nc, climoz.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/AMIP/sst_bc_clim.nc, amipbc_sst_1x1.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/AMIP/sic_bc_clim.nc, amipbc_sic_1x1.nc)
+### Please, add the o2a.nc file to fit with the coupled configuration associated
+#	(${R_INIT}/ATM/IPSLCM5A/ORCA2.3x${RESOL_ATM}/o2a.nc, o2a.nc)
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}_${LMDZ_Physics}, physiq.def), \
+	(${SUBMIT_DIR}/PARAM/gcm.def_${RESOL_ATM_3D},    gcm.def),    \
+        (${SUBMIT_DIR}/PARAM/output.def_low        ,  output.def),    \
+	(${SUBMIT_DIR}/PARAM/config.def_actuel     ,  config.def),    \
+	(${SUBMIT_DIR}/PARAM/run.def, .), \
+	(${SUBMIT_DIR}/PARAM/traceur.def, .)
+
+
+[RestartFiles]
+List= (NONE)
+
+[OutputText]
+List=   (physiq.def, gcm.def, run.def, config.def, traceur.def, output.def)
+
+[OutputFiles]
+List=	( limit.nc,       ${R_OUT_ATM_O}/Boundary/${config_UserChoices_JobName}_clim_limit.nc,   NONE), \
+	( climoz_LMDZ.nc, ${R_OUT_ATM_O}/Boundary/${config_UserChoices_JobName}_climoz_LMDZ.nc,  NONE), \
+	( start.nc,       ${R_OUT_ATM_O}/Restart/${config_UserChoices_JobName}_clim_start.nc,    NONE), \
+	( startphy.nc,    ${R_OUT_ATM_O}/Restart/${config_UserChoices_JobName}_clim_startphy.nc, NONE), \
+	( grilles_gcm.nc, ${R_OUT_ATM_O}/Grid/${config_UserChoices_JobName}_grilles_gcm.nc,      NONE)
Index: CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_clim/config.card
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_clim/config.card	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/CREATE_clim/config.card	(revision 1934)
@@ -0,0 +1,75 @@
+# $Id$
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName= ELC-${RESOL_ATM_3D}
+LongName="create_etat0_limit for climatology run"
+TagName=LMDZ
+#D- Choice of experiment in EXPERIEMENTS directory
+ExpType=LMDZ/CREATE_clim
+#============================
+#D-- leap, noleap, 360d
+CalendarType=360d
+#D-- Begin and end of job 
+#D-- "YYYY-MM-DD"
+DateBegin=1979-01-01
+DateEnd=1979-12-30
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=1Y
+
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ)
+
+#========================================================================
+#D-- Executable -
+[Executable]
+#D- For each component, Real name of executable, Name of executable for oasis
+ATM= (create_etat0_limit.e, create_etat0_limit.e)
+
+#========================================================================
+#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-- ATM -
+[ATM]
+#
+WriteFrequency="1Y"
+# 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=
+# Old component name for restart (if empty, use new name)
+OldName=
+
+
+#========================================================================
+#D-- NO Post -
+[Post]
+RebuildFrequency=NONE
+TimeSeriesFrequency=NONE
+SeasonalFrequency=NONE
+#========================================================================
Index: CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/clim/COMP/lmdz.card
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/clim/COMP/lmdz.card	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/clim/COMP/lmdz.card	(revision 1934)
@@ -0,0 +1,206 @@
+# $Id$
+
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+# Physics package to use : 
+# AP for standard/old physics(defalut), can be used with LMDZ4_AR5 or LMDZ5/trunk sources 
+# NPv3.1 for new physics, to be used with LMDZ5/trunk revision 1554 or later
+LMDZ_Physics=AP
+#LMDZ_Physics=NPv3.1
+
+# You can choose here the name of create_etat0 you want to use
+# If noting is defined here, the job will get the JobName in the file  ../CREATE/config.card
+CREATE=ELC-${RESOL_ATM_3D}
+ByPass_hgardfou_teta=n
+ByPass_hgardfou_mats=n
+LMDZ_NbPeriod_adjust=3
+# LMDZ_NbPeriod_adjust=0
+# To force usage of this Bands file, put LMDZ_NbPeriod_adjust=0 and replace XXXXXXX by Restart Job Name
+# LMDZ_Bands_file_name=${ARCHIVE}/IGCM_OUT/IPSLCM5/XXXXXXX/ATM/Debug/XXXXXXX_Bands_96x95x19_3prc.dat_3
+
+# 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=actuel
+
+# Set OutLevel to choose the level of output files from LMDZ. The file PARAM/output.def_$OutLevel will be used.
+# OutLevel=low/medium/high
+OutLevel=low
+
+# Set COSP activation and Outputs frequency (monthly, daily, HF) = y/n
+LMDZ_COSP_OK=n
+LMDZ_COSP_monthly=y
+LMDZ_COSP_daily=y
+LMDZ_COSP_hf=n
+
+# Set NMC Outputs frequency (monthly, daily, HF) = y/n
+LMDZ_NMC_monthly=n
+LMDZ_NMC_daily=n
+LMDZ_NMC_hf=n
+
+# Nugding : ok_guide=y/n
+ok_guide=n
+
+[InitialStateFiles]
+List=	(${ARCHIVE}/IGCM_OUT/LMDZ/${CREATE}/ATM/Output/Restart/${CREATE}_clim_start.nc,    start.nc),\
+        (${ARCHIVE}/IGCM_OUT/LMDZ/${CREATE}/ATM/Output/Restart/${CREATE}_clim_startphy.nc, startphy.nc) 
+
+[BoundaryFiles]
+List=()
+#List=
+#      (/dmnfs/cont003/p86ghatt/IGCM_OUT/LMDZOR/NUDGING-CORDEX-AFR/${year}${month}/u.nc, u.nc),\
+#      (/dmnfs/cont003/p86ghatt/IGCM_OUT/LMDZOR/NUDGING-CORDEX-AFR/${year}${month}/v.nc, v.nc)
+
+# With AR5 aerosols (including ASSO4M(=SO4), CSSO4M, AIBCM, ASBCM, AIPOMM, ASPOMM, ASSSM, CSSSM, SSSSM, CIDUSTM)
+# With AR5 ozone(tropo + strato from LMDZOR-INCA-REPROBUS simulations)
+ListNonDel= (${ARCHIVE}/IGCM_OUT/LMDZ/${CREATE}/ATM/Output/Boundary/${CREATE}_clim_limit.nc, limit.nc),\
+            (${ARCHIVE}/IGCM_OUT/LMDZ/${CREATE}/ATM/Output/Boundary/${CREATE}_climoz_LMDZ.nc, climoz_LMDZ.nc),\
+            (${R_BC}/ATM/${config_UserChoices_TagName}/${RESOL_ATM}/AR5/HISTORIQUE/aerosols_11YearsClim_1995_v3.nc, aerosols1980.nc),\
+            (${R_BC}/ATM/${config_UserChoices_TagName}/${RESOL_ATM}/AR5/HISTORIQUE/aerosols_11YearsClim_1860_v3.nc, aerosols.nat.nc)
+
+[SmoothFiles]
+# With AR5 aerosols (including ASSO4M(=SO4), CSSO4M, AIBCM, ASBCM, AIPOMM, ASPOMM, ASSSM, CSSSM, SSSSM, CIDUSTM)
+# With AR5 ozone(tropo + strato from LMDZOR-INCA-REPROBUS simulations)
+List=()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}_${LMDZ_Physics}, physiq.def), \
+	(${SUBMIT_DIR}/PARAM/gcm.def_${RESOL_ATM_3D},    gcm.def),    \
+	(${SUBMIT_DIR}/PARAM/cosp_input_nl.txt,                .),    \
+	(${SUBMIT_DIR}/PARAM/cosp_output_nl.txt,               .),    \
+	(${SUBMIT_DIR}/PARAM/output.def_${OutLevel},  output.def),    \
+	(${SUBMIT_DIR}/PARAM/config.def_${ConfType},  config.def),    \
+	(${SUBMIT_DIR}/PARAM/guide.def, .), \
+	(${SUBMIT_DIR}/PARAM/run.def, .), \
+	(${SUBMIT_DIR}/PARAM/traceur.def, .)
+
+[RestartFiles]
+List=	(restart.nc, restart.nc, start.nc), \
+	(restartphy.nc, restartphy.nc,  startphy.nc)
+
+[OutputText]
+List=   (physiq.def, gcm.def, run.def, traceur.def, guide.def, config.def, used_run.def)
+
+[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,  Post_1M_histmthCOSP), \
+	(histmthNMC.nc,   ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmthNMC.nc,   Post_1M_histmthNMC),  \
+	(dynzon.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,       Post_1M_dynzon),      \
+	(histday.nc,      ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,      Post_1D_histday),     \
+	(histdayCOSP.nc,  ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayCOSP.nc,  Post_1D_histdayCOSP), \
+	(histdayNMC.nc,   ${R_OUT_ATM_O_D}/${PREFIX}_1D_histdayNMC.nc,   Post_1D_histdayNMC),  \
+	(histhf.nc,       ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,       Post_HF_histhf),      \
+	(histhfCOSP.nc,   ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfCOSP.nc,   Post_HF_histhfCOSP),  \
+	(histhfNMC.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhfNMC.nc,    Post_HF_histhfNMC),   \
+	(histhf3h.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3h.nc,     Post_3H_histhf3h),    \
+        (histhf3hm.nc,    ${R_OUT_ATM_O_H}/${PREFIX}_3H_histhf3hm.nc,    Post_3H_histhf3hm)
+
+[Post_1M_histmth]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds, aire)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, ftime_th, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, LWupTOA, LWupTOAclr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_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, mc, nettop, paprs, 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, absvisaer, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl)
+ChunckJob2D = NONE
+TimeSeriesVars3D = (upwd, lwcon, iwcon, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc)
+ChunckJob3D = 100Y
+Seasonal=ON
+
+[Post_1M_histmthCOSP]
+Patches = ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds, aire)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE 
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
+Seasonal=OFF
+
+[Post_1M_histmthNMC]
+Patches = ()
+GatherWithInternal = (lon, lat, plev, time_counter, time_counter_bnds)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = (ta, zg, hus, hur, ua, va, wap, psbg, uv, vq, vT, wq, vphi, wT, uxu, vxv, TxT, ozone, ozone_daylight)
+ChunckJob3D = 100Y
+Seasonal=ON
+
+[Post_1M_dynzon]
+Patches = ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
+Seasonal=OFF
+
+[Post_1D_histday]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, iwp, lwp, phi500, precip, pluc, plul, prw, psol, rh2m, rh2m_min, rh2m_max, q2m, sens, sicf, slp, soll, sols, LWdnSFC, SWdnSFC, LWupSFCclr,SWupSFC, t2m, t2m_max, t2m_min, topl, tops, tsol, u10m, v10m, wind10m, wind10max, u850, v850, w850, phi850, q850, t850, u700, v700, w700, phi700, q700, t700, u500, v500, w500, phi500, q500, t500, u200, v200, w200, phi200, q200, t200)
+ChunckJob2D = 200Y
+TimeSeriesVars3D = (rhum, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres)
+ChunckJob3D = 10Y
+Seasonal=OFF
+
+[Post_1D_histdayCOSP]
+Patches = ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds, aire)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE 
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
+Seasonal=OFF
+
+[Post_1D_histdayNMC]
+Patches = ()
+GatherWithInternal = (lon, lat, plev, time_counter, time_counter_bnds)
+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_counter_bnds)
+TimeSeriesVars2D = (cldt, psol, q2m, slp, SWnetOR, SWdownOR, LWdownOR, precip, pluc, plul, t2m, tsol, u10m, v10m, t850, t700, t500, t200, t50, t10, u850, u700, u500, u200, u50, u10, v850, v700, v500, v200, v50, v10, phi850, phi700, phi200, phi50, phi10, w850, w700, w500, w200, w50, w10, q850, q700, q500, q200, q50, q10)
+ChunckJob2D = 50Y
+TimeSeriesVars3D = (temp, theta, ovap, vitu, vitv)
+ChunckJob3D = 10Y
+Seasonal=OFF
+
+[Post_HF_histhfCOSP]
+Patches = ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds, aire)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE 
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
+Seasonal=OFF
+
+[Post_HF_histhfNMC]
+Patches = ()
+GatherWithInternal = (lon, lat, plev, time_counter, time_counter_bnds)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = (ta, zg, hus, hur, ua, va, wap, psbg, uv, vq, vT, wq, vphi, wT, uxu, vxv, TxT)
+ChunckJob3D = 100Y
+Seasonal=OFF
+
+[Post_3H_histhf3h]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds)
+TimeSeriesVars2D = (phis, tsol, t2m, q2m, u10m, v10m, psol, topl, bils, cdrh)
+ChunckJob2D = 50Y
+TimeSeriesVars3D = ()
+ChunckJob3D = OFF
+Seasonal=OFF
+
+[Post_3H_histhf3hm]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds)
+TimeSeriesVars2D = (precip, flat, sens, LWdnSFC, LWupSFC, SWdnSFC, SWupSFC, pluc, snow, LWdnSFCclr, SWdnSFCclr, LWupSFCclr, cldt)
+ChunckJob2D = 50Y
+TimeSeriesVars3D = ()
+ChunckJob3D = OFF
+Seasonal=OFF
Index: CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/clim/config.card
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/clim/config.card	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/EXPERIMENTS/LMDZ/clim/config.card	(revision 1934)
@@ -0,0 +1,101 @@
+# $Id$
+# This is config.card file for restart with an NEW libIGCM simulation tree.
+#
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#D-- UserChoices -
+[UserChoices]
+#============================
+#-- (<8 chars MAX for JobName)
+JobName=tlmdz
+#----- Short Name of Experiment
+ExperimentName=clim
+#----- DEVT TEST PROD
+SpaceName=TEST
+LongName="LMDZ configuration"
+TagName=LMDZ
+#D- Choice of experiment in EXPERIEMENTS directory
+ExpType=LMDZ/clim
+#============================
+#-- leap, noleap, 360d
+CalendarType=360d
+#-- Begin and end of job 
+#-- "YYYY-MM-DD
+DateBegin=1980-01-01
+DateEnd=1980-01-30
+#============================
+#-- 1Y, 1M, 5D, 1D Period Length of one trunk of simulation
+PeriodLength=1M
+#============================
+#-- Total Number of Processors
+JobNumProcTot=32
+#============================
+# config.card configuration options :
+#R_INIT=/u/rech/psl/rpsl035/IGCM/INIT
+#R_INIT=/dmnfs/cont003/p86ipsl/IGCM/INIT
+#R_BC=/u/rech/psl/rpsl035/IGCM/BC
+#R_BC=/dmnfs/cont003/p86ipsl/IGCM/BC
+#============================
+
+#========================================================================
+#D-- Restarts -
+[Restarts]
+#D- by default: config.card describes no restart for all components 
+#D-             ie start from Levitus or limit files
+#D- If you want to restart all components from the same simulation,
+#D-      put OveRule flag to 'y' and set the next 3 parameters
+OverRule=n
+#D- Last day of the experience used as restart for all components
+RestartDate=2000-01-31
+#D- Define restart simulation name for all components
+RestartJobName=EXP00
+#D- Path Server Group Login
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZOR/DEVT/amip
+
+#========================================================================
+#D-- ATM -
+[ATM]
+#
+WriteFrequency="1M 1D HF"
+# 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=2000-01-31
+# Define restart simulation name for this component
+RestartJobName=LMDZOR01
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZOR/DEVT/amip
+# Old component name for restart (if empty, use new name)
+
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ4_AR5)
+
+#========================================================================
+#D-- Executable -
+[Executable]
+#D- For each component, Real name of executable, Name of executable for oasis
+ATM= (gcm.e, lmdz.x)
+
+#========================================================================
+#D-- Post -
+[Post]
+#D- Do we rebuild parallel output, this flag determines
+#D- frequency of rebuild submission (use NONE for DRYRUN=3)
+RebuildFrequency=1Y
+#D- Do we rebuild parallel output from archive (use NONE to use SCRATCHDIR as buffer)
+RebuildFromArchive=NONE
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission (NONE if you don't want)
+TimeSeriesFrequency=10Y
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average (NONE if you don't want)
+SeasonalFrequency=10Y
+#D- Offset for seasonal average first start dates ; same unit as SeasonalFrequency
+#D- Usefull if you do not want to consider the first X simulation's years
+SeasonalFrequencyOffset=0
+
+#========================================================================
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/DRIVER/lmdz.driver
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/DRIVER/lmdz.driver	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/DRIVER/lmdz.driver	(revision 1934)
@@ -0,0 +1,655 @@
+#!/bin/ksh
+#-----------------------------------------------------------------
+function LMDZ_sed
+{
+# Usage : LMDZ_sed filedef var_name myvalue
+#         In file filedef modify the line containing var_name=xxx into var_name=myvalue
+    IGCM_debug_PushStack "LMDZ_sed"
+    
+    # Test if the fichier exist
+    if [ ! -f ${1} ] ; then
+	echo "WARNING : ${1} file does not exist. Following will not be done : LMDZ_sed : ${1} ${2} ${3}"
+	return
+    fi
+
+    sed -e "s/^${2}\ *=.*/${2}= ${3}/" ${1} > ${1}.tmp
+    RET=$?
+    echo "LMDZ_sed : ${1} ${2} ${3}"
+    \mv ${1}.tmp ${1}
+
+    IGCM_debug_PopStack "LMDZ_sed"
+    return $RET
+}
+function LMDZ_sed_default
+{
+# Usage : LMDZ_sed_default filedef var_name
+#         In file filedef modify the line  "var_name=xxx DEFAULT= myvalue"
+#         into "var_name=myvalue"
+    IGCM_debug_PushStack "LMDZ_sed_default"
+
+    DefValue=$( grep ${2} ${1}  | awk -F"DEFAULT.=*" '{print $2}')
+
+    LMDZ_sed  ${1} ${2} $DefValue
+    IGCM_debug_PopStack "LMDZ_sed_default"
+    return $RET
+}
+
+
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' | awk "-F-" '{print $1}')
+
+    [ -f ${SUBMIT_DIR}/../.resol ] && eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x19
+
+    RESOL_ATM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' )
+    RESOL_ATM_X=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $1}' )
+    RESOL_ATM_Y=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $2}' )
+    RESOL_ATM_XY="${RESOL_ATM_X}x${RESOL_ATM_Y}"
+
+    ##-- Calendar type for LMDZ and create_etat0_limit
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=earth_366d
+	    CalendarTypeForCreate=gregorian;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d
+            CalendarTypeForCreate=${CalendarTypeForLmdz};;
+        360d)
+            CalendarTypeForLmdz=earth_360d
+            CalendarTypeForCreate=${CalendarTypeForLmdz};;
+        *)
+            CalendarTypeForLmdz=earth_360d
+            CalendarTypeForCreate=${CalendarTypeForLmdz}
+    esac
+
+    ##- LMDZ physics version
+    ##  Read LMDZ_Physics option in lmdz.card, if not present take default value AP (old physics) 
+    if [ ! X${lmdz_UserChoices_LMDZ_Physics} = X ] ; then
+	LMDZ_Physics=${lmdz_UserChoices_LMDZ_Physics}
+    else
+	LMDZ_Physics=AP
+    fi
+    echo LMDZ physics version : ${LMDZ_Physics}
+
+    ##- Create_etat0_limit version
+    ##  Define variable CREATE only if it is set in lmdz.card section UserChoices.
+    ##  This variable is only used in lmdz.card to choose input files.
+    if [ ! X${lmdz_UserChoices_CREATE} = X ] ; then
+	IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices CREATE
+	CREATE=${lmdz_UserChoices_CREATE}
+	echo create_etat0_limit version : ${CREATE}
+    fi
+
+
+    ## - LMDZ choice of config.def file
+    ##   ConfType must be set in lmdz.card
+    ConfType=${lmdz_UserChoices_ConfType}
+
+    ##- LMDZ output level, to be set in lmdz.card
+    if [ X${lmdz_UserChoices_OutLevel} = X ] ; then
+	# OutLevel is not set. Take default value low.
+	OutLevel=low
+    else
+	OutLevel=${lmdz_UserChoices_OutLevel}
+    fi
+
+    ##-- Output frequency purpose ....
+    ##--  Initialisation  ....
+    ok_mensuel=n
+    ok_journe=n
+    ok_hf=n
+    ok_hf3h=n
+    ok_hf3hm=n
+    ok_stn=n
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) ok_mensuel=y ;;
+	5D|5d|1D|1d) ok_journe=y ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    1M|1m) ok_mensuel=y ;;
+	esac
+	case ${frequency} in
+	    5D|5d|1D|1d) ok_journe=y ;;
+	esac
+	case ${frequency} in
+	    HF|hf) ok_hf=y ;;
+	esac
+	case ${frequency} in 
+	    HF3h|hf3h) ok_hf3h=y ;;
+	esac
+	case ${frequency} in 
+	    HF3hm|hf3hm) ok_hf3hm=y ;;
+	esac
+	case ${frequency} in
+	    STN|stn) ok_stn=y ;;
+	esac
+    done
+
+    ## Read LMDZ_NbPeriod_adjust option in lmdz.card 
+    if [ X"${lmdz_UserChoices_LMDZ_NbPeriod_adjust}" = X"" ] ; then
+	# The variable is not in lmdz.card, set default value
+	LMDZ_NbPeriod_adjust=0
+    else	
+	LMDZ_NbPeriod_adjust=${lmdz_UserChoices_LMDZ_NbPeriod_adjust}
+    fi
+
+    if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+	LMDZ_Bands_file_name=${lmdz_UserChoices_LMDZ_Bands_file_name}
+    fi
+
+
+    ##  Read LMDZ_COSP_OK in lmdz.card
+    if [ ! X${lmdz_UserChoices_LMDZ_COSP_OK} = X ] ; then
+	LMDZ_COSP_OK=${lmdz_UserChoices_LMDZ_COSP_OK}
+    else
+	LMDZ_COSP_OK=n
+    fi
+    ##  Read LMDZ_COSP_monthly in lmdz.card
+    if [ ! X${lmdz_UserChoices_LMDZ_COSP_monthly} = X ] ; then
+	LMDZ_COSP_monthly=${lmdz_UserChoices_LMDZ_COSP_monthly}
+    else
+	LMDZ_COSP_monthly=n
+    fi
+    ##  Read LMDZ_COSP_daily in lmdz.card
+    if [ ! X${lmdz_UserChoices_LMDZ_COSP_daily} = X ] ; then
+	LMDZ_COSP_daily=${lmdz_UserChoices_LMDZ_COSP_daily}
+    else
+	LMDZ_COSP_daily=n
+    fi
+    ##  Read LMDZ_COSP_hf in lmdz.card
+    if [ ! X${lmdz_UserChoices_LMDZ_COSP_hf} = X ] ; then
+	LMDZ_COSP_hf=${lmdz_UserChoices_LMDZ_COSP_hf}
+    else
+	LMDZ_COSP_hf=n
+    fi
+
+    ##  Read LMDZ_NMC_monthly in lmdz.card
+    if [ ! X${lmdz_UserChoices_LMDZ_NMC_monthly} = X ] ; then
+	LMDZ_NMC_monthly=${lmdz_UserChoices_LMDZ_NMC_monthly}
+    else
+	LMDZ_NMC_monthly=n
+    fi
+
+    ## Read LMDZ_NMC_daily in lmdz.card
+    if [ ! X${lmdz_UserChoices_LMDZ_NMC_daily} = X ] ; then
+	LMDZ_NMC_daily=${lmdz_UserChoices_LMDZ_NMC_daily}
+    else
+	LMDZ_NMC_daily=n
+    fi
+    
+    ## Read LMDZ_NMC_hf in lmdz.card
+    if [ ! X${lmdz_UserChoices_LMDZ_NMC_hf} = X ] ; then
+	LMDZ_NMC_hf=${lmdz_UserChoices_LMDZ_NMC_hf}
+    else
+	LMDZ_NMC_hf=n
+    fi
+
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+
+    case ${config_UserChoices_PeriodLength} in
+        *Y|*y) 
+               LMDZ_ecrit_ISCCP=30.
+               LMDZ_periodav=30.
+               if [ "${config_UserChoices_CalendarType}" != "360d" ] ; then
+                  echo Do not consider following warning if your running create_etat0_limit :
+                  echo WARNING !!! For lmdz : calendartype in config.card. PeriodLength=1Y allowed only for CalendarType=360d 
+	       fi
+               ;;
+        *)
+               LMDZ_ecrit_ISCCP=${PeriodLengthInDays}.
+               LMDZ_periodav=${PeriodLengthInDays}.
+               ;;
+    esac
+  
+    ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def
+    if [ ${CumulPeriod} -eq 1 ] ; then
+	RAZ_DATE=1
+    else
+	RAZ_DATE=0
+    fi
+
+    ## Algorithme for special treatment for Bands_xxx file
+    ## For CumulPeriod=1 ; 
+    ##        IF NbPeriod_Adjust = 0 ; set LMDZ_adjust=n
+    ##                                IF LMDZ_Bands_file_name was given in lmdz.card THEN Get Bands file directly from server. Store it later with _0 suffix.
+    ##                                IF no LMDZ_Bands_file_name was given, start without Bands file. Store it later with _0 suffix.
+    ##        ELSE set LMDZ_adjust=y ; start without Bands file ; Store it later in PARAM/ directory in submit directory with suffix _1 ;
+    ##
+    ## For CumulPeriod=2 to LMDZ_NbPeriod_adjust ; 
+    ##       Get Bands file from PARAM/ in submit directory (CumulPeriod-1) ; set LMDZ_adjust=y ; Store Bands file in PARAM/ ;
+    ##
+    ## For CumulPeriod > LMDZ_NbPeriod_adjust ; LMDZ_adjust=n ; Get Bands file from PARAM/ ; Do not store ;
+
+    if [ ${CumulPeriod} -eq 1 ] ; then
+        if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+	    if [ ! X${LMDZ_Bands_file_name} = X ] ; then
+		IGCM_sys_Get ${LMDZ_Bands_file_name} Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat 
+		IGCM_sys_Chmod u+w Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat
+		IGCM_debug_Print 1 "Bands file forced to ${LMDZ_Bands_file_name} and stored in ${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_0"
+	    fi
+        fi
+    fi
+
+    LMDZ_adjust=n
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && LMDZ_adjust=y
+
+    RefPeriod=${LMDZ_NbPeriod_adjust} #0 to use a Bands file from Restart
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && (( RefPeriod = ${CumulPeriod} - 1 )) 
+
+    # Get Bands file from PARAM directory if file exist
+    if [ ${CumulPeriod} -gt 1 && -f ${SUBMIT_DIR}/PARAM/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_${RefPeriod} ] ; then
+        IGCM_sys_Cp ${SUBMIT_DIR}/PARAM/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat
+        IGCM_sys_Chmod u+w Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat
+    fi
+
+    ##-- GHG forcing :
+    ##   If forcing file exist in run directory, read values for the current year and modify config.def
+    if [ -f SOLARANDVOLCANOES.txt ] ; then
+	IPCC_SOLAR=`grep Annee_${year} SOLARANDVOLCANOES.txt | awk -F= '{print $2}'`
+	LMDZ_sed config.def solaire     ${IPCC_SOLAR}
+    else
+	LMDZ_sed_default config.def solaire
+    fi
+    if [ -f CO2.txt ] ; then
+	IPCC_CO2=`grep Annee_${year} CO2.txt | awk -F= '{print $2}'`
+	LMDZ_sed config.def co2_ppm     ${IPCC_CO2}
+    else
+	LMDZ_sed_default config.def co2_ppm
+    fi
+
+
+    if [ -f CH4.txt ] ; then
+	IPCC_CH4=`grep Annee_${year} CH4.txt | awk -F= '{print $2}'`
+	LMDZ_sed config.def CH4_ppb     ${IPCC_CH4}
+    else
+	LMDZ_sed_default config.def CH4_ppb
+    fi
+    if [ -f N2O.txt ] ; then
+	IPCC_N2O=`grep Annee_${year} N2O.txt | awk -F= '{print $2}'`
+	LMDZ_sed config.def N2O_ppb     ${IPCC_N2O}
+    else
+	LMDZ_sed_default config.def N2O_ppb
+    fi
+    if [ -f CFC11.txt ] ; then
+	IPCC_CFC11=`grep Annee_${year} CFC11.txt | awk -F= '{print $2}'`
+	LMDZ_sed config.def CFC11_ppt   ${IPCC_CFC11}
+    else
+	LMDZ_sed_default config.def CFC11_ppt
+    fi
+    if [ -f CFC12.txt ] ; then
+	IPCC_CFC12=`grep Annee_${year} CFC12.txt | awk -F= '{print $2}'`
+	LMDZ_sed config.def CFC12_ppt   ${IPCC_CFC12}
+    else
+	LMDZ_sed_default config.def CFC12_ppt
+    fi
+
+
+    ##-- Add special treatement for CARBON CYCLE
+    if [ X"${lmdz_UserChoices_CARBON_CYCLE}" = X"control" ] || [ X"${lmdz_UserChoices_CARBON_CYCLE}" = X"historical" ] ; then
+	ATM_Carbon_Update
+    fi
+
+
+    ##-- Set LMDZ_COSP_daily1979=y in section UserChoices in lmdz.card to activate 
+    ##   COSP daily output for years 1979 and later (mandatory for historical CMIP5 run).
+    if ([ X"${lmdz_UserChoices_LMDZ_COSP_daily1979}" = X"y" ] &&  [ ${year} -ge 1979 ]) ; then
+	LMDZ_COSP_daily=y
+    fi
+
+
+    ## output.def parameters
+    #  columns in phys_output_filekeys refer to the filenames given in phys_out_filenames in output.def :  
+    #                                            histmth       histday        histhf    histfh3h   histhf3hm   histstn
+    LMDZ_sed output.def phys_out_filekeys       "${ok_mensuel} ${ok_journe}   ${ok_hf}  ${ok_hf3h} ${ok_hf3hm} ${ok_stn}"
+    LMDZ_sed output.def ecrit_ISCCP ${LMDZ_ecrit_ISCCP}
+    LMDZ_sed output.def ok_cosp     ${LMDZ_COSP_OK}
+    LMDZ_sed output.def ok_mensuelCOSP  ${LMDZ_COSP_monthly}
+    LMDZ_sed output.def ok_journeCOSP   ${LMDZ_COSP_daily}
+    LMDZ_sed output.def ok_hfCOSP   ${LMDZ_COSP_hf}
+    LMDZ_sed output.def ok_histNMC  "${LMDZ_NMC_monthly} ${LMDZ_NMC_daily} ${LMDZ_NMC_hf}"
+
+    ## gcm.def parameters : 
+    # Modification only for new physics
+    if [ X${LMDZ_Physics} = X"NPv3.1" ] ; then
+	LMDZ_sed gcm.def iphysiq     5
+    else
+	LMDZ_sed_default gcm.def iphysiq
+    fi
+
+
+    ## run.def parameters
+    LMDZ_sed run.def dayref   ${InitDay}
+    LMDZ_sed run.def nday     ${PeriodLengthInDays}
+    LMDZ_sed run.def raz_date ${RAZ_DATE}
+    LMDZ_sed run.def periodav ${LMDZ_periodav}
+    LMDZ_sed run.def adjust   ${LMDZ_adjust}
+
+    # Set anneeref different for gcm and ce0l
+    # Variable calend is different for a run with create_etat0_limit and the gcm
+    # Test if executable create_etat0_limit is present and lmdz.x is not present
+    if [ -f create_etat0_limit* ] || [ -f ce0l* ] && [ ! -f lmdz.x ] ; then
+	# for case ce0l : always take current year
+	LMDZ_sed run.def anneeref ${year}
+	LMDZ_sed run.def calend   ${CalendarTypeForCreate}
+    else
+	# for case gcm : take first year of simulation
+	LMDZ_sed run.def anneeref ${InitYear}
+	LMDZ_sed run.def calend   ${CalendarTypeForLmdz}
+    fi
+
+    # Activate creation of file grilles_gcm.nc only at first period
+    if [ ${CumulPeriod} -eq 1 ] ; then
+	LMDZ_sed run.def grilles_gcm_netcdf y
+    else
+	LMDZ_sed run.def grilles_gcm_netcdf n
+    fi
+
+    ## Determine from the variable ListOfComponents in config.card coupling to external models 
+    ## and set corresponding parameters in run.def
+    echo ListOfComponents now running : ${config_ListOfComponents[*]}
+
+    if [ X${config_ListOfComponents_SRF} = Xorchidee ] ; then
+	echo "Activate ORCHIDEE, set VEGET=y in run.def"
+	LMDZ_sed run.def VEGET y
+    else
+	echo "No ORCHIDEE, set VEGET=n in run.def"
+	LMDZ_sed run.def VEGET n
+    fi
+
+    if [ X${config_ListOfComponents_CPL} = Xoasis ] ; then
+	echo "Activate coupling to ocean, set type_ocean=couple in run.def"
+	LMDZ_sed run.def type_ocean couple
+    else
+	echo "LMDZ is running in forced mode without ocean model, set type_ocean=force in run.def"
+	LMDZ_sed run.def type_ocean force
+    fi
+
+    if [ X${config_ListOfComponents_CHM} = Xinca ] ; then
+	echo "Activate coupling to INCA, set type_trac=inca in run.def"
+	LMDZ_sed run.def type_trac inca
+    elif [ X${config_ListOfComponents_CHM} = Xreprobus ] ; then
+	echo "Activate coupling to REPROBUS, set type_trac=repr in run.def"
+	LMDZ_sed run.def type_trac repr
+        LMDZ_sed run.def config_inca none
+    else
+	echo "No coupling to chemistry model, set type_trac=lmdz in run.def"
+	LMDZ_sed run.def type_trac lmdz
+        LMDZ_sed run.def config_inca none
+    fi
+
+
+    # config.def : If running with Reprobus, desactivate ok_ade, ok_aie and read_climoz
+    if [ X${config_ListOfComponents_CHM} = Xreprobus ]  ; then
+	LMDZ_sed config.def ok_ade n
+	LMDZ_sed config.def ok_aie n
+	LMDZ_sed config.def aerosol_couple n
+	LMDZ_sed config.def read_climoz 0
+    elif [ X${config_ListOfComponents_CHM} = Xinca ]; then 
+	case ${CHEM} in 
+	    AER)
+		LMDZ_sed config.def ok_ade y
+		LMDZ_sed config.def ok_aie n
+		LMDZ_sed config.def aerosol_couple y
+		LMDZ_sed config.def read_climoz 0
+		;;
+	    NMHC_AER)
+		LMDZ_sed config.def ok_ade y
+		LMDZ_sed config.def ok_aie n
+		LMDZ_sed config.def aerosol_couple y
+		LMDZ_sed config.def read_climoz 0
+		;;
+	    GES)
+		LMDZ_sed config.def ok_ade n
+		LMDZ_sed config.def ok_aie n
+		LMDZ_sed config.def aerosol_couple n
+		LMDZ_sed config.def read_climoz 0
+	esac
+
+    else
+	LMDZ_sed config.def ok_ade y
+	LMDZ_sed config.def ok_aie y
+	LMDZ_sed config.def aerosol_couple n
+	LMDZ_sed config.def read_climoz 2
+    fi
+
+    # guide.def : Activate nudging if ok_guide=y set in lmdz.card
+    if [ ! X${lmdz_UserChoices_ok_guide} = X ] ; then
+	LMDZ_sed guide.def ok_guide  ${lmdz_UserChoices_ok_guide}
+    else
+	LMDZ_sed guide.def ok_guide n
+    fi
+
+    ## Read ByPass_hgardfou_teta option in lmdz.card --> divide teta* by 2 if [ $ByPass_hgardfou_teta = 1 ]
+    ByPass_hgardfou_teta=${lmdz_UserChoices_ByPass_hgardfou_teta}
+    if [ X"${ByPass_hgardfou_teta}" = X"y" ] ; then
+	awk '{ if ($0 ~ /^teta.*=/) {split($0,a,"=") ; print a[1]"="a[2]/2"."} else print $0}' gcm.def > gcm.def.tmp
+	IGCM_sys_Mv gcm.def.tmp gcm.def
+	echo
+	IGCM_debug_Print 1 "ByPass_hgardfou_teta : ^teta*/2 in gcm.def"
+	echo
+	cat gcm.def
+	ByPass_hgardfou_teta=n
+	IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices ByPass_hgardfou_teta "${ByPass_hgardfou_teta}"
+    fi
+
+    ## Read ByPass_hgardfou_mats option in lmdz.card --> purmats=y 2 if [ $ByPass_hgardfou_mats = 1 ]
+    ByPass_hgardfou_mats=${lmdz_UserChoices_ByPass_hgardfou_mats}
+    if [ X"${ByPass_hgardfou_mats}" = X"y" ] ; then
+	sed -e "s/^purmats=.*/purmats=y/" gcm.def > gcm.def.tmp
+        IGCM_sys_Mv gcm.def.tmp gcm.def
+        echo
+        IGCM_debug_Print 1 "ByPass_hgardfou_mats : purmats=y in gcm.def"
+        echo
+        cat gcm.def
+        ByPass_hgardfou_mats=n
+        IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices ByPass_hgardfou_mats "${ByPass_hgardfou_mats}"
+    fi
+
+    if ( ${FirstInitialize} ) ; then
+	if [ -f start.nc ] && [ -f startphy.nc ] ; then
+	    # The restart files already exist. They have been copied from a previous run or from a CREATE job
+	    echo "Restart files for LMDZ have been copied from a previous run"
+	elif [ -f create_etat0_limit* ] || [ -f ce0l* ] && [ ! -f lmdz.x ] ; then
+	    # Do nothing because it is the program for creating restart files that is running and not LMDZ
+	    echo "create_etat0_limit.e will create initial start files"
+	else
+	    # Restart files have to be created
+	    echo "Run create_etat0_limit.e to create restart files start.nc and startphy.nc for LMDZ"
+	    if [ ${DRYRUN} -le 2 ] ; then
+                ##-- suppression of ozone file construction --
+		cp config.def config.def.save
+		LMDZ_sed config.def read_climoz 0
+		IGCM_sys_Cp ${R_EXE}/create_etat0_limit.e ${RUN_DIR}/.
+		if [ ${DRYRUN} -le 1 ] ; then
+		    ${HOST_MPIRUN_COMMAND} -np 1 ./create_etat0_limit.e
+		    IGCM_sys_Put_Out start.nc    ${R_OUT_ATM_R}/${config_UserChoices_JobName}_${PeriodDateBegin}_start.nc
+		    IGCM_sys_Put_Out startphy.nc ${R_OUT_ATM_R}/${config_UserChoices_JobName}_${PeriodDateBegin}_startphy.nc
+		fi
+		IGCM_sys_Mv config.def.save config.def
+	    else
+		echo "EXECUTION of ${HOST_MPIRUN_COMMAND} -np 1 ./create_etat0_limit.e simulated"
+		echo "EXECUTION of ${HOST_MPIRUN_COMMAND} -np 1 ./create_etat0_limit.e simulated for DRYRUN = " $DRYRUN >> stack
+	    fi
+	fi
+    fi
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    # If the file exist, copy Bands_ file to PARAM/ in submit directory. This file will be used for the simulation.
+    # Copy also the same file into ATM/Restart at ARCHIVE directory for backup.
+    if [ -f Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ] ; then 
+      if [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] ; then
+        IGCM_sys_Cp Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${SUBMIT_DIR}/PARAM/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+        IGCM_sys_Put_Out Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_R}/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+      elif [ ${CumulPeriod} -eq 1 ] && [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+        # Special case : first period and no adjust => Save bands file with suffix _0
+        IGCM_sys_Cp Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${SUBMIT_DIR}/PARAM/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_0
+        IGCM_sys_Put_Out Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_R}/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_0
+      fi
+    fi
+
+    # Add special treatement for CARBON CYCLE
+    if [ X"${lmdz_UserChoices_CARBON_CYCLE}" = X"control" ] || [ X"${lmdz_UserChoices_CARBON_CYCLE}" = X"historical" ] ; then
+	ATM_Carbon_Finalize
+    fi
+
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
+
+
+
+function ATM_Carbon_Update
+{
+    # This fuction will be called only if CARBON_CYCLE is set to control or historical in lmdz.card
+
+    IGCM_debug_PushStack "ATM_Carbon_Update"
+
+    typeset InitPeriodCo2 InitDateBeginCo2 InitDateEndCo2
+    typeset CO2SBG CO2MBG CO2LU CO2_ppm CO2_ppm_prec fCO2_ff
+    typeset PathCO2log 
+    typeset PreviousLinelog LastPeriodDateBegin LastPeriodDateEnd LastDatesPeriod LastPREFIX
+    
+    if ( ${FirstInitialize} && [ X"${config_ATM_Restart}" = X"n" ] ) ; then
+        ##--Initialization of fluxes to an undefined value at the first run
+	UndefinedValueCo2=-9999
+	InitPeriodCo2=0
+	InitDateBeginCo2=${UndefinedValueCo2}
+	InitYearCo2=$(( ${year} - 1 )) 
+	InitDateEndCo2=${InitYearCo2}1231
+	fCO2_ff=${UndefinedValueCo2} 
+	CO2SBG=${UndefinedValueCo2}
+	CO2LU=${UndefinedValueCo2}
+	CO2MBG=${UndefinedValueCo2}
+	CO2_ppm_prec=${UndefinedValueCo2}
+        # Get the initial value of atmosp. pco2 
+	CO2_ppm=${lmdz_UserChoices_co2_init}
+	
+	IGCM_sys_Cp ${SUBMIT_DIR}/co2.log.init ${SUBMIT_DIR}/co2.log
+	IGCM_sys_Chmod u+w ${SUBMIT_DIR}/co2.log
+	
+        # Save CO2 values in ExeCpuLog variable contents 5 fields
+	echo "${InitPeriodCo2} ${InitDateBeginCo2} ${InitDateEndCo2} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SBG} ${CO2LU}" |   \
+	    gawk '{printf("  %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \
+                              $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log
+	
+        # Get the path where the log file co2.log is
+	PathCO2log=${SUBMIT_DIR} 
+        # Get Previous line in co2.log
+	PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log )
+	
+    else
+	
+	# NEW parameter in config.card (even for Overule all restarts, you must modify the line after this comment ):
+        # !! NEW for Carbone cycle !! Path where to find old co2.log file for restart CO2.
+	#SubmitRestartPath=
+	if ( ${FirstInitialize} && [ X"${config_ATM_Restart}" = X"y" ] ) ; then
+	    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/config.card ATM SubmitRestartPath
+	    PathCO2log=${config_ATM_SubmitRestartPath}
+	    Date_r=$( IGCM_date_ConvertFormatToGregorian ${config_ATM_RestartDate} )
+	    PreviousLinelog=$( grep ${Date_r} ${PathCO2log}/co2.log )
+	    
+	    IGCM_sys_Cp ${SUBMIT_DIR}/co2.log.init ${SUBMIT_DIR}/co2.log
+	    IGCM_sys_Chmod u+w ${SUBMIT_DIR}/co2.log
+	    
+            # Save Last Line of control in new co2.log file
+	    InitPeriodCo2=0   
+	    InitDateBeginCo2=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $2}' )
+	    InitDateEndCo2=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $3}' )
+	    fCO2_ff=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $6}' )
+	    CO2SBG=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $8}' )
+	    CO2LU=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $9}' )
+	    CO2MBG=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $7}' )
+	    CO2_ppm_prec=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $5}' )
+	    CO2_ppm=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' )
+	    
+	    echo "${InitPeriodCo2} ${InitDateBeginCo2} ${InitDateEndCo2} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SBG} ${CO2LU}" |   \
+		gawk '{printf("  %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \
+                              $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log
+	else
+           # Get restart line in co2.log
+	    PathCO2log=${SUBMIT_DIR}
+	    PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log )
+           # Get the value of atmosp. pco2 in co2.log
+	    CO2_ppm=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' )
+	fi
+	
+    fi
+    
+    if [ ${lmdz_UserChoices_CARBON_CYCLE} = historical ] ; then
+	    # Modify co2_ppm in config.def only for historical run
+	    # If not historical run, keep co2_ppm value set earlier by LMDZ_sed
+	LMDZ_sed config.def co2_ppm     ${CO2_ppm}
+    fi
+    IGCM_debug_PopStack "ATM_Carbon_Update"
+}
+
+
+function ATM_Carbon_Finalize
+{
+    # This fuction will be called only if CARBON_CYCLE is set to control or historical in lmdz.card
+
+    IGCM_debug_PushStack "ATM_Carbon_Finalize"
+
+    typeset CO2SBG_3V CO2SBG CO2MBG CO2SRF CO2LU CO2_ppm CO2_ppm_prec fCO2_ff
+    
+	# Fossil fuel emission
+    if [ ${lmdz_UserChoices_CARBON_CYCLE} = control ] ; then
+            # For Control, no fossil fuel emission :
+	fCO2_ff=0
+	echo "Fossil Fuel for control run  :" ${fCO2_ff}
+    else
+            # For historical, get the fossile fuel file :
+	fCO2_ff=$( gawk "{if (match(\$0,\"${year}\ *${month}\")) {print \$3}}" ${SUBMIT_DIR}/PARAM/CMIP5_gridcar_CO2_emissions_fossil_fuel_Andres_1751-2007_monthly_SC.txt )
+	echo "Fossil Fuel for historical run  :" ${fCO2_ff}
+    fi
+    
+        # Get the value of ocean carbon flux
+    CO2MBG=$( ${SUBMIT_DIR}/COMP/lmdz_analyse_pisces_out.awk ocean.output )
+    CO2MBG=$( echo ${CO2MBG} | sed -e "s/ *//" )
+    echo "Ocean carbon flux  :" ${CO2MBG}
+    
+        # Get the value of land fluxes
+    set +A CO2SBG_3V -- $( ${SUBMIT_DIR}/COMP/lmdz_analyse_stomate_out.awk out_orchidee_0000 )
+    CO2SRF=$( echo ${CO2SBG_3V[0]} | sed -e "s/ *//" )
+    CO2LU=$( echo ${CO2SBG_3V[1]} | sed -e "s/ *//" )
+    CO2SBG=$( echo ${CO2SBG_3V[2]} | sed -e "s/ *//" )
+    echo "Land carbon flux (NEP)         :" ${CO2SRF}
+    echo "LU flux          (FLUC)        :" ${CO2LU}
+    echo "Total Land carbon flux (NBP)   :" ${CO2SBG}
+    
+        # Previous CO2
+    CO2_ppm_prec=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' )
+    echo "Previous CO2  :" ${CO2_ppm_prec}
+    
+        # Formula
+    CO2_ppm=$( echo "${CO2_ppm_prec} + (${fCO2_ff} + ${CO2MBG} + ${CO2SBG}) / 2.12" | bc -l )
+    echo "New CO2  :" ${CO2_ppm}
+    
+        # Save CO2 values in ExeCpuLog variable contents 5 fields
+    echo "${CumulPeriod} ${PeriodDateBegin} ${PeriodDateEnd} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SBG} ${CO2LU}" |   \
+	gawk '{printf("  %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \
+                           $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log
+    
+    IGCM_debug_PopStack "ATM_Carbon_Finalize"
+}
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_L19
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_L19	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_L19	(revision 1934)
@@ -0,0 +1,60 @@
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+
+## Nombre  d'appels des routines de rayonnements ( par jour)                 
+nbapp_rad=12
+
+##  Facteur additif pour l'albedo
+pmagic=0.008
+
+#
+# Parametres "orbitaux/ ere geologique"
+#
+### R_ecc =      Excentricite
+R_ecc = 0.016715
+### R_peri =     Equinoxe
+R_peri = 102.7
+### R_incl =     Inclinaison
+R_incl = 23.441
+### solaire =    Constante solaire
+solaire = _AUTO_: DEFAULT = 1361.
+#
+# Taux gaz a effet de serre
+#
+### co2_ppm =    taux CO2 en ppm
+co2_ppm = _AUTO_: DEFAULT = 367.
+### CH4_ppb =    taux CH4 en ppb
+CH4_ppb = _AUTO_: DEFAULT = 1760.
+### N2O_ppb =    taux N2O en ppb
+N2O_ppb = _AUTO_: DEFAULT = 316.
+### CFC11_ppt =  taux CFC11 en ppt
+CFC11_ppt = _AUTO_: DEFAULT = 741.2
+### CFC12_ppt =  taux CFC12 en ppt
+CFC12_ppt = _AUTO_: DEFAULT = 535.
+#
+# Parametres effets directs/indirects des "aerosols"
+#
+### ok_ade=y/n   flag Aerosol direct effect
+ok_ade=_AUTO_
+### ok_aie=y/n   flag Aerosol indirect effect
+ok_aie=_AUTO_
+### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
+aer_type=actuel
+### aerosol coupling = online (y) or offline (n)
+aerosol_couple=_AUTO_
+###  If aerosol offline : type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
+flag_aerosol=1
+### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b0=1.7
+### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b1=0.2
+#
+# Parametre de lecture de l'ozone
+#
+# Allowed values are 0, 1 and 2
+# 0: do not read an ozone climatology
+# 1: read a single ozone climatology that will be used day and night
+# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
+read_climoz=_AUTO_
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_actuel
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_actuel	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_actuel	(revision 1934)
@@ -0,0 +1,60 @@
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+
+## Nombre  d'appels des routines de rayonnements ( par jour)                 
+nbapp_rad=24
+
+##  Facteur additif pour l'albedo
+pmagic=0.008
+
+#
+# Parametres "orbitaux/ ere geologique"
+#
+### R_ecc =      Excentricite
+R_ecc = 0.016715
+### R_peri =     Equinoxe
+R_peri = 102.7
+### R_incl =     Inclinaison
+R_incl = 23.441
+### solaire =    Constante solaire
+solaire = _AUTO_ : DEFAULT = 1366.0896
+#
+# Taux gaz a effet de serre
+#
+### co2_ppm =    taux CO2 en ppm
+co2_ppm = _AUTO_: DEFAULT = 0.36886500E+03
+### CH4_ppb =    taux CH4 en ppb
+CH4_ppb = _AUTO_ : DEFAULT = 0.17510225E+04
+### N2O_ppb =    taux N2O en ppb
+N2O_ppb = _AUTO_ : DEFAULT = 0.31585000E+03
+### CFC11_ppt =  taux CFC11 en ppt
+CFC11_ppt = _AUTO_ : DEFAULT = 5.18015181E+01
+### CFC12_ppt =  taux CFC12 en ppt
+CFC12_ppt = _AUTO_ : DEFAULT = 0.99862742E+03
+#
+# Parametres effets directs/indirects des "aerosols"
+#
+### ok_ade=y/n   flag Aerosol direct effect
+ok_ade=_AUTO_
+### ok_aie=y/n   flag Aerosol indirect effect
+ok_aie=_AUTO_
+### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
+aer_type=actuel
+### aerosol coupling = online (y) or offline (n)
+aerosol_couple=_AUTO_
+###  If aerosol offline : type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
+flag_aerosol=6 
+### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b0=1.7
+### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b1=0.2
+#
+# Parametre de lecture de l'ozone
+#
+# Allowed values are 0, 1 and 2
+# 0: do not read an ozone climatology
+# 1: read a single ozone climatology that will be used day and night
+# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
+read_climoz=_AUTO_
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_annuel
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_annuel	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_annuel	(revision 1934)
@@ -0,0 +1,60 @@
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+
+## Nombre  d'appels des routines de rayonnements ( par jour)                 
+nbapp_rad=24
+
+##  Facteur additif pour l'albedo
+pmagic=0.0
+
+#
+# Parametres "orbitaux/ ere geologique"
+#
+### R_ecc =      Excentricite
+R_ecc = 0.016715
+### R_peri =     Equinoxe
+R_peri = 102.7
+### R_incl =     Inclinaison
+R_incl = 23.441
+### solaire =    Constante solaire
+solaire = _AUTO_ 
+#
+# Taux gaz a effet de serre
+#
+### co2_ppm =    taux CO2 en ppm
+co2_ppm = _AUTO_
+### CH4_ppb =    taux CH4 en ppb
+CH4_ppb = _AUTO_
+### N2O_ppb =    taux N2O en ppb
+N2O_ppb = _AUTO_
+### CFC11_ppt =  taux CFC11 en ppt
+CFC11_ppt = _AUTO_
+### CFC12_ppt =  taux CFC12 en ppt
+CFC12_ppt = _AUTO_
+#
+# Parametres effets directs/indirects des "aerosols"
+#
+### ok_ade=y/n   flag Aerosol direct effect
+ok_ade=_AUTO_
+### ok_aie=y/n   flag Aerosol indirect effect
+ok_aie=_AUTO_
+### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
+aer_type=annuel
+### aerosol coupling = online (y) or offline (n)
+aerosol_couple=_AUTO_
+###  If aerosol offline : type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
+flag_aerosol=6 
+### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b0=1.7
+### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b1=0.2
+#
+# Parametre de lecture de l'ozone
+#
+# Allowed values are 0, 1 and 2
+# 0: do not read an ozone climatology
+# 1: read a single ozone climatology that will be used day and night
+# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
+read_climoz=_AUTO_
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_preind
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_preind	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/config.def_preind	(revision 1934)
@@ -0,0 +1,60 @@
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+
+## Nombre  d'appels des routines de rayonnements ( par jour)                 
+nbapp_rad=24
+
+##  Facteur additif pour l'albedo
+pmagic=0.0
+
+#
+# Parametres "orbitaux/ ere geologique"
+#
+### R_ecc =      Excentricite
+R_ecc = 0.016715
+### R_peri =     Equinoxe
+R_peri = 102.7
+### R_incl =     Inclinaison
+R_incl = 23.441
+### solaire =    Constante solaire
+solaire = _AUTO_: DEFAULT = 1365.6537
+#
+# Taux gaz a effet de serre
+#
+### co2_ppm =    taux CO2 en ppm
+co2_ppm = _AUTO_: DEFAULT = 0.28472500E+03
+### CH4_ppb =    taux CH4 en ppb
+CH4_ppb = _AUTO_: DEFAULT = 0.79097924E+03
+### N2O_ppb =    taux N2O en ppb
+N2O_ppb = _AUTO_: DEFAULT = 0.27542506E+03
+### CFC11_ppt =  taux CFC11 en ppt
+CFC11_ppt = _AUTO_: DEFAULT = 0.
+### CFC12_ppt =  taux CFC12 en ppt
+CFC12_ppt = _AUTO_: DEFAULT = 0.
+#
+# Parametres effets directs/indirects des "aerosols"
+#
+### ok_ade=y/n   flag Aerosol direct effect
+ok_ade=_AUTO_
+### ok_aie=y/n   flag Aerosol indirect effect
+ok_aie=_AUTO_
+### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
+aer_type=preind
+### aerosol coupling = online (y) or offline (n)
+aerosol_couple=_AUTO_
+###  If aerosol offline : type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
+flag_aerosol=6 
+### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b0=1.7
+### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b1=0.2
+#
+# Parametre de lecture de l'ozone
+#
+# Allowed values are 0, 1 and 2
+# 0: do not read an ozone climatology
+# 1: read a single ozone climatology that will be used day and night
+# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
+read_climoz=_AUTO_
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/cosp_input_nl.txt
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/cosp_input_nl.txt	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/cosp_input_nl.txt	(revision 1934)
@@ -0,0 +1,103 @@
+! (c) British Crown Copyright 2008, the Met Office.
+! All rights reserved.
+! 
+! Redistribution and use in source and binary forms, with or without modification, are permitted 
+! provided that the following conditions are met:
+! 
+!     * Redistributions of source code must retain the above copyright notice, this list 
+!       of conditions and the following disclaimer.
+!     * Redistributions in binary form must reproduce the above copyright notice, this list
+!       of conditions and the following disclaimer in the documentation and/or other materials 
+!       provided with the distribution.
+!     * Neither the name of the Met Office nor the names of its contributors may be used 
+!       to endorse or promote products derived from this software without specific prior written 
+!       permission.
+! 
+! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
+! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
+! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR 
+! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 
+! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, 
+! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER 
+! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT 
+! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+
+! Namelist that sets up the main COSP options
+&COSP_INPUT
+  CMOR_NL='./cmor/cosp_cmor_nl.txt', ! CMOR namelist
+  NPOINTS_IT=10000,! Max number of gridpoints to be processed in one iteration
+  NCOLUMNS=20,  ! Number of subcolumns
+  NLEVELS=39,   ! Number of model levels 
+  USE_VGRID=.true., ! Use fixed vertical grid for outputs? (if .true. then you need to define number of levels with Nlr)
+  NLR=40,       ! Number of levels in statistical outputs (only used if USE_VGRID=.true.)
+  CSAT_VGRID=.true., ! CloudSat vertical grid? (if .true. then the CloudSat standard grid is used for the outputs.
+                     !  USE_VGRID needs also be .true.)
+  FINPUT='histday.nc', ! NetCDF file with 1D inputs
+!  FINPUT='cosp_input_um_2d.nc', ! NetCDF file with 2D inputs
+  !----------------------------------------------------------------------------------
+  !--------------- Inputs related to radar simulations
+  !----------------------------------------------------------------------------------
+  RADAR_FREQ=94.0, ! CloudSat radar frequency (GHz)
+  SURFACE_RADAR=0, ! surface=1, spaceborne=0
+  use_mie_tables=0,! use a precomputed lookup table? yes=1,no=0
+  use_gas_abs=1,   ! include gaseous absorption? yes=1,no=0
+  do_ray=0,        ! calculate/output Rayleigh refl=1, not=0
+  melt_lay=0,      ! melting layer model off=0, on=1
+  k2=-1,           ! |K|^2, -1=use frequency dependent default
+  use_reff=.true., ! True if you want effective radius to be used by radar simulator (always used by lidar)
+  use_precipitation_fluxes=.true.,  ! True if precipitation fluxes are input to the algorithm 
+  !----------------------------------------------------------------------------------
+  !---------------- Inputs related to lidar simulations
+  !----------------------------------------------------------------------------------
+  Nprmts_max_hydro=12, ! Max number of parameters for hydrometeor size distributions
+  Naero=1,             ! Number of aerosol species (Not used)
+  Nprmts_max_aero=1,   ! Max number of parameters for aerosol size distributions (Not used)
+  lidar_ice_type=0,    ! Ice particle shape in lidar calculations (0=ice-spheres ; 1=ice-non-spherical)
+  OVERLAP=3,   !  overlap type: 1=max, 2=rand, 3=max/rand
+  !----------------------------------------------------------------------------------
+  !---------------- Inputs related to ISCCP simulator
+  !----------------------------------------------------------------------------------
+  ISCCP_TOPHEIGHT=1,  !  1 = adjust top height using both a computed
+                       !  infrared brightness temperature and the visible
+                       !  optical depth to adjust cloud top pressure. Note
+                       !  that this calculation is most appropriate to compare
+                       !  to ISCCP data during sunlit hours.
+                      !  2 = do not adjust top height, that is cloud top
+                       !  pressure is the actual cloud top pressure
+                       !  in the model
+                      !  3 = adjust top height using only the computed
+                       !  infrared brightness temperature. Note that this
+                       !  calculation is most appropriate to compare to ISCCP
+                       !  IR only algortihm (i.e. you can compare to nighttime
+                       !  ISCCP data with this option)
+  ISCCP_TOPHEIGHT_DIRECTION=2,   ! direction for finding atmosphere pressure level
+                                 ! with interpolated temperature equal to the radiance
+                                 ! determined cloud-top temperature
+                                 ! 1 = find the *lowest* altitude (highest pressure) level
+                                 ! with interpolated temperature equal to the radiance
+                                 ! determined cloud-top temperature
+                                 ! 2 = find the *highest* altitude (lowest pressure) level
+                                 ! with interpolated temperature equal to the radiance 
+                                 ! determined cloud-top temperature
+                                 ! ONLY APPLICABLE IF top_height EQUALS 1 or 3
+                                 ! 1 = default setting, and matches all versions of 
+                                 ! ISCCP simulator with versions numbers 3.5.1 and lower
+                                 ! 2 = experimental setting  
+  !----------------------------------------------------------------------------------
+  !-------------- RTTOV inputs
+  !----------------------------------------------------------------------------------
+  Platform=1,    ! satellite platform
+  Satellite=15,  ! satellite
+  Instrument=0,  ! instrument
+  Nchannels=8,   ! Number of channels to be computed
+  Channels=1,3,5,6,8,10,11,13,        ! Channel numbers (please be sure that you supply Nchannels)
+  Surfem=0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,  ! Surface emissivity (please be sure that you supply Nchannels)
+  ZenAng=50.0, ! Satellite Zenith Angle
+  CO2=5.241e-04, ! Mixing ratios of trace gases
+  CH4=9.139e-07,
+  N2O=4.665e-07,
+  CO=2.098e-07
+/
+
+
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/cosp_output_nl.txt
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/cosp_output_nl.txt	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/cosp_output_nl.txt	(revision 1934)
@@ -0,0 +1,63 @@
+! (c) British Crown Copyright 2008, the Met Office.
+! All rights reserved.
+! 
+! Redistribution and use in source and binary forms, with or without modification, are permitted 
+! provided that the following conditions are met:
+! 
+!     * Redistributions of source code must retain the above copyright notice, this list 
+!       of conditions and the following disclaimer.
+!     * Redistributions in binary form must reproduce the above copyright notice, this list
+!       of conditions and the following disclaimer in the documentation and/or other materials 
+!       provided with the distribution.
+!     * Neither the name of the Met Office nor the names of its contributors may be used 
+!       to endorse or promote products derived from this software without specific prior written 
+!       permission.
+! 
+! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
+! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
+! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR 
+! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 
+! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, 
+! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER 
+! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT 
+! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+! Namelist that sets up output-related variables. It controls 
+! the instrument simulators to run and the list of variables 
+! to be written to file
+&COSP_OUTPUT
+  ! Simulator flags
+  Lradar_sim=.false.,
+  Llidar_sim=.true.,
+  Lisccp_sim=.true.,
+  Lmisr_sim=.false.,
+  ! Output variables
+  Lalbisccp=.true.,
+  Latb532=.false.,
+  Lboxptopisccp=.true.,
+  Lboxtauisccp=.true.,
+  Lcfad_dbze94=.false.,
+  Lcfad_lidarsr532=.true.,
+  Lclcalipso=.true.,
+  Lclhcalipso=.true.,
+  Lclisccp2=.true.,
+  Lcllcalipso=.true.,
+  Lclmcalipso=.true.,
+  Lcltcalipso=.true.,
+  Lctpisccp=.true.,
+  Ldbze94=.false.,
+  Ltauisccp=.true.,
+  Ltclisccp=.true.,
+  Llongitude=.false.,
+  Llatitude=.false.,
+  Lparasol_refl=.true.,
+  LclMISR=.false.,
+  Lmeantbisccp=.true.,
+  Lmeantbclrisccp=.true.,
+  ! Use lidar and radar
+  Lclcalipso2=.false.,
+  Lcltlidarradar=.false.,
+  ! These are provided for debugging or special purposes
+  Lfrac_out=.false.,
+  Lbeta_mol532=.false.,  
+/
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_144x142x19
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_144x142x19	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_144x142x19	(revision 1934)
@@ -0,0 +1,54 @@
+## $Id$
+## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
+day_step=720
+## periode pour le pas Matsuno (en pas)
+iperiod=5
+## choix de l'operateur de dissipation (star ou  non star )
+lstardis=y
+## nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+## nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+## nombre d'iterations de l'operateur de dissipation  divgrad            
+niterh=2
+## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+tetagdiv=3600.
+## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=5400.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
+tetatemp=5400.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
+purmats=n
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+## iphysiq est modifie pour la nouvelle physique par le lmdz.driver pour les autres cas le valeur par defaut est utilise
+iphysiq=_AUTO_ : DEFAULT = 15
+##  Avec ou sans strato
+ok_strato=n
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=1
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=1.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom                                  
+clat=0.
+## facteur de grossissement du zoom,selon longitude                      
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude                      
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.         
+ysinus=y
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_144x142x39
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_144x142x39	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_144x142x39	(revision 1934)
@@ -0,0 +1,54 @@
+## $Id$
+## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
+day_step=720
+## periode pour le pas Matsuno (en pas)
+iperiod=5
+## choix de l'operateur de dissipation (star ou  non star )
+lstardis=y
+## nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+## nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+## nombre d'iterations de l'operateur de dissipation  divgrad            
+niterh=2
+## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+tetagdiv=3600.
+## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=5400.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
+tetatemp=5400.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
+purmats=n
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+## iphysiq est modifie pour la nouvelle physique par le lmdz.driver pour les autres cas le valeur par defaut est utilise
+iphysiq=_AUTO_ : DEFAULT = 15
+##  Avec ou sans strato
+ok_strato=y
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=2
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=5.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom                                  
+clat=0.
+## facteur de grossissement du zoom,selon longitude                      
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude                      
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.         
+ysinus=y
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_56x55x19
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_56x55x19	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_56x55x19	(revision 1934)
@@ -0,0 +1,72 @@
+## $Header: /home/ioipsl/CVSROOT/IPSLCM4_v1_OASIS3/EXP00/gcm.def_ORCA4xLMD7245,v 1.1.1.1 2005/03/23 15:44:29 adm Exp $
+## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
+day_step=240
+## periode pour le pas Matsuno (en pas)
+iperiod=5
+## choix de l'operateur de dissipation (star ou  non star )
+lstardis=y
+## nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+## nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+## nombre d'iterations de l'operateur de dissipation  divgrad            
+niterh=2
+## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+tetagdiv=36000.
+## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=18000.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
+tetatemp=18000.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
+purmats=n
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+## iphysiq est modifie pour la nouvelle physique par le lmdz.driver pour les autres cas le valeur par defaut est utilise
+iphysiq=_AUTO_ : DEFAULT = 5
+##  Activation ou non de la parametrisation de Hines pour la strato
+ok_strato=n
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=1
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=1.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom                                  
+clat=0.
+## facteur de grossissement du zoom,selon longitude                      
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude                      
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.         
+ysinus=y
+##  [FC] Transport par les ondes stationnaires selon Francis Codron
+##  Choix possibles :
+##Â  * dissuv=y, diss_uvwin=n, dissuvz=n (local)
+##  * dissuv=n, diss_uvwin=n, dissuvz=y (moyenne zonale)
+##  * dissuv=y, diss_uvwin=y, dissuvz=n (lissage sur 120Â° de longitude)
+##Â  Dans le permier cas, mettre fact_uv Ã  0.2 ou 0.3 dans dissip_p.F,
+##  sinon c'est un peu fort...
+dissuv=n
+diss_uvwin=n
+dissuvz=n
+## [FC] fichiers instantanes de controle une fois par jour
+out_diss=n
+## Orosm pour le lissage de l'orog.
+ok_orosm=n
+# GKLIFT : Gravity Wave Drag calculation parameter
+GKLIFT=0.5
+## Fract pour le fractionnement des mailles
+ok_fract=.true.
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x71x19
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x71x19	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x71x19	(revision 1934)
@@ -0,0 +1,54 @@
+## $Id$
+## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
+day_step=480
+## periode pour le pas Matsuno (en pas)
+iperiod=5
+## choix de l'operateur de dissipation (star ou  non star )
+lstardis=y
+## nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+## nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+## nombre d'iterations de l'operateur de dissipation  divgrad            
+niterh=2
+## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+tetagdiv=10800.
+## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=18000.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
+tetatemp=18000.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
+purmats=n
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+## iphysiq est modifie pour la nouvelle physique par le lmdz.driver pour les autres cas le valeur par defaut est utilise
+iphysiq=_AUTO_ : DEFAULT = 10
+##  Avec ou sans strato
+ok_strato=n
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=1
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=1.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom                                  
+clat=0.
+## facteur de grossissement du zoom,selon longitude                      
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude                      
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.         
+ysinus=y
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x95x19
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x95x19	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x95x19	(revision 1934)
@@ -0,0 +1,54 @@
+## $Id$
+## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
+day_step=480
+## periode pour le pas Matsuno (en pas)
+iperiod=5
+## choix de l'operateur de dissipation (star ou  non star )
+lstardis=y
+## nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+## nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+## nombre d'iterations de l'operateur de dissipation  divgrad            
+niterh=2
+## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+tetagdiv=5400.
+## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=7200.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
+tetatemp=7200.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
+purmats=n
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+## iphysiq est modifie pour la nouvelle physique par le lmdz.driver pour les autres cas le valeur par defaut est utilise
+iphysiq=_AUTO_ : DEFAULT = 10
+##  Avec ou sans strato
+ok_strato=n
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=1
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=1.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom                                  
+clat=0.
+## facteur de grossissement du zoom,selon longitude                      
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude                      
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.         
+ysinus=y
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x95x39
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x95x39	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/gcm.def_96x95x39	(revision 1934)
@@ -0,0 +1,54 @@
+## $Id$
+## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
+day_step=480
+## periode pour le pas Matsuno (en pas)
+iperiod=5
+## choix de l'operateur de dissipation (star ou  non star )
+lstardis=y
+## nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+## nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+## nombre d'iterations de l'operateur de dissipation  divgrad            
+niterh=2
+## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+tetagdiv=5400.
+## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=5400.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
+tetatemp=5400.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
+purmats=n
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+## iphysiq est modifie pour la nouvelle physique par le lmdz.driver pour les autres cas le valeur par defaut est utilise
+iphysiq=_AUTO_ : DEFAULT = 10
+##  Avec ou sans strato
+ok_strato=y 
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=2
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=5.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom                                  
+clat=0.
+## facteur de grossissement du zoom,selon longitude                      
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude                      
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.         
+ysinus=y
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/guide.def
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/guide.def	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/guide.def	(revision 1934)
@@ -0,0 +1,16 @@
+# ok_guide=y : guidage active
+ok_guide=_AUTO_ : DEFFAULT = n
+#
+guide_u= y
+guide_v= y
+guide_T= n
+guide_P= n
+guide_Q= n
+
+ini_anal=n
+tau_min_u=0.0208333
+tau_max_u=0.125
+tau_min_v=0.0208333
+tau_max_v=0.125
+#tau_min_T=0.0208333
+#tau_max_T=10
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_high
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_high	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_high	(revision 1934)
@@ -0,0 +1,131 @@
+#
+# Parametres fichiers de sortie
+#
+# Noms des fichiers
+phys_out_filenames=      histmth histday histhf   histhf3h  histhf3hm histstn
+# Sortir ou non les fichiers
+phys_out_filekeys= _AUTO_
+# Niveaux de sorties
+phys_out_filelevels=     5       5       5        5         5         5
+### Type de fichier : global (n) ou stations (y)
+phys_out_filestations =  n       n       n        n         n         y
+### Frequences des sorties
+phys_out_filetimesteps = 1.mth,  1.day,  0.25day, 0.125day, 0.125day, 1800.s
+#
+#
+### variables niveaux pressions 100, 50, 10 hPa a 3h-inst - F.Lott
+flag_t100       = 1, 7, 7, 5, 10, 10
+flag_u100       = 1, 7, 7, 5, 10, 10
+flag_v100       = 1, 7, 7, 5, 10, 10
+flag_w100       = 1, 7, 7, 5, 10, 10
+flag_t50        = 1, 7, 7, 5, 10, 10
+flag_u50        = 1, 7, 7, 5, 10, 10
+flag_v50        = 1, 7, 7, 5, 10, 10
+flag_w50        = 1, 7, 7, 5, 10, 10
+flag_t10        = 1, 7, 7, 5, 10, 10
+flag_u10        = 1, 7, 7, 5, 10, 10
+flag_v10        = 1, 7, 7, 5, 10, 10
+flag_w10        = 1, 7, 7, 5, 10, 10
+# more variables for level 5 histday/month: vitw and pres
+# Vertical wind
+flag_vitw         =  2, 3, 7, 6, 6, 6
+# Air pressure
+flag_pres         =  2, 3, 7, 6, 6, 6
+# w500 for histday
+flag_w500         =  1, 5, 3, 10, 10, 10
+# t700 for histday
+flag_t700         =  1, 5, 3, 10, 10, 10
+# Cloud fraction
+flag_rneb         =  2, 5, 10, 5, 10, 10
+# Aga sortie avec niveau 2
+flag_Aga          =  2, 10, 10, 10, 10, 10
+# niveaux maximum 5 mensuel pour variables CFMIP
+flag_rlu          =  5, 10, 10, 10, 10, 10
+flag_rsu          =  5, 10, 10, 10, 10, 10
+flag_rld          =  5, 10, 10, 10, 10, 10
+flag_rsd          =  5, 10, 10, 10, 10, 10
+flag_rlucs        =  5, 10, 10, 10, 10, 10
+flag_rsucs        =  5, 10, 10, 10, 10, 10
+flag_rldcs        =  5, 10, 10, 10, 10, 10
+flag_rsdcs        =  5, 10, 10, 10, 10, 10
+flag_tnt          =  5, 10, 10, 10, 10, 10
+flag_dtdyn        =  4, 10, 10, 10, 10, 10
+flag_dtphy        =  4, 10, 10, 10, 10, 10
+flag_dtlschr      =  5, 10, 10, 10, 10, 10
+flag_tntr         =  5, 10, 10, 10, 10, 10
+flag_tntc         =  5, 10, 10, 10, 10, 10
+flag_tntscpbl     =  5, 10, 10, 10, 10, 10
+flag_tnhus        =  5, 10, 10, 10, 10, 10
+flag_dqdyn        =  4, 10, 10, 10, 10, 10
+flag_tnhusc       =  5, 10, 10, 10, 10, 10
+flag_dqvdf        =  5, 10, 10, 10, 10, 10
+flag_tnhusscpbl   =  5, 10, 10, 10, 10, 10
+flag_dqphy        =  2, 10, 10, 10, 10, 10
+flag_evu          =  5, 10, 10, 10, 10, 10
+flag_kz           =  4, 10, 10, 10, 10, 10
+###
+# more variables for t tendency monthly
+###
+flag_dtvdf        =  5, 10, 10, 10, 10, 10
+flag_dteva        =  5, 10, 10, 10, 10, 10
+flag_dtcon        =  4, 10, 10, 10, 10, 10
+flag_dtajs        =  4, 10, 10, 10, 10, 10
+flag_dtlsc        =  4, 10, 10, 10, 10, 10
+flag_dtoro        =  4, 10, 10, 10, 10, 10
+flag_dtlif        =  4, 10, 10, 10, 10, 10
+flag_dthin        =  4, 10, 10, 10, 10, 10
+flag_dtec         =  4, 10, 10, 10, 10, 10
+### more variables for q tendency monthly
+flag_dqeva        =  5, 10, 10, 10, 10, 10
+flag_dqcon        =  4, 10, 10, 10, 10, 10
+flag_dqajs        =  4, 10, 10, 10, 10, 10
+flag_dqlsc        =  5, 10, 10, 10, 10, 10
+### nouvelle physique monthly
+flag_dtthe        =  4, 10, 10, 10, 10, 10
+flag_dtwak        =  4, 5, 10, 10, 10, 10
+flag_dqthe        =  4, 10, 10, 10, 10, 10
+flag_dqwak        =  4, 5, 10, 10, 10, 10
+flag_ftime_th     =  4, 10, 10, 5, 10, 10
+#
+# Parametres niveau de sorties differents fichiers
+#
+### parametres pour hist*NMC.nc
+#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
+#                histdayNMC: 8 ou 17
+lev_histdayNMC=8
+### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
+ok_histNMC = _AUTO_
+
+#
+# Parametres simulateur COSP (CFMIP Observational Simulator Package)
+#
+### ok_cosp=y/n flag simulateur COSP
+ok_cosp= _AUTO_
+### freq_COSP = frequence d'appel de COSP en secondes
+freq_COSP=10800.
+### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
+ok_mensuelCOSP= _AUTO_
+### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
+ok_journeCOSP= _AUTO_
+### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
+ok_hfCOSP= _AUTO_
+#
+# Parametres simulateur ISCCP
+#
+### ok_isccp=y/n flag simulateur ISCCP
+ok_isccp=n
+### freqin_isccp = frequence input en secondes du simulateur ISCCP
+freq_ISCCP=10800.
+### freqout_isccp = frequence output en jours du simulateur ISCCP
+ecrit_ISCCP= _AUTO_
+### top_height = flag choix calcul nuages par le simulateur en utilisant 
+# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
+# - top_height = 1 -> algo IR-VIS
+# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
+# - top_height = 3 -> algo IR
+top_height = 1
+### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
+# - overlap=1    Max overlap
+# - overlap=2    Random overlap
+# - overlap=3    Max/Random overlap
+overlap = 3
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_low
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_low	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_low	(revision 1934)
@@ -0,0 +1,60 @@
+#
+# Parametres fichiers de sortie
+#
+# Noms des fichiers
+phys_out_filenames=      histmth histday histhf   histhf3h  histhf3hm histstn
+# Sortir ou non les fichiers
+phys_out_filekeys=  _AUTO_
+# Niveaux de sorties
+phys_out_filelevels=     5       2       2        5         5         5
+### Type de fichier : global (n) ou stations (y)
+phys_out_filestations =  n       n       n        n         n         y
+### Frequences des sorties
+phys_out_filetimesteps = 1.mth,  1.day,  0.25day, 0.125day, 0.125day, 1800.s
+#
+#
+# Parametres niveau de sorties differents fichiers
+#
+### parametres pour hist*NMC.nc
+#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
+#                histdayNMC: 8 ou 17
+lev_histdayNMC=8
+### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
+ok_histNMC = _AUTO_
+
+#
+# Parametres simulateur COSP (CFMIP Observational Simulator Package)
+#
+### ok_cosp=y/n flag simulateur COSP
+ok_cosp= _AUTO_
+### freq_COSP = frequence d'appel de COSP en secondes
+freq_COSP=10800.
+### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
+ok_mensuelCOSP= _AUTO_
+### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
+ok_journeCOSP= _AUTO_
+### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
+ok_hfCOSP= _AUTO_
+
+#
+# Parametres simulateur ISCCP
+#
+### ok_isccp=y/n flag simulateur ISCCP
+ok_isccp=n
+### freqin_isccp = frequence input en secondes du simulateur ISCCP
+freq_ISCCP=10800.
+### freqout_isccp = frequence output en jours du simulateur ISCCP
+ecrit_ISCCP= _AUTO_
+### top_height = flag choix calcul nuages par le simulateur en utilisant 
+# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
+# - top_height = 1 -> algo IR-VIS
+# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
+# - top_height = 3 -> algo IR
+top_height = 1
+### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
+# - overlap=1    Max overlap
+# - overlap=2    Random overlap
+# - overlap=3    Max/Random overlap
+overlap = 3
+
+
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_medium
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_medium	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/output.def_medium	(revision 1934)
@@ -0,0 +1,131 @@
+#
+# Parametres fichiers de sortie
+#
+# Noms des fichiers
+phys_out_filenames=      histmth histday histhf   histhf3h  histhf3hm histstn
+# Sortir ou non les fichiers
+phys_out_filekeys= _AUTO_
+# Niveaux de sorties
+phys_out_filelevels=     5       5       2        5         5         5
+### Type de fichier : global (n) ou stations (y)
+phys_out_filestations =  n       n       n        n         n         y
+### Frequences des sorties
+phys_out_filetimesteps = 1.mth,  1.day,  0.25day, 0.125day, 0.125day, 1800.s
+#
+#
+### variables niveaux pressions 100, 50, 10 hPa a 3h-inst - F.Lott
+flag_t100       = 1, 7, 7, 5, 10, 10
+flag_u100       = 1, 7, 7, 5, 10, 10
+flag_v100       = 1, 7, 7, 5, 10, 10
+flag_w100       = 1, 7, 7, 5, 10, 10
+flag_t50        = 1, 7, 7, 5, 10, 10
+flag_u50        = 1, 7, 7, 5, 10, 10
+flag_v50        = 1, 7, 7, 5, 10, 10
+flag_w50        = 1, 7, 7, 5, 10, 10
+flag_t10        = 1, 7, 7, 5, 10, 10
+flag_u10        = 1, 7, 7, 5, 10, 10
+flag_v10        = 1, 7, 7, 5, 10, 10
+flag_w10        = 1, 7, 7, 5, 10, 10
+#
+# more variables for level 5 : vitw and pres
+# Vertical wind
+flag_vitw         =  2, 3, 7, 6, 6, 6
+# Air pressure
+flag_pres         =  2, 3, 7, 6, 6, 6
+# w500 for histday
+flag_w500         =  1, 5, 3, 10, 10, 10
+# t700 for histday
+flag_t700         =  1, 5, 3, 10, 10, 10
+# Cloud fraction
+flag_rneb         =  2, 5, 10, 5, 10, 10
+# Aga sortie avec niveau 2
+flag_Aga          =  2, 10, 10, 10, 10, 10
+# niveaux maximum 5 mensuel pour variables CFMIP
+flag_rlu          =  5, 10, 10, 10, 10, 10
+flag_rsu          =  5, 10, 10, 10, 10, 10
+flag_rld          =  5, 10, 10, 10, 10, 10
+flag_rsd          =  5, 10, 10, 10, 10, 10
+flag_rlucs        =  5, 10, 10, 10, 10, 10
+flag_rsucs        =  5, 10, 10, 10, 10, 10
+flag_rldcs        =  5, 10, 10, 10, 10, 10
+flag_rsdcs        =  5, 10, 10, 10, 10, 10
+flag_tnt          =  5, 10, 10, 10, 10, 10
+flag_dtdyn        =  4, 10, 10, 10, 10, 10
+flag_dtphy        =  4, 10, 10, 10, 10, 10
+flag_dtlschr      =  5, 10, 10, 10, 10, 10
+flag_tntr         =  5, 10, 10, 10, 10, 10
+flag_tntc         =  5, 10, 10, 10, 10, 10
+flag_tntscpbl     =  5, 10, 10, 10, 10, 10
+flag_tnhus        =  5, 10, 10, 10, 10, 10
+flag_dqdyn        =  4, 10, 10, 10, 10, 10
+flag_tnhusc       =  5, 10, 10, 10, 10, 10
+flag_dqvdf        =  5, 10, 10, 10, 10, 10
+flag_tnhusscpbl   =  5, 10, 10, 10, 10, 10
+flag_dqphy        =  2, 10, 10, 10, 10, 10
+flag_evu          =  5, 10, 10, 10, 10, 10
+flag_kz           =  4, 10, 10, 10, 10, 10
+###
+# more variables for t tendency monthly
+###
+flag_dtvdf        =  5, 10, 10, 10, 10, 10
+flag_dteva        =  5, 10, 10, 10, 10, 10
+flag_dtcon        =  4, 10, 10, 10, 10, 10
+flag_dtajs        =  4, 10, 10, 10, 10, 10
+flag_dtlsc        =  4, 10, 10, 10, 10, 10
+flag_dtoro        =  4, 10, 10, 10, 10, 10
+flag_dtlif        =  4, 10, 10, 10, 10, 10
+flag_dthin        =  4, 10, 10, 10, 10, 10
+flag_dtec         =  4, 10, 10, 10, 10, 10
+### more variables for q tendency monthly
+flag_dqeva        =  5, 10, 10, 10, 10, 10
+flag_dqcon        =  4, 10, 10, 10, 10, 10
+flag_dqajs        =  4, 10, 10, 10, 10, 10
+flag_dqlsc        =  5, 10, 10, 10, 10, 10
+### nouvelle physique monthly
+flag_dtthe        =  4, 10, 10, 10, 10, 10
+flag_dtwak        =  4, 5, 10, 10, 10, 10
+flag_dqthe        =  4, 10, 10, 10, 10, 10
+flag_dqwak        =  4, 5, 10, 10, 10, 10
+flag_ftime_th     =  4, 10, 10, 5, 10, 10
+#
+# Parametres niveau de sorties differents fichiers
+#
+### parametres pour hist*NMC.nc
+#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
+#                histdayNMC: 8 ou 17
+lev_histdayNMC=8
+### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
+ok_histNMC = _AUTO_
+#
+# Parametres simulateur COSP (CFMIP Observational Simulator Package)
+#
+### ok_cosp=y/n flag simulateur COSP
+ok_cosp= _AUTO_
+### freq_COSP = frequence d'appel de COSP en secondes
+freq_COSP=10800.
+### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
+ok_mensuelCOSP= _AUTO_
+### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
+ok_journeCOSP= _AUTO_
+### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
+ok_hfCOSP= _AUTO_
+#
+# Parametres simulateur ISCCP
+#
+### ok_isccp=y/n flag simulateur ISCCP
+ok_isccp=n
+### freqin_isccp = frequence input en secondes du simulateur ISCCP
+freq_ISCCP=10800.
+### freqout_isccp = frequence output en jours du simulateur ISCCP
+ecrit_ISCCP= _AUTO_
+### top_height = flag choix calcul nuages par le simulateur en utilisant 
+# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
+# - top_height = 1 -> algo IR-VIS
+# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
+# - top_height = 3 -> algo IR
+top_height = 1
+### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
+# - overlap=1    Max overlap
+# - overlap=2    Random overlap
+# - overlap=3    Max/Random overlap
+overlap = 3
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L19_AP
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L19_AP	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L19_AP	(revision 1934)
@@ -0,0 +1,97 @@
+## $Id$
+#
+##  Choix ou non  de  New oliq               
+new_oliq=y
+##  Activation ou non de la parametrisation de Hines pour la strato
+ok_hines=n
+##  Orodr  ou  non   pour l orographie              
+ok_orodr=y
+##  Orolf  ou  non   pour l orographie              
+ok_orolf=y
+##   Si = .T. ,  lecture du fichier limit avec la bonne annee             
+ok_limitvrai=n
+##  Flag  pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect
+iflag_con=30
+## Facteur multiplication des precip convectives dans KE
+cvl_corr=1.0
+#
+#
+### niveau du diagnostique de conservation d energie
+if_ebil=0
+#
+# parametres KE
+#
+### epmax = Efficacite precipitation maximale
+epmax = .99
+### ok_adj_ema = ?? pas utilise
+ok_adj_ema = n
+### iflag_clw Flag calcul eau liquide
+# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
+# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
+# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
+iflag_clw = 1
+# 
+# parametres nuages
+#
+### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp) 
+cld_lc_lsc = 2.6e-4
+### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp) 
+cld_lc_con = 2.6e-4
+### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) 
+cld_tau_lsc = 3600.
+### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) 
+cld_tau_con = 3600.
+### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) 
+ffallv_lsc = 1.
+### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) 
+ffallv_con = 1.
+### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
+coef_eva = 2.e-5
+### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) 
+reevap_ice = y
+### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) 
+# - iflag_cldcon<=-1 diag. rain_Tiedtke 
+# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
+# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
+# - iflag_cldcon=3,   ratqs=ratqss
+iflag_cldcon = 3
+### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages 
+# - iflag_pdf=0, version ratqs, 
+# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
+# -              des PDFS de Sandrine Bony 
+iflag_pdf = 1
+### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) 
+fact_cldcon = 1.
+### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
+facttemps = 0.
+## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) 
+ok_newmicro = y
+### iflag_ratqs=0 correspond a la version IPCC AR4 
+iflag_ratqs=0
+### ratqsbas     ratqs en bas si iflag_cldcon=1 
+ratqsbas = 0.005
+### ratqshaut    ratqs en haut pour ratqss "stables"  
+ratqshaut = 0.33
+### rad_froid    rayon cristaux des nuages de glace "froids" 
+rad_froid = 35
+### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev 
+rad_chau1=12
+### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3  
+rad_chau2=11
+# 
+# Coefficient et parametres sur les drags
+#
+f_cdrag_ter=1.
+f_cdrag_stable=1.
+f_cdrag_oce=0.8
+f_rugoro=0.
+### ksta_ter      coef.diffusion minimale sur terre/sic/lic 
+ksta_ter=1.e-7
+### cdmmax =     cdrag maximum pour le moment
+cdmmax = 2.5E-3
+### cdhmax =     cdrag maximum pour l energie
+cdhmax = 2.0E-3
+
+### nsplit_phys (D=1) 
+nsplit_phys=1
+
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L39_AP
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L39_AP	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L39_AP	(revision 1934)
@@ -0,0 +1,97 @@
+## $Id$
+#
+##  Choix ou non  de  New oliq               
+new_oliq=y
+##  Activation ou non de la parametrisation de Hines pour la strato
+ok_hines=y
+##  Orodr  ou  non   pour l orographie              
+ok_orodr=y
+##  Orolf  ou  non   pour l orographie              
+ok_orolf=y
+##   Si = .T. ,  lecture du fichier limit avec la bonne annee             
+ok_limitvrai=n
+##  Flag  pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect
+iflag_con=30
+## Facteur multiplication des precip convectives dans KE
+cvl_corr=1.0
+#
+#
+
+### niveau du diagnostique de conservation d energie
+if_ebil=0
+#
+# parametres KE
+#
+### epmax = Efficacite precipitation maximale
+epmax = .999
+### ok_adj_ema = ?? pas utilise
+ok_adj_ema = n
+### iflag_clw Flag calcul eau liquide
+# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
+# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
+# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
+iflag_clw = 1
+# 
+# parametres nuages
+#
+### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp) 
+cld_lc_lsc = 4.16e-4
+### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp) 
+cld_lc_con = 4.16e-4
+### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) 
+cld_tau_lsc = 1800.
+### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) 
+cld_tau_con = 1800.
+### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) 
+ffallv_lsc = 0.5
+### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) 
+ffallv_con = 0.5
+### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
+coef_eva = 2.e-5
+### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) 
+reevap_ice = y
+### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) 
+# - iflag_cldcon<=-1 diag. rain_Tiedtke 
+# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
+# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
+# - iflag_cldcon=3,   ratqs=ratqss
+iflag_cldcon = 3
+### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages 
+# - iflag_pdf=0, version ratqs, 
+# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
+# -              des PDFS de Sandrine Bony 
+iflag_pdf = 1
+### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) 
+fact_cldcon = 1.
+### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
+facttemps = 0.
+## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) 
+ok_newmicro = y
+### iflag_ratqs=0 correspond a la version IPCC AR4 
+iflag_ratqs=0
+### ratqsbas     ratqs en bas si iflag_cldcon=1 
+ratqsbas = 0.005
+### ratqshaut    ratqs en haut pour ratqss "stables"  
+ratqshaut = 0.33
+### rad_froid    rayon cristaux des nuages de glace "froids" 
+rad_froid = 35
+### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev 
+rad_chau1=12
+### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3  
+rad_chau2=11
+# 
+# Coefficient et parametres sur les drags
+#
+f_cdrag_ter=1.
+f_cdrag_stable=1.
+f_cdrag_oce=0.8
+f_rugoro=0.
+### ksta_ter      coef.diffusion minimale sur terre/sic/lic 
+ksta_ter=1.e-7
+### cdmmax =     cdrag maximum pour le moment
+cdmmax = 2.5E-3
+### cdhmax =     cdrag maximum pour l energie
+cdhmax = 2.0E-3
+
+### nsplit_phys (D=1) 
+nsplit_phys=1
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L39_NPv3.1
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L39_NPv3.1	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/physiq.def_L39_NPv3.1	(revision 1934)
@@ -0,0 +1,132 @@
+##############################################################################
+# PARAMETRES ET CLES DE CONFIGURATION LMDZ 
+# Cette parie sera bientot extrait dans un fichier config.def
+##############################################################################
+# D: Valeur par default
+#
+##############################################################################
+# PARAMETRES ET CLES DE LA PHYSIQUE LMDZ
+##############################################################################
+# D: Valeur par default
+#
+#Parametres Orographiques et cdrags
+#**********************************
+##   Si=.T. ,  lecture du fichier limit avec la bonne annee 
+ok_limitvrai=n
+#Cdrags 
+f_cdrag_stable=1.
+f_cdrag_ter=1.
+f_cdrag_oce=0.7
+#cdmmax
+cdmmax=2.5E-3
+#cdhmax
+cdhmax=2.0E-3
+##  Orodr  ou  non   pour l orographie
+ok_orodr=y
+##  Orolf  ou  non   pour l orographie              
+ok_orolf=y
+#Rugoro
+f_rugoro=0.
+#
+# Rayonnement
+#*********************
+# activation du nouveau code de rayonnement RRTM
+# 0 : Ancien code et 1 : RRTM (D=0)
+iflag_rrtm=0
+#
+##  Activation ou non de la parametrisation de Hines pour la strato
+#******************************************************************
+ok_hines=y
+#
+# Parametres nuages
+#******************
+#seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
+cld_lc_lsc=6.e-4
+cld_lc_con=6.e-4
+#constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.)
+cld_tau_lsc=1800.
+cld_tau_con=1800.
+#facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
+ffallv_lsc=1.35
+ffallv_con=1.35
+#coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
+coef_eva=1e-4
+reevap_ice=y
+#calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4) 
+iflag_cldcon=6
+fact_cldcon=1.
+facttemps=0.
+#calcul eau condensee et fraction nuageuse a partir des PDF (D:0, 0:version avec ratqs sinon nouvelles PDFs)  
+iflag_pdf=1
+#calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
+ok_newmicro=y
+iflag_ratqs=2
+ratqsbas=0.002
+ratqshaut=0.25
+#rayon effectif pour les nuages de glace (D:35)
+rad_froid=35
+#rayons effectifs pour les nuages eau liq (tailles des gouttes d eau nuageuse, D: 13 9)
+rad_chau1=12
+rad_chau2=11
+##  Choix ou non  de  New oliq               
+new_oliq=y
+#
+rei_min=20.
+rei_max=61.29
+#
+t_glace_min=258.
+t_glace_max=273.13
+#
+# Flags Convection 
+#*****************
+#flag  pour la convection (D:2, 1:LMD, 2:Tiedtke, 3:KE nouvelle physique, 30:KE IPCC)
+iflag_con=3
+#niveau de sorties des diagnos sur la conservation energie
+if_ebil=0
+#efficacite de precepitation maxim (D:.993)
+epmax=0.997
+#ajustement convectif sec au debut de la convection (D: n 0)
+ok_adj_ema=n
+iflag_clw=0
+#fermeture convective (D:1, 0:AR4, 2:ALE and ALP)
+iflag_clos=2
+#loi de melange a l entrainement  (D:1, 0:AR4=PDF plate, 1=PDF)
+iflag_mix=1
+#poids des PDFs plate et en cloche (D: 1 0)
+qqa1=1.
+qqa2=0.
+# Facteur multiplication des precip convectives dans KE
+cvl_corr=1.0
+#
+Fmax=0.65
+#
+# Flags Couche limite
+#********************
+#shema de CL des thermiques (D:0, 0:ajustement sec,=>1:versions thermiques) 
+iflag_thermals=15
+#pas de temps du spliting pour les thermiques
+nsplit_thermals=1
+tau_thermals=0.
+iflag_thermals_ed=10
+iflag_thermals_optflux=0
+#shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada) 
+iflag_pbl=8
+#diffusion turbulente
+ksta_ter=1.e-7
+ksta=1.e-10
+#ok_kzmin : calcul Kzmin dans la CL de surface
+ok_kzmin=y
+#couplage ou non avec la convection
+iflag_coupl=5 (0:AR4 et 1:nouvelle physique)
+seuil_inversion=-0.08
+#rendement sur la puissance fournie par les thermiques a la convection
+alp_bl_k=0.5
+#
+# Flags wakes
+#************
+#activation ou non des wakes (D:0, 0:AR4, 1:nouvelle physique) 
+iflag_wake=1
+alp_offset=0.
+
+nsplit_phys=2
+
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/run.def
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/run.def	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/run.def	(revision 1934)
@@ -0,0 +1,55 @@
+# $Id$
+#
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+INCLUDEDEF=output.def
+INCLUDEDEF=config.def
+INCLUDEDEF=reprobus.def
+INCLUDEDEF=guide.def
+## Type de calendrier utilise
+## valeur possible: earth_360d (defaut), earth_365d, earth_366d
+calend= _AUTO_
+## Jour de l etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+dayref= _AUTO_
+##  Annee de l etat  initial (   avec  4  chiffres   )
+anneeref= _AUTO_
+## Nombre de jours d integration
+nday= _AUTO_
+## Remise a zero de la date initiale
+raz_date= _AUTO_
+## periode de sortie des variables de controle (en pas)
+iconser=240
+## periode d ecriture du fichier histoire (en jour)
+iecri=1
+## flag de sortie dynzon
+ok_dynzon=n
+## periode de stockage fichier dynzon (en jour)
+periodav= _AUTO_
+## activation du calcul d equilibrage de charge
+adjust= _AUTO_
+## activation du filtre fft
+use_filtre_fft=y
+## niveau d impression de controle
+prt_level=0
+###
+### Informations sur la configuration utilisee
+###
+### type_ocean = force / slab  /couple
+### type_ocean is modified automatically by the lmdz.driver 
+type_ocean=_AUTO_
+### version_ocean = nemo / opa8
+version_ocean=nemo
+### activation passage courants surface ocean-glace
+cpl_current=y
+### VEGET= y si ORCHIDEE, =n si bucket
+### VEGET is modified automatically by the lmdz.driver 
+VEGET=_AUTO_
+### Choice of tracers type_trac=lmdz/inca/repr
+### type_trac is set automatically by the lmdz.driver
+type_trac=_AUTO_
+### config_inca=aero/chem needed only if type_trac=inca
+config_inca=_AUTO_ 
+### Only concerning the program ce0l : for creating file grilles_gcm.nc
+grilles_gcm_netcdf=_AUTO_
+
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/traceur.def
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/traceur.def	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/PARAM/traceur.def	(revision 1934)
@@ -0,0 +1,4 @@
+3
+14 14 H2Ov
+10 10 H2Ol
+10 10 Aga
Index: CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/POST/monitoring01_lmdz.cfg
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/POST/monitoring01_lmdz.cfg	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/GENERAL/POST/monitoring01_lmdz.cfg	(revision 1934)
@@ -0,0 +1,84 @@
+#**************************************************************
+# Author: Patrick Brockmann
+# Contact: Patrick.Brockmann@cea.fr
+# $Date: 2009/05/14 15:09:35 $
+# $Name: ATLAS_608_1_5 $
+# $Revision: 1.1.1.2 $
+# History:
+# Modification:
+#**************************************************************
+
+##################################################################################################################
+#-----------------------------------------------------------------------------------------------------------------
+#  field | files patterns | files additionnal | operations | title | units | calcul of area
+#-----------------------------------------------------------------------------------------------------------------
+#
+#   field = name of the field to monitore 
+#   files patterns = pattern of the files used to monitore the field (must be discriminant) 
+#   files additionnal = additionnal files accessible by $FER_DATA
+#   operations = operations to calculate the field 
+#   title = variable title
+#   units = variable units
+#   calcul of area = expression to use for the weight average
+#
+# Notes:
+#   - Comment lines begin by # character.
+#   - Separator between fields is | character.
+#   - Operations must use the ferret syntax, if several files are used precise with the syntax d=x to refer to the x dataset.
+#     They must be enclosed by parenthesis if more than one variable is used.
+#   - files patterns,files additionnal,operations,title,units must be enclosed with character ". 
+#     Use "" if empty.
+#   - Fields will be presented through an html page with thumbnails global, north, south, land, ocean.
+#     Use a field name with one of these pattern to display it under the corresponding thumbnail.
+#
+##################################################################################################################
+
+color=AECDFF
+smooth=12
+#-----------------------------------------------------------------------------------------------------------------
+#  field | files patterns | files additionnal | operations | title | units | calcul of area
+#-----------------------------------------------------------------------------------------------------------------
+ bils_global    | "bils"                       | "" | "bils[d=1]"               | "Surf. total heat flux (GLOBAL)"         | "W/m^2"     | "aire[d=1]"   
+ nettop_global  | "tops topl"                  | "" | "(tops[d=1]-topl[d=2])"   | "TOA. total heat flux (GLOBAL)"          | "W/m^2"     | "aire[d=1]"
+ nettop0_global | "tops0 topl0"                | "" | "(tops0[d=1]-topl0[d=2])" | "TOA. total heat flux ClearSky (GLOBAL)" | "W/m^2"     | "aire[d=1]"
+ t2m_global     | "t2m"                        | "" | "(t2m[d=1]-273.15)"       | "Temperature at 2 meters (GLOBAL)"       | "degrees C" | "aire[d=1]"
+ t2m_land       | "t2m pourc_ter pourc_lic"    | "" | "(t2m[d=1]-273.15)"       | "Temperature at 2 meters (LAND)"         | "degrees C" | "aire[d=1]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)"
+ tsol_global    | "tsol"                       | "" | "(tsol[d=1]-273.15)"      | "Surface Temperature (GLOBAL)"           | "degrees C" | "aire[d=1]"
+ tsol_land      | "tsol pourc_ter pourc_lic"   | "" | "(tsol[d=1]-273.15)"      | "Surface Temperature (LAND)"             | "degrees C" | "aire[d=1]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)"
+ precip_global  | "precip"                     | "" | "precip[d=1]*86400"       | "Precipitation Totale liq+sol (GLOBAL)"  | "mm/d"      | "aire[d=1]"
+ precip_land    | "precip pourc_ter pourc_lic" | "" | "precip[d=1]*86400"       | "Precipitation Totale liq+sol (LAND)"    | "mm/d"      | "aire[d=1]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)"
+ evap_global    | "evap"                       | "" | "evap[d=1]"               | "Evaporation (GLOBAL)"                   | "kg/(s*m2)" | "aire[d=1]"
+ evap_land      | "evap pourc_ter pourc_lic"   | "" | "evap[d=1]"               | "Evaporation (LAND)"                     | "kg/(s*m2)" | "aire[d=1]*((pourc_ter[d=2]+pourc_lic[d=3])/100.)"
+ flat_global    | "lat_ter pourc_ter lat_sic pourc_sic lat_lic pourc_lic lat_oce pourc_oce" | "" | "(lat_ter[d=1]*pourc_ter[d=2]+lat_sic[d=3]*pourc_sic[d=4]+lat_lic[d=5]*pourc_lic[d=6]+lat_oce[d=7]*pourc_oce[d=8])" | "Latent Heat Flux (GLOBAL)"      | "W/m^2" | "aire[d=1]"
+ od550aer_forcing       | "od550aer"              | "" | "od550aer[d=1]"           | "Total aerosol optical depth at 550nm"   | "-"         | "aire[d=1]"   
+ od550lt1aer_forcing    | "od550lt1aer"           | "" | "od550lt1aer[d=1]"        | "Fine mode optical depth"                | "-"         | "aire[d=1]"   
+ absvisaer_forcing      | "absvisaer"             | "" | "absvisaer[d=1]"          | "Absorption aerosol visible optical depth" | "-"       | "aire[d=1]"   
+ swtoaas_ant_forcing    | "swtoaas_ant"           | "" | "swtoaas_ant[d=1]"        | "Anthropogenic aerosol radiative forcing all-sky at TOA"         | "W/m^2"   | "aire[d=1]"   
+ swtoacs_ant_forcing    | "swtoacs_ant"           | "" | "swtoacs_ant[d=1]"        | "Anthropogenic aerosol radiative forcing clear-sky at TOA"       | "W/m^2"   | "aire[d=1]" 
+ swtoacf_ant_forcing    | "swtoacf_ant"           | "" | "swtoacf_ant[d=1]"        | "Anthropogenic aerosol impact on cloud radiative forcing at TOA" | "W/m^2"   | "aire[d=1]" 
+ swsrfas_ant_forcing    | "swsrfas_ant"           | "" | "swsrfas_ant[d=1]"        | "Anthropogenic aerosol radiative forcing all-sky at SRF"         | "W/m^2"   | "aire[d=1]"  
+ swsrfcs_ant_forcing    | "swsrfcs_ant"           | "" | "swsrfcs_ant[d=1]"        | "Anthropogenic aerosol radiative forcing clear-sky at SRF"       | "W/m^2"   | "aire[d=1]" 
+ swsrfcf_ant_forcing    | "swsrfcf_ant"           | "" | "swsrfcf_ant[d=1]"        | "Anthropogenic aerosol impact on cloud radiative forcing at SRF" | "W/m^2"   | "aire[d=1]"
+ swtoaas_nat_forcing    | "swtoaas_nat"           | "" | "swtoaas_nat[d=1]"        | "Natural aerosol radiative forcing all-sky at TOA"               | "W/m^2"   | "aire[d=1]"   
+ swtoacs_nat_forcing    | "swtoacs_nat"           | "" | "swtoacs_nat[d=1]"        | "Natural aerosol radiative forcing clear-sky at TOA"             | "W/m^2"   | "aire[d=1]" 
+ swtoacf_nat_forcing    | "swtoacf_nat"           | "" | "swtoacf_nat[d=1]"        | "Natural aerosol impact on cloud radiative forcing at TOA"       | "W/m^2"   | "aire[d=1]" 
+ swsrfas_nat_forcing    | "swsrfas_nat"           | "" | "swsrfas_nat[d=1]"        | "Natural aerosol radiative forcing all-sky at SRF"               | "W/m^2"   | "aire[d=1]"  
+ swsrfcs_nat_forcing    | "swsrfcs_nat"           | "" | "swsrfcs_nat[d=1]"        | "Natural aerosol radiative forcing clear-sky at SRF"             | "W/m^2"   | "aire[d=1]" 
+ swsrfcf_nat_forcing    | "swsrfcf_nat"           | "" | "swsrfcf_nat[d=1]"        | "Natural aerosol impact on cloud radiative forcing at SRF"       | "W/m^2"   | "aire[d=1]"
+ topswai_forcing        | "topswai"               | "" | "topswai[d=1]"            | "AIE at TOA"                             | "W/m^2" | "aire[d=1]"
+ solswai_forcing        | "solswai"               | "" | "solswai[d=1]"            | "AIE at SRF"                             | "W/m^2" | "aire[d=1]"
+ swtoacf_zero_forcing   | "swtoacf_zero"          | "" | "swtoacf_zero[d=1]"       | "Cloud radiative forcing (allsky-clearsky fluxes) at TOA" | "W/m^2" | "aire[d=1]"
+ swsrfcf_zero_forcing   | "swsrfcf_zero"          | "" | "swsrfcf_zero[d=1]"       | "Cloud radiative forcing (allsky-clearsky fluxes) at SRF" | "W/m^2" | "aire[d=1]"
+ tro3_strato_forcing    | "tro3"                  | "" | "tro3[x=@din,y=@din,k=1:7@sum,d=1]"  | "ozone strato k=1:7@sum" | "Mole" | ""
+ tro3_tropo_forcing     | "tro3"                  | "" | "tro3[x=@din,y=@din,k=8:17@sum,d=1]" | "ozone trop k=8:17@sum"  | "Mole" | ""
+ solaire_forcing        | "solaire"               | "" | "solaire[d=1]"   | "Constante solaire"      | "W/m^2" | "1"
+ co2_ppm_forcing        | "co2_ppm"               | "" | "co2_ppm[d=1]"   | "Concentration du CO2"   | "ppm"   | "1"
+ CH4_ppb_forcing        | "CH4_ppb"               | "" | "CH4_ppb[d=1]"   | "Concentration du CH4"   | "ppb"   | "1"
+ N2O_ppb_forcing        | "N2O_ppb"               | "" | "N2O_ppb[d=1]"   | "Concentration du N2O"   | "ppb"   | "1"
+ CFC11_ppt_forcing      | "CFC11_ppt"             | "" | "CFC11_ppt[d=1]" | "Concentration du CFC11" | "ppt"   | "1"
+ CFC12_ppt_forcing      | "CFC12_ppt"             | "" | "CFC12_ppt[d=1]" | "Concentration du CFC12" | "ppt"   | "1"
+## R_ecc_forcing      | "R_ecc"                  | "" | "R_ecc[d=1]"   | "Excentricite"      | "-"     | "1"
+## R_peri_forcing     | "R_peri"                 | "" | "R_peri[d=1]"  | "Equinoxe"          | "-"     | "1"
+## R_incl_forcing     | "R_incl"                 | "" | "R_incl[d=1]"  | "Inclinaison"       | "deg"   | "1"
+#-----------------------------------------------------------------------------------------------------------------
+
Index: CONFIG/UNIFORM/v5/LMDZ_v5/README
===================================================================
--- CONFIG/UNIFORM/v5/LMDZ_v5/README	(revision 1934)
+++ CONFIG/UNIFORM/v5/LMDZ_v5/README	(revision 1934)
@@ -0,0 +1,13 @@
+How to create a job and run directory :
+
+Choose in EXPERIMENTS the experiment to be done and copy the corresponding config.card to modipsl/config/LMDZ_v5/, beside this README file. Modify the config.card as usual, for at least the JobName. Apply ../../util/ins_job and the creation of a new run directory will be done :
+  - the run directory will be created with the same name as JobName in config.card. If a directory with the same name already exists, nothing will be done.
+  - the PARAM, DRIVER and POST directories from GENERAL/ will be copied into the new run directory.
+  - the COMP directory from EXPERIMENTS/ will be copied into the new run directory, according to the variable ExpType in config.card.
+  - the config.card will be moved into the run directory.
+  - a job will be created in the run directory.
+
+Once the run directory has been created with ins_job, this directory is used as a normal run directory. More run directories can be created using ins_job several times. ins_job will never write in an existing run directory.
+
+
+For details see http://forge.ipsl.jussieu.fr/igcmg/wiki/Modipsl_exec
