Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/AA_make
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/AA_make	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/AA_make	(revision 1242)
@@ -0,0 +1,124 @@
+#- $Id$
+ifneq ($(MAKECMDGOALS),$(shell cat .resol))
+ifneq ($(MAKECMDGOALS),)
+$(shell rm -rf ../../modeles/INCA3/config)
+$(shell rm -f ../../modeles/INCA3/src/INCA_PP/*)
+endif
+endif
+
+
+all :
+	if [ -s ./.resol ] ; then gmake `cat .resol` ; else gmake AERxLMD9695 ; fi
+
+
+CH4_AERxLMD9672 : libioipsl inca3_lmdz9672_ch4aer
+	echo "CH4_AERxLMD9672" > .resol
+
+CH4xLMD9672 : libioipsl inca3_lmdz9672_ch4
+	echo "CH4xLMD9672" > .resol
+
+NMHCxLMD9672 : libioipsl inca3_lmdz9672_nmhc
+	echo "NMHCxLMD9672" > .resol
+
+NMHC_AERxLMD9672 : libioipsl inca3_lmdz9672_nmhcaer
+	echo "NMHC_AERxLMD9672" > .resol
+
+AERxLMD9672 : libioipsl inca3_lmdz9672_aer
+	echo "AERxLMD9672" > .resol
+
+GESxLMD9672 : libioipsl inca3_lmdz9672_ges
+	echo "GESxLMD9672" > .resol
+
+CH4_AERxLMD9695 : libioipsl inca3_lmdz9695_ch4aer
+	echo "CH4_AERxLMD9695" > .resol
+
+CH4xLMD9695 : libioipsl inca3_lmdz9695_ch4
+	echo "CH4xLMD9695" > .resol
+
+NMHCxLMD9695 : libioipsl inca3_lmdz9695_nmhc
+	echo "NMHCxLMD9695" > .resol
+
+NMHC_AERxLMD9695 : libioipsl inca3_lmdz9695_nmhcaer
+	echo "NMHC_AERxLMD9695" > .resol
+
+AERxLMD9695 : libioipsl inca3_lmdz9695_aer
+	echo "AERxLMD9695" > .resol
+
+GESxLMD9695 : libioipsl inca3_lmdz9695_ges
+	echo "GESxLMD9695" > .resol
+
+
+inca3_lmdz9672_aer : 
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x72x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie AER -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/AER/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4;  ./makelmdz_fcm -d 96x72x19 -v false -chimie INCA -parallel mpi -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x72x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9672_ch4aer :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x72x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie CH4_AER -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/CH4_AER/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x72x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x72x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+inca3_lmdz9672_ch4 :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x72x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie CH4 -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/CH4/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x72x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x72x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9672_nmhcaer :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x72x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie NMHC_AER -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/NMHC_AER/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x72x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x72x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9672_nmhc :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x72x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie NMHC -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/NMHC/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x72x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x72x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9672_ges :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x72x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie GES -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/GES/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x72x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x72x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9695_aer : 
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie AER -parallel mpi  -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/AER/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4;  ./makelmdz_fcm -d 96x95x19 -v false -chimie INCA -parallel mpi -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x95x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9695_ch4aer :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie CH4_AER -parallel mpi -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/CH4_AER/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x95x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x95x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+inca3_lmdz9695_ch4 :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie CH4 -parallel mpi -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/CH4/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x95x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x95x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9695_nmhcaer :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie NMHC_AER -parallel mpi -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/NMHC_AER/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x95x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x95x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9695_nmhc :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie NMHC -parallel mpi -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/NMHC/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x95x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x95x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+
+inca3_lmdz9695_ges :
+#	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_seq.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie GES -parallel mpi -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/GES/inca.dat ../../bin/inca.dat ; )
+	(cd ../../modeles/LMDZ4; ./makelmdz_fcm -d 96x95x19  -v false -chimie INCA -parallel mpi  -arch $(FCM_ARCH) -g reg gcm ; cp bin/gcm_96x95x19_phylmd_para_inca.e ../../bin/gcm.e ; )
+
+libioipsl :
+	(cd ../../modeles/IOIPSL/src ; $(M_K) -f Makefile)
+
+clean :
+	(rm -rf $(LIBDIR)/* ; rm -f $(BINDIR)/* ; rm -rf ../../modeles/LMDZ4/libo/* ../../modeles/LMDZ4/.lock ; )
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/AA_make.ldef
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/AA_make.ldef	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/AA_make.ldef	(revision 1242)
@@ -0,0 +1,14 @@
+#- $Id$
+#---------------------------------------------------------------------
+#- Creation des elements relatifs a LMDZ4INCA :
+#-      LMDZ4 (lib, executable : gcm.e)
+#-      INCA.2 (lib)
+#-      IOIPSL (lib + modules F90)
+#---------------------------------------------------------------------
+SHELL = /bin/sh
+#---------------------------------------------------------------------
+LIBDIR = ../../../../lib
+BINDIR = ../../../../bin
+MODDIR = $(LIBDIR)
+#---------------------------------------------------------------------
+#-
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/COMP/lmdz.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/COMP/lmdz.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/COMP/lmdz.card	(revision 1242)
@@ -0,0 +1,37 @@
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List=	()
+
+[BoundaryFiles]
+List=   (${R_INIT}/ATM/${config_UserChoices_TagName}/Amip6/amipbc_sst_360x180_${year}.nc, amipbc_sst_1x1.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Amip6/amipbc_sic_360x180_${year}.nc, amipbc_sic_1x1.nc)
+ListNonDel= (${R_INIT}/ATM/${config_UserChoices_TagName}/Albedo.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/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, .)
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def, .), \
+	(${SUBMIT_DIR}/PARAM/gcm.def, gcm.def), \
+	(${SUBMIT_DIR}/PARAM/run.def, .), \
+	(${SUBMIT_DIR}/PARAM/traceur.def, .) 
+
+[RestartFiles]
+List=	(start.nc, start.nc, start.nc), \
+	(startphy.nc, startphy.nc,  startphy.nc), \
+	(limit.nc, limit.nc,  limit.nc)
+
+[OutputText]
+List=   (ftrace.out.1.0)
+
+[OutputFiles]
+List=	()
+
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/COMP/lmdz.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/COMP/lmdz.driver	(revision 1242)
@@ -0,0 +1,133 @@
+#- $Id: lmdz.driver 868 2010-01-08 13:20:18Z acosce $
+#!/bin/ksh
+
+#-----------------------------------------------------------------
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' )
+    case ${RESOL_ATM} in
+	LMD7245)   BandsResol=72x45x19   ;;
+	LMD9671)   BandsResol=96x71x19   ;;
+	LMD9672)   BandsResol=96x72x19   ;;
+	LMD9695)   BandsResol=96x95x19   ;;
+	LMD14496)  BandsResol=144x96x19  ;;
+	LMD144142) BandsResol=144x142x19 ;;
+    esac
+
+
+    ##-- Calendar type for LMDZ
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=earth_366d;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d;;
+        360d)
+            CalendarTypeForLmdz=earth_360d;;
+        *)
+            CalendarTypeForLmdz=earth_360d
+    esac
+
+    ##- Default number of processor for lmdz
+    NUM_PROC_ATM=1
+    if [ X"${BATCH_NUM_PROC_TOT}" != X ] ; then
+	NUM_PROC_ATM=${BATCH_NUM_PROC_TOT}
+    fi
+
+    ##--Frequency purpose ....
+    ##--  Initialisation  ....
+    OK_instan=n
+    ##-- 
+    OK_journe=n
+    OK_mensuel=y
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) OK_mensuel=y ;;
+	5D|5d|1D|1d) OK_journe=n ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    5D|5d|1D|1d) OK_journe=n ;;
+	esac
+    done
+
+    ##--Variables used by LMDZ --
+    PAT_INST=$(     grep 'OK_instan'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_JOUR=$(     grep 'OK_journe'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_MOIS=$(     grep 'OK_mensuel'  ${SUBMIT_DIR}/PARAM/physiq.def )
+
+    PAT_iphysiq=$(  grep 'iphysiq'     ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_iperiod=$(  grep 'iperiod'     ${SUBMIT_DIR}/PARAM/gcm.def | tail -1)
+    PAT_day_step=$( grep 'day_step'    ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_ecritphy=$( grep 'ecritphy'    ${SUBMIT_DIR}/PARAM/gcm.def )
+
+
+
+    ##-- This could be define in lmdz.card, inside section [UserChoices]
+    ##-- Otherwise we get the value in *.def
+    iperiod=$(grep iperiod ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}' | tail -1)
+    iphysiq=$(grep iphysiq ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}')
+    # day_step : number of steps per day (multiple of iperiod)
+    (( day_step = 48 * iphysiq ))
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    ecritphy=${PeriodLengthInDays}
+
+    ##-- 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
+
+    ## Mise en forme du fichier physiq.def
+    sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   	\
+	-e "s/OK_journe=.*/OK_journe=${OK_journe}/"   	\
+	-e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" 	\
+	physiq.def > physiq.def.tmp
+    IGCM_sys_Mv physiq.def.tmp physiq.def
+
+    ## Mise en forme du fichier gcm.def
+    sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" 	\
+	-e "s/${PAT_day_step}/day_step=${day_step}/" 	\
+	-e "s/${PAT_iperiod}/iperiod=${iperiod}/" 	\
+	-e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" 	\
+	gcm.def > gcm.def.tmp
+    IGCM_sys_Mv gcm.def.tmp gcm.def
+
+    greg_dat=$( IGCM_date_ConvertFormatToGregorian $PeriodDateBegin ) 
+    jul_dat=$( IGCM_date_ConvertGregorianDateToJulian $greg_dat )
+    nbjour=$( expr \( $jul_dat \% 1000 \) ) 
+    yractu=$( expr $greg_dat / 10000  )
+    echo  "dayref = " $nbjour  "year actu = " $yractu
+
+    ## Mise en forme du fichier run.def
+    sed -e "s/_dayref_/${nbjour}/" 			\
+	-e "s/_anneeref_/${yractu}/" 			\
+        -e "s/_calend_/${CalendarTypeForLmdz}/"        	\
+	-e "s/_nday_/${PeriodLengthInDays}/" 		\
+	-e "s/_raz_date_/${RAZ_DATE}/" 			\
+	run.def > run.def.tmp
+    IGCM_sys_Mv run.def.tmp run.def
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/gcm.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/gcm.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/gcm.def	(revision 1242)
@@ -0,0 +1,74 @@
+## $Id: gcm.def 626 2009-05-15 15:04:37Z acosce $
+## 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
+## periode de la dissipation (en pas)
+idissip=10
+## 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=3600.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=3600.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## guide ou pas
+ok_guide=.false.
+## utilisation de restart de la version LMDZ3
+old_restart=
+## avec ou sans physique
+##physic=n
+iflag_phys=1
+## periode de la physique (en pas)
+iphysiq=10
+## frequence (en  jours ) de l'ecriture du fichier histphy               
+ecritphy=1
+## frequence (en jours ) de l'ecriture du fichier histmth
+ecrit_mth=1
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=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
+## Nombre  d'appels des routines de rayonnements ( par jour)
+nbapp_rad=12
+##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)
+iflag_con=30
+## 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/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/physiq.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/physiq.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/physiq.def	(revision 1242)
@@ -0,0 +1,92 @@
+##- $Id: physiq.def 607 2009-04-21 11:08:53Z acosce $
+type_ocean=force 
+VEGET=n
+OK_journe=n
+OK_mensuel=y
+OK_instan=n
+if_ebil=0
+#
+# parametres KE
+#
+epmax = .99
+ok_adj_ema = n
+iflag_clw = 1
+# 
+# parametres nuages
+#
+cld_lc_lsc = 0.00026
+cld_lc_con = 0.00026
+cld_tau_lsc = 3600.
+cld_tau_con = 3600.
+ffallv_lsc = 1.
+ffallv_con = 1.
+coef_eva = 0.00002
+reevap_ice = y
+iflag_cldcon = 3
+iflag_pdf = 1
+fact_cldcon = 1.
+#facttemps = 1.e-4
+facttemps = 0.0001
+ok_newmicro = y
+ratqsbas = 0.005
+ratqshaut = 0.33
+rad_froid = 35
+rad_chau1=12
+rad_chau2=11
+#ksta_ter=1.e-7
+ksta_ter=0.0000001
+#
+# parametres outputs
+#
+#niveau de sortie "hf" lev_histhf avec
+# - lev_histhf=3 => defaut
+# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele
+lev_histhf=3
+#niveau de sortie "day" lev_histday
+# - lev_histday=2 => defaut
+# - lev_histday=3 => + champs 3D => F. Lott
+# - lev_histday=4 => + champs sous-surfaces
+lev_histday=2
+#niveau de sortie "mth" lev_histmth avec
+# - lev_histmth=2 => defaut
+# - lev_histmth=3 => albedo, rugosite sous-surfaces
+# - lev_histmth=4 => champs tendances 3d
+lev_histmth=4
+#
+# parametres climatique
+#
+R_ecc = 0.016715
+R_peri = 102.7
+R_incl = 23.441
+solaire = 1365.
+co2_ppm = 348.
+CH4_ppb = 1650.
+N2O_ppb = 306.
+CFC11_ppt = 280.
+CFC12_ppt = 484.
+#
+# parametres simulateur ISCCP
+#
+top_height = 3
+#overlap = 1, 2 ou 3
+overlap = 3
+#cdmmax
+#cdmmax = 2.5E-3
+cdmmax = 0.0025
+#cdhmax
+#cdhmax = 2.0E-3
+cdhmax = 0.002
+#autres flags
+iflag_pbl = 1
+iflag_thermals = 0
+iflag_clos=0
+iflag_ratqs=0
+#
+# parametre chimie/aerosols
+#
+config_inca=none
+aerosol_couple=n
+flag_aerosol=1
+new_aod=.true.
+ok_ade=n
+ok_aie=n
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/run.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/run.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/run.def	(revision 1242)
@@ -0,0 +1,23 @@
+## $Id: run.def 859 2009-12-18 14:26:57Z acosce $
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+## Type de calendrier utilise
+## valeur possible: earth_360d (defaut), earth_365d, earth_366d
+calend=_calend_
+## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+dayref=_dayref_
+##  Annee de l'etat  initial (   avec  4  chiffres   )
+anneeref=_anneeref_
+## raz_date
+raz_date=_raz_date_
+## Nombre de jours d'integration
+nday=_nday_
+## periode de sortie des variables de controle (en pas)
+iconser= 480
+## periode d'ecriture du fichier histoire (en jour)
+iecri=10
+## periode de stockage fichier histmoy (en jour)
+periodav=1.
+## unite de sortie des impressions
+lunout=6
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/traceur.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/PARAM/traceur.def	(revision 1242)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/config.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/config.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/CREATE/config.card	(revision 1242)
@@ -0,0 +1,76 @@
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=CREATE_L
+LongName='Call for create_etat0_limit for AMIP run.'
+TagName=LMDZINCA
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Début et fin de Job
+#D-- "YYYY-MM-DD"
+DateBegin=2000-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, LMDZ4_V3_1)
+
+#========================================================================
+#D-- Executable -
+[Executable]
+Name=create_etat0_limit.e
+#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 (=> JOB_OS)
+RestartJobName=
+#D- Path Server Group Login (=> PSGL)
+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=
+MonitoringFrequency=
+TimeSeriesFrequency=
+SeasonalFrequency=
+#========================================================================
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/inca.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/inca.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/inca.card	(revision 1242)
@@ -0,0 +1,55 @@
+#- $Id$
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+	
+
+[BoundaryFiles]
+List=	    (/dmnfs/cont003/p24data/ECMWF320x160/AN${year}/165_${year}${month}.nc      , u10mec.nc)\
+	    (/dmnfs/cont003/p24data/ECMWF320x160/AN${year}/166_${year}${month}.nc      , v10mec.nc)\
+	    (/dmnfs/cont003/p24data/ECMWF96x72/AN${year}/u10m_ecmwf_${year}${month}.nc , u10m.nc  )\
+	    (/dmnfs/cont003/p24data/ECMWF96x72/AN${year}/v10m_ecmwf_${year}${month}.nc , v10m.nc  )
+
+ListNonDel= (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/oxydants.nc             ,   oxydants.nc    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/npp.nc                  ,   npp.nc         )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/aircraft.nc             ,   aircraft.nc    )\
+	    (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/sflx_lmdz_phy.present.nc,   sflx.nc	    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/phototable.dat_V4.2	                 ,   phototable.dat )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/landuse.nc              ,   landuse.nc     )\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/rhvlmdz.txt             ,   rhv.dat        )\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/wthlmdz.txt             ,   wth.dat	    )\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/clylmdz.txt             ,   cly.dat	    )\
+            (${R_INIT}/CHM/LMDZINCA/AER/rhvEC11.09092002.txt	                 ,   rhvEC.txt	    )\
+            (${R_INIT}/CHM/LMDZINCA/AER/wthEC11.09092002.txt	                 ,   wthEC.txt	    )\
+            (${R_INIT}/CHM/LMDZINCA/AER/clyEC11.09092002.txt	                 ,   clyEC.txt	    )
+
+[ParametersFiles]
+List=	
+
+[RestartFiles]
+# List restart that have to be saved/restored each loop (file out, saved, and in) :
+List= (precipinsoil.dat, precipinsoil.dat, precipinsoil.dat)
+
+
+
+[OutputText]
+List=  
+# avec la // : out_orchidee_*
+
+[OutputFiles]
+List=	(inca_avgr.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc,  Post_1M_inca_avgr  ),\
+	(forcage.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_forcage.nc,      NONE)
+
+
+[Post_1M_inca_avgr]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter)
+TimeSeriesVars2D = (OD550_ASBCM, OD550_ASPOMM, OD550_ASSO4M, OD550_CSSO4M, OD550_SSSSM, OD550_ASSSM, OD550_CSSSM, OD550_CIDUSTM, OD550_AIBCM, OD550_AIPOMM)
+ChunckJob2D = NONE
+TimeSeriesVars3D = (CH4,NO,NO2,CO,O3,HNO3)
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/inca.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/inca.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/inca.driver	(revision 1242)
@@ -0,0 +1,40 @@
+#- $Id$
+#!/bin/ksh
+
+#D- Driver du script pour INCA
+
+#-----------------------------------------------------------------
+function CHM_Initialize
+{
+    IGCM_debug_PushStack "CHM_Initialize"
+
+    RESOL_CHM=$( echo $RESOL_ATM | awk "-FD" '{print $2}' )
+    echo 'resol chimie', $RESOL_CHM
+    echo "INCA Tag : LMDZ-INCA" ${config_UserChoices_TagName}
+
+    IGCM_debug_PopStack "CHM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function CHM_Update
+{
+    IGCM_debug_PushStack "CHM_Update"
+    IGCM_debug_PopStack "CHM_Update"
+}
+
+#-----------------------------------------------------------------
+function CHM_Finalize
+{
+    IGCM_debug_PushStack "CHM_Finalize"
+
+    echo FINALIZE CHM !!!
+#
+#    IGCM_sys_Mkdir ${R_OUT_EXE}
+#    cd ${MODIPSL}
+#    tar cf src_${config_UserChoices_JobName}.tar  modeles
+#    IGCM_sys_Mv src_${config_UserChoices_JobName}.tar ${R_OUT_EXE}/.
+#    cd -
+#
+
+    IGCM_debug_PopStack "CHM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/lmdz.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/lmdz.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/lmdz.card	(revision 1242)
@@ -0,0 +1,46 @@
+#- $Id$
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+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
+
+[InitialStateFiles]
+# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job
+List= 	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/start.nc, start.nc) \
+	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/startphy.nc, startphy.nc) \
+	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/precipinsoil.dat, precipinsoil.dat) 
+
+
+[BoundaryFiles]
+List=   (${R_BC}/CHM/LMDZORINCA/${RESOL_ATM}_leap/EL-${RESOL_ATM}_${year}1231_limit.nc , limit.nc) 
+ListNonDel= ()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def, .), \
+	(${SUBMIT_DIR}/PARAM/gcm.def, gcm.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, Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat)
+
+
+[OutputFiles]
+List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,      Post_1M_histmth), \
+
+[Post_1M_histmth]
+Patches= (Patch_20091030_histcom_time_axis)
+GatherWithInternal = (lon, lat, presnivs, time_counter, aire)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux_ter, taux_sic, taux_lic, taux_oce, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, 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)
+ChunckJob2D = NONE
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/lmdz.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/COMP/lmdz.driver	(revision 1242)
@@ -0,0 +1,160 @@
+#- $Id$
+#!/bin/ksh
+
+#-----------------------------------------------------------------
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' )
+    case ${RESOL_ATM} in
+	LMD7245)   BandsResol=72x45x19   ;;
+	LMD9671)   BandsResol=96x71x19   ;;
+	LMD9672)   BandsResol=96x72x19   ;;
+	LMD9695)   BandsResol=96x95x19   ;;
+	LMD14496)  BandsResol=144x96x19  ;;
+	LMD144142) BandsResol=144x142x19 ;;
+    esac
+
+
+    ##-- Calendar type for LMDZ
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=earth_366d;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d;;
+        360d)
+            CalendarTypeForLmdz=earth_360d;;
+        *)
+            CalendarTypeForLmdz=earth_360d
+    esac
+
+    ##- Default number of processor for lmdz
+    NUM_PROC_ATM=1
+    if [ X"${BATCH_NUM_PROC_TOT}" != X ] ; then
+	NUM_PROC_ATM=${BATCH_NUM_PROC_TOT}
+    fi
+
+    ##--Frequency purpose ....
+    ##--  Initialisation  ....
+    OK_instan=n
+    ##-- 
+    OK_journe=n
+    OK_mensuel=y
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) OK_mensuel=y ;;
+	5D|5d|1D|1d) OK_journe=n ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    5D|5d|1D|1d) OK_journe=n ;;
+	esac
+    done
+
+    ##--Variables used by LMDZ --
+    PAT_INST=$(     grep 'OK_instan'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_JOUR=$(     grep 'OK_journe'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_MOIS=$(     grep 'OK_mensuel'  ${SUBMIT_DIR}/PARAM/physiq.def )
+
+    PAT_iphysiq=$(  grep 'iphysiq'     ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_iperiod=$(  grep 'iperiod'     ${SUBMIT_DIR}/PARAM/gcm.def | tail -1)
+    PAT_day_step=$( grep 'day_step'    ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_ecritphy=$( grep 'ecritphy'    ${SUBMIT_DIR}/PARAM/gcm.def )
+
+
+    ## Read LMDZ_NbPeriod_adjust option in lmdz.card 
+    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_NbPeriod_adjust
+    eval LMDZ_NbPeriod_adjust=\${lmdz_UserChoices_LMDZ_NbPeriod_adjust} > /dev/null 2>&1
+
+    if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+        IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_Bands_file_name
+        eval LMDZ_Bands_file_name=\${lmdz_UserChoices_LMDZ_Bands_file_name} > /dev/null 2>&1
+    fi
+
+
+    ##-- This could be define in lmdz.card, inside section [UserChoices]
+    ##-- Otherwise we get the value in *.def
+    iperiod=$(grep iperiod ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}' | tail -1)
+    iphysiq=$(grep iphysiq ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}')
+    # day_step : number of steps per day (multiple of iperiod)
+    (( day_step = 48 * iphysiq ))
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    ecritphy=${PeriodLengthInDays}
+
+    ##-- 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
+
+    if [ ${CumulPeriod} -eq 1 ] ; then
+        if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+            IGCM_sys_Get ${LMDZ_Bands_file_name} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+            IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0
+            IGCM_debug_Print 1 "Bands file forced to ${LMDZ_Bands_file_name} and stored in ${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0"
+        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 )) 
+
+    [ ${CumulPeriod} -gt 1 ] && IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+
+    ## Mise en forme du fichier physiq.def
+    sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   	\
+	-e "s/OK_journe=.*/OK_journe=${OK_journe}/"   	\
+	-e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" 	\
+	physiq.def > physiq.def.tmp
+    IGCM_sys_Mv physiq.def.tmp physiq.def
+
+    ## Mise en forme du fichier gcm.def
+    sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" 	\
+	-e "s/${PAT_day_step}/day_step=${day_step}/" 	\
+	-e "s/${PAT_iperiod}/iperiod=${iperiod}/" 	\
+	-e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" 	\
+	gcm.def > gcm.def.tmp
+    IGCM_sys_Mv gcm.def.tmp gcm.def
+
+    greg_dat=$( IGCM_date_ConvertFormatToGregorian $PeriodDateBegin ) 
+    jul_dat=$( IGCM_date_ConvertGregorianDateToJulian $greg_dat )
+    nbjour=$( expr \( $jul_dat \% 1000 \) ) 
+    yractu=$( expr $greg_dat / 10000  )
+    echo  "dayref = " $nbjour  "year actu = " $yractu
+
+    ## Mise en forme du fichier run.def
+    sed -e "s/_dayref_/${nbjour}/" 			\
+	-e "s/_anneeref_/${yractu}/" 			\
+        -e "s/_calend_/${CalendarTypeForLmdz}/"        	\
+	-e "s/_nday_/${PeriodLengthInDays}/" 		\
+	-e "s/_raz_date_/${RAZ_DATE}/" 			\
+	-e "s/_adjust_/${LMDZ_adjust}/"                	\
+	run.def > run.def.tmp
+    IGCM_sys_Mv run.def.tmp run.def
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/Bands_96x72x19_2prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
@@ -0,0 +1,8 @@
+ 0  37
+ 1  36
+ 0  37
+ 1  36
+ 0  36
+ 1  37
+ 0  3410
+ 1  3408
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/Bands_96x72x19_4prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
@@ -0,0 +1,16 @@
+ 0  8
+ 1  28
+ 2  28
+ 3  9
+ 0  19
+ 1  18
+ 2  18
+ 3  18
+ 0  10
+ 1  52
+ 2  8
+ 3  3
+ 0  1701
+ 1  1716
+ 2  1705
+ 3  1696
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/gcm.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/gcm.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/gcm.def	(revision 1242)
@@ -0,0 +1,74 @@
+## $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
+## periode de la dissipation (en pas)
+idissip=10
+## 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=3600.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=3600.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## guide ou pas
+ok_guide=.false.
+## utilisation de restart de la version LMDZ3
+old_restart=
+## avec ou sans physique
+##physic=n
+iflag_phys=1
+## periode de la physique (en pas)
+iphysiq=10
+## frequence (en  jours ) de l'ecriture du fichier histphy               
+ecritphy=1
+## frequence (en jours ) de l'ecriture du fichier histmth
+ecrit_mth=1
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=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
+## Nombre  d'appels des routines de rayonnements ( par jour)
+nbapp_rad=12
+##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)
+iflag_con=30
+## 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/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/offline.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/offline.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/offline.def	(revision 1242)
@@ -0,0 +1,10 @@
+## $Id$
+T
+4
+T
+-2.
+48.1
+1
+T
+6
+2
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/orchidee.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/orchidee.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/orchidee.def	(revision 1242)
@@ -0,0 +1,44 @@
+## $Id$
+#
+# SECHIBA
+#
+STOMATE_OK_CO2=FALSE
+# STOMATE_OK_STOMATE is not set
+# STOMATE_OK_DGVM is not set
+# STOMATE_WATCHOUT is not set
+#SECHIBA_restart_in=default
+SECHIBA_restart_in=start_sech.nc
+SECHIBA_rest_out=sechiba_rest.nc
+SECHIBA_reset_time=y
+SECHIBA_reset_time is not set
+OUTPUT_FILE=sechiba_out.nc
+WRITE_STEP=2592000
+SECHIBA_HISTLEVEL=10
+STOMATE_OUTPUT_FILE=stomate_history.nc
+STOMATE_HIST_DT=10.
+STOMATE_HISTLEVEL=0
+SECHIBA_DAY=0.0
+SECHIBA_ZCANOP=0.5
+DT_SLOW=86400.
+# IMPOSE_VEG is not set
+VEGETATION_FILE=carteveg5km.nc
+DIFFUCO_LEAFCI=233.
+CONDVEG_SNOWA=default
+# IMPOSE_AZE is not set
+SOILALB_FILE=soils_param.nc
+SOILTYPE_FILE=soils_param.nc 
+ENERBIL_TSURF=280.
+HYDROL_SNOW=0.0
+HYDROL_SNOWAGE=0.0
+HYDROL_SNOWICE=0.0
+HYDROL_SNOWICEAGE=0.0
+HYDROL_HDRY=1.0
+HYDROL_HUMR=1.0
+HYDROL_BQSB=default
+HYDROL_GQSB=0.0
+HYDROL_DSG=0.0
+HYDROL_DSP=default
+HYDROL_QSV=0.0
+THERMOSOIL_TPRO=280.
+RIVER_ROUTING=y
+ROUTING_FILE=routing.nc
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/physiq.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/physiq.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/physiq.def	(revision 1242)
@@ -0,0 +1,92 @@
+##- $Id$
+type_ocean=force 
+VEGET=n
+OK_journe=n
+OK_mensuel=y
+OK_instan=n
+if_ebil=0
+#
+# parametres KE
+#
+epmax = .99
+ok_adj_ema = n
+iflag_clw = 1
+# 
+# parametres nuages
+#
+cld_lc_lsc = 0.00026
+cld_lc_con = 0.00026
+cld_tau_lsc = 3600.
+cld_tau_con = 3600.
+ffallv_lsc = 1.
+ffallv_con = 1.
+coef_eva = 0.00002
+reevap_ice = y
+iflag_cldcon = 3
+iflag_pdf = 1
+fact_cldcon = 1.
+#facttemps = 1.e-4
+facttemps = 0.0001
+ok_newmicro = y
+ratqsbas = 0.005
+ratqshaut = 0.33
+rad_froid = 35
+rad_chau1=12
+rad_chau2=11
+#ksta_ter=1.e-7
+ksta_ter=0.0000001
+#
+# parametres outputs
+#
+#niveau de sortie "hf" lev_histhf avec
+# - lev_histhf=3 => defaut
+# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele
+lev_histhf=3
+#niveau de sortie "day" lev_histday
+# - lev_histday=2 => defaut
+# - lev_histday=3 => + champs 3D => F. Lott
+# - lev_histday=4 => + champs sous-surfaces
+lev_histday=2
+#niveau de sortie "mth" lev_histmth avec
+# - lev_histmth=2 => defaut
+# - lev_histmth=3 => albedo, rugosite sous-surfaces
+# - lev_histmth=4 => champs tendances 3d
+lev_histmth=4
+#
+# parametres climatique
+#
+R_ecc = 0.016715
+R_peri = 102.7
+R_incl = 23.441
+solaire = 1365.
+co2_ppm = 348.
+CH4_ppb = 1650.
+N2O_ppb = 306.
+CFC11_ppt = 280.
+CFC12_ppt = 484.
+#
+# parametres simulateur ISCCP
+#
+top_height = 3
+#overlap = 1, 2 ou 3
+overlap = 3
+#cdmmax
+#cdmmax = 2.5E-3
+cdmmax = 0.0025
+#cdhmax
+#cdhmax = 2.0E-3
+cdhmax = 0.002
+#autres flags
+iflag_pbl = 1
+iflag_thermals = 0
+iflag_clos=0
+iflag_ratqs=0
+#
+# parametre chimie/aerosols
+#
+config_inca=aero
+aerosol_couple=n
+flag_aerosol=1
+new_aod=.true.
+ok_ade=n
+ok_aie=n
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/run.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/run.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/run.def	(revision 1242)
@@ -0,0 +1,25 @@
+## $Id$
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+## activation du calcul d equilibrage de charge
+adjust=_adjust_
+## Type de calendrier utilise
+## valeur possible: earth_360d (defaut), earth_365d, earth_366d
+calend=_calend_
+## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+dayref=_dayref_
+##  Annee de l'etat  initial (   avec  4  chiffres   )
+anneeref=_anneeref_
+## raz_date
+raz_date=_raz_date_
+## Nombre de jours d'integration
+nday=_nday_
+## periode de sortie des variables de controle (en pas)
+iconser= 480
+## periode d'ecriture du fichier histoire (en jour)
+iecri=10
+## periode de stockage fichier histmoy (en jour)
+periodav=1.
+## unite de sortie des impressions
+lunout=6
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/traceur.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/PARAM/traceur.def	(revision 1242)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/POST/monitoring01_inca.cfg
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/POST/monitoring01_inca.cfg	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/POST/monitoring01_inca.cfg	(revision 1242)
@@ -0,0 +1,51 @@
+#**************************************************************
+# 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
+#   - files patterns, operations,title,units must be enclosed with character '"'
+#   - 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=F0D5F4
+smooth=12
+#-----------------------------------------------------------------------------------------------------------------
+#  field | files patterns | files additionnal | operations | title | units | calcul of area
+#-----------------------------------------------------------------------------------------------------------------
+ OD550_AERO_global | "OD550_ASBCM OD550_ASPOMM OD550_ASSO4M OD550_CSSO4M OD550_SSSSM OD550_ASSSM OD550_CSSSM OD550_CIDUSTM OD550_AIBCM OD550_AIPOMM" | LMDZ4.0_9695_grid.nc | "(OD550_ASBCM[d=1]+ OD550_ASPOMM[d=2]+ OD550_ASSO4M[d=3]+ OD550_CSSO4M[d=4]+ OD550_SSSSM[d=5]+ OD550_ASSSM[d=6]+ OD550_CSSSM[d=7]+ OD550_CIDUSTM[d=8]+ OD550_AIBCM[d=9]+ OD550_AIPOMM[d=10])" | "total  aerosol optical depth @550nm"  |"(-)" | "aire[d=11]"
+NOX_surf_global		| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=19]+NO2[d=2,k=19])" | " NOX a la surface"   | "VMR" | "aire[d=3]"
+CH4_surf_global		| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=19]"                | " CH4 a la surface"   | "VMR" | "aire[d=2]"
+CO_surf_global		| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=19]"                 | " CO a la surface"    | "VMR" | "aire[d=2]"
+O3_surf_global		| "_O3"   	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=19]"                 | " O3 a la surface"    | "VMR" | "aire[d=2]"
+HNO3_surf_global	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=19]"               | " HNO3 a la surface"  | "?" 	 | "aire[d=2]"
+NOX_545HPa_global    	| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=12]+NO2[d=2,k=12])" | " NOX a 545 HPa"      | "VMR" | "aire[d=3]"
+CH4_545HPa_global    	| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=12]"                | " CH4 a 545 HPa"      | "VMR" | "aire[d=2]"
+CO_545HPa_global     	| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=12]"                 | " CO a 545 HPa"       | "VMR" | "aire[d=2]"
+O3_545HPa_global     	| "_O3"    	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=12]"                 | " O3 a 545 HPa"       | "VMR" | "aire[d=2]"
+HNO3_545HPa_global   	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=12]"               | " HNO3 a 545 HPa"     | "?"   | "aire[d=2]"
+
+#-----------------------------------------------------------------------------------------------------------------
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/config.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/config.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_AER/config.card	(revision 1242)
@@ -0,0 +1,106 @@
+#- $Id$
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=LI_A_00
+LongName='LMDz4INCA with New scripts'
+TagName=LMDZINCA/AER
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Début et fin de Job
+#D-- "YYYY-MM-DD"
+DateBegin=2000-01-01
+DateEnd=2000-12-31
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=1M
+#============================
+#-- PBS Class (required on brodie mono/multi)
+#JobClass=mono
+JobClass=multi
+#============================
+#-- Total Number of Processors (minimum is 2 for a coupled configuration)
+#JobNumProcTot=2
+JobNumProcTot=4
+#============================
+#-- Run Options : here for oasis coupling
+JobRunOptions='"-np "${config_UserChoices_JobNumProcTot}'
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ4)
+CHM= (inca, INCA3) 
+
+#========================================================================
+#D-- Executable -
+[Executable]
+Name=gcm.e
+#D- For each component, Real name of executable, Name of executable for oasis
+ATM= (gcm.e, gcm.e)
+CHM= (inca.dat, inca.dat)
+
+#========================================================================
+#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=1869-12-30
+#D- Define restart simulation name (=> JOB_OS)
+RestartJobName=JobName
+#D- Path Server Group Login (=> PSGL)
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/AER/
+# Attention login depend de la machine
+
+#========================================================================
+#D-- Post -
+[Post]
+#D- Do we rebuild parallel output, this flag determines
+#D- frequency of rebuild submission (use NONE for DRYRUN=3)
+RebuildFrequency=1M
+#D- Do we rebuild parallel output from archive
+RebuildFromArchive=true
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+TimeSeriesFrequency=1M
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average
+SeasonalFrequency=5Y
+
+#========================================================================
+#D-- ATM -
+[ATM]
+#
+WriteFrequency="1M 1D HF"
+# If config_Restarts_OverRule == 'n' all params are read
+Restart=n
+# Last day of the experience used as restart
+RestartDate=2000-01-01
+# Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/AER/
+# Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
+#D-- CHM -
+[CHM]
+WriteFrequency="1M"
+Restart=n
+#D-- Last day of the experience used as restart
+RestartDate=2000-01-01
+#D-- Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/AER/
+#D-- Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/inca.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/inca.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/inca.card	(revision 1242)
@@ -0,0 +1,44 @@
+# $Id: inca.card 73 2007-06-21 14:42:36Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+
+[BoundaryFiles]
+List=	() 
+
+ListNonDel= (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/o3clim.nc               ,   o3clim.nc      )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/npp.nc	                 ,   npp.nc	    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/aircraft.nc             ,   aircraft.nc    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/so4.nc	                 ,   so4.nc	    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/phototable.dat_V4.2	                 ,   phototable.dat )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/landuse.nc              ,   landuse.nc     )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/sflx_Present.nc         ,   sflx.nc        ) 
+
+[ParametersFiles]
+List=	
+
+[RestartFiles]
+# List restart that have to be saved/restored each loop (file out, saved, and in) :
+List= (NONE)
+
+
+
+[OutputText]
+List=  
+# avec la // : out_orchidee_*
+
+[OutputFiles]
+List=	(inca_avgr.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc,    Post_1M_inca_avgr)
+
+
+[Post_1M_inca_avgr]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = (CH4,NO,NO2,CO,O3,HNO3)
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/inca.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/inca.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/inca.driver	(revision 1242)
@@ -0,0 +1,40 @@
+#- $Id: inca.driver 119 2007-07-31 13:56:33Z acosce $
+#!/bin/ksh
+
+#D- Driver du script pour INCA
+
+#-----------------------------------------------------------------
+function CHM_Initialize
+{
+    IGCM_debug_PushStack "CHM_Initialize"
+
+    RESOL_CHM=$( echo $RESOL_ATM | awk "-FD" '{print $2}' )
+    echo 'resol chimie', $RESOL_CHM
+    echo "INCA Tag : LMDZ-INCA" ${config_UserChoices_TagName}
+
+    IGCM_debug_PopStack "CHM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function CHM_Update
+{
+    IGCM_debug_PushStack "CHM_Update"
+    IGCM_debug_PopStack "CHM_Update"
+}
+
+#-----------------------------------------------------------------
+function CHM_Finalize
+{
+    IGCM_debug_PushStack "CHM_Finalize"
+
+    echo FINALIZE CHM !!!
+#
+#    IGCM_sys_Mkdir ${R_OUT_EXE}
+#    cd ${MODIPSL}
+#    tar cf src_${config_UserChoices_JobName}.tar  modeles
+#    IGCM_sys_Mv src_${config_UserChoices_JobName}.tar ${R_OUT_EXE}/.
+#    cd -
+#
+
+    IGCM_debug_PopStack "CHM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/lmdz.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/lmdz.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/lmdz.card	(revision 1242)
@@ -0,0 +1,44 @@
+# $Id: lmdz.card 73 2007-06-21 14:42:36Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+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
+
+[InitialStateFiles]
+# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job
+List= 	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/start.nc, start.nc) \
+	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/startphy.nc, startphy.nc)
+
+[BoundaryFiles]
+List=   (${R_BC}/CHM/LMDZORINCA/${RESOL_ATM}_leap/EL-${RESOL_ATM}_${year}1231_limit.nc , limit.nc) 
+ListNonDel= ()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def, .), \
+	(${SUBMIT_DIR}/PARAM/gcm.def, gcm.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, Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat)
+
+
+[OutputFiles]
+List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,  Post_1M_histmth   )	
+
+[Post_1M_histmth]
+Patches= (Patch_20091030_histcom_time_axis)
+GatherWithInternal = (lon, lat, presnivs, time_counter, aire)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux_ter, taux_sic, taux_lic, taux_oce, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, 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)
+ChunckJob2D = NONE
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/lmdz.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/COMP/lmdz.driver	(revision 1242)
@@ -0,0 +1,160 @@
+#- $Id: lmdz.driver 841 2009-12-08 09:48:13Z acosce $
+#!/bin/ksh
+
+#-----------------------------------------------------------------
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' )
+    case ${RESOL_ATM} in
+	LMD7245)   BandsResol=72x45x19   ;;
+	LMD9671)   BandsResol=96x71x19   ;;
+	LMD9672)   BandsResol=96x72x19   ;;
+	LMD9695)   BandsResol=96x95x19   ;;
+	LMD14496)  BandsResol=144x96x19  ;;
+	LMD144142) BandsResol=144x142x19 ;;
+    esac
+
+
+    ##-- Calendar type for LMDZ
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=earth_366d;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d;;
+        360d)
+            CalendarTypeForLmdz=earth_360d;;
+        *)
+            CalendarTypeForLmdz=earth_360d
+    esac
+
+    ##- Default number of processor for lmdz
+    NUM_PROC_ATM=1
+    if [ X"${BATCH_NUM_PROC_TOT}" != X ] ; then
+	NUM_PROC_ATM=${BATCH_NUM_PROC_TOT}
+    fi
+
+    ##--Frequency purpose ....
+    ##--  Initialisation  ....
+    OK_instan=n
+    ##-- 
+    OK_journe=n
+    OK_mensuel=y
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) OK_mensuel=y ;;
+	5D|5d|1D|1d) OK_journe=n ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    5D|5d|1D|1d) OK_journe=n ;;
+	esac
+    done
+
+    ##--Variables used by LMDZ --
+    PAT_INST=$(     grep 'OK_instan'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_JOUR=$(     grep 'OK_journe'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_MOIS=$(     grep 'OK_mensuel'  ${SUBMIT_DIR}/PARAM/physiq.def )
+
+    PAT_iphysiq=$(  grep 'iphysiq'     ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_iperiod=$(  grep 'iperiod'     ${SUBMIT_DIR}/PARAM/gcm.def | tail -1)
+    PAT_day_step=$( grep 'day_step'    ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_ecritphy=$( grep 'ecritphy'    ${SUBMIT_DIR}/PARAM/gcm.def )
+
+
+    ## Read LMDZ_NbPeriod_adjust option in lmdz.card 
+    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_NbPeriod_adjust
+    eval LMDZ_NbPeriod_adjust=\${lmdz_UserChoices_LMDZ_NbPeriod_adjust} > /dev/null 2>&1
+
+    if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+        IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_Bands_file_name
+        eval LMDZ_Bands_file_name=\${lmdz_UserChoices_LMDZ_Bands_file_name} > /dev/null 2>&1
+    fi
+
+
+    ##-- This could be define in lmdz.card, inside section [UserChoices]
+    ##-- Otherwise we get the value in *.def
+    iperiod=$(grep iperiod ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}' | tail -1)
+    iphysiq=$(grep iphysiq ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}')
+    # day_step : number of steps per day (multiple of iperiod)
+    (( day_step = 48 * iphysiq ))
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    ecritphy=${PeriodLengthInDays}
+
+    ##-- 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
+
+    if [ ${CumulPeriod} -eq 1 ] ; then
+        if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+            IGCM_sys_Get ${LMDZ_Bands_file_name} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+            IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0
+            IGCM_debug_Print 1 "Bands file forced to ${LMDZ_Bands_file_name} and stored in ${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0"
+        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 )) 
+
+    [ ${CumulPeriod} -gt 1 ] && IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+
+    ## Mise en forme du fichier physiq.def
+    sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   	\
+	-e "s/OK_journe=.*/OK_journe=${OK_journe}/"   	\
+	-e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" 	\
+	physiq.def > physiq.def.tmp
+    IGCM_sys_Mv physiq.def.tmp physiq.def
+
+    ## Mise en forme du fichier gcm.def
+    sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" 	\
+	-e "s/${PAT_day_step}/day_step=${day_step}/" 	\
+	-e "s/${PAT_iperiod}/iperiod=${iperiod}/" 	\
+	-e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" 	\
+	gcm.def > gcm.def.tmp
+    IGCM_sys_Mv gcm.def.tmp gcm.def
+
+    greg_dat=$( IGCM_date_ConvertFormatToGregorian $PeriodDateBegin ) 
+    jul_dat=$( IGCM_date_ConvertGregorianDateToJulian $greg_dat )
+    nbjour=$( expr \( $jul_dat \% 1000 \) ) 
+    yractu=$( expr $greg_dat / 10000  )
+    echo  "dayref = " $nbjour  "year actu = " $yractu
+
+    ## Mise en forme du fichier run.def
+    sed -e "s/_dayref_/${nbjour}/" 			\
+	-e "s/_anneeref_/${yractu}/" 			\
+        -e "s/_calend_/${CalendarTypeForLmdz}/"        	\
+	-e "s/_nday_/${PeriodLengthInDays}/" 		\
+	-e "s/_raz_date_/${RAZ_DATE}/" 			\
+	-e "s/_adjust_/${LMDZ_adjust}/"                	\
+	run.def > run.def.tmp
+    IGCM_sys_Mv run.def.tmp run.def
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/Bands_96x72x19_2prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
@@ -0,0 +1,8 @@
+ 0  37
+ 1  36
+ 0  37
+ 1  36
+ 0  36
+ 1  37
+ 0  3539
+ 1  3279
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/Bands_96x72x19_4prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
@@ -0,0 +1,16 @@
+ 0  8
+ 1  29
+ 2  27
+ 3  9
+ 0  19
+ 1  18
+ 2  18
+ 3  18
+ 0  13
+ 1  34
+ 2  23
+ 3  3
+ 0  1848
+ 1  1682
+ 2  1670
+ 3  1618
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/gcm.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/gcm.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/gcm.def	(revision 1242)
@@ -0,0 +1,74 @@
+## $Id: gcm.def 74 2007-06-21 14:46:35Z acosce $
+## 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
+## periode de la dissipation (en pas)
+idissip=10
+## 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=3600.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=3600.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## guide ou pas
+ok_guide=.false.
+## utilisation de restart de la version LMDZ3
+old_restart=
+## avec ou sans physique
+##physic=n
+iflag_phys=1
+## periode de la physique (en pas)
+iphysiq=10
+## frequence (en  jours ) de l'ecriture du fichier histphy               
+ecritphy=1
+## frequence (en jours) de l'ecriture du fichier histmth
+ecrit_mth=1
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=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
+## Nombre  d'appels des routines de rayonnements ( par jour)
+nbapp_rad=12
+##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)
+iflag_con=30
+## 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/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/offline.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/offline.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/offline.def	(revision 1242)
@@ -0,0 +1,10 @@
+## $Id: offline.def 74 2007-06-21 14:46:35Z acosce $
+T
+4
+T
+-2.
+48.1
+1
+T
+6
+2
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/orchidee.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/orchidee.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/orchidee.def	(revision 1242)
@@ -0,0 +1,44 @@
+## $Id: orchidee.def 74 2007-06-21 14:46:35Z acosce $
+#
+# SECHIBA
+#
+STOMATE_OK_CO2=FALSE
+# STOMATE_OK_STOMATE is not set
+# STOMATE_OK_DGVM is not set
+# STOMATE_WATCHOUT is not set
+#SECHIBA_restart_in=default
+SECHIBA_restart_in=start_sech.nc
+SECHIBA_rest_out=sechiba_rest.nc
+SECHIBA_reset_time=y
+SECHIBA_reset_time is not set
+OUTPUT_FILE=sechiba_out.nc
+WRITE_STEP=2592000
+SECHIBA_HISTLEVEL=10
+STOMATE_OUTPUT_FILE=stomate_history.nc
+STOMATE_HIST_DT=10.
+STOMATE_HISTLEVEL=0
+SECHIBA_DAY=0.0
+SECHIBA_ZCANOP=0.5
+DT_SLOW=86400.
+# IMPOSE_VEG is not set
+VEGETATION_FILE=carteveg5km.nc
+DIFFUCO_LEAFCI=233.
+CONDVEG_SNOWA=default
+# IMPOSE_AZE is not set
+SOILALB_FILE=soils_param.nc
+SOILTYPE_FILE=soils_param.nc 
+ENERBIL_TSURF=280.
+HYDROL_SNOW=0.0
+HYDROL_SNOWAGE=0.0
+HYDROL_SNOWICE=0.0
+HYDROL_SNOWICEAGE=0.0
+HYDROL_HDRY=1.0
+HYDROL_HUMR=1.0
+HYDROL_BQSB=default
+HYDROL_GQSB=0.0
+HYDROL_DSG=0.0
+HYDROL_DSP=default
+HYDROL_QSV=0.0
+THERMOSOIL_TPRO=280.
+RIVER_ROUTING=y
+ROUTING_FILE=routing.nc
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/physiq.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/physiq.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/physiq.def	(revision 1242)
@@ -0,0 +1,87 @@
+## $Id: physiq.def 74 2007-06-21 14:46:35Z acosce $
+type_ocean=force 
+VEGET=n
+OK_journe=n
+OK_mensuel=y
+OK_instan=n
+if_ebil=0
+#
+# parametres KE
+#
+epmax = .99
+ok_adj_ema = n
+iflag_clw = 1
+# 
+# parametres nuages
+#
+cld_lc_lsc = 0.00026
+cld_lc_con = 0.00026
+cld_tau_lsc = 3600.
+cld_tau_con = 3600.
+ffallv_lsc = 1.
+ffallv_con = 1.
+coef_eva = 0.00002
+reevap_ice = y
+iflag_cldcon = -1
+iflag_pdf = 1
+fact_cldcon = 1.
+#facttemps = 1.e-4
+facttemps = 0.0001
+ok_newmicro = y
+ratqsbas = 0.005
+ratqshaut = 0.33
+rad_froid = 35
+rad_chau1=12
+rad_chau2=11
+#ksta_ter=1.e-7
+ksta_ter=0.0000001
+#
+# parametres outputs
+#
+#niveau de sortie "hf" lev_histhf avec
+# - lev_histhf=3 => defaut
+# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele
+lev_histhf=3
+#niveau de sortie "day" lev_histday
+# - lev_histday=2 => defaut
+# - lev_histday=3 => + champs 3D => F. Lott
+# - lev_histday=4 => + champs sous-surfaces
+lev_histday=2
+#niveau de sortie "mth" lev_histmth avec
+# - lev_histmth=2 => defaut
+# - lev_histmth=3 => albedo, rugosite sous-surfaces
+# - lev_histmth=4 => champs tendances 3d
+lev_histmth=4
+#
+# parametres climatique
+#
+R_ecc = 0.016715
+R_peri = 102.7
+R_incl = 23.441
+solaire = 1365.
+co2_ppm = 348.
+CH4_ppb = 1650.
+N2O_ppb = 306.
+CFC11_ppt = 280.
+CFC12_ppt = 484.
+#
+# parametres simulateur ISCCP
+#
+top_height = 3
+#overlap = 1, 2 ou 3
+overlap = 3
+#cdmmax
+#cdmmax = 2.5E-3
+cdmmax = 0.0025
+#cdhmax
+#cdhmax = 2.0E-3
+cdhmax = 0.002
+#autres flags
+iflag_pbl = 1
+iflag_thermals = 0
+iflag_clos=0
+iflag_ratqs=0
+#
+# parametre chimie/aerosols
+#
+config_inca=chem
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/run.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/run.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/run.def	(revision 1242)
@@ -0,0 +1,25 @@
+## $Id: run.def 74 2007-06-21 14:46:35Z acosce $
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+## activation du calcul d equilibrage de charge
+adjust=_adjust_
+## Type de calendrier utilise
+## valeur possible: earth_360d (defaut), earth_365d, earth_366d
+calend=_calend_
+## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+dayref=_dayref_
+##  Annee de l'etat  initial (   avec  4  chiffres   )
+anneeref=_anneeref_
+## raz_date
+raz_date=_raz_date_
+## Nombre de jours d'integration
+nday=_nday_
+## periode de sortie des variables de controle (en pas)
+iconser= 480
+## periode d'ecriture du fichier histoire (en jour)
+iecri=10
+## periode de stockage fichier histmoy (en jour)
+periodav=1.
+## unite de sortie des impressions
+lunout=6
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/traceur.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/PARAM/traceur.def	(revision 1242)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/POST/monitoring01_inca.cfg
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/POST/monitoring01_inca.cfg	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/POST/monitoring01_inca.cfg	(revision 1242)
@@ -0,0 +1,50 @@
+#**************************************************************
+# 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
+#   - files patterns, operations,title,units must be enclosed with character '"'
+#   - 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=F0D5F4
+smooth=12
+#-----------------------------------------------------------------------------------------------------------------
+#  field | files patterns | files additionnal | operations | title | units | calcul of area
+#-----------------------------------------------------------------------------------------------------------------
+NOX_surf_global		| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=19]+NO2[d=2,k=19])" | " NOX a la surface"   | "VMR" | "aire[d=3]"
+CH4_surf_global		| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=19]"                | " CH4 a la surface"   | "VMR" | "aire[d=2]"
+CO_surf_global		| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=19]"                 | " CO a la surface"    | "VMR" | "aire[d=2]"
+O3_surf_global		| "_O3"   	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=19]"                 | " O3 a la surface"    | "VMR" | "aire[d=2]"
+HNO3_surf_global	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=19]"               | " HNO3 a la surface"  | "?" 	 | "aire[d=2]"
+NOX_545HPa_global    	| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=12]+NO2[d=2,k=12])" | " NOX a 545 HPa"      | "VMR" | "aire[d=3]"
+CH4_545HPa_global    	| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=12]"                | " CH4 a 545 HPa"      | "VMR" | "aire[d=2]"
+CO_545HPa_global     	| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=12]"                 | " CO a 545 HPa"       | "VMR" | "aire[d=2]"
+O3_545HPa_global     	| "_O3"    	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=12]"                 | " O3 a 545 HPa"       | "VMR" | "aire[d=2]"
+HNO3_545HPa_global   	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=12]"               | " HNO3 a 545 HPa"     | "?"   | "aire[d=2]"
+
+#-----------------------------------------------------------------------------------------------------------------
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/config.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/config.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4/config.card	(revision 1242)
@@ -0,0 +1,106 @@
+# $Id: config.card 72 2007-06-21 14:40:18Z acosce $
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=LI_C_00
+LongName='LMDz4INCA with New scripts'
+TagName=LMDZINCA/CH4
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Début et fin de Job
+#D-- "YYYY-MM-DD"
+DateBegin=2000-01-01
+DateEnd=2000-12-31
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=1M
+#============================
+#-- PBS Class (required on brodie mono/multi)
+#JobClass=mono
+JobClass=multi
+#============================
+#-- Total Number of Processors (minimum is 2 for a coupled configuration)
+#JobNumProcTot=2
+JobNumProcTot=4
+#============================
+#-- Run Options : here for oasis coupling
+JobRunOptions='"-np "${config_UserChoices_JobNumProcTot}'
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ4)
+CHM= (inca, INCA3) 
+
+#========================================================================
+#D-- Executable -
+[Executable]
+Name=gcm.e
+#D- For each component, Real name of executable, Name of executable for oasis
+ATM= (gcm.e, gcm.e)
+CHM= (inca.dat, inca.dat)
+
+#========================================================================
+#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=1869-12-30
+#D- Define restart simulation name (=> JOB_OS)
+RestartJobName=JobName
+#D- Path Server Group Login (=> PSGL)
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/CH4/
+# Attention login depend de la machine
+
+#========================================================================
+#D-- Post -
+[Post]
+#D- Do we rebuild parallel output, this flag determines
+#D- frequency of rebuild submission (use NONE for DRYRUN=3)
+RebuildFrequency=1M
+#D- Do we rebuild parallel output from archive
+RebuildFromArchive=true
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+TimeSeriesFrequency=1M
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average
+SeasonalFrequency=5Y
+
+#========================================================================
+#D-- ATM -
+[ATM]
+#
+WriteFrequency="1M 1D HF"
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=2000-01-01
+# Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/CH4/
+# Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
+#D-- CHM -
+[CHM]
+WriteFrequency="1M"
+Restart=n
+#D-- Last day of the experience used as restart
+RestartDate=2000-01-01
+#D-- Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/CH4/
+#D-- Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/inca.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/inca.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/inca.card	(revision 1242)
@@ -0,0 +1,56 @@
+# $Id$
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+ 
+[InitialStateFiles]
+List= ()
+	
+	
+
+[BoundaryFiles]
+List=	    (/dmnfs/cont003/p24data/ECMWF320x160/AN${year}/165_${year}${month}.nc      , u10mec.nc)\
+	    (/dmnfs/cont003/p24data/ECMWF320x160/AN${year}/166_${year}${month}.nc      , v10mec.nc)\
+	    (/dmnfs/cont003/p24data/ECMWF96x72/AN${year}/u10m_ecmwf_${year}${month}.nc , u10m.nc  )\
+	    (/dmnfs/cont003/p24data/ECMWF96x72/AN${year}/v10m_ecmwf_${year}${month}.nc , v10m.nc  )
+
+ListNonDel= (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/o3clim.nc           ,   o3clim.nc      )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/npp.nc	             ,   npp.nc	        )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/aircraft.nc         ,   aircraft.nc    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/so4.nc	             ,   so4.nc	        )\
+	    (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/sflx_lmdz_phy.present.nc,   sflx.nc	)\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/phototable.dat_V4.2	             ,   phototable.dat )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/landuse.nc          ,   landuse.nc     )\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/rhvlmdz.txt         ,	 rhv.dat        )\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/wthlmdz.txt         ,   wth.dat	)\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/clylmdz.txt         ,	 cly.dat	)\
+            (${R_INIT}/CHM/LMDZINCA/AER/rhvEC11.09092002.txt	             ,   rhvEC.txt	)\
+            (${R_INIT}/CHM/LMDZINCA/AER/wthEC11.09092002.txt	             ,   wthEC.txt	)\
+            (${R_INIT}/CHM/LMDZINCA/AER/clyEC11.09092002.txt	             ,   clyEC.txt	)
+
+[ParametersFiles]
+List=	
+
+[RestartFiles]
+# List restart that have to be saved/restored each loop (file out, saved, and in) :
+List= (precipinsoil.dat, precipinsoil.dat, precipinsoil.dat)
+
+
+
+[OutputText]
+List=  
+# avec la // : out_orchidee_*
+
+[OutputFiles]
+List=	(inca_avgr.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc,   Post_1M_inca_avgr  ),\
+	(forcage.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_forcage.nc,      NONE)
+
+
+[Post_1M_inca_avgr]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter)
+TimeSeriesVars2D = (OD550_ASBCM, OD550_ASPOMM, OD550_ASSO4M, OD550_CSSO4M, OD550_SSSSM, OD550_ASSSM, OD550_CSSSM, OD550_CIDUSTM, OD550_AIBCM, OD550_AIPOMM)
+ChunckJob2D = NONE
+TimeSeriesVars3D = (CH4,NO,NO2,CO,O3,HNO3)
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/inca.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/inca.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/inca.driver	(revision 1242)
@@ -0,0 +1,40 @@
+#- $Id$
+#!/bin/ksh
+
+#D- Driver du script pour INCA
+
+#-----------------------------------------------------------------
+function CHM_Initialize
+{
+    IGCM_debug_PushStack "CHM_Initialize"
+
+    RESOL_CHM=$( echo $RESOL_ATM | awk "-FD" '{print $2}' )
+    echo 'resol chimie', $RESOL_CHM
+    echo "INCA Tag : LMDZ-INCA" ${config_UserChoices_TagName}
+
+    IGCM_debug_PopStack "CHM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function CHM_Update
+{
+    IGCM_debug_PushStack "CHM_Update"
+    IGCM_debug_PopStack "CHM_Update"
+}
+
+#-----------------------------------------------------------------
+function CHM_Finalize
+{
+    IGCM_debug_PushStack "CHM_Finalize"
+
+    echo FINALIZE CHM !!!
+#
+#    IGCM_sys_Mkdir ${R_OUT_EXE}
+#    cd ${MODIPSL}
+#    tar cf src_${config_UserChoices_JobName}.tar  modeles
+#    IGCM_sys_Mv src_${config_UserChoices_JobName}.tar ${R_OUT_EXE}/.
+#    cd -
+#
+
+    IGCM_debug_PopStack "CHM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/lmdz.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/lmdz.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/lmdz.card	(revision 1242)
@@ -0,0 +1,45 @@
+# $Id$
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+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
+
+[InitialStateFiles]
+# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job
+List= 	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/start.nc, start.nc) \
+	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/startphy.nc, startphy.nc)\
+	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/precipinsoil.dat, precipinsoil.dat)
+
+[BoundaryFiles]
+List=   (${R_BC}/CHM/LMDZORINCA/${RESOL_ATM}_leap/EL-${RESOL_ATM}_${year}1231_limit.nc , limit.nc) 
+ListNonDel= ()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def, .), \
+	(${SUBMIT_DIR}/PARAM/gcm.def, gcm.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, Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat)
+
+
+[OutputFiles]
+List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,  Post_1M_histmth  )	
+
+[Post_1M_histmth]
+Patches= (Patch_20091030_histcom_time_axis)
+GatherWithInternal = (lon, lat, presnivs, time_counter, aire)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux_ter, taux_sic, taux_lic, taux_oce, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, 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)
+ChunckJob2D = NONE
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/lmdz.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/COMP/lmdz.driver	(revision 1242)
@@ -0,0 +1,160 @@
+#- $Id$
+#!/bin/ksh
+
+#-----------------------------------------------------------------
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' )
+    case ${RESOL_ATM} in
+	LMD7245)   BandsResol=72x45x19   ;;
+	LMD9671)   BandsResol=96x71x19   ;;
+	LMD9672)   BandsResol=96x72x19   ;;
+	LMD9695)   BandsResol=96x95x19   ;;
+	LMD14496)  BandsResol=144x96x19  ;;
+	LMD144142) BandsResol=144x142x19 ;;
+    esac
+
+
+    ##-- Calendar type for LMDZ
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=earth_366d;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d;;
+        360d)
+            CalendarTypeForLmdz=earth_360d;;
+        *)
+            CalendarTypeForLmdz=earth_360d
+    esac
+
+    ##- Default number of processor for lmdz
+    NUM_PROC_ATM=1
+    if [ X"${BATCH_NUM_PROC_TOT}" != X ] ; then
+	NUM_PROC_ATM=${BATCH_NUM_PROC_TOT}
+    fi
+
+    ##--Frequency purpose ....
+    ##--  Initialisation  ....
+    OK_instan=n
+    ##-- 
+    OK_journe=n
+    OK_mensuel=y
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) OK_mensuel=y ;;
+	5D|5d|1D|1d) OK_journe=n ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    5D|5d|1D|1d) OK_journe=n ;;
+	esac
+    done
+
+    ##--Variables used by LMDZ --
+    PAT_INST=$(     grep 'OK_instan'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_JOUR=$(     grep 'OK_journe'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_MOIS=$(     grep 'OK_mensuel'  ${SUBMIT_DIR}/PARAM/physiq.def )
+
+    PAT_iphysiq=$(  grep 'iphysiq'     ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_iperiod=$(  grep 'iperiod'     ${SUBMIT_DIR}/PARAM/gcm.def | tail -1)
+    PAT_day_step=$( grep 'day_step'    ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_ecritphy=$( grep 'ecritphy'    ${SUBMIT_DIR}/PARAM/gcm.def )
+
+
+    ## Read LMDZ_NbPeriod_adjust option in lmdz.card 
+    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_NbPeriod_adjust
+    eval LMDZ_NbPeriod_adjust=\${lmdz_UserChoices_LMDZ_NbPeriod_adjust} > /dev/null 2>&1
+
+    if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+        IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_Bands_file_name
+        eval LMDZ_Bands_file_name=\${lmdz_UserChoices_LMDZ_Bands_file_name} > /dev/null 2>&1
+    fi
+
+
+    ##-- This could be define in lmdz.card, inside section [UserChoices]
+    ##-- Otherwise we get the value in *.def
+    iperiod=$(grep iperiod ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}' | tail -1)
+    iphysiq=$(grep iphysiq ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}')
+    # day_step : number of steps per day (multiple of iperiod)
+    (( day_step = 48 * iphysiq ))
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    ecritphy=${PeriodLengthInDays}
+
+    ##-- 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
+
+    if [ ${CumulPeriod} -eq 1 ] ; then
+        if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+            IGCM_sys_Get ${LMDZ_Bands_file_name} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+            IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0
+            IGCM_debug_Print 1 "Bands file forced to ${LMDZ_Bands_file_name} and stored in ${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0"
+        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 )) 
+
+    [ ${CumulPeriod} -gt 1 ] && IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+
+    ## Mise en forme du fichier physiq.def
+    sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   	\
+	-e "s/OK_journe=.*/OK_journe=${OK_journe}/"   	\
+	-e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" 	\
+	physiq.def > physiq.def.tmp
+    IGCM_sys_Mv physiq.def.tmp physiq.def
+
+    ## Mise en forme du fichier gcm.def
+    sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" 	\
+	-e "s/${PAT_day_step}/day_step=${day_step}/" 	\
+	-e "s/${PAT_iperiod}/iperiod=${iperiod}/" 	\
+	-e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" 	\
+	gcm.def > gcm.def.tmp
+    IGCM_sys_Mv gcm.def.tmp gcm.def
+
+    greg_dat=$( IGCM_date_ConvertFormatToGregorian $PeriodDateBegin ) 
+    jul_dat=$( IGCM_date_ConvertGregorianDateToJulian $greg_dat )
+    nbjour=$( expr \( $jul_dat \% 1000 \) ) 
+    yractu=$( expr $greg_dat / 10000  )
+    echo  "dayref = " $nbjour  "year actu = " $yractu
+
+    ## Mise en forme du fichier run.def
+    sed -e "s/_dayref_/${nbjour}/" 			\
+	-e "s/_anneeref_/${yractu}/" 			\
+        -e "s/_calend_/${CalendarTypeForLmdz}/"        	\
+	-e "s/_nday_/${PeriodLengthInDays}/" 		\
+	-e "s/_raz_date_/${RAZ_DATE}/" 			\
+	-e "s/_adjust_/${LMDZ_adjust}/"                	\
+	run.def > run.def.tmp
+    IGCM_sys_Mv run.def.tmp run.def
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/Bands_96x72x19_2prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
@@ -0,0 +1,8 @@
+ 0  37
+ 1  36
+ 0  37
+ 1  36
+ 0  36
+ 1  37
+ 0  3479
+ 1  3339
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/Bands_96x72x19_4prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
@@ -0,0 +1,16 @@
+ 0  8
+ 1  29
+ 2  27
+ 3  9
+ 0  19
+ 1  18
+ 2  18
+ 3  18
+ 0  11
+ 1  27
+ 2  32
+ 3  3
+ 0  1758
+ 1  1719
+ 2  1681
+ 3  1660
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/gcm.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/gcm.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/gcm.def	(revision 1242)
@@ -0,0 +1,74 @@
+## $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
+## periode de la dissipation (en pas)
+idissip=10
+## 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=3600.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=3600.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## guide ou pas
+ok_guide=.false.
+## utilisation de restart de la version LMDZ3
+old_restart=
+## avec ou sans physique
+##physic=n
+iflag_phys=1
+## periode de la physique (en pas)
+iphysiq=10
+## frequence (en  jours ) de l'ecriture du fichier histphy               
+ecritphy=1
+## frequence (en jours) de l'ecriture du fichier histmth
+ecrit_mth=1
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=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
+## Nombre  d'appels des routines de rayonnements ( par jour)
+nbapp_rad=12
+##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)
+iflag_con=30
+## 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/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/offline.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/offline.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/offline.def	(revision 1242)
@@ -0,0 +1,10 @@
+## $Id$
+T
+4
+T
+-2.
+48.1
+1
+T
+6
+2
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/orchidee.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/orchidee.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/orchidee.def	(revision 1242)
@@ -0,0 +1,44 @@
+## $Id$
+#
+# SECHIBA
+#
+STOMATE_OK_CO2=FALSE
+# STOMATE_OK_STOMATE is not set
+# STOMATE_OK_DGVM is not set
+# STOMATE_WATCHOUT is not set
+#SECHIBA_restart_in=default
+SECHIBA_restart_in=start_sech.nc
+SECHIBA_rest_out=sechiba_rest.nc
+SECHIBA_reset_time=y
+SECHIBA_reset_time is not set
+OUTPUT_FILE=sechiba_out.nc
+WRITE_STEP=2592000
+SECHIBA_HISTLEVEL=10
+STOMATE_OUTPUT_FILE=stomate_history.nc
+STOMATE_HIST_DT=10.
+STOMATE_HISTLEVEL=0
+SECHIBA_DAY=0.0
+SECHIBA_ZCANOP=0.5
+DT_SLOW=86400.
+# IMPOSE_VEG is not set
+VEGETATION_FILE=carteveg5km.nc
+DIFFUCO_LEAFCI=233.
+CONDVEG_SNOWA=default
+# IMPOSE_AZE is not set
+SOILALB_FILE=soils_param.nc
+SOILTYPE_FILE=soils_param.nc 
+ENERBIL_TSURF=280.
+HYDROL_SNOW=0.0
+HYDROL_SNOWAGE=0.0
+HYDROL_SNOWICE=0.0
+HYDROL_SNOWICEAGE=0.0
+HYDROL_HDRY=1.0
+HYDROL_HUMR=1.0
+HYDROL_BQSB=default
+HYDROL_GQSB=0.0
+HYDROL_DSG=0.0
+HYDROL_DSP=default
+HYDROL_QSV=0.0
+THERMOSOIL_TPRO=280.
+RIVER_ROUTING=y
+ROUTING_FILE=routing.nc
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/physiq.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/physiq.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/physiq.def	(revision 1242)
@@ -0,0 +1,93 @@
+## $Id$
+type_ocean=force 
+VEGET=n
+OK_journe=n
+OK_mensuel=y
+OK_instan=n
+if_ebil=0
+#
+# parametres KE
+#
+epmax = .99
+ok_adj_ema = n
+iflag_clw = 1
+# 
+# parametres nuages
+#
+cld_lc_lsc = 0.00026
+cld_lc_con = 0.00026
+cld_tau_lsc = 3600.
+cld_tau_con = 3600.
+ffallv_lsc = 1.
+ffallv_con = 1.
+coef_eva = 0.00002
+reevap_ice = y
+iflag_cldcon = -1
+iflag_pdf = 1
+fact_cldcon = 1.
+#facttemps = 1.e-4
+facttemps = 0.0001
+ok_newmicro = y
+ratqsbas = 0.005
+ratqshaut = 0.33
+rad_froid = 35
+rad_chau1=12
+rad_chau2=11
+#ksta_ter=1.e-7
+ksta_ter=0.0000001
+#
+# parametres outputs
+#
+#niveau de sortie "hf" lev_histhf avec
+# - lev_histhf=3 => defaut
+# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele
+lev_histhf=3
+#niveau de sortie "day" lev_histday
+# - lev_histday=2 => defaut
+# - lev_histday=3 => + champs 3D => F. Lott
+# - lev_histday=4 => + champs sous-surfaces
+lev_histday=2
+#niveau de sortie "mth" lev_histmth avec
+# - lev_histmth=2 => defaut
+# - lev_histmth=3 => albedo, rugosite sous-surfaces
+# - lev_histmth=4 => champs tendances 3d
+lev_histmth=4
+#
+# parametres climatique
+#
+R_ecc = 0.016715
+R_peri = 102.7
+R_incl = 23.441
+solaire = 1365.
+co2_ppm = 348.
+CH4_ppb = 1650.
+N2O_ppb = 306.
+CFC11_ppt = 280.
+CFC12_ppt = 484.
+#
+# parametres simulateur ISCCP
+#
+top_height = 3
+#overlap = 1, 2 ou 3
+overlap = 3
+#cdmmax
+#cdmmax = 2.5E-3
+cdmmax = 0.0025
+#cdhmax
+#cdhmax = 2.0E-3
+cdhmax = 0.002
+#autres flags
+iflag_pbl = 1
+iflag_thermals = 0
+iflag_clos=0
+iflag_ratqs=0
+#
+# parametre chimie/aerosols
+#
+config_inca=aero
+aerosol_couple=n
+flag_aerosol=1
+new_aod=.true.
+ok_ade=n
+ok_aie=n
+
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/run.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/run.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/run.def	(revision 1242)
@@ -0,0 +1,25 @@
+## $Id$
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+## activation du calcul d equilibrage de charge
+adjust=_adjust_
+## Type de calendrier utilise
+## valeur possible: earth_360d (defaut), earth_365d, earth_366d
+calend=_calend_
+## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+dayref=_dayref_
+##  Annee de l'etat  initial (   avec  4  chiffres   )
+anneeref=_anneeref_
+## raz_date
+raz_date=_raz_date_
+## Nombre de jours d'integration
+nday=_nday_
+## periode de sortie des variables de controle (en pas)
+iconser= 480
+## periode d'ecriture du fichier histoire (en jour)
+iecri=10
+## periode de stockage fichier histmoy (en jour)
+periodav=1.
+## unite de sortie des impressions
+lunout=6
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/traceur.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/PARAM/traceur.def	(revision 1242)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/POST/monitoring01_inca.cfg
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/POST/monitoring01_inca.cfg	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/POST/monitoring01_inca.cfg	(revision 1242)
@@ -0,0 +1,51 @@
+#**************************************************************
+# 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
+#   - files patterns, operations,title,units must be enclosed with character '"'
+#   - 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=F0D5F4
+smooth=12
+#-----------------------------------------------------------------------------------------------------------------
+#  field | files patterns | files additionnal | operations | title | units | calcul of area
+#-----------------------------------------------------------------------------------------------------------------
+ OD550_AERO_global | "OD550_ASBCM OD550_ASPOMM OD550_ASSO4M OD550_CSSO4M OD550_SSSSM OD550_ASSSM OD550_CSSSM OD550_CIDUSTM OD550_AIBCM OD550_AIPOMM" | LMDZ4.0_9695_grid.nc | "(OD550_ASBCM[d=1]+ OD550_ASPOMM[d=2]+ OD550_ASSO4M[d=3]+ OD550_CSSO4M[d=4]+ OD550_SSSSM[d=5]+ OD550_ASSSM[d=6]+ OD550_CSSSM[d=7]+ OD550_CIDUSTM[d=8]+ OD550_AIBCM[d=9]+ OD550_AIPOMM[d=10])" | "total  aerosol optical depth @550nm"  |"(-)" | "aire[d=11]"
+NOX_surf_global		| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=19]+NO2[d=2,k=19])" | " NOX a la surface"   | "VMR" | "aire[d=3]"
+CH4_surf_global		| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=19]"                | " CH4 a la surface"   | "VMR" | "aire[d=2]"
+CO_surf_global		| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=19]"                 | " CO a la surface"    | "VMR" | "aire[d=2]"
+O3_surf_global		| "_O3"   	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=19]"                 | " O3 a la surface"    | "VMR" | "aire[d=2]"
+HNO3_surf_global	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=19]"               | " HNO3 a la surface"  | "?" 	 | "aire[d=2]"
+NOX_545HPa_global    	| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=12]+NO2[d=2,k=12])" | " NOX a 545 HPa"      | "VMR" | "aire[d=3]"
+CH4_545HPa_global    	| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=12]"                | " CH4 a 545 HPa"      | "VMR" | "aire[d=2]"
+CO_545HPa_global     	| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=12]"                 | " CO a 545 HPa"       | "VMR" | "aire[d=2]"
+O3_545HPa_global     	| "_O3"    	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=12]"                 | " O3 a 545 HPa"       | "VMR" | "aire[d=2]"
+HNO3_545HPa_global   	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=12]"               | " HNO3 a 545 HPa"     | "?"   | "aire[d=2]"
+
+#-----------------------------------------------------------------------------------------------------------------
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/config.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/config.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_CH4_AER/config.card	(revision 1242)
@@ -0,0 +1,107 @@
+# $Id$
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=LI_CA_00
+LongName='LMDz4INCA with New scripts'
+TagName=LMDZINCA/CH4_AER
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Début et fin de Job
+#D-- "YYYY-MM-DD"
+DateBegin=2000-01-01
+DateEnd=2000-12-31
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=1M
+#============================
+#-- PBS Class (required on brodie mono/multi)
+#JobClass=mono
+JobClass=multi
+#============================
+#-- Total Number of Processors (minimum is 2 for a coupled configuration)
+#JobNumProcTot=2
+JobNumProcTot=4
+#============================
+#-- Run Options : here for oasis coupling
+JobRunOptions='"-np "${config_UserChoices_JobNumProcTot}'
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ4)
+CHM= (inca, INCA3) 
+
+#========================================================================
+#D-- Executable -
+[Executable]
+Name=gcm.e
+#D- For each component, Real name of executable, Name of executable for oasis
+#ATM= (gcm_${config_UserChoices_TagName}.e, gcm.e)
+ATM= (gcm.e, gcm.e)
+CHM= (inca.dat, inca.dat)
+
+#========================================================================
+#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=1869-12-30
+#D- Define restart simulation name (=> JOB_OS)
+RestartJobName=JobName
+#D- Path Server Group Login (=> PSGL)
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/CH4_AER/
+# Attention login depend de la machine
+
+#========================================================================
+#D-- Post -
+[Post]
+#D- Do we rebuild parallel output, this flag determines
+#D- frequency of rebuild submission (use NONE for DRYRUN=3)
+RebuildFrequency=1M
+#D- Do we rebuild parallel output from archive
+RebuildFromArchive=true
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+TimeSeriesFrequency=1M
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average
+SeasonalFrequency=5Y
+
+#========================================================================
+#D-- ATM -
+[ATM]
+#
+WriteFrequency="1M 1D HF"
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=2000-01-01
+# Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/CH4_AER/
+# Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
+#D-- CHM -
+[CHM]
+WriteFrequency="1M"
+Restart=n
+#D-- Last day of the experience used as restart
+RestartDate=2000-01-01
+#D-- Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/CH4_AER/
+#D-- Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/inca.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/inca.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/inca.card	(revision 1242)
@@ -0,0 +1,48 @@
+#- $Id: inca.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+[BoundaryFiles]
+List=
+
+ListNonDel= (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/oxydants.nc             ,   oxydants.nc    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/npp.nc	                 ,   npp.nc	    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/aircraft.nc             ,   aircraft.nc    )\
+	    (${R_INIT}/CHM/LMDZINCA/NMHC/INCA${RESOL_CHM}/sflx_Present.nc        ,   sflx.nc        )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/phototable.dat_V4.2                      ,   phototable.dat )\
+	    (${R_INIT}/CHM/LMDZINCA/NMHC/INCA${RESOL_CHM}/landuse.nc             ,   landuse.nc     )\
+	    (${R_INIT}/CHM/LMDZINCA/GES/INCA${RESOL_CHM}/fco2.andres.fos.glo.fy95.m.lmdz96.nc            ,andres.nc        )\
+	    (${R_INIT}/CHM/LMDZINCA/GES/INCA${RESOL_CHM}/fco2.casa.nep.glo.fy90.m.lmdz96.nc              ,casa_m.nc        )\
+	    (${R_INIT}/CHM/LMDZINCA/GES/INCA${RESOL_CHM}/fco2.taka.oce.glo.fy95.m.lmdz96.nc              ,taka.nc          )\
+	    (${R_INIT}/CHM/LMDZINCA/GES/INCA${RESOL_CHM}/fco2.casa.nep.glo.fy90.3hr.lmdz96.${month}.nc   ,casa_h.nc        )
+
+
+[ParametersFiles]
+List=	
+
+[RestartFiles]
+# List restart that have to be saved/restored each loop (file out, saved, and in) :
+List= (NONE) 
+
+
+
+[OutputText]
+List=  
+# avec la // : out_orchidee_*
+
+[OutputFiles]
+List=	(inca_avgr.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc,  Post_1M_inca_avgr  )
+
+
+[Post_1M_inca_avgr]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = (CH4,NO,NO2,CO,O3,HNO3)
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/inca.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/inca.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/inca.driver	(revision 1242)
@@ -0,0 +1,40 @@
+#- $Id: inca.driver 119 2007-07-31 13:56:33Z acosce $
+#!/bin/ksh
+
+#D- Driver du script pour INCA
+
+#-----------------------------------------------------------------
+function CHM_Initialize
+{
+    IGCM_debug_PushStack "CHM_Initialize"
+
+    RESOL_CHM=$( echo $RESOL_ATM | awk "-FD" '{print $2}' )
+    echo 'resol chimie', $RESOL_CHM
+    echo "INCA Tag : LMDZ-INCA" ${config_UserChoices_TagName}
+
+    IGCM_debug_PopStack "CHM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function CHM_Update
+{
+    IGCM_debug_PushStack "CHM_Update"
+    IGCM_debug_PopStack "CHM_Update"
+}
+
+#-----------------------------------------------------------------
+function CHM_Finalize
+{
+    IGCM_debug_PushStack "CHM_Finalize"
+
+    echo FINALIZE CHM !!!
+#
+#    IGCM_sys_Mkdir ${R_OUT_EXE}
+#    cd ${MODIPSL}
+#    tar cf src_${config_UserChoices_JobName}.tar  modeles
+#    IGCM_sys_Mv src_${config_UserChoices_JobName}.tar ${R_OUT_EXE}/.
+#    cd -
+#
+
+    IGCM_debug_PopStack "CHM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/lmdz.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/lmdz.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/lmdz.card	(revision 1242)
@@ -0,0 +1,44 @@
+#- $Id: lmdz.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+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
+
+[InitialStateFiles]
+# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job
+List= 	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/start.nc, start.nc) \
+	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/startphy.nc, startphy.nc)
+
+
+[BoundaryFiles]
+List=   (${R_BC}/CHM/LMDZORINCA/${RESOL_ATM}_leap/EL-${RESOL_ATM}_${year}1231_limit.nc , limit.nc) 
+ListNonDel= ()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def, .), \
+	(${SUBMIT_DIR}/PARAM/gcm.def, gcm.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, Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat)
+
+
+[OutputFiles]
+List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,    Post_1M_histmth)	
+
+[Post_1M_histmth]
+Patches= (Patch_20091030_histcom_time_axis)
+GatherWithInternal = (lon, lat, presnivs, time_counter, aire)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux_ter, taux_sic, taux_lic, taux_oce, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, 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)
+ChunckJob2D = NONE
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/lmdz.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/COMP/lmdz.driver	(revision 1242)
@@ -0,0 +1,160 @@
+#- $Id: lmdz.driver 841 2009-12-08 09:48:13Z acosce $
+#!/bin/ksh
+
+#-----------------------------------------------------------------
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' )
+    case ${RESOL_ATM} in
+	LMD7245)   BandsResol=72x45x19   ;;
+	LMD9671)   BandsResol=96x71x19   ;;
+	LMD9672)   BandsResol=96x72x19   ;;
+	LMD9695)   BandsResol=96x95x19   ;;
+	LMD14496)  BandsResol=144x96x19  ;;
+	LMD144142) BandsResol=144x142x19 ;;
+    esac
+
+
+    ##-- Calendar type for LMDZ
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=earth_366d;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d;;
+        360d)
+            CalendarTypeForLmdz=earth_360d;;
+        *)
+            CalendarTypeForLmdz=earth_360d
+    esac
+
+    ##- Default number of processor for lmdz
+    NUM_PROC_ATM=1
+    if [ X"${BATCH_NUM_PROC_TOT}" != X ] ; then
+	NUM_PROC_ATM=${BATCH_NUM_PROC_TOT}
+    fi
+
+    ##--Frequency purpose ....
+    ##--  Initialisation  ....
+    OK_instan=n
+    ##-- 
+    OK_journe=n
+    OK_mensuel=y
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) OK_mensuel=y ;;
+	5D|5d|1D|1d) OK_journe=n ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    5D|5d|1D|1d) OK_journe=n ;;
+	esac
+    done
+
+    ##--Variables used by LMDZ --
+    PAT_INST=$(     grep 'OK_instan'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_JOUR=$(     grep 'OK_journe'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_MOIS=$(     grep 'OK_mensuel'  ${SUBMIT_DIR}/PARAM/physiq.def )
+
+    PAT_iphysiq=$(  grep 'iphysiq'     ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_iperiod=$(  grep 'iperiod'     ${SUBMIT_DIR}/PARAM/gcm.def | tail -1)
+    PAT_day_step=$( grep 'day_step'    ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_ecritphy=$( grep 'ecritphy'    ${SUBMIT_DIR}/PARAM/gcm.def )
+
+
+    ## Read LMDZ_NbPeriod_adjust option in lmdz.card 
+    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_NbPeriod_adjust
+    eval LMDZ_NbPeriod_adjust=\${lmdz_UserChoices_LMDZ_NbPeriod_adjust} > /dev/null 2>&1
+
+    if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+        IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_Bands_file_name
+        eval LMDZ_Bands_file_name=\${lmdz_UserChoices_LMDZ_Bands_file_name} > /dev/null 2>&1
+    fi
+
+
+    ##-- This could be define in lmdz.card, inside section [UserChoices]
+    ##-- Otherwise we get the value in *.def
+    iperiod=$(grep iperiod ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}' | tail -1)
+    iphysiq=$(grep iphysiq ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}')
+    # day_step : number of steps per day (multiple of iperiod)
+    (( day_step = 48 * iphysiq ))
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    ecritphy=${PeriodLengthInDays}
+
+    ##-- 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
+
+    if [ ${CumulPeriod} -eq 1 ] ; then
+        if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+            IGCM_sys_Get ${LMDZ_Bands_file_name} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+            IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0
+            IGCM_debug_Print 1 "Bands file forced to ${LMDZ_Bands_file_name} and stored in ${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0"
+        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 )) 
+
+    [ ${CumulPeriod} -gt 1 ] && IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+
+    ## Mise en forme du fichier physiq.def
+    sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   	\
+	-e "s/OK_journe=.*/OK_journe=${OK_journe}/"   	\
+	-e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" 	\
+	physiq.def > physiq.def.tmp
+    IGCM_sys_Mv physiq.def.tmp physiq.def
+
+    ## Mise en forme du fichier gcm.def
+    sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" 	\
+	-e "s/${PAT_day_step}/day_step=${day_step}/" 	\
+	-e "s/${PAT_iperiod}/iperiod=${iperiod}/" 	\
+	-e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" 	\
+	gcm.def > gcm.def.tmp
+    IGCM_sys_Mv gcm.def.tmp gcm.def
+
+    greg_dat=$( IGCM_date_ConvertFormatToGregorian $PeriodDateBegin ) 
+    jul_dat=$( IGCM_date_ConvertGregorianDateToJulian $greg_dat )
+    nbjour=$( expr \( $jul_dat \% 1000 \) ) 
+    yractu=$( expr $greg_dat / 10000  )
+    echo  "dayref = " $nbjour  "year actu = " $yractu
+
+    ## Mise en forme du fichier run.def
+    sed -e "s/_dayref_/${nbjour}/" 			\
+	-e "s/_anneeref_/${yractu}/" 			\
+        -e "s/_calend_/${CalendarTypeForLmdz}/"        	\
+	-e "s/_nday_/${PeriodLengthInDays}/" 		\
+	-e "s/_raz_date_/${RAZ_DATE}/" 			\
+	-e "s/_adjust_/${LMDZ_adjust}/"                	\
+	run.def > run.def.tmp
+    IGCM_sys_Mv run.def.tmp run.def
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/Bands_96x72x19_2prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
@@ -0,0 +1,8 @@
+ 0  37
+ 1  36
+ 0  37
+ 1  36
+ 0  40
+ 1  33
+ 0  3660
+ 1  3158
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/Bands_96x72x19_4prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
@@ -0,0 +1,16 @@
+ 0  8
+ 1  29
+ 2  27
+ 3  9
+ 0  19
+ 1  18
+ 2  18
+ 3  18
+ 0  14
+ 1  34
+ 2  22
+ 3  3
+ 0  1898
+ 1  1727
+ 2  1664
+ 3  1529
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/gcm.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/gcm.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/gcm.def	(revision 1242)
@@ -0,0 +1,74 @@
+## $Id: gcm.def 71 2007-06-21 14:26:03Z acosce $
+## 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
+## periode de la dissipation (en pas)
+idissip=10
+## 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=3600.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=3600.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## guide ou pas
+ok_guide=.false.
+## utilisation de restart de la version LMDZ3
+old_restart=
+## avec ou sans physique
+##physic=n
+iflag_phys=1
+## periode de la physique (en pas)
+iphysiq=10
+## frequence (en  jours ) de l'ecriture du fichier histphy               
+ecritphy=1
+## frequence (en jours ) de l'ecriture du fichier histmth
+ecrit_mth=1
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=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
+## Nombre  d'appels des routines de rayonnements ( par jour)
+nbapp_rad=12
+##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)
+iflag_con=30
+## 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/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/offline.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/offline.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/offline.def	(revision 1242)
@@ -0,0 +1,10 @@
+## $Id: offline.def 71 2007-06-21 14:26:03Z acosce $
+T
+4
+T
+-2.
+48.1
+1
+T
+6
+2
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/orchidee.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/orchidee.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/orchidee.def	(revision 1242)
@@ -0,0 +1,44 @@
+## $Id: orchidee.def 71 2007-06-21 14:26:03Z acosce $
+#
+# SECHIBA
+#
+STOMATE_OK_CO2=FALSE
+# STOMATE_OK_STOMATE is not set
+# STOMATE_OK_DGVM is not set
+# STOMATE_WATCHOUT is not set
+#SECHIBA_restart_in=default
+SECHIBA_restart_in=start_sech.nc
+SECHIBA_rest_out=sechiba_rest.nc
+SECHIBA_reset_time=y
+SECHIBA_reset_time is not set
+OUTPUT_FILE=sechiba_out.nc
+WRITE_STEP=2592000
+SECHIBA_HISTLEVEL=10
+STOMATE_OUTPUT_FILE=stomate_history.nc
+STOMATE_HIST_DT=10.
+STOMATE_HISTLEVEL=0
+SECHIBA_DAY=0.0
+SECHIBA_ZCANOP=0.5
+DT_SLOW=86400.
+# IMPOSE_VEG is not set
+VEGETATION_FILE=carteveg5km.nc
+DIFFUCO_LEAFCI=233.
+CONDVEG_SNOWA=default
+# IMPOSE_AZE is not set
+SOILALB_FILE=soils_param.nc
+SOILTYPE_FILE=soils_param.nc 
+ENERBIL_TSURF=280.
+HYDROL_SNOW=0.0
+HYDROL_SNOWAGE=0.0
+HYDROL_SNOWICE=0.0
+HYDROL_SNOWICEAGE=0.0
+HYDROL_HDRY=1.0
+HYDROL_HUMR=1.0
+HYDROL_BQSB=default
+HYDROL_GQSB=0.0
+HYDROL_DSG=0.0
+HYDROL_DSP=default
+HYDROL_QSV=0.0
+THERMOSOIL_TPRO=280.
+RIVER_ROUTING=y
+ROUTING_FILE=routing.nc
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/physiq.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/physiq.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/physiq.def	(revision 1242)
@@ -0,0 +1,87 @@
+##- $Id: physiq.def 71 2007-06-21 14:26:03Z acosce $
+type_ocean=force 
+VEGET=n
+OK_journe=n
+OK_mensuel=y
+OK_instan=n
+if_ebil=0
+#
+# parametres KE
+#
+epmax = .99
+ok_adj_ema = n
+iflag_clw = 1
+# 
+# parametres nuages
+#
+cld_lc_lsc = 0.00026
+cld_lc_con = 0.00026
+cld_tau_lsc = 3600.
+cld_tau_con = 3600.
+ffallv_lsc = 1.
+ffallv_con = 1.
+coef_eva = 0.00002
+reevap_ice = y
+iflag_cldcon = 3
+iflag_pdf = 1
+fact_cldcon = 1.
+#facttemps = 1.e-4
+facttemps = 0.0001
+ok_newmicro = y
+ratqsbas = 0.005
+ratqshaut = 0.33
+rad_froid = 35
+rad_chau1=12
+rad_chau2=11
+#ksta_ter=1.e-7
+ksta_ter=0.0000001
+#
+# parametres outputs
+#
+#niveau de sortie "hf" lev_histhf avec
+# - lev_histhf=3 => defaut
+# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele
+lev_histhf=3
+#niveau de sortie "day" lev_histday
+# - lev_histday=2 => defaut
+# - lev_histday=3 => + champs 3D => F. Lott
+# - lev_histday=4 => + champs sous-surfaces
+lev_histday=2
+#niveau de sortie "mth" lev_histmth avec
+# - lev_histmth=2 => defaut
+# - lev_histmth=3 => albedo, rugosite sous-surfaces
+# - lev_histmth=4 => champs tendances 3d
+lev_histmth=4
+#
+# parametres climatique
+#
+R_ecc = 0.016715
+R_peri = 102.7
+R_incl = 23.441
+solaire = 1365.
+co2_ppm = 348.
+CH4_ppb = 1650.
+N2O_ppb = 306.
+CFC11_ppt = 280.
+CFC12_ppt = 484.
+#
+# parametres simulateur ISCCP
+#
+top_height = 3
+#overlap = 1, 2 ou 3
+overlap = 3
+#cdmmax
+#cdmmax = 2.5E-3
+cdmmax = 0.0025
+#cdhmax
+#cdhmax = 2.0E-3
+cdhmax = 0.002
+#autres flags
+iflag_pbl = 1
+iflag_thermals = 0
+iflag_clos=0
+iflag_ratqs=0
+#
+# parametre chimie/aerosols
+#
+config_inca=chem
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/run.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/run.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/run.def	(revision 1242)
@@ -0,0 +1,25 @@
+## $Id: run.def 71 2007-06-21 14:26:03Z acosce $
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+## activation du calcul d equilibrage de charge
+adjust=_adjust_
+## Type de calendrier utilise
+## valeur possible: earth_360d (defaut), earth_365d, earth_366d
+calend=_calend_
+## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+dayref=_dayref_
+##  Annee de l'etat  initial (   avec  4  chiffres   )
+anneeref=_anneeref_
+## raz_date
+raz_date=_raz_date_
+## Nombre de jours d'integration
+nday=_nday_
+## periode de sortie des variables de controle (en pas)
+iconser= 480
+## periode d'ecriture du fichier histoire (en jour)
+iecri=10
+## periode de stockage fichier histmoy (en jour)
+periodav=1.
+## unite de sortie des impressions
+lunout=6
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/traceur.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/PARAM/traceur.def	(revision 1242)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/config.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/config.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_GES/config.card	(revision 1242)
@@ -0,0 +1,106 @@
+#- $Id: config.card 69 2007-06-21 14:08:22Z acosce $
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=LI_G_00
+LongName='LMDz4INCA with New scripts'
+TagName=LMDZINCA/GES
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Début et fin de Job
+#D-- "YYYY-MM-DD"
+DateBegin=2000-01-01
+DateEnd=2000-12-31
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=1M
+#============================
+#-- PBS Class (required on brodie mono/multi)
+#JobClass=mono
+JobClass=multi
+#============================
+#-- Total Number of Processors (minimum is 2 for a coupled configuration)
+#JobNumProcTot=2
+JobNumProcTot=4
+#============================
+#-- Run Options : here for oasis coupling
+JobRunOptions='"-np "${config_UserChoices_JobNumProcTot}'
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ4)
+CHM= (inca, INCA3) 
+
+#========================================================================
+#D-- Executable -
+[Executable]
+Name=gcm.e
+#D- For each component, Real name of executable, Name of executable for oasis
+ATM= (gcm.e, gcm.e)
+CHM= (inca.dat, inca.dat)
+
+#========================================================================
+#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=1869-12-30
+#D- Define restart simulation name (=> JOB_OS)
+RestartJobName=JobName
+#D- Path Server Group Login (=> PSGL)
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/GES/
+# Attention login depend de la machine
+
+#========================================================================
+#D-- Post -
+[Post]
+#D- Do we rebuild parallel output, this flag determines
+#D- frequency of rebuild submission (use NONE for DRYRUN=3)
+RebuildFrequency=1M
+#D- Do we rebuild parallel output from archive
+RebuildFromArchive=true
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+TimeSeriesFrequency=1M
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average
+SeasonalFrequency=5Y
+
+#========================================================================
+#D-- ATM -
+[ATM]
+#
+WriteFrequency="1M 1D HF"
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=2000-01-01
+# Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/GES/
+# Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
+#D-- CHM -
+[CHM]
+WriteFrequency="1M"
+Restart=n
+#D-- Last day of the experience used as restart
+RestartDate=2000-01-01
+#D-- Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/GES/
+#D-- Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/inca.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/inca.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/inca.card	(revision 1242)
@@ -0,0 +1,44 @@
+#- $Id: inca.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+[BoundaryFiles]
+List=
+
+ListNonDel= (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/o3clim.nc               ,   o3clim.nc      )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/npp.nc	                 ,   npp.nc	    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/aircraft.nc             ,   aircraft.nc    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/so4.nc	                 ,   so4.nc	    )\
+	    (${R_INIT}/CHM/LMDZINCA/NMHC/INCA${RESOL_CHM}/sflx_Present.nc        ,   sflx.nc        )\
+	    (${R_INIT}/CHM/LMDZINCA/NMHC/phototable.dat_V4.1	                 ,   phototable.dat )\
+	    (${R_INIT}/CHM/LMDZINCA/NMHC/INCA${RESOL_CHM}/landuse.nc             ,   landuse.nc     )
+
+[ParametersFiles]
+List=	
+
+[RestartFiles]
+# List restart that have to be saved/restored each loop (file out, saved, and in) :
+List= (NONE) 
+
+
+
+[OutputText]
+List=  
+# avec la // : out_orchidee_*
+
+[OutputFiles]
+List=	(inca_avgr.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc,  Post_1M_inca_avgr   )
+
+
+[Post_1M_inca_avgr]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter)
+TimeSeriesVars2D = ()
+ChunckJob2D = NONE
+TimeSeriesVars3D = (CH4,NO,NO2,CO,O3,HNO3)
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/inca.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/inca.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/inca.driver	(revision 1242)
@@ -0,0 +1,40 @@
+#- $Id: inca.driver 119 2007-07-31 13:56:33Z acosce $
+#!/bin/ksh
+
+#D- Driver du script pour INCA
+
+#-----------------------------------------------------------------
+function CHM_Initialize
+{
+    IGCM_debug_PushStack "CHM_Initialize"
+
+    RESOL_CHM=$( echo $RESOL_ATM | awk "-FD" '{print $2}' )
+    echo 'resol chimie', $RESOL_CHM
+    echo "INCA Tag : LMDZ-INCA" ${config_UserChoices_TagName}
+
+    IGCM_debug_PopStack "CHM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function CHM_Update
+{
+    IGCM_debug_PushStack "CHM_Update"
+    IGCM_debug_PopStack "CHM_Update"
+}
+
+#-----------------------------------------------------------------
+function CHM_Finalize
+{
+    IGCM_debug_PushStack "CHM_Finalize"
+
+    echo FINALIZE CHM !!!
+#
+#    IGCM_sys_Mkdir ${R_OUT_EXE}
+#    cd ${MODIPSL}
+#    tar cf src_${config_UserChoices_JobName}.tar  modeles
+#    IGCM_sys_Mv src_${config_UserChoices_JobName}.tar ${R_OUT_EXE}/.
+#    cd -
+#
+
+    IGCM_debug_PopStack "CHM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/lmdz.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/lmdz.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/lmdz.card	(revision 1242)
@@ -0,0 +1,44 @@
+#- $Id: lmdz.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+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
+
+[InitialStateFiles]
+# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job
+List= 	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/start.nc, start.nc) \
+	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/startphy.nc, startphy.nc)
+
+
+[BoundaryFiles]
+List=   (${R_BC}/CHM/LMDZORINCA/${RESOL_ATM}_leap/EL-${RESOL_ATM}_${year}1231_limit.nc , limit.nc) 
+ListNonDel= ()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def, .), \
+	(${SUBMIT_DIR}/PARAM/gcm.def, gcm.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, Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat)
+
+
+[OutputFiles]
+List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,    Post_1M_histmth  )	
+
+[Post_1M_histmth]
+Patches= (Patch_20091030_histcom_time_axis)
+GatherWithInternal = (lon, lat, presnivs, time_counter, aire)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux_ter, taux_sic, taux_lic, taux_oce, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, 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)
+ChunckJob2D = NONE
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/lmdz.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/COMP/lmdz.driver	(revision 1242)
@@ -0,0 +1,160 @@
+#- $Id: lmdz.driver 841 2009-12-08 09:48:13Z acosce $
+#!/bin/ksh
+
+#-----------------------------------------------------------------
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' )
+    case ${RESOL_ATM} in
+	LMD7245)   BandsResol=72x45x19   ;;
+	LMD9671)   BandsResol=96x71x19   ;;
+	LMD9672)   BandsResol=96x72x19   ;;
+	LMD9695)   BandsResol=96x95x19   ;;
+	LMD14496)  BandsResol=144x96x19  ;;
+	LMD144142) BandsResol=144x142x19 ;;
+    esac
+
+
+    ##-- Calendar type for LMDZ
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=earth_366d;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d;;
+        360d)
+            CalendarTypeForLmdz=earth_360d;;
+        *)
+            CalendarTypeForLmdz=earth_360d
+    esac
+
+    ##- Default number of processor for lmdz
+    NUM_PROC_ATM=1
+    if [ X"${BATCH_NUM_PROC_TOT}" != X ] ; then
+	NUM_PROC_ATM=${BATCH_NUM_PROC_TOT}
+    fi
+
+    ##--Frequency purpose ....
+    ##--  Initialisation  ....
+    OK_instan=n
+    ##-- 
+    OK_journe=n
+    OK_mensuel=y
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) OK_mensuel=y ;;
+	5D|5d|1D|1d) OK_journe=n ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    5D|5d|1D|1d) OK_journe=n ;;
+	esac
+    done
+
+    ##--Variables used by LMDZ --
+    PAT_INST=$(     grep 'OK_instan'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_JOUR=$(     grep 'OK_journe'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_MOIS=$(     grep 'OK_mensuel'  ${SUBMIT_DIR}/PARAM/physiq.def )
+
+    PAT_iphysiq=$(  grep 'iphysiq'     ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_iperiod=$(  grep 'iperiod'     ${SUBMIT_DIR}/PARAM/gcm.def | tail -1)
+    PAT_day_step=$( grep 'day_step'    ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_ecritphy=$( grep 'ecritphy'    ${SUBMIT_DIR}/PARAM/gcm.def )
+
+
+    ## Read LMDZ_NbPeriod_adjust option in lmdz.card 
+    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_NbPeriod_adjust
+    eval LMDZ_NbPeriod_adjust=\${lmdz_UserChoices_LMDZ_NbPeriod_adjust} > /dev/null 2>&1
+
+    if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+        IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_Bands_file_name
+        eval LMDZ_Bands_file_name=\${lmdz_UserChoices_LMDZ_Bands_file_name} > /dev/null 2>&1
+    fi
+
+
+    ##-- This could be define in lmdz.card, inside section [UserChoices]
+    ##-- Otherwise we get the value in *.def
+    iperiod=$(grep iperiod ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}' | tail -1)
+    iphysiq=$(grep iphysiq ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}')
+    # day_step : number of steps per day (multiple of iperiod)
+    (( day_step = 48 * iphysiq ))
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    ecritphy=${PeriodLengthInDays}
+
+    ##-- 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
+
+    if [ ${CumulPeriod} -eq 1 ] ; then
+        if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+            IGCM_sys_Get ${LMDZ_Bands_file_name} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+            IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0
+            IGCM_debug_Print 1 "Bands file forced to ${LMDZ_Bands_file_name} and stored in ${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0"
+        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 )) 
+
+    [ ${CumulPeriod} -gt 1 ] && IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+
+    ## Mise en forme du fichier physiq.def
+    sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   	\
+	-e "s/OK_journe=.*/OK_journe=${OK_journe}/"   	\
+	-e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" 	\
+	physiq.def > physiq.def.tmp
+    IGCM_sys_Mv physiq.def.tmp physiq.def
+
+    ## Mise en forme du fichier gcm.def
+    sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" 	\
+	-e "s/${PAT_day_step}/day_step=${day_step}/" 	\
+	-e "s/${PAT_iperiod}/iperiod=${iperiod}/" 	\
+	-e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" 	\
+	gcm.def > gcm.def.tmp
+    IGCM_sys_Mv gcm.def.tmp gcm.def
+
+    greg_dat=$( IGCM_date_ConvertFormatToGregorian $PeriodDateBegin ) 
+    jul_dat=$( IGCM_date_ConvertGregorianDateToJulian $greg_dat )
+    nbjour=$( expr \( $jul_dat \% 1000 \) ) 
+    yractu=$( expr $greg_dat / 10000  )
+    echo  "dayref = " $nbjour  "year actu = " $yractu
+
+    ## Mise en forme du fichier run.def
+    sed -e "s/_dayref_/${nbjour}/" 			\
+	-e "s/_anneeref_/${yractu}/" 			\
+        -e "s/_calend_/${CalendarTypeForLmdz}/"        	\
+	-e "s/_nday_/${PeriodLengthInDays}/" 		\
+	-e "s/_raz_date_/${RAZ_DATE}/" 			\
+	-e "s/_adjust_/${LMDZ_adjust}/"                	\
+	run.def > run.def.tmp
+    IGCM_sys_Mv run.def.tmp run.def
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/Bands_96x72x19_2prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
@@ -0,0 +1,8 @@
+ 0  37
+ 1  36
+ 0  37
+ 1  36
+ 0  36
+ 1  37
+ 0  3460
+ 1  3358
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/Bands_96x72x19_4prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
@@ -0,0 +1,16 @@
+ 0  8
+ 1  29
+ 2  27
+ 3  9
+ 0  19
+ 1  18
+ 2  18
+ 3  18
+ 0  3
+ 1  29
+ 2  30
+ 3  11
+ 0  1772
+ 1  1670
+ 2  1644
+ 3  1732
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/gcm.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/gcm.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/gcm.def	(revision 1242)
@@ -0,0 +1,74 @@
+## $Id: gcm.def 71 2007-06-21 14:26:03Z acosce $
+## 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
+## periode de la dissipation (en pas)
+idissip=10
+## 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=3600.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=3600.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## guide ou pas
+ok_guide=.false.
+## utilisation de restart de la version LMDZ3
+old_restart=
+## avec ou sans physique
+##physic=n
+iflag_phys=1
+## periode de la physique (en pas)
+iphysiq=10
+## frequence (en  jours ) de l'ecriture du fichier histphy               
+ecritphy=1
+## frequence (en jours ) de l'ecriture du fichier histmth
+ecrit_mth=1
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=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
+## Nombre  d'appels des routines de rayonnements ( par jour)
+nbapp_rad=12
+##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)
+iflag_con=30
+## 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/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/offline.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/offline.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/offline.def	(revision 1242)
@@ -0,0 +1,10 @@
+## $Id: offline.def 71 2007-06-21 14:26:03Z acosce $
+T
+4
+T
+-2.
+48.1
+1
+T
+6
+2
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/orchidee.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/orchidee.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/orchidee.def	(revision 1242)
@@ -0,0 +1,44 @@
+## $Id: orchidee.def 71 2007-06-21 14:26:03Z acosce $
+#
+# SECHIBA
+#
+STOMATE_OK_CO2=FALSE
+# STOMATE_OK_STOMATE is not set
+# STOMATE_OK_DGVM is not set
+# STOMATE_WATCHOUT is not set
+#SECHIBA_restart_in=default
+SECHIBA_restart_in=start_sech.nc
+SECHIBA_rest_out=sechiba_rest.nc
+SECHIBA_reset_time=y
+SECHIBA_reset_time is not set
+OUTPUT_FILE=sechiba_out.nc
+WRITE_STEP=2592000
+SECHIBA_HISTLEVEL=10
+STOMATE_OUTPUT_FILE=stomate_history.nc
+STOMATE_HIST_DT=10.
+STOMATE_HISTLEVEL=0
+SECHIBA_DAY=0.0
+SECHIBA_ZCANOP=0.5
+DT_SLOW=86400.
+# IMPOSE_VEG is not set
+VEGETATION_FILE=carteveg5km.nc
+DIFFUCO_LEAFCI=233.
+CONDVEG_SNOWA=default
+# IMPOSE_AZE is not set
+SOILALB_FILE=soils_param.nc
+SOILTYPE_FILE=soils_param.nc 
+ENERBIL_TSURF=280.
+HYDROL_SNOW=0.0
+HYDROL_SNOWAGE=0.0
+HYDROL_SNOWICE=0.0
+HYDROL_SNOWICEAGE=0.0
+HYDROL_HDRY=1.0
+HYDROL_HUMR=1.0
+HYDROL_BQSB=default
+HYDROL_GQSB=0.0
+HYDROL_DSG=0.0
+HYDROL_DSP=default
+HYDROL_QSV=0.0
+THERMOSOIL_TPRO=280.
+RIVER_ROUTING=y
+ROUTING_FILE=routing.nc
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/physiq.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/physiq.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/physiq.def	(revision 1242)
@@ -0,0 +1,87 @@
+##- $Id: physiq.def 71 2007-06-21 14:26:03Z acosce $
+type_ocean=force 
+VEGET=n
+OK_journe=n
+OK_mensuel=y
+OK_instan=n
+if_ebil=0
+#
+# parametres KE
+#
+epmax = .99
+ok_adj_ema = n
+iflag_clw = 1
+# 
+# parametres nuages
+#
+cld_lc_lsc = 0.00026
+cld_lc_con = 0.00026
+cld_tau_lsc = 3600.
+cld_tau_con = 3600.
+ffallv_lsc = 1.
+ffallv_con = 1.
+coef_eva = 0.00002
+reevap_ice = y
+iflag_cldcon = 3
+iflag_pdf = 1
+fact_cldcon = 1.
+#facttemps = 1.e-4
+facttemps = 0.0001
+ok_newmicro = y
+ratqsbas = 0.005
+ratqshaut = 0.33
+rad_froid = 35
+rad_chau1=12
+rad_chau2=11
+#ksta_ter=1.e-7
+ksta_ter=0.0000001
+#
+# parametres outputs
+#
+#niveau de sortie "hf" lev_histhf avec
+# - lev_histhf=3 => defaut
+# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele
+lev_histhf=3
+#niveau de sortie "day" lev_histday
+# - lev_histday=2 => defaut
+# - lev_histday=3 => + champs 3D => F. Lott
+# - lev_histday=4 => + champs sous-surfaces
+lev_histday=2
+#niveau de sortie "mth" lev_histmth avec
+# - lev_histmth=2 => defaut
+# - lev_histmth=3 => albedo, rugosite sous-surfaces
+# - lev_histmth=4 => champs tendances 3d
+lev_histmth=4
+#
+# parametres climatique
+#
+R_ecc = 0.016715
+R_peri = 102.7
+R_incl = 23.441
+solaire = 1365.
+co2_ppm = 348.
+CH4_ppb = 1650.
+N2O_ppb = 306.
+CFC11_ppt = 280.
+CFC12_ppt = 484.
+#
+# parametres simulateur ISCCP
+#
+top_height = 3
+#overlap = 1, 2 ou 3
+overlap = 3
+#cdmmax
+#cdmmax = 2.5E-3
+cdmmax = 0.0025
+#cdhmax
+#cdhmax = 2.0E-3
+cdhmax = 0.002
+#autres flags
+iflag_pbl = 1
+iflag_thermals = 0
+iflag_ratqs=0
+iflag_clos=0
+#
+# parametre chimie/aerosols
+#
+config_inca=chem
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/run.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/run.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/run.def	(revision 1242)
@@ -0,0 +1,25 @@
+## $Id: run.def 71 2007-06-21 14:26:03Z acosce $
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+## activation du calcul d equilibrage de charge
+adjust=_adjust_
+## Type de calendrier utilise
+## valeur possible: earth_360d (defaut), earth_365d, earth_366d
+calend=_calend_
+## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+dayref=_dayref_
+##  Annee de l'etat  initial (   avec  4  chiffres   )
+anneeref=_anneeref_
+## raz_date
+raz_date=_raz_date_
+## Nombre de jours d'integration
+nday=_nday_
+## periode de sortie des variables de controle (en pas)
+iconser= 480
+## periode d'ecriture du fichier histoire (en jour)
+iecri=10
+## periode de stockage fichier histmoy (en jour)
+periodav=1.
+## unite de sortie des impressions
+lunout=6
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/traceur.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/PARAM/traceur.def	(revision 1242)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/POST/monitoring01_inca.cfg
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/POST/monitoring01_inca.cfg	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/POST/monitoring01_inca.cfg	(revision 1242)
@@ -0,0 +1,50 @@
+#**************************************************************
+# 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
+#   - files patterns, operations,title,units must be enclosed with character '"'
+#   - 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=F0D5F4
+smooth=12
+#-----------------------------------------------------------------------------------------------------------------
+#  field | files patterns | files additionnal | operations | title | units | calcul of area
+#-----------------------------------------------------------------------------------------------------------------
+NOX_surf_global		| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=19]+NO2[d=2,k=19])" | " NOX a la surface"   | "VMR" | "aire[d=3]"
+CH4_surf_global		| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=19]"                | " CH4 a la surface"   | "VMR" | "aire[d=2]"
+CO_surf_global		| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=19]"                 | " CO a la surface"    | "VMR" | "aire[d=2]"
+O3_surf_global		| "_O3"   	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=19]"                 | " O3 a la surface"    | "VMR" | "aire[d=2]"
+HNO3_surf_global	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=19]"               | " HNO3 a la surface"  | "?" 	 | "aire[d=2]"
+NOX_545HPa_global    	| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=12]+NO2[d=2,k=12])" | " NOX a 545 HPa"      | "VMR" | "aire[d=3]"
+CH4_545HPa_global    	| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=12]"                | " CH4 a 545 HPa"      | "VMR" | "aire[d=2]"
+CO_545HPa_global     	| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=12]"                 | " CO a 545 HPa"       | "VMR" | "aire[d=2]"
+O3_545HPa_global     	| "_O3"    	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=12]"                 | " O3 a 545 HPa"       | "VMR" | "aire[d=2]"
+HNO3_545HPa_global   	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=12]"               | " HNO3 a 545 HPa"     | "?"   | "aire[d=2]"
+
+#-----------------------------------------------------------------------------------------------------------------
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/config.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/config.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC/config.card	(revision 1242)
@@ -0,0 +1,106 @@
+#- $Id: config.card 69 2007-06-21 14:08:22Z acosce $
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=LI_N_00
+LongName='LMDz4INCA with New scripts'
+TagName=LMDZINCA/NMHC
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Début et fin de Job
+#D-- "YYYY-MM-DD"
+DateBegin=2000-01-01
+DateEnd=2000-12-31
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=1M
+#============================
+#-- PBS Class (required on brodie mono/multi)
+#JobClass=mono
+JobClass=multi
+#============================
+#-- Total Number of Processors (minimum is 2 for a coupled configuration)
+#JobNumProcTot=2
+JobNumProcTot=4
+#============================
+#-- Run Options : here for oasis coupling
+JobRunOptions='"-np "${config_UserChoices_JobNumProcTot}'
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ4)
+CHM= (inca, INCA3) 
+
+#========================================================================
+#D-- Executable -
+[Executable]
+Name=gcm.e
+#D- For each component, Real name of executable, Name of executable for oasis
+ATM= (gcm.e, gcm.e)
+CHM= (inca.dat, inca.dat)
+
+#========================================================================
+#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=1869-12-30
+#D- Define restart simulation name (=> JOB_OS)
+RestartJobName=JobName
+#D- Path Server Group Login (=> PSGL)
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/NMHC/
+# Attention login depend de la machine
+
+#========================================================================
+#D-- Post -
+[Post]
+#D- Do we rebuild parallel output, this flag determines
+#D- frequency of rebuild submission (use NONE for DRYRUN=3)
+RebuildFrequency=1M
+#D- Do we rebuild parallel output from archive
+RebuildFromArchive=true
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+TimeSeriesFrequency=1M
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average
+SeasonalFrequency=5Y
+
+#========================================================================
+#D-- ATM -
+[ATM]
+#
+WriteFrequency="1M 1D HF"
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=2000-01-01
+# Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/NMHC/
+# Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
+#D-- CHM -
+[CHM]
+WriteFrequency="1M"
+Restart=n
+#D-- Last day of the experience used as restart
+RestartDate=2000-01-01
+#D-- Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/NMHC/
+#D-- Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/inca.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/inca.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/inca.card	(revision 1242)
@@ -0,0 +1,54 @@
+#- $Id: inca.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+[BoundaryFiles]
+List=	    (/dmnfs/cont003/p24data/ECMWF320x160/AN${year}/165_${year}${month}.nc      , u10mec.nc)\
+	    (/dmnfs/cont003/p24data/ECMWF320x160/AN${year}/166_${year}${month}.nc      , v10mec.nc)\
+	    (/dmnfs/cont003/p24data/ECMWF96x72/AN${year}/u10m_ecmwf_${year}${month}.nc , u10m.nc  )\
+	    (/dmnfs/cont003/p24data/ECMWF96x72/AN${year}/v10m_ecmwf_${year}${month}.nc , v10m.nc  )
+
+ListNonDel= (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/o3clim.nc               ,   o3clim.nc      )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/npp.nc	                 ,   npp.nc	    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/aircraft.nc             ,   aircraft.nc    )\
+	    (${R_INIT}/CHM/LMDZINCA/CH4/INCA${RESOL_CHM}/so4.nc	                 ,   so4.nc	    )\
+	    (${R_INIT}/CHM/LMDZINCA/NMHC_AER/INCA${RESOL_CHM}/sflx_lmdz_phy.present.nc,   sflx.nc        )\
+	    (${R_INIT}/CHM/LMDZINCA/NMHC/phototable.dat_V4.1	                 ,   phototable.dat )\
+	    (${R_INIT}/CHM/LMDZINCA/NMHC/INCA${RESOL_CHM}/landuse.nc             ,   landuse.nc     )\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/rhvlmdz.txt             ,   rhv.dat        )\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/wthlmdz.txt             ,   wth.dat	    )\
+            (${R_INIT}/CHM/LMDZINCA/AER/INCA${RESOL_CHM}/clylmdz.txt             ,   cly.dat	    )\
+            (${R_INIT}/CHM/LMDZINCA/AER/rhvEC11.09092002.txt	                 ,   rhvEC.txt	    )\
+            (${R_INIT}/CHM/LMDZINCA/AER/wthEC11.09092002.txt	                 ,   wthEC.txt	    )\
+            (${R_INIT}/CHM/LMDZINCA/AER/clyEC11.09092002.txt	                 ,   clyEC.txt	    )
+
+[ParametersFiles]
+List=	
+
+[RestartFiles]
+# List restart that have to be saved/restored each loop (file out, saved, and in) :
+List= (precipinsoil.dat, precipinsoil.dat, precipinsoil.dat) 
+
+
+
+[OutputText]
+List=  
+# avec la // : out_orchidee_*
+
+[OutputFiles]
+List=	(inca_avgr.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_inca_avgr.nc,     Post_1M_inca_avgr ),\
+	(forcage.nc,      ${R_OUT_CHM_O_M}/${PREFIX}_1M_forcage.nc,      NONE)
+
+
+[Post_1M_inca_avgr]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter)
+TimeSeriesVars2D = (OD550_ASBCM, OD550_ASPOMM, OD550_ASSO4M, OD550_CSSO4M, OD550_SSSSM, OD550_ASSSM, OD550_CSSSM, OD550_CIDUSTM, OD550_AIBCM, OD550_AIPOMM)
+ChunckJob2D = NONE
+TimeSeriesVars3D = (CH4,NO,NO2,CO,O3,HNO3)
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/inca.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/inca.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/inca.driver	(revision 1242)
@@ -0,0 +1,40 @@
+#- $Id: inca.driver 119 2007-07-31 13:56:33Z acosce $
+#!/bin/ksh
+
+#D- Driver du script pour INCA
+
+#-----------------------------------------------------------------
+function CHM_Initialize
+{
+    IGCM_debug_PushStack "CHM_Initialize"
+
+    RESOL_CHM=$( echo $RESOL_ATM | awk "-FD" '{print $2}' )
+    echo 'resol chimie', $RESOL_CHM
+    echo "INCA Tag : LMDZ-INCA" ${config_UserChoices_TagName}
+
+    IGCM_debug_PopStack "CHM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function CHM_Update
+{
+    IGCM_debug_PushStack "CHM_Update"
+    IGCM_debug_PopStack "CHM_Update"
+}
+
+#-----------------------------------------------------------------
+function CHM_Finalize
+{
+    IGCM_debug_PushStack "CHM_Finalize"
+
+    echo FINALIZE CHM !!!
+#
+#    IGCM_sys_Mkdir ${R_OUT_EXE}
+#    cd ${MODIPSL}
+#    tar cf src_${config_UserChoices_JobName}.tar  modeles
+#    IGCM_sys_Mv src_${config_UserChoices_JobName}.tar ${R_OUT_EXE}/.
+#    cd -
+#
+
+    IGCM_debug_PopStack "CHM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/lmdz.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/lmdz.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/lmdz.card	(revision 1242)
@@ -0,0 +1,45 @@
+#- $Id: lmdz.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+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
+
+[InitialStateFiles]
+# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job
+List= 	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/start.nc, start.nc) \
+	(${R_INIT}/CHM/${config_UserChoices_TagName}/INCA${RESOL_CHM}/startphy.nc, startphy.nc)
+
+
+[BoundaryFiles]
+List=   (${R_BC}/CHM/LMDZORINCA/${RESOL_ATM}_leap/EL-${RESOL_ATM}_${year}1231_limit.nc , limit.nc) 
+ListNonDel= ()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def, .), \
+	(${SUBMIT_DIR}/PARAM/gcm.def, gcm.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, Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat)
+
+
+[OutputFiles]
+List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,   Post_1M_histmth)	
+
+[Post_1M_histmth]
+Patches= (Patch_20091030_histcom_time_axis)
+GatherWithInternal = (lon, lat, presnivs, time_counter, aire)
+TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux_ter, taux_sic, taux_lic, taux_oce, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, 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)
+ChunckJob2D = NONE
+TimeSeriesVars3D = ()
+ChunckJob3D = NONE
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/lmdz.driver	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/COMP/lmdz.driver	(revision 1242)
@@ -0,0 +1,160 @@
+#- $Id: lmdz.driver 841 2009-12-08 09:48:13Z acosce $
+#!/bin/ksh
+
+#-----------------------------------------------------------------
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' )
+    case ${RESOL_ATM} in
+	LMD7245)   BandsResol=72x45x19   ;;
+	LMD9671)   BandsResol=96x71x19   ;;
+	LMD9672)   BandsResol=96x72x19   ;;
+	LMD9695)   BandsResol=96x95x19   ;;
+	LMD14496)  BandsResol=144x96x19  ;;
+	LMD144142) BandsResol=144x142x19 ;;
+    esac
+
+
+    ##-- Calendar type for LMDZ
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=earth_366d;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d;;
+        360d)
+            CalendarTypeForLmdz=earth_360d;;
+        *)
+            CalendarTypeForLmdz=earth_360d
+    esac
+
+    ##- Default number of processor for lmdz
+    NUM_PROC_ATM=1
+    if [ X"${BATCH_NUM_PROC_TOT}" != X ] ; then
+	NUM_PROC_ATM=${BATCH_NUM_PROC_TOT}
+    fi
+
+    ##--Frequency purpose ....
+    ##--  Initialisation  ....
+    OK_instan=n
+    ##-- 
+    OK_journe=n
+    OK_mensuel=y
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) OK_mensuel=y ;;
+	5D|5d|1D|1d) OK_journe=n ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    5D|5d|1D|1d) OK_journe=n ;;
+	esac
+    done
+
+    ##--Variables used by LMDZ --
+    PAT_INST=$(     grep 'OK_instan'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_JOUR=$(     grep 'OK_journe'   ${SUBMIT_DIR}/PARAM/physiq.def )
+    PAT_MOIS=$(     grep 'OK_mensuel'  ${SUBMIT_DIR}/PARAM/physiq.def )
+
+    PAT_iphysiq=$(  grep 'iphysiq'     ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_iperiod=$(  grep 'iperiod'     ${SUBMIT_DIR}/PARAM/gcm.def | tail -1)
+    PAT_day_step=$( grep 'day_step'    ${SUBMIT_DIR}/PARAM/gcm.def )
+    PAT_ecritphy=$( grep 'ecritphy'    ${SUBMIT_DIR}/PARAM/gcm.def )
+
+
+    ## Read LMDZ_NbPeriod_adjust option in lmdz.card 
+    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_NbPeriod_adjust
+    eval LMDZ_NbPeriod_adjust=\${lmdz_UserChoices_LMDZ_NbPeriod_adjust} > /dev/null 2>&1
+
+    if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+        IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices LMDZ_Bands_file_name
+        eval LMDZ_Bands_file_name=\${lmdz_UserChoices_LMDZ_Bands_file_name} > /dev/null 2>&1
+    fi
+
+
+    ##-- This could be define in lmdz.card, inside section [UserChoices]
+    ##-- Otherwise we get the value in *.def
+    iperiod=$(grep iperiod ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}' | tail -1)
+    iphysiq=$(grep iphysiq ${SUBMIT_DIR}/PARAM/gcm.def | awk -F= '{print $2}')
+    # day_step : number of steps per day (multiple of iperiod)
+    (( day_step = 48 * iphysiq ))
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    ecritphy=${PeriodLengthInDays}
+
+    ##-- 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
+
+    if [ ${CumulPeriod} -eq 1 ] ; then
+        if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+            IGCM_sys_Get ${LMDZ_Bands_file_name} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+            IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0
+            IGCM_debug_Print 1 "Bands file forced to ${LMDZ_Bands_file_name} and stored in ${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_0"
+        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 )) 
+
+    [ ${CumulPeriod} -gt 1 ] && IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat
+
+    ## Mise en forme du fichier physiq.def
+    sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   	\
+	-e "s/OK_journe=.*/OK_journe=${OK_journe}/"   	\
+	-e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" 	\
+	physiq.def > physiq.def.tmp
+    IGCM_sys_Mv physiq.def.tmp physiq.def
+
+    ## Mise en forme du fichier gcm.def
+    sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" 	\
+	-e "s/${PAT_day_step}/day_step=${day_step}/" 	\
+	-e "s/${PAT_iperiod}/iperiod=${iperiod}/" 	\
+	-e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" 	\
+	gcm.def > gcm.def.tmp
+    IGCM_sys_Mv gcm.def.tmp gcm.def
+
+    greg_dat=$( IGCM_date_ConvertFormatToGregorian $PeriodDateBegin ) 
+    jul_dat=$( IGCM_date_ConvertGregorianDateToJulian $greg_dat )
+    nbjour=$( expr \( $jul_dat \% 1000 \) ) 
+    yractu=$( expr $greg_dat / 10000  )
+    echo  "dayref = " $nbjour  "year actu = " $yractu
+
+    ## Mise en forme du fichier run.def
+    sed -e "s/_dayref_/${nbjour}/" 			\
+	-e "s/_anneeref_/${yractu}/" 			\
+        -e "s/_calend_/${CalendarTypeForLmdz}/"        	\
+	-e "s/_nday_/${PeriodLengthInDays}/" 		\
+	-e "s/_raz_date_/${RAZ_DATE}/" 			\
+	-e "s/_adjust_/${LMDZ_adjust}/"                	\
+	run.def > run.def.tmp
+    IGCM_sys_Mv run.def.tmp run.def
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && IGCM_sys_Put_Out Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/Bands_96x72x19_2prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/Bands_96x72x19_2prc.dat	(revision 1242)
@@ -0,0 +1,8 @@
+ 0  37
+ 1  36
+ 0  37
+ 1  36
+ 0  34
+ 1  39
+ 0  3456
+ 1  3362
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/Bands_96x72x19_4prc.dat
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/Bands_96x72x19_4prc.dat	(revision 1242)
@@ -0,0 +1,16 @@
+ 0  8
+ 1  28
+ 2  28
+ 3  9
+ 0  19
+ 1  18
+ 2  18
+ 3  18
+ 0  9
+ 1  31
+ 2  25
+ 3  8
+ 0  1762
+ 1  1674
+ 2  1653
+ 3  1729
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/gcm.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/gcm.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/gcm.def	(revision 1242)
@@ -0,0 +1,74 @@
+## $Id: gcm.def 71 2007-06-21 14:26:03Z acosce $
+## 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
+## periode de la dissipation (en pas)
+idissip=10
+## 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=3600.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=3600.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## guide ou pas
+ok_guide=.false.
+## utilisation de restart de la version LMDZ3
+old_restart=
+## avec ou sans physique
+##physic=n
+iflag_phys=1
+## periode de la physique (en pas)
+iphysiq=10
+## frequence (en  jours ) de l'ecriture du fichier histphy               
+ecritphy=1
+## frequence (en jours ) de l'ecriture du fichier histmth
+ecrit_mth=1
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=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
+## Nombre  d'appels des routines de rayonnements ( par jour)
+nbapp_rad=12
+##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect)
+iflag_con=30
+## 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/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/offline.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/offline.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/offline.def	(revision 1242)
@@ -0,0 +1,10 @@
+## $Id: offline.def 71 2007-06-21 14:26:03Z acosce $
+T
+4
+T
+-2.
+48.1
+1
+T
+6
+2
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/orchidee.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/orchidee.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/orchidee.def	(revision 1242)
@@ -0,0 +1,44 @@
+## $Id: orchidee.def 71 2007-06-21 14:26:03Z acosce $
+#
+# SECHIBA
+#
+STOMATE_OK_CO2=FALSE
+# STOMATE_OK_STOMATE is not set
+# STOMATE_OK_DGVM is not set
+# STOMATE_WATCHOUT is not set
+#SECHIBA_restart_in=default
+SECHIBA_restart_in=start_sech.nc
+SECHIBA_rest_out=sechiba_rest.nc
+SECHIBA_reset_time=y
+SECHIBA_reset_time is not set
+OUTPUT_FILE=sechiba_out.nc
+WRITE_STEP=2592000
+SECHIBA_HISTLEVEL=10
+STOMATE_OUTPUT_FILE=stomate_history.nc
+STOMATE_HIST_DT=10.
+STOMATE_HISTLEVEL=0
+SECHIBA_DAY=0.0
+SECHIBA_ZCANOP=0.5
+DT_SLOW=86400.
+# IMPOSE_VEG is not set
+VEGETATION_FILE=carteveg5km.nc
+DIFFUCO_LEAFCI=233.
+CONDVEG_SNOWA=default
+# IMPOSE_AZE is not set
+SOILALB_FILE=soils_param.nc
+SOILTYPE_FILE=soils_param.nc 
+ENERBIL_TSURF=280.
+HYDROL_SNOW=0.0
+HYDROL_SNOWAGE=0.0
+HYDROL_SNOWICE=0.0
+HYDROL_SNOWICEAGE=0.0
+HYDROL_HDRY=1.0
+HYDROL_HUMR=1.0
+HYDROL_BQSB=default
+HYDROL_GQSB=0.0
+HYDROL_DSG=0.0
+HYDROL_DSP=default
+HYDROL_QSV=0.0
+THERMOSOIL_TPRO=280.
+RIVER_ROUTING=y
+ROUTING_FILE=routing.nc
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/physiq.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/physiq.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/physiq.def	(revision 1242)
@@ -0,0 +1,93 @@
+##- $Id: physiq.def 71 2007-06-21 14:26:03Z acosce $
+type_ocean=force 
+VEGET=n
+OK_journe=n
+OK_mensuel=y
+OK_instan=n
+if_ebil=0
+#
+# parametres KE
+#
+epmax = .99
+ok_adj_ema = n
+iflag_clw = 1
+# 
+# parametres nuages
+#
+cld_lc_lsc = 0.00026
+cld_lc_con = 0.00026
+cld_tau_lsc = 3600.
+cld_tau_con = 3600.
+ffallv_lsc = 1.
+ffallv_con = 1.
+coef_eva = 0.00002
+reevap_ice = y
+iflag_cldcon = 3
+iflag_pdf = 1
+fact_cldcon = 1.
+#facttemps = 1.e-4
+facttemps = 0.0001
+ok_newmicro = y
+ratqsbas = 0.005
+ratqshaut = 0.33
+rad_froid = 35
+rad_chau1=12
+rad_chau2=11
+#ksta_ter=1.e-7
+ksta_ter=0.0000001
+#
+# parametres outputs
+#
+#niveau de sortie "hf" lev_histhf avec
+# - lev_histhf=3 => defaut
+# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele
+lev_histhf=3
+#niveau de sortie "day" lev_histday
+# - lev_histday=2 => defaut
+# - lev_histday=3 => + champs 3D => F. Lott
+# - lev_histday=4 => + champs sous-surfaces
+lev_histday=2
+#niveau de sortie "mth" lev_histmth avec
+# - lev_histmth=2 => defaut
+# - lev_histmth=3 => albedo, rugosite sous-surfaces
+# - lev_histmth=4 => champs tendances 3d
+lev_histmth=4
+#
+# parametres climatique
+#
+R_ecc = 0.016715
+R_peri = 102.7
+R_incl = 23.441
+solaire = 1365.
+co2_ppm = 348.
+CH4_ppb = 1650.
+N2O_ppb = 306.
+CFC11_ppt = 280.
+CFC12_ppt = 484.
+#
+# parametres simulateur ISCCP
+#
+top_height = 3
+#overlap = 1, 2 ou 3
+overlap = 3
+#cdmmax
+#cdmmax = 2.5E-3
+cdmmax = 0.0025
+#cdhmax
+#cdhmax = 2.0E-3
+cdhmax = 0.002
+#autres flags
+iflag_pbl = 1
+iflag_thermals = 0
+iflag_ratqs=0
+iflag_clos=0
+#
+# parametre chimie/aerosols
+#
+config_inca=aero
+aerosol_couple=n
+flag_aerosol=1
+new_aod=.true.
+ok_ade=n
+ok_aie=n
+
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/run.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/run.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/run.def	(revision 1242)
@@ -0,0 +1,25 @@
+## $Id: run.def 71 2007-06-21 14:26:03Z acosce $
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+## activation du calcul d equilibrage de charge
+adjust=_adjust_
+## Type de calendrier utilise
+## valeur possible: earth_360d (defaut), earth_365d, earth_366d
+calend=_calend_
+## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+dayref=_dayref_
+##  Annee de l'etat  initial (   avec  4  chiffres   )
+anneeref=_anneeref_
+## raz_date
+raz_date=_raz_date_
+## Nombre de jours d'integration
+nday=_nday_
+## periode de sortie des variables de controle (en pas)
+iconser= 480
+## periode d'ecriture du fichier histoire (en jour)
+iecri=10
+## periode de stockage fichier histmoy (en jour)
+periodav=1.
+## unite de sortie des impressions
+lunout=6
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/traceur.def	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/PARAM/traceur.def	(revision 1242)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/POST/monitoring01_inca.cfg
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/POST/monitoring01_inca.cfg	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/POST/monitoring01_inca.cfg	(revision 1242)
@@ -0,0 +1,51 @@
+#**************************************************************
+# 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
+#   - files patterns, operations,title,units must be enclosed with character '"'
+#   - 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=F0D5F4
+smooth=12
+#-----------------------------------------------------------------------------------------------------------------
+#  field | files patterns | files additionnal | operations | title | units | calcul of area
+#-----------------------------------------------------------------------------------------------------------------
+ OD550_AERO_global | "OD550_ASBCM OD550_ASPOMM OD550_ASSO4M OD550_CSSO4M OD550_SSSSM OD550_ASSSM OD550_CSSSM OD550_CIDUSTM OD550_AIBCM OD550_AIPOMM" | LMDZ4.0_9695_grid.nc | "(OD550_ASBCM[d=1]+ OD550_ASPOMM[d=2]+ OD550_ASSO4M[d=3]+ OD550_CSSO4M[d=4]+ OD550_SSSSM[d=5]+ OD550_ASSSM[d=6]+ OD550_CSSSM[d=7]+ OD550_CIDUSTM[d=8]+ OD550_AIBCM[d=9]+ OD550_AIPOMM[d=10])" | "total  aerosol optical depth @550nm"  |"(-)" | "aire[d=11]"
+NOX_surf_global		| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=19]+NO2[d=2,k=19])" | " NOX a la surface"   | "VMR" | "aire[d=3]"
+CH4_surf_global		| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=19]"                | " CH4 a la surface"   | "VMR" | "aire[d=2]"
+CO_surf_global		| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=19]"                 | " CO a la surface"    | "VMR" | "aire[d=2]"
+O3_surf_global		| "_O3"   	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=19]"                 | " O3 a la surface"    | "VMR" | "aire[d=2]"
+HNO3_surf_global	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=19]"               | " HNO3 a la surface"  | "?" 	 | "aire[d=2]"
+NOX_545HPa_global    	| "NO NO2" 	|  LMDZ4.0_9695_grid.nc | "(NO[d=1,k=12]+NO2[d=2,k=12])" | " NOX a 545 HPa"      | "VMR" | "aire[d=3]"
+CH4_545HPa_global    	| "CH4"    	|  LMDZ4.0_9695_grid.nc | "CH4[d=1,k=12]"                | " CH4 a 545 HPa"      | "VMR" | "aire[d=2]"
+CO_545HPa_global     	| "CO"     	|  LMDZ4.0_9695_grid.nc | "CO[d=1,k=12]"                 | " CO a 545 HPa"       | "VMR" | "aire[d=2]"
+O3_545HPa_global     	| "_O3"    	|  LMDZ4.0_9695_grid.nc | "O3[d=1,k=12]"                 | " O3 a 545 HPa"       | "VMR" | "aire[d=2]"
+HNO3_545HPa_global   	| "HNO3"   	|  LMDZ4.0_9695_grid.nc | "HNO3[d=1,k=12]"               | " HNO3 a 545 HPa"     | "?"   | "aire[d=2]"
+
+#-----------------------------------------------------------------------------------------------------------------
Index: CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/config.card
===================================================================
--- CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/config.card	(revision 1242)
+++ CONFIG/LMDZINCA/tags/LMDZINCA_v3_1/EXP_NMHC_AER/config.card	(revision 1242)
@@ -0,0 +1,106 @@
+#- $Id: config.card 69 2007-06-21 14:08:22Z acosce $
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=LI_NA_00
+LongName='LMDz4INCA with New scripts'
+TagName=LMDZINCA/NMHC_AER
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Début et fin de Job
+#D-- "YYYY-MM-DD"
+DateBegin=2000-01-01
+DateEnd=2000-12-31
+#============================
+#D-- 1Y, 1M, 5D, 1D
+PeriodLength=1M
+#============================
+#-- PBS Class (required on brodie mono/multi)
+#JobClass=mono
+JobClass=multi
+#============================
+#-- Total Number of Processors (minimum is 2 for a coupled configuration)
+#JobNumProcTot=2
+JobNumProcTot=4
+#============================
+#-- Run Options : here for oasis coupling
+JobRunOptions='"-np "${config_UserChoices_JobNumProcTot}'
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ4)
+CHM= (inca, INCA3) 
+
+#========================================================================
+#D-- Executable -
+[Executable]
+Name=gcm.e
+#D- For each component, Real name of executable, Name of executable for oasis
+ATM= (gcm.e, gcm.e)
+CHM= (inca.dat, inca.dat)
+
+#========================================================================
+#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=1869-12-30
+#D- Define restart simulation name (=> JOB_OS)
+RestartJobName=JobName
+#D- Path Server Group Login (=> PSGL)
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/NMHC_AER/
+# Attention login depend de la machine
+
+#========================================================================
+#D-- Post -
+[Post]
+#D- Do we rebuild parallel output, this flag determines
+#D- frequency of rebuild submission (use NONE for DRYRUN=3)
+RebuildFrequency=1M
+#D- Do we rebuild parallel output from archive
+RebuildFromArchive=true
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+TimeSeriesFrequency=1M
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average
+SeasonalFrequency=5Y
+
+#========================================================================
+#D-- ATM -
+[ATM]
+#
+WriteFrequency="1M 1D HF"
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=2000-01-01
+# Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/NMHC_AER/
+# Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
+#D-- CHM -
+[CHM]
+WriteFrequency="1M"
+Restart=n
+#D-- Last day of the experience used as restart
+RestartDate=2000-01-01
+#D-- Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZINCA/NMHC_AER/
+#D-- Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
