Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/AA_make
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/AA_make	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/AA_make	(revision 1738)
@@ -0,0 +1,207 @@
+#- $Id$
+
+all : 
+	if [ -s ./.resol ] ; then gmake `cat .resol` ; else gmake AERxLMD9695 ; fi
+
+AERxLMD9671 : libioipsl liborchidee aer_lmdz9671
+	echo "AERxLMD9671" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=AER" > .chimie
+
+AERxLMD9672 : libioipsl liborchidee aer_lmdz9672
+	echo "AERxLMD9672" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=AER" > .chimie
+
+CH4_AERxLMD9672 : libioipsl liborchidee ch4aer_lmdz9672
+	echo "CH4_AERxLMD9672" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=CH4_AER" > .chimie
+
+CH4xLMD9672 : libioipsl liborchidee ch4_lmdz9672
+	echo "CH4xLMD9672" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=CH4" > .chimie
+
+NMHCxLMD9672 : libioipsl liborchidee nmhc_lmdz9672
+	echo "NMHCxLMD9672" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=NMHC" > .chimie
+
+NMHC_AERxLMD9672 : libioipsl liborchidee nmhcaer_lmdz9672
+	echo "NMHC_AERxLMD9672" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=NMHC_AER" > .chimie
+
+GESxLMD9672 : libioipsl liborchidee ges_lmdz9672
+	echo "GESxLMD9672" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=GES" > .chimie
+
+AERxLMD9695 : libioipsl liborchidee aer_lmdz9695
+	echo "AERxLMD9695" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=AER" > .chimie
+
+CH4_AERxLMD9695 : libioipsl liborchidee ch4aer_lmdz9695
+	echo "CH4_AERxLMD9695" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=CH4_AER" > .chimie
+
+CH4xLMD9695 : libioipsl liborchidee ch4_lmdz9695
+	echo "CH4xLMD9695" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=CH4" > .chimie
+
+NMHCxLMD9695 : libioipsl liborchidee nmhc_lmdz9695
+	echo "NMHCxLMD9695" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=NMHC" > .chimie
+
+NMHC_AERxLMD9695 : libioipsl liborchidee nmhcaer_lmdz9695
+	echo "NMHC_AERxLMD9695" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=NMHC_AER" > .chimie
+
+GESxLMD9695 : libioipsl liborchidee ges_lmdz9695
+	echo "GESxLMD9695" >.resol
+	echo "RESOL_ATM_3D=96x95x19" >>.resol
+	echo "CHEM=GES" > .chimie
+
+AERxLMD144142 : libioipsl liborchidee aer_lmdz144142
+	echo "AERxLMD144142" >.resol
+	echo "RESOL_ATM_3D=144x142x19" >>.resol
+	echo "CHEM=AER" > .chimie
+
+
+NMHCxLMD9695-L39 : libioipsl liborchidee nmhc_lmdz969539
+	echo "NMHCxLMD9695-L39" >.resol
+	echo "RESOL_ATM_3D=96x95x39" >>.resol
+	echo "CHEM=NMHC" > .chimie
+
+AERxLMD9695-L39 : libioipsl liborchidee aer_lmdz969539
+	echo "AERxLMD9695-L39" >.resol
+	echo "RESOL_ATM_3D=96x95x39" >>.resol
+	echo "CHEM=AER" > .chimie
+
+NMHC_AERxLMD9695-L39 : libioipsl liborchidee nmhcaer_lmdz969539
+	echo "NMHC_AERxLMD9695-L39" >.resol
+	echo "RESOL_ATM_3D=96x95x39" >>.resol
+	echo "CHEM=NMHC_AER" > .chimie
+
+GESxLMD9695-L39 : libioipsl liborchidee ges_lmdz969539
+	echo "GESxLMD9695-L39" >.resol
+	echo "RESOL_ATM_3D=96x95x39" >>.resol
+	echo "CHEM=GES" > .chimie
+
+
+libioipsl :
+	(cd ../../modeles/IOIPSL/src ; $(M_K) -f Makefile)
+
+
+liborchidee :
+	(cd ../../modeles/ORCHIDEE/ ; $(M_K) -f Makefile)
+#-Q- platine    (rm -f ../../lib/parallel.mod)
+
+aer_lmdz9671:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie AER -parallel mpi -resol 96x71 -arch $(FCM_ARCH); cp SIMULATIONS/AER/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x71x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x71x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+aer_lmdz9672:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie AER -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/AER/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x72x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+ch4aer_lmdz9672:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie CH4_AER -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/CH4_AER/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x72x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+ch4_lmdz9672:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie CH4 -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/CH4/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x72x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+nmhc_lmdz9672:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie NMHC -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/NMHC/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x72x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+nmhcaer_lmdz9672:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie NMHC_AER -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/NMHC_AER/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x72x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+ges_lmdz9672:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie GES -parallel mpi -arch $(FCM_ARCH); cp SIMULATIONS/GES/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x72x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x72x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+aer_lmdz9695:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie AER -parallel mpi -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/AER/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+ch4aer_lmdz9695:
+	(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/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+ch4_lmdz9695:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie CH4 -parallel mpi -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/CH4/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+nmhc_lmdz9695:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie NMHC -parallel mpi  -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/NMHC/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+nmhcaer_lmdz9695:
+	(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/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+ges_lmdz9695:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie GES -parallel mpi -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/GES/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+aer_lmdz144142:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie AER -parallel mpi -resol 144x142 -arch $(FCM_ARCH); cp SIMULATIONS/AER/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x19 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x19_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 144x142x19 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_144x142x19_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+nmhc_lmdz969539:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie NMHC -parallel mpi  -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/NMHC/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x39 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x39_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x39 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x39_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+aer_lmdz969539:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie AER -parallel mpi  -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/AER/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x39 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x39_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x39 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x39_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+nmhcaer_lmdz969539:
+	(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/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x39 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x39_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x39 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x39_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+ges_lmdz969539:
+	(cd ../../modeles/INCA3; ./makeinca_fcm -chimie GES -parallel mpi  -resol 96x95 -arch $(FCM_ARCH); cp SIMULATIONS/GES/inca.dat ../../bin/inca.dat ; )	
+#	(cd ../../modeles/CREATE_E0; ./makelmdz_fcm -chimie INCA -d 96x95x39 -parallel mpi -arch $(FCM_ARCH) ce0l ; cp bin/ce0l_96x95x39_phylmd_para_inca.e ../../bin/create_etat0_limit.e ; )
+	(cd ../../modeles/LMDZ; ./makelmdz_fcm -cpp ORCHIDEE_NOOPENMP -d 96x95x39 -chimie INCA -v true -parallel mpi  -arch $(FCM_ARCH) gcm ; cp bin/gcm_96x95x39_phylmd_para_orch_inca.e ../../bin/gcm.e ; )
+
+
+clean :
+	(rm -rf $(LIBDIR)/* ; rm -f $(BINDIR)/* ; rm -rf ../../modeles/LMDZ/libo/* ../../modeles/LMDZ/.lock ; rm -rf ../../modeles/INCA3/config ; )
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/AA_make.ldef
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/AA_make.ldef	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/AA_make.ldef	(revision 1738)
@@ -0,0 +1,17 @@
+#- $Id$
+#---------------------------------------------------------------------
+#- Creation des elements relatifs a CPL+ORCA4+LMDZ.3.3 :
+#-      LMDZ.3.3 (lib, executable : gcm.e)
+#-      ORCHIDEE (lib + modules F90)
+#-      ORCA 4 (lib, executable : opa)
+#-      UCL (lib opa)
+#-      IOIPSL (lib + modules F90)
+#-      CPL (lib, libmpi2, executable : oasis)
+#---------------------------------------------------------------------
+SHELL = /bin/sh
+#---------------------------------------------------------------------
+LIBDIR = ../../../../lib
+BINDIR = ../../../../bin
+MODDIR = $(LIBDIR)
+#---------------------------------------------------------------------
+#-
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/inca.card
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/inca.card	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/inca.card	(revision 1738)
@@ -0,0 +1,30 @@
+#- $Id: inca.card 1246 2011-01-17 16:03:35Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+
+[BoundaryFiles]
+List=	    ()
+ListNonDel= ()
+
+[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=	()
+
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/inca.driver
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/inca.driver	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/inca.driver	(revision 1738)
@@ -0,0 +1,40 @@
+#- $Id: inca.driver 184 2007-10-19 13:38:59Z 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/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/lmdz.card
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/lmdz.card	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/lmdz.card	(revision 1738)
@@ -0,0 +1,44 @@
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+
+[BoundaryFiles]
+List=   (${R_INIT}/ATM/${config_UserChoices_TagName}/AMIP/amipbc_sst_360x180_${year}.nc, amipbc_sst_1x1.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/AMIP/amipbc_sic_360x180_${year}.nc, amipbc_sic_1x1.nc), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Ozone/HYBRIDE/v2.clim/tro3_${year}.new.nc, climoz.nc)
+
+ListNonDel= (${R_INIT}/ATM/${config_UserChoices_TagName}/Albedo.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/ECDYN.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/ECPHY.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Relief.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/Rugos.nc, .), \
+	(${R_INIT}/ATM/${config_UserChoices_TagName}/landiceref.nc, .)
+### Please, add the o2a.nc file to fit with the coupled configuration associated
+#	(${R_INIT}/ATM/IPSLCM5A/ORCA2.3xLMD9695/o2a.nc, o2a.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= (NONE)
+
+[OutputText]
+List=   (physiq.def, gcm.def, run.def, traceur.def, guide.def)
+
+[OutputFiles]
+List=	( limit.nc,       ${R_OUT_ATM_O}/Boundary/${config_UserChoices_JobName}_${year}_limit.nc,       NONE), \
+	( climoz_LMDZ.nc, ${R_OUT_ATM_O}/Boundary/${config_UserChoices_JobName}_${year}_climoz_LMDZ.nc, NONE), \
+	( start.nc,       ${R_OUT_ATM_O}/Restart/${config_UserChoices_JobName}_${year}_start.nc,        NONE), \
+	( startphy.nc,    ${R_OUT_ATM_O}/Restart/${config_UserChoices_JobName}_${year}_startphy.nc,     NONE), \
+        ( grilles_gcm.nc, ${R_OUT_ATM_O}/Grid/${config_UserChoices_JobName}_grilles_gcm.nc,      NONE)
+
+
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/lmdz.driver	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/COMP/lmdz.driver	(revision 1738)
@@ -0,0 +1,133 @@
+#!/bin/ksh
+#-----------------------------------------------------------------
+function LMDZ_sed
+{
+    IGCM_debug_PushStack "LMDZ_sed"
+
+    sed -e "s/^${2}\ *=.*/${2}= ${3}/" ${1} > ${1}.tmp
+    RET=$?
+    echo "LMDZ_sed : ${1} ${2} ${3}"
+    \mv ${1}.tmp ${1}
+
+    IGCM_debug_PopStack "LMDZ_sed"
+    return $RET
+}
+
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    [ -f ${SUBMIT_DIR}/../.resol ] && eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x19
+
+    RESOL_ATM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' | awk "-F-" '{print $1}' )
+
+    ##-- Calendar type for LMDZ
+    case ${config_UserChoices_CalendarType} in
+        leap|gregorian)
+            CalendarTypeForLmdz=gregorian;;
+        noleap)
+            CalendarTypeForLmdz=earth_365d;;
+        360d)
+            CalendarTypeForLmdz=earth_360d;;
+        *)
+            CalendarTypeForLmdz=earth_360d
+    esac
+
+    ##- PARAM directory
+#    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices PARAM_DIR
+    PARAM_DIR=${lmdz_UserChoices_PARAM_DIR}
+    echo parameter file directory : ${PARAM_DIR}
+
+    ##--Frequency purpose ....
+    ##--  Initialisation  ....
+    OK_instan=n
+    OK_journe=n
+    OK_mensuel=n
+    ok_hf=n
+
+    case ${config_UserChoices_PeriodLength} in
+	1Y|1y|1M|1m) OK_mensuel=y ;;
+	5D|5d|1D|1d) OK_journe=y ;;
+    esac
+
+    for frequency in ${config_ATM_WriteFrequency} ; do
+	case ${frequency} in
+	    5D|5d|1D|1d) OK_journe=y ;;
+	esac
+	case ${frequency} in
+	    HF|hf) ok_hf=y ;;
+	esac
+    done
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    case ${config_UserChoices_PeriodLength} in
+        *Y|*y) 
+               LMDZ_ecrit_mth=30. 
+               LMDZ_ecrit_ISCCP=30.
+               LMDZ_periodav=30.
+               ;;
+        *)
+               LMDZ_ecrit_mth=${PeriodLengthInDays}.
+               LMDZ_ecrit_ISCCP=${PeriodLengthInDays}.
+               LMDZ_periodav=${PeriodLengthInDays}.
+               ;;
+    esac
+
+    ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def
+    if [ ${CumulPeriod} -eq 1 ] ; then
+	RAZ_DATE=1
+    else
+	RAZ_DATE=0
+    fi
+
+    LMDZ_adjust=n
+
+    ## physiq.def parameters modified from initial physiq_L${RESOL_ATM_Z}
+    LMDZ_sed physiq.def OK_instan   ${OK_instan}
+    LMDZ_sed physiq.def OK_journe   ${OK_journe}
+    LMDZ_sed physiq.def OK_mensuel  ${OK_mensuel}
+    LMDZ_sed physiq.def ok_hf       ${ok_hf}
+    LMDZ_sed physiq.def ecrit_mth   ${LMDZ_ecrit_mth}
+    LMDZ_sed physiq.def ecrit_ISCCP ${LMDZ_ecrit_ISCCP}
+#    LMDZ_sed physiq.def ok_mensuelCOSP  ${LMDZ_COSP_monthly}
+#    LMDZ_sed physiq.def ok_journeCOSP   ${LMDZ_COSP_daily}
+#    LMDZ_sed physiq.def ok_hfCOSP   ${LMDZ_COSP_hf}
+#    LMDZ_sed physiq.def ok_histNMC  "${LMDZ_NMC_monthly}, ${LMDZ_NMC_daily}, ${LMDZ_NMC_hf}"
+
+    ## gcm.def parameters : no change since gcm.def_${RESOL_ATM_3D} is used and already modified
+
+    ## run.def parameters
+    LMDZ_sed run.def dayref   ${InitDay}
+    LMDZ_sed run.def anneeref ${year}
+    LMDZ_sed run.def calend   ${CalendarTypeForLmdz}
+    LMDZ_sed run.def nday     ${PeriodLengthInDays}
+    LMDZ_sed run.def raz_date ${RAZ_DATE}
+    LMDZ_sed run.def periodav ${LMDZ_periodav}
+    LMDZ_sed run.def adjust   ${LMDZ_adjust}
+    LMDZ_sed run.def use_filtre_fft n
+    if ( ${FirstInitialize} ) ; then
+
+    echo cat run
+    cat run.def
+
+    fi 
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/gcm.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/gcm.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/gcm.def	(revision 1738)
@@ -0,0 +1,74 @@
+## $Id: gcm.def 414 2008-10-03 08:59:31Z 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/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/physiq.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/physiq.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/physiq.def	(revision 1738)
@@ -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=aero
+aerosol_couple=n
+flag_aerosol=1
+new_aod=.true.
+ok_ade=n
+ok_aie=n
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/run.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/run.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/run.def	(revision 1738)
@@ -0,0 +1,24 @@
+## $Id: run.def 71 2007-06-21 14:26:03Z 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
+grilles_gcm_netcdf=y
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/traceur.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/PARAM/traceur.def	(revision 1738)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/create_sst_anomaly.ksh
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/create_sst_anomaly.ksh	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/create_sst_anomaly.ksh	(revision 1738)
@@ -0,0 +1,100 @@
+#!/bin/ksh
+set -vx
+#
+# create_sst_anomaly.ksh
+#
+# This script will correct the sea-surface temperature (SST) from a coupled model simulation 
+# using the difference between an AMIP reference SST field and an historical coupled model 
+# reference SST field. The two reference SST fields are both monthly mean over the period 1979-2005.
+# 
+# The output is stored in the directory out_anomaly, set below, and is directly usable with 
+# the program ce0l (create_etat0_limit) in LMDZ5 revision 1508 or later. For each year there 
+# are two output files : one at model original grid and one at a higer resolution grid(with 
+# suffix _HR.nc). If the SST are to be used at a higer resolution grid or with a zoom it is 
+# better to use the high resolution SST (_HR.nc). 
+#
+# 1 - Create reference file
+# 2 - Add reference file to the scenario simulation
+# 3 - Interpolate to higer regular resolution
+# 
+#
+# Author and contact Josefine Ghattas, IPSL, 2011
+#
+########## User definitions
+# Three time series are needed : amip_TS, hist_TS and simul_TS. They should all contain 
+# the same variable tsol_oce at the same model grid resolution with monthly values.
+
+# amip_TS : AMIP SST containg at least the reference period 1979-2005
+amip_TS=/dmnfs/cont003/p86musat/IGCM_OUT/LMDZOR/PROD/amip/v3.amip1/ATM/Analyse/TS_MO/v3.amip1_19790101_20091231_1M_tsol_oce.nc
+
+# hist_TS : model output SST from reference simulation containg at least the reference period 1979-2005
+hist_TS=/dmnfs/cont003/p86denv/IGCM_OUT/IPSLCM5A/PROD/historical/v3.historical1/ATM/Analyse/TS_MO/v3.historical1_18500101_20051231_1M_tsol_oce.nc
+
+# simul_TS : model output SST from the simulation to be corrected containing at least the years [year_begin,year_end] to be set below
+simul_TS=/dmnfs/cont003/p86denv/IGCM_OUT/IPSLCM5A/PROD/rcp45/v3.rcp45.1/ATM/Analyse/TS_MO/v3.rcp45.1_20060101_21651231_1M_tsol_oce.nc
+
+# year_begin and year_end : interval for simul_TS
+year_begin=2006
+year_end=2095
+
+# outname : basename for output files
+outname=v3.rcp45.1_tsol_oce_
+
+# out_simul : directory where the original simulation will be stored splitted in years
+out_simul=$SCRATCHDIR/OUTTEST/simul
+
+# out_anomaly : directory where the final simulation with anomaly will be created and stored
+out_anomaly=$SCRATCHDIR/OUTTEST/anomaly
+
+# rundir : directory where temporary files will be created and stored
+rundir=$SCRATCHDIR/OUTTEST/rundir
+
+########## End User definitions
+
+########## 0 Create output and run directories
+mkdir -p $out_simul
+mkdir -p $out_anomaly
+
+mkdir -p $rundir 
+cd $rundir
+
+
+########## 1 Create reference file
+# extract reference period from AMIP SST and create reference file
+cdo seldate,1979-01-01,2005-12-31 $amip_TS amip_ts_ref.nc
+cdo ymonmean amip_ts_ref.nc amip_ref.nc
+
+# extract reference period from histroical simulation SST and create reference file
+cdo seldate,1979-01-01,2005-12-31 $hist_TS hist_ts_ref.nc
+cdo ymonmean hist_ts_ref.nc hist_ref.nc
+
+# create difference file AMIP - HIST
+cdo sub amip_ref.nc hist_ref.nc refdiff.nc
+
+
+########## 2 Add reference file to the scenario simulation
+# split time serie into yearly files
+cdo splityear $simul_TS $out_simul/${outname}
+
+# add refdiff anomalies to each year of RCP4.5
+year=$year_begin
+while [[ ${year} -le ${year_end} ]] ; do
+   cdo add $out_simul/${outname}$year.nc2 refdiff.nc $out_anomaly/${outname}anomaly_amip_${year}.nc
+
+   # next year
+   year=`expr $year + 1`
+done
+
+########## 3 Interpolate to higer regular resolution
+year=$year_begin
+while [[ ${year} -le ${year_end} ]] ; do
+   cdo remapbic,r720x360 $out_anomaly/${outname}anomaly_amip_$year.nc $out_anomaly/${outname}anomaly_amip_${year}_HR.nc
+
+   # next year
+   year=`expr $year + 1`
+done
+
+
+exit
+
+
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/era2gcm.ksh
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/era2gcm.ksh	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/era2gcm.ksh	(revision 1738)
@@ -0,0 +1,78 @@
+#!/bin/ksh
+#set -xv
+# Function that does a linear interpolation of the variable in 
+# filein towards the destionation grid in gridfile. Results in are 
+# stored in fileout.
+#
+# Exemple for use :
+# ./era2gcm.ksh grilles_gcm.nc u_filein.nc u_fileout.nc u
+#
+###################################################################
+# O - Input arguments
+###################################################################
+# gridfile : File containing the target grid (grilles_gcm.nc)
+    gridfile=$1
+# filein : File containing the variable at source grid
+    filein=$2
+# fileout : Output file with the variable at target grid
+    fileout=$3
+# varin : Variable name in source file
+    varin=$4
+
+###################################################################
+# 1 - Define variables
+###################################################################
+# Find number of time step
+    tmax=`ncdump -h ${filein} | grep time | head -1 | awk ' { print $6 } '`
+    tmax=$( echo ${tmax} | awk '-F(' '{print $2}' )
+    
+# Choose grid and output variable name 
+    if [ ${varin} = 'u' ] || [ ${varin} = 'vitu' ] ; then 
+	varout=uwnd
+	grille='grille_u'
+    elif [ ${varin} = 'v' ] || [ ${varin} = 'vitv' ] ; then 
+	varout=vwnd
+	grille='grille_v'
+    elif [ ${varin} = 'ta' ] || [ ${varin} = 'temp' ] ; then 
+	varout=air
+	grille='grille_t'
+    elif [ ${varin} = 'r' ] || [ ${varin} = 'hur' ] ; then 
+	varout=rh
+	grille='grille_t'
+    elif [ ${varin} = 'msl' ] ; then 
+	varout=sp
+	grille='grille_t'
+    else 
+	echo Error : ${varin} unknown!!!!
+	exit
+    fi 
+
+###################################################################
+# 2 - Create ferret script for interpolation
+###################################################################
+    cat << eod >  ${varin}.jnl
+set memory/size=50
+
+use "${gridfile}"
+use "${filein}"
+let ${varout}=${varin}
+
+define grid/like=${varout}[d=2]/x=grille_u[d=1]/y=grille_u[d=1] grille_u
+define grid/like=${varout}[d=2]/x=grille_v[d=1]/y=grille_v[d=1] grille_v
+define grid/like=${varout}[d=2]/x=grille_v[d=1]/y=grille_u[d=1] grille_t
+
+save/clobber/file="${fileout}" ${varout}[d=2,g=${grille},l=1]
+repeat/l=1:${tmax} save/file="${fileout}"/append ${varout}[d=2,g=${grille}]
+
+exit
+eod
+    
+###################################################################
+# 3 - Launch interpolation with ferret
+###################################################################
+    ferret -nojnl <<eod >  /dev/null
+go ${varin}.jnl
+quit
+eod
+
+#}
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/interp_from_TS.ksh
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/interp_from_TS.ksh	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/interp_from_TS.ksh	(revision 1738)
@@ -0,0 +1,170 @@
+#!/bin/ksh
+set -xv
+set verbose echo
+#
+# This script will interpolate horizontaly time series from IPSLCM5 model 
+# output to the desired model grid 
+#
+# Contact Josefine Ghattas
+###############################################################################
+###############################################################################
+###############################################################################
+# BEGIN USER DEFINITION
+###############################################################################
+
+# gridfile : file grilles_gcm.nc containing longitudes and latitudes of the destination grid
+gridfile=/work/cont003/p86ghatt/GUIDAGE/INTERP_FROM_TS/LMD96x95_grilles_gcm.nc
+
+
+# indir : path for directory containing the files to interpolate
+indir=/dmnfs/cont003/p86denv/IGCM_OUT/IPSLCM5A/PROD/rcp45/v3.rcp45.1/ATM/Analyse/TS_HF
+
+# varlist : list of variables to interpolate
+#varlist="vitu vitv temp hur"
+varlist="vitu"
+
+# filelist_var : list of files to interpolate for variable var
+filelist_vitu="v3.rcp45.1_20060101_20151231_HF_vitu.nc  v3.rcp45.1_20360101_20451231_HF_vitu.nc  v3.rcp45.1_20660101_20751231_HF_vitu.nc v3.rcp45.1_20160101_20251231_HF_vitu.nc  v3.rcp45.1_20460101_20551231_HF_vitu.nc  v3.rcp45.1_20760101_20851231_HF_vitu.nc v3.rcp45.1_20260101_20351231_HF_vitu.nc  v3.rcp45.1_20560101_20651231_HF_vitu.nc  v3.rcp45.1_20860101_20951231_HF_vitu.nc"
+
+filelist_vitv="v3.rcp45.1_20060101_20151231_HF_vitv.nc  v3.rcp45.1_20360101_20451231_HF_vitv.nc  v3.rcp45.1_20660101_20751231_HF_vitv.nc v3.rcp45.1_20160101_20251231_HF_vitv.nc  v3.rcp45.1_20460101_20551231_HF_vitv.nc  v3.rcp45.1_20760101_20851231_HF_vitv.nc v3.rcp45.1_20260101_20351231_HF_vitv.nc  v3.rcp45.1_20560101_20651231_HF_vitv.nc  v3.rcp45.1_20860101_20951231_HF_vitv.nc"
+
+filelist_temp="v3.rcp45.1_20060101_20151231_HF_temp.nc  v3.rcp45.1_20360101_20451231_HF_temp.nc  v3.rcp45.1_20660101_20751231_HF_temp.nc v3.rcp45.1_20160101_20251231_HF_temp.nc  v3.rcp45.1_20460101_20551231_HF_temp.nc  v3.rcp45.1_20760101_20851231_HF_temp.nc v3.rcp45.1_20260101_20351231_HF_temp.nc  v3.rcp45.1_20560101_20651231_HF_temp.nc  v3.rcp45.1_20860101_20951231_HF_temp.nc"
+
+filelist_hur="v3.rcp45.1_20060101_20151231_HF_hur.nc  v3.rcp45.1_20360101_20451231_HF_hur.nc  v3.rcp45.1_20660101_20751231_HF_hur.nc v3.rcp45.1_20160101_20251231_HF_hur.nc  v3.rcp45.1_20460101_20551231_HF_hur.nc  v3.rcp45.1_20760101_20851231_HF_hur.nc v3.rcp45.1_20260101_20351231_HF_hur.nc  v3.rcp45.1_20560101_20651231_HF_hur.nc  v3.rcp45.1_20860101_20951231_HF_hur.nc"
+
+# Define variables for output
+# Complet file name for output : $outdir/$year/$outprefix_$var_$year_$month.nc
+# outdir : directory to store output files
+outdir=$DMFDIR/IGCM_OUT/LMDZOR/NUDGEFILES/LMD9695/RCP45
+
+# outprefix : prefix to output files
+outprefix=v3.rcp45.1_HF_
+# outsuffix : suffix to output files
+outsuffix=.nc
+
+# Time variables
+# first_year : first year to interpolate
+first_year=2006
+# last_year : last year to interpolate
+last_year=2095
+
+# rundir : temporary run directory
+rundir=$CSCRATCHDIR/RUNDIR/interprun
+###############################################################################
+# END USER DEFINITION
+###############################################################################
+###############################################################################
+###############################################################################
+# Create and move to run directory
+if [ ! -d ${rundir} ] ; then
+    mkdir -p ${rundir}
+fi
+cp -f era2gcm.ksh $rundir/.
+cd $rundir
+
+# Create output directory if it does not exist
+if [ ! -d $outdir ] ; then
+    mkdir -p $outdir
+fi
+
+
+###############################################################################
+# Loop over variable list
+#   * Loop over corresponding file list
+#     1 - Interpolate the whole time-serie to destination grid
+#     2 - Devide files into monthly files
+#   * Loop over all years for each variable
+#     3 - Add a time step last in each file coming from next month 
+###############################################################################
+for var in $varlist ; do
+    case $var in
+	vitu)
+	    filelist=$filelist_vitu ;;
+	vitv)
+	    filelist=$filelist_vitv ;;
+        temp)
+            filelist=$fileslist_temp ;;
+        hur)
+            filelist=$filelist_hur ;;
+    esac
+    
+    echo $var : $filelist; echo
+ 
+    for file in $filelist ; do
+# 1 - Interpolate time serie to destiation grid
+	filein=$indir/$file
+	file2=$outdir/$file
+	cp -f $filein file1.nc
+	echo Interpolate file $filein and save in file $file2
+	./era2gcm.ksh $gridfile file1.nc file2.nc $var
+	
+# 2 - Devide files into monthly files
+	cdo splityear file2.nc outyear_
+    done
+    ls -lrt
+
+    year=$first_year
+    while [[ ${year} -le ${last_year} ]] ; do
+	# split in month
+	cdo splitmon outyear_$year.nc outmonth_
+
+	# move monthly files to output directory
+	mkdir -p $outdir/$year
+	month=01
+	for month in 01 02 03 04 05 06 07 08 09 10 11 12 ; do
+	    mv outmonth_$month.nc $outdir/$year/${outprefix}${var}_${year}_${month}${outsuffix}
+	done
+	
+	# update year
+	year=`expr $year + 1`
+    done
+    
+# Loop over all years for one variables
+# 3 - Add first time step last in file
+    echo Before add step
+    ls -lrt $outdir
+    ls -lrt $outdir/*
+
+    year=$first_year
+    month=01
+    while [[ ${year} -le ${last_year} ]] ; do
+	case $month in
+	    12)
+		nextyear=`expr $year + 1`
+		nextmonth=01
+		;;
+	    *)
+		nextyear=$year
+		nextmonth=`expr $month + 1`
+		if [[ $nextmonth -le 9 ]] ; then
+		    nextmonth=0$nextmonth
+		fi
+		;;
+	esac
+	
+	file=$outdir/$year/${outprefix}${var}_${year}_${month}${outsuffix}
+	nextfile=$outdir/$year/${outprefix}${var}_${nextyear}_${nextmonth}${outsuffix}
+	if [ -f ${nextfile} ] ; then
+	    echo First time step from file ${nextfile} will be added last in file ${file}
+	    cp -f $file file.nc
+	    cp -f $nextfile next.nc
+	    ncks -O -d time,0,0 next.nc tmp.nc
+	    ncrcat -O file.nc tmp.nc tmptot.nc
+	    mv tmptot.nc $file
+	else
+	    echo Next file ${nextfile} does not exist. Nothing will be done for ${file}
+	fi
+	
+	# Update month and year
+	year=$nextyear
+	month=$nextmonth
+    
+    done
+	
+# End of loop variable list
+done
+
+### Copy gridfile in output dir
+cp -f $gridfile $outdir/grilles_gcm.nc
+	    
+exit
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/interp_from_era.ksh
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/interp_from_era.ksh	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/SCRIPT/interp_from_era.ksh	(revision 1738)
@@ -0,0 +1,147 @@
+#!/bin/ksh
+set verbose echo
+#
+# This script will interpolate horizontaly ERA-I
+# output to the desired model grid 
+#
+# Contact Josefine Ghattas
+###############################################################################
+###############################################################################
+###############################################################################
+# BEGIN USER DEFINITION
+###############################################################################
+
+# gridfile : file grilles_gcm.nc containing longitudes and latitudes of the destination grid
+gridfile=/work/cont003/p86ghatt/GUIDAGE/INTERP_FROM_TS/LMD96x95_grilles_gcm.nc
+
+# indir : path for directory containing the files to interpolate
+#indir=/u/rech/psl/rpsl376/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL
+indir=/dmnfs/cont003/p24data/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL
+# insuffix : suffix in input files
+insuffix=.aphei.GLOBAL_075.nc
+
+# varlist : list of variables to interpolate
+varlist="u v ta r"
+
+# Variables for output :
+# outdir : directory to store output files
+outdir=/dmnfs/cont003/p86ghatt/IGCM_OUT/LMDZOR/NUDGEFILES/LMD9695/ERAI
+# outsuffix : suffix to output files
+outsuffix=_erai.nc
+
+# Time variables :
+# first_year : first year to interpolate
+first_year=2000
+# last_year : last year to interpolate
+last_year=2000
+
+# rundir : temporary run directory
+rundir=$CSCRATCHDIR/RUNDIR/interprun
+###############################################################################
+# END USER DEFINITION
+###############################################################################
+###############################################################################
+###############################################################################
+if [ ! -d ${rundir} ] ; then
+    mkdir -p ${rundir}
+fi
+cp -f era2gcm.ksh $rundir/.
+cd $rundir
+
+###############################################################################
+# For all variables do :
+#   1) Interpolate all month and years for each variable
+#   2) Add a time step last in each file coming from next month 
+###############################################################################
+# Loop over all variables
+for var in $varlist ; do
+    year=$first_year
+    month=01
+
+    # 1 - Loop over all files and interpolate the variable
+    while [[ ${year} -le ${last_year} ]] ; do
+        # Create output directory if it does not exist
+	if [ ! -d $outdir/${year} ] ; then
+	    mkdir -p $outdir/${year}
+	fi
+
+	# Define source and target file
+	filein=${indir}/${year}/${var}.${year}${month}${insuffix}
+	fileout=${outdir}/${year}/${var}_${year}${month}${outsuffix}
+
+	# Interpolate filein and save in fileout
+	echo Interpolate file $filein and save in file $fileout
+	./era2gcm.ksh $gridfile $filein tmp.nc $var
+	mv tmp.nc $fileout
+	echo
+
+
+	# Find next month and year
+	case $month in
+	    12)
+		nextyear=`expr $year + 1`
+		nextmonth=01
+		;;
+	    *)
+		nextyear=$year
+		nextmonth=`expr $month + 1`
+		if [[ $nextmonth -le 9 ]] ; then
+		    nextmonth=0$nextmonth
+		fi
+		;;
+	esac
+	# Update month and year
+	year=$nextyear
+	month=$nextmonth
+    done
+
+
+    # 2 - Loop over all files and add a time step : 
+    #     take the first time step from next coming month and add last in current month
+    echo Before add step
+
+    year=$first_year
+    month=01
+    while [[ ${year} -le ${last_year} ]] ; do
+	# Find next month and year
+	case $month in
+	    12)
+		nextyear=`expr $year + 1`
+		nextmonth=01
+		;;
+	    *)
+		nextyear=$year
+		nextmonth=`expr $month + 1`
+		if [[ $nextmonth -le 9 ]] ; then
+		    nextmonth=0$nextmonth
+		fi
+		;;
+	esac
+	
+	file=${outdir}/${year}/${var}_${year}${month}${outsuffix}
+	nextfile=${outdir}/${year}/${var}_${nextyear}${nextmonth}${outsuffix}
+	if [ -f ${nextfile} ] ; then
+	    echo First time step from file ${nextfile} will be added last in file ${file}
+	    cp -f $file file.nc
+	    cp -f $nextfile next.nc
+	    ncks -O -d TIME,0,0 next.nc tmp.nc
+	    ncrcat -O file.nc tmp.nc tmptot.nc
+	    mv tmptot.nc $file
+	else
+	    echo Next file ${nextfile} does not exist. Nothing will be done for ${file}
+	fi
+
+	# Update month and year
+	year=$nextyear
+	month=$nextmonth
+    done
+
+
+    
+# End of loop variable list
+done
+
+### Copy gridfile in output dir
+cp -f $gridfile $outdir/grilles_gcm.nc
+    
+exit
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/config.card
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/config.card	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/CREATE/config.card	(revision 1738)
@@ -0,0 +1,93 @@
+# $Id$
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=CE0L
+LongName="create_etat0_limit for interannuel run"
+TagName=LMDZORINCA
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Begin and end of job 
+#D-- "YYYY-MM-DD"
+DateBegin=1979-01-01
+DateEnd=1979-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, CREATE_E0)
+CHM=(inca, INCA3)
+
+#========================================================================
+#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)
+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=
+#D- Define restart simulation name
+RestartJobName=
+#D- Path Server Group Login
+RestartPath=
+# Attention login depend de la machine
+
+#========================================================================
+#D-- ATM -
+[ATM]
+#
+WriteFrequency="1Y"
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=
+# Define restart simulation name
+RestartJobName=
+RestartPath=
+# Old component name for restart (if empty, use new name)
+OldName=
+
+#========================================================================
+#D-- CHM -
+[CHM]
+#
+WriteFrequency="1Y"
+# If config_Restarts_OverRule == 'n' all params are read
+Restart= n
+# Last day of the experience used as restart
+RestartDate=
+# Define restart simulation name
+RestartJobName=
+RestartPath=
+# Old component name for restart (if empty, use new name)
+OldName=
+
+
+#========================================================================
+#D-- NO Post -
+[Post]
+RebuildFrequency=NONE
+TimeSeriesFrequency=NONE
+SeasonalFrequency=NONE
+#========================================================================
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_AER
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_AER	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_AER	(revision 1738)
@@ -0,0 +1,51 @@
+#- $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/ECMWF96x95/AN${year}/u10m_ecmwf_${year}${month}.nc , u10m.nc  )\
+	    (/dmnfs/cont003/p24data/ECMWF96x95/AN${year}/v10m_ecmwf_${year}${month}.nc , v10m.nc  )
+
+ListNonDel= (${R_INIT}/CHM/LMDZORINCA/AER/INCA${RESOL_CHM}/oxydants_AR5WERF_AN2006.nc				,   oxydants.nc     )\
+	    (${R_INIT}/CHM/LMDZORINCA/CH4/INCA${RESOL_CHM}/npp.nc                  				,   npp.nc          )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA${RESOL_CHM}/sflx_lmdz_phy_AR5werf_2006_BioLathpres.nc	,   sflx.nc	    )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA${RESOL_CHM}/aircraft_ar5_2000_lmdz_phy.nc  			,   aircraft.nc     )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC/phototable.dat_V4.1	                 				,   phototable.dat )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC/INCA${RESOL_CHM}/landuse.nc              				,   landuse.nc      )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/rhvEC11.May2009v2.txt           					,   rhvEC.txt	    )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/wthEC11.09092002.txt	                 				,   wthEC.txt	    )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/clyEC11.09092002.txt	                 				,   clyEC.txt	    )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA${RESOL_CHM}/wthlmdz.txt	 					,   wth.dat	    )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA${RESOL_CHM}/clylmdz.txt  					,   cly.dat	    )\
+ 	    (${R_INIT}/CHM/LMDZORINCA/AER/INCA${RESOL_CHM}/rhvlmdz.txt  					,   rhv.dat         )
+
+[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=  
+
+[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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_AER_L39
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_AER_L39	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_AER_L39	(revision 1738)
@@ -0,0 +1,51 @@
+#- $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/ECMWF96x95/AN${year}/u10m_ecmwf_${year}${month}.nc , u10m.nc  )\
+	    (/dmnfs/cont003/p24data/ECMWF96x95/AN${year}/v10m_ecmwf_${year}${month}.nc , v10m.nc  )
+
+ListNonDel= (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/oxydants_LOINA39_02_AN2006.nc	,   oxydants.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/npp.nc     			,   npp.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/aircraft.nc			,   aircraft.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/sflx.nc      		,   sflx.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/phototable.dat      	,   phototable.dat	)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/landuse.nc      		,   landuse.nc		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/rhvEC11.May2009v2.txt		,   rhvEC.txt		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/wthEC11.09092002.txt		,   wthEC.txt		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/clyEC11.09092002.txt		,   clyEC.txt		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/wthlmdz.txt			,   wth.dat		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/clylmdz.txt			,   cly.dat		)\
+ 	    (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/rhvlmdz.txt			,   rhv.dat		)
+
+[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=  
+
+[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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_GES
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_GES	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_GES	(revision 1738)
@@ -0,0 +1,45 @@
+#- $Id: inca.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+[BoundaryFiles]
+List=	    ()
+
+ListNonDel= (${R_INIT}/CHM/LMDZORINCA/AER/INCA${RESOL_CHM}/oxydants_AR5WERF_AN2006.nc				,   oxydants.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/CH4/INCA${RESOL_CHM}/npp.nc						,   npp.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA${RESOL_CHM}/sflx_lmdz_phy_AR5werf_2006_BioLathpres.nc	,   sflx.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA${RESOL_CHM}/aircraft_ar5_2000_lmdz_phy.nc			,   aircraft.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/CH4/phototable.dat_V4.2							,   phototable.dat	)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC/INCA${RESOL_CHM}/landuse.nc						,   landuse.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/GES/INCA${RESOL_CHM}/FOS_LMDZ9696.nc								,   andres.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/GES/INCA${RESOL_CHM}/BIO_m_LMDZ9696.nc								,   casa_m.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/GES/INCA${RESOL_CHM}/OCE_LMDZ9696.nc								,   taka.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/GES/INCA${RESOL_CHM}/BIO_3h_LMDZ9696.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=  
+
+[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 = (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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_GES_L39
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_GES_L39	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_GES_L39	(revision 1738)
@@ -0,0 +1,46 @@
+#- $Id: inca.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+[BoundaryFiles]
+List=
+
+ListNonDel= (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/oxydants_LOINA39_02_AN2006.nc		,   oxydants.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/npp.nc				,   npp.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/aircraft.nc				,   aircraft.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/sflx.nc				,   sflx.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/phototable.dat			,   phototable.dat	)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/landuse.nc				,   landuse.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/GES/INCA9695_L39/FOS_LMDZ9696.nc				,   andres.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/GES/INCA9695_L39/BIO_m_LMDZ9696.nc			,   casa_m.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/GES/INCA9695_L39/OCE_LMDZ9696.nc				,   taka.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/GES/INCA9695_L39/BIO_3h_LMDZ9696.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=  
+
+[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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC	(revision 1738)
@@ -0,0 +1,41 @@
+#- $Id: inca.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+[BoundaryFiles]
+List=	    ()
+
+ListNonDel= (${R_INIT}/CHM/LMDZORINCA/CH4/INCA${RESOL_CHM}/o3clim.nc             				,   o3clim.nc       )\
+	    (${R_INIT}/CHM/LMDZORINCA/CH4/INCA${RESOL_CHM}/so4.nc				         	,   so4.nc	    )\
+	    (${R_INIT}/CHM/LMDZORINCA/CH4/INCA${RESOL_CHM}/npp.nc                  				,   npp.nc          )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA${RESOL_CHM}/sflx_lmdz_phy_AR5werf_2006_BioLathpres.nc	,   sflx.nc	    )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA${RESOL_CHM}/aircraft_ar5_2000_lmdz_phy.nc  			,   aircraft.nc     )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC/phototable.dat_V4.1	                 				,   phototable.dat  )\
+	    (${R_INIT}/CHM/LMDZORINCA/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=  
+
+[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 = (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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_AER
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_AER	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_AER	(revision 1738)
@@ -0,0 +1,52 @@
+#- $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/ECMWF96x95/AN${year}/u10m_ecmwf_${year}${month}.nc , u10m.nc  )\
+	    (/dmnfs/cont003/p24data/ECMWF96x95/AN${year}/v10m_ecmwf_${year}${month}.nc , v10m.nc  )
+
+ListNonDel= (${R_INIT}/CHM/LMDZORINCA/CH4/INCA${RESOL_CHM}/o3clim.nc             				,   o3clim.nc      )\
+	    (${R_INIT}/CHM/LMDZORINCA/CH4/INCA${RESOL_CHM}/so4.nc				         	,   so4.nc	    )\
+	    (${R_INIT}/CHM/LMDZORINCA/CH4/INCA${RESOL_CHM}/npp.nc                  				,   npp.nc          )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA${RESOL_CHM}/sflx_lmdz_phy_AR5werf_2006_BioLathpres.nc	,   sflx.nc	    )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA${RESOL_CHM}/aircraft_ar5_2000_lmdz_phy.nc  			,   aircraft.nc     )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC/phototable.dat_V4.1	                 				,   phototable.dat )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC/INCA${RESOL_CHM}/landuse.nc              				,   landuse.nc      )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/rhvEC11.May2009v2.txt           					,   rhvEC.txt	    )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/wthEC11.09092002.txt	                 				,   wthEC.txt	    )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/clyEC11.09092002.txt	                 				,   clyEC.txt	    )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA${RESOL_CHM}/wthlmdz.txt	 					,   wth.dat	    )\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA${RESOL_CHM}/clylmdz.txt  					,   cly.dat	    )\
+ 	    (${R_INIT}/CHM/LMDZORINCA/AER/INCA${RESOL_CHM}/rhvlmdz.txt  					,   rhv.dat         )
+
+[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=  
+
+[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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_AER_L39
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_AER_L39	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_AER_L39	(revision 1738)
@@ -0,0 +1,52 @@
+#- $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/ECMWF96x95/AN${year}/u10m_ecmwf_${year}${month}.nc , u10m.nc  )\
+	    (/dmnfs/cont003/p24data/ECMWF96x95/AN${year}/v10m_ecmwf_${year}${month}.nc , v10m.nc  )
+
+ListNonDel= (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/o3clim.nc      	,   o3clim.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/npp.nc      	,   npp.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/aircraft.nc      	,   aircraft.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/so4.nc      	,   so4.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/sflx.nc      	,   sflx.nc		)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/phototable.dat      ,   phototable.dat	)\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/landuse.nc      	,   landuse.nc		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/rhvEC11.May2009v2.txt	,   rhvEC.txt		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/wthEC11.09092002.txt	,   wthEC.txt		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/clyEC11.09092002.txt	,   clyEC.txt		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/wthlmdz.txt		,   wth.dat		)\
+            (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/clylmdz.txt		,   cly.dat		)\
+ 	    (${R_INIT}/CHM/LMDZORINCA/AER/INCA9695_L39/rhvlmdz.txt		,   rhv.dat		)
+
+[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=  
+
+[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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_L39
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_L39	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.card_NMHC_L39	(revision 1738)
@@ -0,0 +1,45 @@
+#- $Id: inca.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+
+[InitialStateFiles]
+List= ()
+	
+[BoundaryFiles]
+List=
+
+
+ListNonDel= (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/o3clim.nc        	,   o3clim.nc       )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/so4.nc		,   so4.nc	    )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/npp.nc           	,   npp.nc          )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/sflx.nc		,   sflx.nc	    )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/aircraft.nc 	,   aircraft.nc     )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/phototable.dat  	,   phototable.dat  )\
+	    (${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/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=  
+
+[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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.driver
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.driver	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/inca.driver	(revision 1738)
@@ -0,0 +1,92 @@
+#- $Id: inca.driver 119 2007-07-31 13:56:33Z acosce $
+#!/bin/ksh
+
+#D- Driver du script pour INCA
+#-----------------------------------------------------------------
+function INCA_sed_replace
+{
+    IGCM_debug_PushStack "INCA_sed_replace"
+
+    sed -e "s/${2}/${3}/" ${1} > ${1}.tmp
+    RET=$?
+    echo "INCA_sed : ${1} ${2} ${3}"
+    \mv ${1}.tmp ${1}
+
+    IGCM_debug_PopStack "INCA_sed_replace"
+    return $RET
+}
+function INCA_sed_complete
+{
+    IGCM_debug_PushStack "INCA_sed_complete"
+
+    sed -e "s/^${2}\ *=.*/${2}= ${3}/" ${1} > ${1}.tmp
+    RET=$?
+    echo "INCA_sed : ${1} ${2} ${3}"
+    \mv ${1}.tmp ${1}
+
+    IGCM_debug_PopStack "INCA_sed_complete"
+    return $RET
+}
+
+#-----------------------------------------------------------------
+function CHM_Initialize
+{
+    IGCM_debug_PushStack "CHM_Initialize"
+
+    RESOL_CHM=$( echo $RESOL_ATM | awk "-FD" '{print $2}' )
+
+    [ -f ${SUBMIT_DIR}/../.resol ] && eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x19
+    RESOL_CHM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' | awk "-F-" '{print $1}')
+    echo 'resol chimie', $RESOL_CHM
+
+    [ -f ${SUBMIT_DIR}/../.chimie ] && eval $(grep CHEM ${SUBMIT_DIR}/../.chimie) 
+    INCA_sed_replace ${SUBMIT_DIR}/config.card \${CHEM} ${CHEM}
+    echo "INCA Tag : LMDZ-INCA" ${config_UserChoices_TagName}
+
+    IGCM_debug_PopStack "CHM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function CHM_Update
+{
+    IGCM_debug_PushStack "CHM_Update"
+
+    case ${CHEM} in 
+	AER)
+	INCA_sed_complete physiq.def config_inca aero
+	;;
+	NMHC) 
+	INCA_sed_complete physiq.def config_inca chem
+	;;
+	NMHC_AER)
+	INCA_sed_complete physiq.def config_inca aero
+	;;
+	CH4)
+	INCA_sed_complete physiq.def config_inca chem
+	;;
+	CH4_AER)
+	INCA_sed_complete physiq.def config_inca aero
+	;;
+	GES)
+	INCA_sed_complete physiq.def config_inca chem
+	esac
+
+    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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.card_L19
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.card_L19	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.card_L19	(revision 1738)
@@ -0,0 +1,65 @@
+#- $Id: lmdz.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+# Physics package to use :
+# AP for standard/old physics(defalut), can be used with LMDZ4_AR5 or LMDZ5/trunk sources with 19 or 39 levels
+# NPv3.1 for new physics, to be used with LMDZ5/trunk revision 1554 or later with 39 levels
+LMDZ_Physics=AP
+#LMDZ_Physics=NPv3.1
+
+
+
+LMDZ_NbPeriod_adjust=3
+#LMDZ_NbPeriod_adjust=0
+# To force usage of this Bands file, put LMDZ_NbPeriod_adjust=0 and replace XXXXXXX by Name/Resolution/NB Prco Job Restart
+# LMDZ_Bands_file_name=${ARCHIVE}/IGCM_OUT/LMDZORINCA/XXXXXXX/ATM/Debug/XXXXXXX_Bands_XXXXXXX_XXXprc.dat_3
+# LMDZ_Bands_file_name=/dmnfs11/cont003/p86cozic/IGCM_OUT/LMDZORINCA/AER/LOI3A_BR/ATM/Debug/LOI3A_BR_Bands_96x95x19_32prc.dat_1
+# LMDZ_Bands_file_name=/dmnfs11/cont003/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER/LOI3NA_CB/ATM/Debug/LOI3NA_CB_Bands_96x95x19_32prc.dat_1
+
+# Other specific parameters to this experience 
+aer_type=scenario
+pmagic=0.0
+
+# Nugding : ok_guide=y/n
+ok_guide=n
+
+
+[InitialStateFiles]
+# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job
+List= 	(/dmnfs/cont003/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER/LOI_AR5WERF_E0_2006/ATM/Restart/LOI_AR5WERF_E0_2006_20051231_restart.nc, start.nc),\
+	(/dmnfs/cont003/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER/LOI_AR5WERF_E0_2006/ATM/Restart/LOI_AR5WERF_E0_2006_20051231_restartphy.nc, startphy.nc)
+
+[BoundaryFiles]
+List=   (${R_BC}/CHM/LMDZORINCA/${RESOL_ATM}_leap/EL-${RESOL_ATM}_${year}1231_limit.nc 	, limit.nc	),\
+	(${R_BC}/ATM/LMDZORINCA/${RESOL_ATM}/NUDGE_FILES/u_ecmwf_${year}${month}.nc	, u.nc		),\
+	(${R_BC}/ATM/LMDZORINCA/${RESOL_ATM}/NUDGE_FILES/v_ecmwf_${year}${month}.nc	, v.nc		)
+
+ListNonDel= ()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}_${LMDZ_Physics}	, physiq.def		), \
+	(${SUBMIT_DIR}/PARAM/gcm.def_L${RESOL_ATM_Z}_${LMDZ_Physics}	, gcm.def		), \
+	(${SUBMIT_DIR}/PARAM/run.def					, .			), \
+        (${SUBMIT_DIR}/PARAM/guide.def					, .			)
+ 
+
+[RestartFiles]
+List=	(restart.nc, restart.nc, start.nc), \
+	(restartphy.nc, restartphy.nc,  startphy.nc) 
+
+[OutputText]
+List=   (physiq.def, gcm.def, run.def, guide.def, traceur.def)
+
+
+[OutputFiles]
+List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,   Post_1M_histmth   )	
+
+[Post_1M_histmth]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds, 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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.card_L39
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.card_L39	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.card_L39	(revision 1738)
@@ -0,0 +1,64 @@
+#- $Id: lmdz.card 70 2007-06-21 14:20:42Z acosce $
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+# Physics package to use :
+# AP for standard/old physics(defalut), can be used with LMDZ4_AR5 or LMDZ5/trunk sources with 19 or 39 levels
+# NPv3.1 for new physics, to be used with LMDZ5/trunk revision 1554 or later with 39 levels
+LMDZ_Physics=AP
+#LMDZ_Physics=NPv3.1
+
+
+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=/ccc/store/cont003/dsm/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER/L4OI3NA_39_Bands/ATM/Debug/L4OI3NA_39_Bands_Bands_96x95x39_32prc.dat_1
+
+# Other specific parameters to this experience 
+aer_type=scenario
+pmagic=0.0
+
+# Nugding : ok_guide=y/n
+ok_guide=n
+
+
+[InitialStateFiles]
+# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job
+List= 	(${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/start.nc	, start.nc	) \
+	(${R_INIT}/CHM/LMDZORINCA/NMHC_AER/INCA9695_L39/startphy.nc	, startphy.nc	)
+
+
+[BoundaryFiles]
+List=   (${R_BC}/CHM/LMDZORINCA/LMD9695_leap/EL-LMD9695_${year}1231_limit.nc , limit.nc)
+
+#	(${R_BC}/ATM/LMDZORINCA/LMD9695/NUDGE_FILES/u_ecmwf_${year}${month}.nc	, u.nc		),\
+#	(${R_BC}/ATM/LMDZORINCA/LMD9695/NUDGE_FILES/v_ecmwf_${year}${month}.nc	, v.nc		)
+ 
+ListNonDel= ()
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}_${LMDZ_Physics}	, physiq.def	), \
+	(${SUBMIT_DIR}/PARAM/gcm.def_L${RESOL_ATM_Z}_${LMDZ_Physics}	, gcm.def	), \
+	(${SUBMIT_DIR}/PARAM/run.def					, .		), \
+	(${SUBMIT_DIR}/PARAM/traceur.def				, .		), \
+	(${SUBMIT_DIR}/PARAM/guide.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)
+
+
+[OutputFiles]
+List=	(histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,   Post_1M_histmth   )	
+
+[Post_1M_histmth]
+Patches= ()
+GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds, 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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.driver
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.driver	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/lmdz.driver	(revision 1738)
@@ -0,0 +1,166 @@
+#- $Id: lmdz.driver 841 2009-12-08 09:48:13Z acosce $
+#!/bin/ksh
+
+#-----------------------------------------------------------------
+function LMDZ_sed
+{
+    IGCM_debug_PushStack "LMDZ_sed"
+
+    sed -e "s/^${2}\ *=.*/${2}= ${3}/" ${1} > ${1}.tmp
+    RET=$?
+    echo "LMDZ_sed : ${1} ${2} ${3}"
+    \mv ${1}.tmp ${1}
+
+    IGCM_debug_PopStack "LMDZ_sed"
+    return $RET
+}
+
+#-----------------------------------------------------------------
+
+function ATM_Initialize
+{
+    IGCM_debug_PushStack "ATM_Initialize"
+
+    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' )
+
+    [ -f ${SUBMIT_DIR}/../.resol ] && eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x19
+
+    RESOL_ATM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' | awk "-F-" '{print $1}')
+
+
+    ##-- 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
+
+
+   ##- LMDZ physics version
+   ##  Read LMDZ_Physics option in lmdz.card, if not present take default value AP (old physics)
+   if [ ! X${lmdz_UserChoices_LMDZ_Physics} = X ] ; then
+       LMDZ_Physics=${lmdz_UserChoices_LMDZ_Physics}
+   else
+       LMDZ_Physics=AP
+   fi
+   echo LMDZ physics version : ${LMDZ_Physics}
+
+
+    ##- 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
+
+
+    ##-- Output frequency purpose ....
+    ##--  Initialisation  ....
+    ok_mensuel=y
+    ok_journe=n
+    ok_hf=n
+    ok_hf3h=n
+    ok_hf3hm=n
+    ok_stn=n
+
+
+
+    ## 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
+
+
+    IGCM_debug_PopStack "ATM_Initialize"
+}
+
+#-----------------------------------------------------------------
+function ATM_Update
+{
+    IGCM_debug_PushStack "ATM_Update"
+
+    ecritphy=${PeriodLengthInDays}
+    LMDZ_ecrit_ISCCP=${PeriodLengthInDays}.
+    LMDZ_periodav=${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
+
+
+    ## algo : For CumulPeriod=1 ; 
+    ##        IF NbPeriod_Adjust = 0, Get Bands file directly from server, store it with _0 suffix
+    ##        ELSE LMDZ_adjust=y ; save Bands file on file server ;
+    ##        For CumulPeriod=2 to LMDZ_NbPeriod_adjust ; Get Bands file from server (CumulPeriod-1) ; LMDZ_adjust=y ; save Bands file on file server ;
+    ##        For CumulPeriod > LMDZ_NbPeriod_adjust ; LMDZ_adjust=n ; Get Bands file from server (CumulPeriod=LMDZ_NbPeriod_adjust 0 to force it) ; Use Bands ;
+
+    if [ ${CumulPeriod} -eq 1 ] ; then
+        if [ ${LMDZ_NbPeriod_adjust} -eq 0 ] ; then
+            if [ ! X${LMDZ_Bands_file_name} = X] ; then
+               IGCM_sys_Get ${LMDZ_Bands_file_name} Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat
+               IGCM_sys_Put_Out Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${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_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_0"
+      	       IGCM_sys_Chmod u+w Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat
+            fi
+        fi
+    fi
+
+    LMDZ_adjust=n
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && LMDZ_adjust=y
+
+    RefPeriod=${LMDZ_NbPeriod_adjust} #0 to use a Bands file from Restart
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && (( RefPeriod = ${CumulPeriod} - 1 )) 
+
+    [ ${CumulPeriod} -gt 1 ] && ( IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat; IGCM_sys_Chmod u+w Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat )
+
+    ## physiq.def parameters modified from initial physiq_L${RESOL_ATM_Z}
+    #  columns in phys_output_filekeys refer to the filenames given in phys_out_filenames in output.def :  
+    #                                            histmth       histday        histhf    histfh3h   histhf3hm   histstn
+    LMDZ_sed physiq.def phys_out_filekeys       "${ok_mensuel} ${ok_journe}   ${ok_hf}  ${ok_hf3h} ${ok_hf3hm} ${ok_stn}"
+    LMDZ_sed physiq.def ecrit_ISCCP ${LMDZ_ecrit_ISCCP}
+
+    if [ ! X${lmdz_UserChoices_ok_guide} = X ] ; then
+	LMDZ_sed guide.def ok_guide  ${lmdz_UserChoices_ok_guide}
+    fi
+
+    ## run.def parameters
+    LMDZ_sed run.def dayref   ${InitDay}
+    LMDZ_sed run.def anneeref ${InitYear}
+    LMDZ_sed run.def calend   ${CalendarTypeForLmdz}
+    LMDZ_sed run.def nday     ${PeriodLengthInDays}
+    LMDZ_sed run.def raz_date ${RAZ_DATE}
+    LMDZ_sed run.def periodav ${LMDZ_periodav}
+    LMDZ_sed run.def adjust   ${LMDZ_adjust}
+
+    ## Other specific parameters to this experience 
+    if [ ! X${lmdz_UserChoices_aer_type} = X ] ; then
+        LMDZ_sed physiq.def aer_type ${lmdz_UserChoices_aer_type}
+    fi
+    if [ ! X${lmdz_UserChoices_pmagic} = X ] ; then
+        LMDZ_sed physiq.def pmagic   ${lmdz_UserChoices_pmagic}
+    fi
+
+    IGCM_debug_PopStack "ATM_Update"
+}
+
+#-----------------------------------
+function ATM_Finalize
+{
+    IGCM_debug_PushStack "ATM_Finalize"
+
+    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && IGCM_sys_Put_Out Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_${CumulPeriod}
+    echo FINALIZE ATM !
+
+    IGCM_debug_PopStack "ATM_Finalize"
+}
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/orchidee.card
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/orchidee.card	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/orchidee.card	(revision 1738)
@@ -0,0 +1,63 @@
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+NEWHYDROL=n
+sechiba_LEVEL=11
+VEGET_UPDATE=1Y
+LAND_COVER_CHANGE=y
+
+[InitialStateFiles]
+List=	(${R_INIT}/SRF/LMDZORINCA/soils_param.nc,      . ), \
+	(${R_INIT}/SRF/LMDZORINCA/routing.nc,          . ), \
+	(${R_BC}/SRF/LMDZORINCA/PFTmap_1850to2005_AR5_LUHa.rc2/PFTmap_IPCC_${year}.nc, PFTmap.nc)
+## Pour les annees 2006-2100 suivant le RCP 
+#(${R_BC}/SRF/LMDZORINCA/PFTmap_2006to2100_AR5_RCP8.5_MESSAGE/PFTmap_IPCC_${year}.nc, PFTmap.nc)
+#(${R_BC}/SRF/LMDZORINCA/PFTmap_2006to2100_AR5_RCP4.5_MiniCAM//PFTmap_IPCC_${year}.nc, PFTmap.nc)
+#(${R_BC}/SRF/LMDZORINCA/PFTmap_2006to2100_AR5_RCP6.0_AIM//PFTmap_IPCC_${year}.nc, PFTmap.nc)
+#(${R_BC}/SRF/LMDZORINCA/PFTmap_2006to2100_AR5_RCP3-PD_IMAGE//PFTmap_IPCC_${year}.nc, PFTmap.nc)
+
+[BoundaryFiles]
+List= ()
+ListNonDel= ()
+
+[SmoothFiles]
+List=(${R_BC}/SRF/LMDZORINCA/PFTmap_1850to2005_AR5_LUHa.rc2/PFTmap_IPCC_${year_p1}.nc, PFTmap.nc, 12:12:)
+## Pour les annees 2006-2100 suivant le RCP 
+#(${R_BC}/SRF/LMDZORINCA/PFTmap_2006to2100_AR5_RCP8.5_MESSAGE/PFTmap_IPCC_${year_p1}.nc, PFTmap.nc, 12:12:)
+#(${R_BC}/SRF/LMDZORINCA/PFTmap_2006to2100_AR5_RCP4.5_MiniCAM//PFTmap_IPCC_${year_p1}.nc, PFTmap.nc, 12:12:)
+#(${R_BC}/SRF/LMDZORINCA/PFTmap_2006to2100_AR5_RCP6.0_AIM//PFTmap_IPCC_${year_p1}.nc, PFTmap.nc, 12:12:)
+#(${R_BC}/SRF/LMDZORINCA/PFTmap_2006to2100_AR5_RCP3-PD_IMAGE//PFTmap_IPCC_${year_p1}.nc, PFTmap.nc, 12:12:)
+
+[ParametersFiles]
+List=	(${SUBMIT_DIR}/PARAM/orchidee.def, .)
+
+[RestartFiles]
+# List restart that have to be saved/restored each loop (file out, saved, and in) :
+List=	(sechiba_rest_out.nc, sechiba_rest.nc, sechiba_rest_in.nc)
+
+[OutputText]
+List=   (out_orchidee, orchidee.def)
+
+[OutputFiles]
+List=	(sechiba_history.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_history.nc, Post_1M_sechiba_history), \
+        (sechiba_out_2.nc, ${R_OUT_SRF_O_H}/${PREFIX}_HF_sechiba_out_2.nc, Post_HF_sechiba_out_2), \
+        (orchidee_watchout.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_watchout.nc, NONE)
+
+[Post_1M_sechiba_history]
+Patches= ()
+GatherWithInternal= (lon, lat, veget, time_counter, time_counter_bnds, Areas, Contfrac)
+TimeSeriesVars2D= (nobiofrac, alb_nir, alb_vis, bqsb, evap, fluxlat, fluxsens, gqsb, netrad, qair, rain, runoff, snow, snownobio, snowf, subli, tair, temp_sol, tsol_max, tsol_min, drainage,mrsos, mrso, mrros, mrro, prveg, evspsblveg, evspsblsoi, tran, treeFrac, grassFrac, cropFrac, baresoilFrac, residualFrac)
+ChunckJob2D= NONE
+TimeSeriesVars3D= (lai, maxvegetfrac, vegetfrac, CO2FLUX, ptn, nee)
+ChunckJob3D= NONE
+Seasonal= ON
+
+[Post_HF_sechiba_out_2]
+Patches= ()
+GatherWithInternal= (lon, lat, veget, time_counter, time_counter_bnds, Areas, Contfrac)
+TimeSeriesVars2D= (mrsos, mrro)
+ChunckJob2D= 50Y
+TimeSeriesVars3D= (ptn)
+ChunckJob3D= 10Y
+Seasonal= OFF
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/orchidee.driver
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/orchidee.driver	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/orchidee.driver	(revision 1738)
@@ -0,0 +1,172 @@
+#!/bin/ksh
+
+#D- Driver du script pour ORCHIDEE
+
+function ORCHIDEE_sed
+{
+    IGCM_debug_PushStack "ORCHIDEE_sed"
+    
+    sed -e "s/^${1}\ *=.*/${1}= ${2}/" \
+	orchidee.def > orchidee.def.tmp
+    RET=$?
+    echo "ORCHIDEE_sed : ${1} ${2}"
+    \mv orchidee.def.tmp orchidee.def
+    
+    IGCM_debug_PopStack "ORCHIDEE_sed"
+    return $RET
+}
+
+#-----------------------------------------------------------------
+function SRF_Initialize
+{
+    IGCM_debug_PushStack "SRF_Initialize"
+
+    RESOL_SRF=ALL
+
+    for frequency in ${config_SRF_WriteFrequency} ; do
+        case ${frequency} in
+            HF|hf) SRF_ok_hf=y ;;
+        esac
+    done
+
+    IGCM_debug_PopStack "SRF_Initialize"
+}
+
+#-----------------------------------------------------------------
+function SRF_PeriodStart
+{
+
+    IGCM_debug_PushStack "SRF_PeriodStart"
+    if [ ${year} -gt 2004 ] ; then
+
+        PFTmap="${R_BC}/SRF/LMDZORINCA/ORCHIDEE_PFTmap_2006to2100_AR5_RCP8.5_MESSAGE.nc"
+        if [ ${year} -eq 2005 ] ; then
+            IGCM_sys_Rm -f PFTmap.nc
+        fi
+
+    else
+        PFTmap="${R_BC}/SRF/LMDZORINCA/PFTmap_1850to2005_AR5_LUHa.rc2.nc"
+    fi
+
+    BC_SRF_PFTmap=${PFTmap}
+
+    IGCM_debug_PopStack "SRF_PeriodStart"
+
+}
+
+#-----------------------------------------------------------------
+function SRF_Update
+{
+    IGCM_debug_PushStack "SRF_Update"
+
+    typeset SECHIBA_WRITE_STEP
+
+    SRF_WriteFrequency=$( echo ${config_SRF_WriteFrequency} | sed -e 's/\([0-9]*[yYmMdDs]\).*/\1/' ) 
+    case ${SRF_WriteFrequency} in
+	*Y|*y) 
+	    WriteInYears=$( echo ${SRF_WriteFrequency} | awk -F '[yY]' '{print $1}' )
+	    PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' )
+	    (( SECHIBA_WRITE_STEP = PeriodLengthInDays * WriteInYears / PeriodLengthInYears * 86400 )) ;;
+	1M)
+	    case ${config_UserChoices_PeriodLength} in
+	    *Y|*y)
+		SECHIBA_WRITE_STEP=-1.
+		;;
+	    *M|*m)
+		SECHIBA_WRITE_STEP=-1.
+		;;
+	    *)
+		(( SECHIBA_WRITE_STEP = $( IGCM_date_DaysInMonth $year $month ) * 86400 ))
+		;;
+	    esac
+	    ;;
+	*M|*m) 
+	    WriteInMonths=$( echo ${SRF_WriteFrequency} | awk -F '[mM]' '{print $1}' )
+	    case ${config_UserChoices_PeriodLength} in
+	    *Y|*y)
+		PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' )
+		(( SECHIBA_WRITE_STEP = PeriodLengthInDays * 86400 / PeriodLengthInYears / 12  ))
+		;;
+	    *M|*m)
+		PeriodLengthInMonths=$( echo ${config_UserChoices_PeriodLength} | awk -F '[mM]' '{print $1}' )
+		(( SECHIBA_WRITE_STEP = PeriodLengthInDays * WriteInMonths  * 86400 / PeriodLengthInMonths  ))
+		;;
+	    *)
+		(( SECHIBA_WRITE_STEP = $( IGCM_date_DaysInMonth $year $month ) * 86400 ))
+		;;
+	    esac
+	    ;;
+	5D|5d) 
+	    (( SECHIBA_WRITE_STEP = 5 * 86400 )) ;;
+	1D|1d) 
+	    (( SECHIBA_WRITE_STEP = 86400 )) ;;
+	*s)
+	    WriteInSeconds=$( echo ${SRF_WriteFrequency} | awk -F '[s]' '{print $1}' )
+	    (( SECHIBA_WRITE_STEP = WriteInSeconds )) ;;
+	*) 
+	    IGCM_debug_Exit "SRF_Update " ${SRF_WriteFrequency} " invalid WriteFrequency : choose in 1Y, 1M, 5D, 1D." 
+	    IGCM_debug_Verif_Exit ;;
+    esac
+
+    # Dans le cas de test :
+#    if [ ${year} -eq 1979 ] ; then
+#	ORCHIDEE_sed VEGET_REINIT y
+#	ORCHIDEE_sed VEGET_YEAR 131
+#    fi
+#    
+#    # Dans le cas historique :
+#    if [ ${year} -eq 2005 ] ; then
+#	ORCHIDEE_sed VEGET_REINIT y
+#	ORCHIDEE_sed VEGET_YEAR 1
+##	ORCHIDEE_sed VEGET_UPDATE '0Y'
+#    fi
+
+#     # Dans le cas sans changement de carte aprÃšs 2005 
+#     if [ ${year} -gt 2005 ] ; then
+# 	ORCHIDEE_sed VEGET_UPDATE '0Y'
+#     fi
+ 
+    ORCHIDEE_sed HYDROL_CWRR ${orchidee_UserChoices_NEWHYDROL}
+
+    ORCHIDEE_sed WRITE_STEP ${SECHIBA_WRITE_STEP}
+    ORCHIDEE_sed SECHIBA_HISTLEVEL ${orchidee_UserChoices_sechiba_LEVEL}
+
+    # Outputs HF in HISTFILE2 if required
+    [ X${SRF_ok_hf} = Xy ] && ORCHIDEE_sed SECHIBA_HISTFILE2 y
+    ORCHIDEE_sed SECHIBA_HISTLEVEL2 1
+    ORCHIDEE_sed WRITE_STEP2 10800.0
+
+    if ( [ ${CumulPeriod} -ne 1 ] || [ "${config_SRF_Restart}" != "n" ] ) ; then
+	ORCHIDEE_sed SECHIBA_restart_in sechiba_rest_in.nc
+    fi
+
+# Desactaivate river routing for high resolution or very low resolution
+    if ( [ ${RESOL_ATM_3D} = "280x280x39-ZAf" ] || [ ${RESOL_ATM_3D} = "360x180x19-ZIn" ] || [ ${RESOL_ATM_3D} = "32x24x11" ] ) ; then
+	ORCHIDEE_sed  RIVER_ROUTING n
+    fi
+
+# Parameters set in orchidee.card
+    if [ ! X${orchidee_UserChoices_VEGET_UPDATE} = X ] ; then
+	ORCHIDEE_sed VEGET_UPDATE   ${orchidee_UserChoices_VEGET_UPDATE}
+    fi
+    if [ ! X${orchidee_UserChoices_LAND_COVER_CHANGE} = X ] ; then
+	ORCHIDEE_sed LAND_COVER_CHANGE ${orchidee_UserChoices_LAND_COVER_CHANGE}
+    fi
+
+    #IGCM_sys_Cp ${RUN_DIR}/orchidee.def ${RUN_DIR}/run.def
+    #IGCM_sys_Put_Out ${RUN_DIR}/run.def ${R_SAVE}/${PREFIX}_run.def
+
+    IGCM_debug_PopStack "SRF_Update"
+}
+
+#-----------------------------------------------------------------
+function SRF_Finalize
+{
+    IGCM_debug_PushStack "SRF_Finalize"
+
+    #IGCM_sys_Put_Out ${RUN_DIR}/used_run.def ${R_SAVE}/${PREFIX}_used_run.def
+
+    echo FINALIZE SRF !!!
+
+    IGCM_debug_PopStack "SRF_Finalize"
+}
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/stomate.card
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/stomate.card	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/stomate.card	(revision 1738)
@@ -0,0 +1,47 @@
+[Compatibility]
+libIGCM=1.0
+
+[UserChoices]
+stomate_LEVEL=10
+
+[InitialStateFiles]
+List=   ()
+
+[BoundaryFiles]
+List=	()
+ListNonDel= (${R_BC}/SRF/LMDZORINCA/reftemp.nc, .)
+
+[SmoothFiles]
+List=	()
+
+[ParametersFiles]
+List=   ()
+
+[RestartFiles]
+# List restart that have to be saved/restored each loop (file out, saved, and in) :
+List=	(stomate_rest_out.nc, stomate_rest.nc, stomate_rest_in.nc)
+
+[OutputText]
+List=   ()
+
+[OutputFiles]
+List=	(stomate_history.nc, ${R_OUT_SBG_O_M}/${PREFIX}_1M_stomate_history.nc, Post_1M_stomate_history) \
+    (stomate_ipcc_history.nc, ${R_OUT_SBG_O_M}/${PREFIX}_1M_stomate_ipcc_history.nc, Post_1M_stomate_ipcc_history)
+
+[Post_1M_stomate_history]
+Patches= ()
+GatherWithInternal= (lon, lat, PFT, time_counter, time_counter_bnds, Areas)
+TimeSeriesVars2D= (T2M_MONTH,CONTFRAC,RESOLUTION_X,RESOLUTION_Y,CONVFLUX,CFLUX_PROD10,CFLUX_PROD100,HARVEST_ABOVE)
+ChunckJob2D = NONE
+TimeSeriesVars3D = (CO2FLUX_MONTHLY,LAI,VEGET,VEGET_MAX,NPP,GPP,HET_RESP,MAINT_RESP,GROWTH_RESP,AGE,HEIGHT,TOTAL_M,TOTAL_BM_LITTER,TOTAL_SOIL_CARB,CO2_FIRE)
+ChunckJob3D = NONE
+Seasonal=ON
+
+[Post_1M_stomate_ipcc_history]
+Patches= ()
+GatherWithInternal= (lon, lat, PFT, time_counter, time_counter_bnds, Areas, CONTFRAC)
+TimeSeriesVars2D=(cVeg, cLitter, cSoil, cProduct, lai, gpp, ra, npp, rh, fFire, fHarvest, fLuc, nbp, fVegLitter, fLitterSoil, cLeaf, cWood, cRoot, cMisc, cLitterAbove, cLitterBelow, cSoilFast, cSoilMedium, cSoilSlow, landCoverFrac, treeFracPrimDec, treeFracPrimEver, c3PftFrac, c4PftFrac, rGrowth, rMaint, nppLeaf, nppWood, nppRoot, nep)
+ChunckJob2D = NONE
+TimeSeriesVars3D=()
+ChunckJob3D = NONE
+Seasonal=ON
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/stomate.driver
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/stomate.driver	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/COMP/stomate.driver	(revision 1738)
@@ -0,0 +1,94 @@
+#!/bin/ksh
+
+#D- Driver du script pour ORCHIDEE
+
+#-----------------------------------------------------------------
+function SBG_Initialize
+{
+    IGCM_debug_PushStack "SBG_Initialize"
+
+    RESOL_SBG=ALL
+    
+    IGCM_debug_PopStack "SBG_Initialize"
+}
+
+#-----------------------------------------------------------------
+function SBG_Update
+{
+    IGCM_debug_PushStack "SBG_Update"
+
+    typeset STOMATE_WRITE_STEP
+
+    case ${config_SBG_WriteFrequency} in
+	*Y|*y) 
+	    WriteInYears=$( echo ${config_SBG_WriteFrequency} | awk -F '[yY]' '{print $1}' )
+	    PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' )
+	    (( STOMATE_WRITE_STEP = PeriodLengthInDays * WriteInYears / PeriodLengthInYears )) ;;
+	1M|1m) 
+	    case ${config_UserChoices_PeriodLength} in
+	    *Y|*y)
+		STOMATE_WRITE_STEP=-1.
+		;;
+	    *M|*m)
+		STOMATE_WRITE_STEP=-1.
+		;;
+	    *)
+		(( STOMATE_WRITE_STEP = $( IGCM_date_DaysInMonth $year $month ) )) ;;
+	    esac
+	    ;;
+	*M|*m) 
+	    WriteInMonths=$( echo ${config_SBG_WriteFrequency} | awk -F '[mM]' '{print $1}' )
+	    case ${config_UserChoices_PeriodLength} in
+	    *Y|*y)
+		PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' )
+		(( STOMATE_WRITE_STEP = 30 ))
+		;;
+	    *M|*m)
+		PeriodLengthInMonths=$( echo ${config_UserChoices_PeriodLength} | awk -F '[mM]' '{print $1}' )
+		(( STOMATE_WRITE_STEP = PeriodLengthInDays * WriteInMonths / PeriodLengthInMonths ))
+		;;
+	    *)
+		(( STOMATE_WRITE_STEP = $( IGCM_date_DaysInMonth $year $month ) )) ;;
+	    esac
+	    ;;
+	5D|5d) 
+	    (( STOMATE_WRITE_STEP = 5 )) ;;
+	1D|1d) 
+	    (( STOMATE_WRITE_STEP = 1 )) ;;
+	*s)
+	    WriteInSeconds=$( echo ${config_SBG_WriteFrequency} | awk -F '[s]' '{print $1}' )
+	    (( STOMATE_WRITE_STEP = 1 )) ;;
+	*) 
+	    IGCM_debug_Exit "SBG_Update " ${config_SBG_WriteFrequency} " invalid WriteFrequency : choose in 1Y, 1M, 5D, 1D." 
+	    IGCM_debug_Verif_Exit ;;
+    esac
+
+    ORCHIDEE_sed STOMATE_OK_STOMATE y
+
+    ORCHIDEE_sed STOMATE_HIST_DT ${STOMATE_WRITE_STEP}
+    ORCHIDEE_sed STOMATE_HISTLEVEL ${stomate_UserChoices_stomate_LEVEL}
+
+    if ( [ ${CumulPeriod} -eq 1 ] && [ "${config_SBG_Restart}" = "n" ] ) ; then
+	echo "STOMATE : without restart"
+    else
+	ORCHIDEE_sed STOMATE_RESTART_FILEIN stomate_rest_in.nc
+    fi
+
+    IGCM_debug_PopStack "SBG_Update"
+}
+
+#-----------------------------------------------------------------
+function SBG_Finalize
+{
+#set -vx
+    IGCM_debug_PushStack "SBG_Finalize"
+
+    if [ $( IGCM_date_DaysBetweenGregorianDate ${PeriodDateEnd} ${year}1230 ) -ge 0 ] ; then
+	IGCM_sys_Put_Out stomate_Cforcing.nc ${R_OUT_SBG_R}/${config_UserChoices_JobName}_${PeriodDateEnd}_stomate_Cforcing.nc
+	IGCM_sys_Put_Out stomate_forcing.nc  ${R_OUT_SBG_R}/${config_UserChoices_JobName}_${PeriodDateEnd}_stomate_forcing.nc
+# 	rm -f stomate_Cforcing.nc
+# 	rm -f stomate_forcing.nc
+    fi
+
+    IGCM_debug_PopStack "SBG_Finalize"
+}
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L19_AP
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L19_AP	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L19_AP	(revision 1738)
@@ -0,0 +1,53 @@
+## $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
+## 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
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+iphysiq=10
+##  Avec ou sans strato
+ok_strato=n
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=1
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=1.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom
+clat=0.
+## facteur de grossissement du zoom,selon longitude
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.
+ysinus=y  
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L39_AP
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L39_AP	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L39_AP	(revision 1738)
@@ -0,0 +1,53 @@
+## $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
+## choix de l'operateur de dissipation (star ou  non star )
+lstardis=y
+## nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+## nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+## nombre d'iterations de l'operateur de dissipation  divgrad            
+niterh=2
+## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+tetagdiv=5400.
+## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=5400.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=5400.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+iphysiq=10
+##  Avec ou sans strato
+ok_strato=y
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=2
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=5.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom
+clat=0.
+## facteur de grossissement du zoom,selon longitude
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.
+ysinus=y  
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L39_NPv3.1
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L39_NPv3.1	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/gcm.def_L39_NPv3.1	(revision 1738)
@@ -0,0 +1,53 @@
+## $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
+## choix de l'operateur de dissipation (star ou  non star )
+lstardis=y
+## nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+## nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+## nombre d'iterations de l'operateur de dissipation  divgrad            
+niterh=2
+## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+tetagdiv=5400.
+## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=5400.
+## temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=5400.
+## coefficient pour gamdissip                                            
+coefdis=0.
+## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+purmats=n
+## avec ou sans physique                                                 
+physic=y
+## periode de la physique (en pas)                                       
+iphysiq=5
+##  Avec ou sans strato
+ok_strato=y
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=2
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=5.e-5
+## longitude en degres du centre du zoom                                 
+clon=0.
+## latitude en degres du centre du zoom
+clat=0.
+## facteur de grossissement du zoom,selon longitude
+grossismx=1.0
+## facteur de grossissement du zoom ,selon latitude
+grossismy=1.0
+##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale
+fxyhypb=y
+## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+dzoomx=0.0
+## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+dzoomy=0.0
+##raideur du zoom en  X
+taux=3.
+##raideur du zoom en  Y
+tauy=3.
+##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.
+ysinus=y  
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/guide.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/guide.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/guide.def	(revision 1738)
@@ -0,0 +1,17 @@
+# ok_guide=y : guidage active
+ok_guide=n
+#
+ncep=.false. # true si ncep (ou era 40) , false si ecmwf
+guide_u= y
+guide_v= y
+guide_T= n
+guide_P= n
+guide_Q= n
+
+ini_anal=n
+tau_min_u=0.0208333
+tau_max_u=0.1
+tau_min_v=0.0208333
+tau_max_v=0.1
+tau_min_T=0.0208333
+tau_max_T=10
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/offline.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/offline.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/offline.def	(revision 1738)
@@ -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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/orchidee.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/orchidee.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/orchidee.def	(revision 1738)
@@ -0,0 +1,887 @@
+#
+#**************************************************************************
+#                    Namelist for ORCHIDEE
+#**************************************************************************
+#
+#
+#**************************************************************************
+#          OPTIONS NOT SET
+#**************************************************************************
+#
+#
+#**************************************************************************
+#          Management of display in the run of ORCHIDEE
+#**************************************************************************
+
+# Model chatting level
+# level of online diagnostics in STOMATE (0-4)
+# With this variable, you can determine how much online information STOMATE
+#  gives during the run. 0 means virtually no info.
+BAVARD = 1
+# default = 1
+
+# Flag for debug information
+# This option allows to switch on the output of debug
+#         information without recompiling the code.
+DEBUG_INFO = n
+#default = n
+
+# ORCHIDEE will print more messages
+# This flag permits to print more debug messages in the run.
+LONGPRINT = n
+#default = n
+
+#---------------------------------------------------------------------
+
+# To reset the time coming from SECHIBA restart file
+# This option allows the model to override the time
+#  found in the restart file of SECHIBA with the time
+#  of the first call. That is the restart time of the GCM.
+SECHIBA_reset_time = y
+# default = n
+
+#**************************************************************************
+#          Files : incoming / forcing / restart /output
+#**************************************************************************
+# Ancillary files :
+#---------------------------------------------------------------------
+
+# Name of file from which the vegetation map is to be read
+# If !IMPOSE_VEG
+# If LAND_USE 
+#   default = pft_new.nc
+#   The name of the file to be opened to read a vegetation
+#   map (in pft) is to be given here. 
+# If !LAND_USE
+#   default = ../surfmap/carteveg5km.nc
+#   The name of the file to be opened to read the vegetation
+#   map is to be given here. Usualy SECHIBA runs with a 5kmx5km
+#   map which is derived from the IGBP one. We assume that we have
+#   a classification in 87 types. This is Olson modified by Viovy.
+VEGETATION_FILE = PFTmap.nc
+
+
+# Name of file from which the bare soil albedo
+# If !IMPOSE_AZE
+# The name of the file to be opened to read the soil types from 
+#  which we derive then the bare soil albedos. This file is 1x1 
+#  deg and based on the soil colors defined by Wilson and Henderson-Seller.
+SOILALB_FILE = soils_param.nc
+# default = ../surfmap/soils_param.nc
+
+# Name of file from which soil types are read
+# If !IMPOSE_VEG
+# The name of the file to be opened to read the soil types. 
+#  The data from this file is then interpolated to the grid of
+#  of the model. The aim is to get fractions for sand loam and
+#  clay in each grid box. This information is used for soil hydrology
+#  and respiration.
+SOILTYPE_FILE = soils_param.nc
+# default = ../surfmap/soils_param.nc
+
+# Name of file from which the reference
+# The name of the file to be opened to read
+#  temperature is read
+#  the reference surface temperature.
+#  The data from this file is then interpolated
+#  to the grid of the model.
+#  The aim is to get a reference temperature either
+#  to initialize the corresponding prognostic model
+#  variable correctly (ok_dgvm = TRUE) or to impose it
+#  as boundary condition (ok_dgvm = FALSE)
+REFTEMP_FILE = reftemp.nc
+# default = reftemp.nc
+
+# Input and output restart file for SECHIBA :
+#---------------------------------------------------------------------
+
+# Name of restart to READ for initial conditions
+# This is the name of the file which will be opened
+#  to extract the initial values of all prognostic
+#  values of the model. This has to be a netCDF file.
+#  Not truly COADS compliant. NONE will mean that
+#  no restart file is to be expected.
+SECHIBA_restart_in = NONE
+# default = NONE
+
+# Name of restart files to be created by SECHIBA
+# This variable give the name for the restart files. 
+#  The restart software within IOIPSL will add .nc if needed.
+SECHIBA_rest_out = sechiba_rest_out.nc
+# default = sechiba_rest_out.nc
+
+# Input and output restart file for STOMATE :
+#---------------------------------------------------------------------
+
+# Name of restart to READ for initial conditions of STOMATE
+# If STOMATE_OK_STOMATE || STOMATE_WATCHOUT
+# This is the name of the file which will be opened of STOMATE
+#   to extract the initial values of all prognostic values of STOMATE.
+STOMATE_RESTART_FILEIN = NONE
+# default = NONE
+
+# Name of restart files to be created by STOMATE
+# If STOMATE_OK_STOMATE || STOMATE_WATCHOUT
+# This is the name of the file which will be opened
+#        to write the final values of all prognostic values
+#        of STOMATE.
+STOMATE_RESTART_FILEOUT = stomate_rest_out.nc
+# default = stomate_restart.nc
+
+# Forcing files for TESTSTOMATE and FORCESOIL
+#---------------------------------------------------------------------
+
+# Name of STOMATE's forcing file
+# Name that will be given to STOMATE's offline forcing file
+#STOMATE_FORCING_NAME = stomate_forcing.nc
+#default = NONE
+
+# Size of STOMATE forcing data in memory (MB)
+# This variable determines how many
+#  forcing states will be kept in memory.
+#  Must be a compromise between memory
+#  use and frequeny of disk access.
+STOMATE_FORCING_MEMSIZE = 50
+# default = 50
+
+# Name of STOMATE's carbon forcing file
+# Name that will be given to STOMATE's carbon offline forcing file
+#STOMATE_CFORCING_NAME = stomate_Cforcing.nc
+# default = NONE
+
+
+# Produced forcing file name (SECHIBA puis STOMATE) :
+#---------------------------------------------------------------------
+
+# ORCHIDEE will write out its forcing to a file
+# This flag allows to write to a file all the variables
+#  which are used to force the land-surface. The file 
+#  has exactly the same format than a normal off-line forcing
+#  and thus this forcing can be used for forcing ORCHIDEE.
+#ORCHIDEE_WATCHOUT = y
+# default = n
+
+# Filenane for the ORCHIDEE forcing file
+# If ORCHIDEE_WATCHOUT
+# This is the name of the file in which the
+#  forcing used here will be written for later use. 
+WATCHOUT_FILE = orchidee_watchout.nc
+# default = orchidee_watchout.nc
+
+# ORCHIDEE will write out with this frequency
+# If ORCHIDEE_WATCHOUT
+# This flag indicates the frequency of the write of the variables. 
+DT_WATCHOUT = 1800
+# default = dt
+
+# STOMATE does minimum service
+# set to TRUE if you want STOMATE to read
+#  and write its start files and keep track
+#  of longer-term biometeorological variables.
+#  This is useful if OK_STOMATE is not set,
+#  but if you intend to activate STOMATE later.
+#  In that case, this run can serve as a 
+#  spinup for longer-term biometeorological
+#  variables.
+#STOMATE_WATCHOUT = y
+# default = n
+
+# Output file name (SECHIBA and STOMATE) :
+#---------------------------------------------------------------------
+# Name of file in which the output is going
+# This file is going to be created by the model
+#  to be written
+#  and will contain the output from the model.
+#  This file is a truly COADS compliant netCDF file.
+#  It will be generated by the hist software from
+#  the IOIPSL package.
+OUTPUT_FILE = sechiba_history.nc
+# default = cabauw_out.nc
+
+# Flag to switch on histfile 2 for SECHIBA (hi-frequency ?)
+# This Flag switch on the second SECHIBA writing for hi (or low) 
+#  frequency writing. This second output is optional and not written
+#  by default.
+SECHIBA_HISTFILE2 = FALSE
+# default  = FALSE
+
+# Name of file in which the output number 2 is going
+#   to be written
+# If SECHIBA_HISTFILE2
+# This file is going to be created by the model
+#   and will contain the output 2 from the model.
+SECHIBA_OUTPUT_FILE2 = sechiba_out_2.nc
+# default  = sechiba_out_2.nc
+
+# Name of file in which STOMATE's output is going to be written
+# This file is going to be created by the model
+#  and will contain the output from the model.
+#  This file is a truly COADS compliant netCDF file.
+#  It will be generated by the hist software from
+#  the IOIPSL package.
+STOMATE_OUTPUT_FILE = stomate_history.nc
+# default = stomate_history.nc
+
+# Write levels for outputs files (number of variables) :
+#---------------------------------------------------------------------
+
+# SECHIBA history output level (0..10)
+# Chooses the list of variables in the history file. 
+#  Values between 0: nothing is written; 10: everything is 
+#  written are available More details can be found on the web under documentation.
+#  web under documentation.
+SECHIBA_HISTLEVEL = 5
+# default = 5
+
+# SECHIBA history 2 output level (0..10)
+# If SECHIBA_HISTFILE2
+# Chooses the list of variables in the history file. 
+#   Values between 0: nothing is written; 10: everything is 
+#   written are available More details can be found on the web under documentation.
+#   web under documentation.
+# First level contains all ORCHIDEE outputs.
+SECHIBA_HISTLEVEL2 = 1
+# default = 1
+
+# STOMATE history output level (0..10)
+#  0: nothing is written; 10: everything is written
+STOMATE_HISTLEVEL = 10
+# default = 10
+
+#--------------------------------------------------------------------
+# STOMATE_IPCC_OUTPUT_FILE
+# This file is going to be created by the model
+#    and will contain the output from the model.
+#    This file is a truly COADS compliant netCDF file.
+#    It will be generated by the hist software from
+#    the IOIPSL package.
+# Name of file in which STOMATE's output is going
+# to be written
+STOMATE_IPCC_OUTPUT_FILE = stomate_ipcc_history.nc
+# default = stomate_ipcc_history.nc
+
+# STOMATE_IPCC_HIST_DT
+# Time step of the STOMATE IPCC history file
+# STOMATE IPCC history time step (d)
+STOMATE_IPCC_HIST_DT = -1
+# default = 0.
+
+# Write frequency for output files (SECHIBA in seconds et
+# STOMATE in days) :
+#---------------------------------------------------------------------
+# Frequency in seconds at which to WRITE output
+# This variables gives the frequency the output of
+#  the model should be written into the netCDF file.
+#  It does not affect the frequency at which the
+#  operations such as averaging are done.
+WRITE_STEP = 86400.0
+# default = 86400.0
+
+# Frequency in seconds at which to WRITE output
+# If SECHIBA_HISTFILE2
+# This variables gives the frequency the output 2 of
+#   the model should be written into the netCDF file.
+#   It does not affect the frequency at which the
+#   operations such as averaging are done.
+#   That is IF the coding of the calls to histdef
+#   are correct !
+WRITE_STEP2 = 1800.0
+# default = 1800.0
+
+# STOMATE history time step (d)
+# Time step of the STOMATE history file
+# Care : this variable must be higher than DT_SLOW
+STOMATE_HIST_DT = 10.
+# default = 10.
+
+#---------------------------------------------------------------------
+# FORCESOIL CARBON spin up parametrization
+#---------------------------------------------------------------------
+
+# Number of time steps per year for carbon spinup.
+FORCESOIL_STEP_PER_YEAR = 12
+# default = 12
+
+# Number of years saved for carbon spinup.
+FORCESOIL_NB_YEAR = 1
+# default = 1
+
+#---------------------------------------------------------------------
+# Parametrization :
+#---------------------------------------------------------------------
+
+# Activate STOMATE?
+# set to TRUE if STOMATE is to be activated
+STOMATE_OK_STOMATE = n
+# default = n
+
+# Activate DGVM?
+# set to TRUE if Dynamic Vegetation DGVM is to be activated
+STOMATE_OK_DGVM = n
+# default = n
+
+# Activate CO2?
+# set to TRUE if photosynthesis is to be activated
+STOMATE_OK_CO2 = y
+# default = n
+
+# Flag to force the value of atmospheric CO2 for vegetation.
+# If this flag is set to true, the ATM_CO2 parameter is used
+#  to prescribe the atmospheric CO2.
+# This Flag is only use in couple mode.
+FORCE_CO2_VEG = FALSE
+# default = FALSE
+
+# Value for atm CO2.
+# If FORCE_CO2_VEG (in not forced mode)
+# Value to prescribe the atm CO2.
+#  For pre-industrial simulations, the value is 286.2 .
+#  348. for 1990 year.
+ATM_CO2 = 350.
+# default = 350.
+
+# constant tree mortality
+# If yes, then a constant mortality is applied to trees. 
+#  Otherwise, mortality is a function of the trees' 
+#  vigour (as in LPJ).
+LPJ_GAP_CONST_MORT = y
+# default = y
+
+# no fire allowed
+# With this variable, you can allow or not
+#  the estimation of CO2 lost by fire
+FIRE_DISABLE = n
+# default = n
+
+# Average method for z0
+# If this flag is set to true (y) then the neutral Cdrag
+#  is averaged instead of the log(z0). This should be
+#  the prefered option. We still wish to keep the other
+#  option so we can come back if needed. If this is
+#  desired then one should set Z0CDRAG_AVE = n
+Z0CDRAG_AVE = y
+# default = y
+
+# parameters describing the surface (vegetation + soil) :
+#---------------------------------------------------------------------
+#
+# Should the vegetation be prescribed
+# This flag allows the user to impose a vegetation distribution
+#  and its characterisitcs. It is espacially interesting for 0D
+#  simulations. On the globe it does not make too much sense as
+#  it imposes the same vegetation everywhere
+IMPOSE_VEG = n
+# default = n
+
+# Flag to use old "interpolation" of vegetation map.
+# IF NOT IMPOSE_VEG and NOT LAND_USE
+#  If you want to recover the old (ie orchidee_1_2 branch) 
+#   "interpolation" of vegetation map.
+SLOWPROC_VEGET_OLD_INTERPOL = n
+# default = n
+
+# Vegetation distribution within the mesh (0-dim mode)
+# If IMPOSE_VEG
+# The fraction of vegetation is read from the restart file. If
+#  it is not found there we will use the values provided here.
+SECHIBA_VEG__01 = 0.2
+SECHIBA_VEG__02 = 0.0
+SECHIBA_VEG__03 = 0.0
+SECHIBA_VEG__04 = 0.0
+SECHIBA_VEG__05 = 0.0
+SECHIBA_VEG__06 = 0.0
+SECHIBA_VEG__07 = 0.0
+SECHIBA_VEG__08 = 0.0
+SECHIBA_VEG__09 = 0.0
+SECHIBA_VEG__10 = 0.8
+SECHIBA_VEG__11 = 0.0
+SECHIBA_VEG__12 = 0.0
+SECHIBA_VEG__13 = 0.0
+# default = 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8, 0.0, 0.0, 0.0
+
+# Maximum vegetation distribution within the mesh (0-dim mode)
+# If IMPOSE_VEG
+# The fraction of vegetation is read from the restart file. If
+#  it is not found there we will use the values provided here.
+SECHIBA_VEGMAX__01 = 0.2
+SECHIBA_VEGMAX__02 = 0.0
+SECHIBA_VEGMAX__03 = 0.0
+SECHIBA_VEGMAX__04 = 0.0
+SECHIBA_VEGMAX__05 = 0.0
+SECHIBA_VEGMAX__06 = 0.0
+SECHIBA_VEGMAX__07 = 0.0
+SECHIBA_VEGMAX__08 = 0.0
+SECHIBA_VEGMAX__09 = 0.0
+SECHIBA_VEGMAX__10 = 0.8
+SECHIBA_VEGMAX__11 = 0.0
+SECHIBA_VEGMAX__12 = 0.0
+SECHIBA_VEGMAX__13 = 0.0
+# default = 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8, 0.0, 0.0, 0.0
+
+# LAI for all vegetation types (0-dim mode)
+# If IMPOSE_VEG
+# The maximum LAI used in the 0dim mode. The values should be found
+#  in the restart file. The new values of LAI will be computed anyway
+#  at the end of the current day. The need for this variable is caused
+#  by the fact that the model may stop during a day and thus we have not
+#  yet been through the routines which compute the new surface conditions.
+SECHIBA_LAI__01 = 0.
+SECHIBA_LAI__02 = 8.
+SECHIBA_LAI__03 = 8.
+SECHIBA_LAI__04 = 4.
+SECHIBA_LAI__05 = 4.5
+SECHIBA_LAI__06 = 4.5
+SECHIBA_LAI__07 = 4.
+SECHIBA_LAI__08 = 4.5
+SECHIBA_LAI__09 = 4.
+SECHIBA_LAI__10 = 2.
+SECHIBA_LAI__11 = 2.
+SECHIBA_LAI__12 = 2.
+SECHIBA_LAI__13 = 2.
+# default = 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2.
+
+# Height for all vegetation types (m)
+# If IMPOSE_VEG
+# The height used in the 0dim mode. The values should be found
+#  in the restart file. The new values of height will be computed anyway
+#  at the end of the current day. The need for this variable is caused
+#  by the fact that the model may stop during a day and thus we have not
+#  yet been through the routines which compute the new surface conditions.
+SLOWPROC_HEIGHT__01 = 0.
+SLOWPROC_HEIGHT__02 = 50.
+SLOWPROC_HEIGHT__03 = 50.
+SLOWPROC_HEIGHT__04 = 30.
+SLOWPROC_HEIGHT__05 = 30.
+SLOWPROC_HEIGHT__06 = 30.
+SLOWPROC_HEIGHT__07 = 20.
+SLOWPROC_HEIGHT__08 = 20.
+SLOWPROC_HEIGHT__09 = 20.
+SLOWPROC_HEIGHT__10 = .2
+SLOWPROC_HEIGHT__11 = .2
+SLOWPROC_HEIGHT__12 = .4
+SLOWPROC_HEIGHT__13 = .4
+# default = 0., 30., 30., 20., 20., 20., 15., 15., 15., .5, .6, 1.0, 1.0
+
+
+# Fraction of the 3 soil types (0-dim mode)
+# If IMPOSE_VEG
+# Determines the fraction for the 3 soil types
+#  in the mesh in the following order : sand loam and clay.
+SOIL_FRACTIONS__01 = 0.28
+SOIL_FRACTIONS__02 = 0.52
+SOIL_FRACTIONS__03 = 0.20
+# default = 0.28, 0.52, 0.20
+
+# Fraction of other surface types within the mesh (0-dim mode)
+# If IMPOSE_VEG
+# The fraction of ice, lakes, etc. is read from the restart file. If
+#  it is not found there we will use the values provided here.
+#  For the moment, there is only ice.
+SECHIBA_FRAC_NOBIO = 0.0
+# default = 0.0
+
+# Fraction of the clay fraction (0-dim mode)
+# If IMPOSE_VEG
+# Determines the fraction of clay in the grid box.
+CLAY_FRACTION = 0.2
+# default = 0.2
+
+# Should the surface parameters be prescribed
+# This flag allows the user to impose the surface parameters
+#  (Albedo Roughness and Emissivity). It is espacially interesting for 0D
+#  simulations. On the globe it does not make too much sense as
+#  it imposes the same vegetation everywhere
+IMPOSE_AZE = n
+# default = n
+
+# Emissivity of the surface for LW radiation
+# If IMPOSE_AZE
+# The surface emissivity used for compution the LE emission
+#  of the surface in a 0-dim version. Values range between 
+#  0.97 and 1.. The GCM uses 0.98.
+CONDVEG_EMIS = 1.0
+# default = 1.0
+
+# SW visible albedo for the surface
+# If IMPOSE_AZE
+# Surface albedo in visible wavelengths to be used 
+#  on the point if a 0-dim version of SECHIBA is used. 
+#  Look at the description of the forcing data for 
+#  the correct value.
+CONDVEG_ALBVIS = 0.25
+# default = 0.25
+
+# SW near infrared albedo for the surface
+# If IMPOSE_AZE
+# Surface albedo in near infrared wavelengths to be used 
+#  on the point if a 0-dim version of SECHIBA is used. 
+#  Look at the description of the forcing data for 
+#  the correct value.
+CONDVEG_ALBNIR = 0.25
+# default = 0.25
+
+# Surface roughness (m)
+# If IMPOSE_AZE
+# Surface rougness to be used on the point if a 0-dim version
+#  of SECHIBA is used. Look at the description of the forcing  
+#  data for the correct value.
+CONDVEG_Z0 = 0.15
+# default = 0.15_stnd
+
+# Height to be added to the height of the first level (m)
+# If IMPOSE_AZE
+# ORCHIDEE assumes that the atmospheric level height is counted
+#  from the zero wind level. Thus to take into account the roughness
+#  of tall vegetation we need to correct this by a certain fraction
+#  of the vegetation height. This is called the roughness height in
+#  ORCHIDEE talk.
+ROUGHHEIGHT = 0.0
+# default = 0.0
+
+# The snow albedo used by SECHIBA
+# This option allows the user to impose a snow albedo.
+#  Default behaviour is to use the model of snow albedo
+#  developed by Chalita (1993).
+CONDVEG_SNOWA = default
+# default = use the model of snow albedo developed by Chalita
+
+# Switch bare soil albedo dependent (if TRUE) on soil wetness
+# If TRUE, the model for bare soil albedo is the old formulation.
+#  Then it depend on the soil dry or wetness. If FALSE, it is the 
+#  new computation that is taken, it is only function of soil color.
+ALB_BARE_MODEL = FALSE
+# default = FALSE
+
+# Initial snow mass if not found in restart
+# The initial value of snow mass if its value is not found
+#   in the restart file. This should only be used if the model is 
+#   started without a restart file.
+HYDROL_SNOW = 0.0
+# default = 0.0
+
+
+# Initial snow age if not found in restart
+# The initial value of snow age if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+HYDROL_SNOWAGE = 0.0
+# default = 0.0
+
+# Initial snow amount on ice, lakes, etc. if not found in restart
+# The initial value of snow if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+HYDROL_SNOW_NOBIO = 0.0
+# default = 0.0
+
+# Initial snow age on ice, lakes, etc. if not found in restart
+# The initial value of snow age if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+HYDROL_SNOW_NOBIO_AGE = 0.0
+# default = 0.0
+
+# Initial soil moisture stress if not found in restart
+# The initial value of soil moisture stress if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+HYDROL_HUMR = 1.0
+# default = 1.0
+
+# Total depth of soil reservoir
+HYDROL_SOIL_DEPTH = 4.
+# default = 2.
+
+# Root profile
+# Default values were defined for 2 meters soil depth.
+# For 4 meters soil depth, you may use those ones :
+# 5., .4, .4, 1., .8, .8, 1., 1., .8, 4., 1., 4., 1.
+HYDROL_HUMCSTE = 5., .4, .4, 1., .8, .8, 1., 1., .8, 4., 1., 4., 1.
+# default =  5., .8, .8, 1., .8, .8, 1., 1., .8, 4., 4., 4., 4.
+
+# Initial restart deep soil moisture if not found in restart
+# The initial value of deep soil moisture if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file. Default behaviour is a saturated soil.
+HYDROL_BQSB = default
+# default = Maximum quantity of water (Kg/M3) * Total depth of soil reservoir = 150. * 2
+
+# Initial upper soil moisture if not found in restart
+# The initial value of upper soil moisture if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+HYDROL_GQSB = 0.0
+# default = 0.0
+
+# Initial upper reservoir depth if not found in restart
+# The initial value of upper reservoir depth if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+HYDROL_DSG = 0.0
+# default = 0.0
+
+# Initial dry soil above upper reservoir if not found in restart
+# The initial value of dry soil above upper reservoir if its value 
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file. The default behaviour
+#  is to compute it from the variables above. Should be OK most of 
+#  the time.
+HYDROL_DSP = default
+# default = Total depth of soil reservoir - HYDROL_BQSB / Maximum quantity of water (Kg/M3) = 0.0
+
+# Initial water on canopy if not found in restart
+# The initial value of moisture on canopy if its value 
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+HYDROL_QSV = 0.0
+# default = 0.0
+
+# Soil moisture on each soil tile and levels
+# The initial value of mc if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+HYDROL_MOISTURE_CONTENT = 0.3
+# default = 0.3
+
+# US_NVM_NSTM_NSLM
+# The initial value of us (relative moisture) if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+US_INIT = 0.0
+# default = 0.0
+
+# Coefficient for free drainage at bottom
+# The initial value of free drainage if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+FREE_DRAIN_COEF = 1.0, 1.0, 1.0
+# default = 1.0, 1.0, 1.0
+
+# Bare soil evap on each soil if not found in restart
+# The initial value of bare soils evap if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+EVAPNU_SOIL = 0.0
+# default = 0.0
+
+# Initial temperature if not found in restart
+# The initial value of surface temperature if its value is not found
+#  in the restart file. This should only be used if the model is 
+#  started without a restart file.
+ENERBIL_TSURF = 280.
+# default = 280.
+
+# Initial Soil Potential Evaporation
+# The initial value of soil potential evaporation if its value 
+#  is not found in the restart file. This should only be used if
+#  the model is started without a restart file. 
+ENERBIL_EVAPOT = 0.0
+# default = 0.0
+
+# Initial soil temperature profile if not found in restart
+# The initial value of the temperature profile in the soil if 
+#   its value is not found in the restart file. This should only 
+#   be used if the model is started without a restart file. Here
+#   we only require one value as we will assume a constant 
+#   throughout the column.
+THERMOSOIL_TPRO = 280.
+# default = 280.
+
+# Initial leaf CO2 level if not found in restart
+# The initial value of leaf_ci if its value is not found
+#  in the restart file. This should only be used if the model is
+#  started without a restart file.
+DIFFUCO_LEAFCI = 233.
+# default = 233.
+
+
+# Keep cdrag coefficient from gcm.
+# Set to .TRUE. if you want q_cdrag coming from GCM.
+#  Keep cdrag coefficient from gcm for latent and sensible heat fluxes.
+#  TRUE if q_cdrag on initialization is non zero (FALSE for off-line runs).
+CDRAG_FROM_GCM = y
+# default =  IF q_cdrag == 0 ldq_cdrag_from_gcm = .FALSE. ELSE .TRUE.
+
+
+# Artificial parameter to increase or decrease canopy resistance
+# Add from Nathalie - the 28 of March 2006 - advice from Fred Hourdin
+# By PFT.
+RVEG_PFT = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.
+# default = 1.
+
+
+# Interception reservoir coefficient.
+# Transforms leaf area index into size of interception reservoir
+#  for slowproc_derivvar or stomate.
+SECHIBA_QSINT = 0.02
+# default = 0.1
+
+#**************************************************************************
+# LAND_USE
+#**************************************************************************
+
+# Read a land_use vegetation map
+# pft values are needed, max time axis is 293
+LAND_USE = y
+# default = n
+
+# Year of the land_use vegetation map readed
+# year off the pft map
+# If LAND_USE (11 = 1860 - 1850 +1 for PFTmap.20C3M.nc, 1 for PFTmap_IPCC_2000.nc) 
+VEGET_YEAR = 1
+# default = 282
+
+# booleen to indicate that a new LAND USE file will be used (since 1.9.5 version).
+# The parameter is used to bypass veget_year count 
+# and reinitialize it with VEGET_YEAR parameter.
+# Then it is possible to change LAND USE file.
+# If LAND_USE
+VEGET_REINIT = n
+# default = n
+
+# Update vegetation frequency (since 2.0 version)
+# The veget datas will be update each this time step.
+# If LAND_USE
+VEGET_UPDATE = XY
+# default = 1Y
+
+# treat land use modifications
+# With this variable, you can use a Land Use map
+# to simulate anthropic modifications such as   
+# deforestation.                                
+# If LAND_USE
+LAND_COVER_CHANGE = x
+# default = y
+
+#**************************************************************************
+
+# agriculture allowed?
+# With this variable, you can determine
+#  whether agriculture is allowed
+AGRICULTURE = y
+# default = y
+
+# Harvert model for agricol PFTs.
+# Compute harvest above ground biomass for agriculture.
+# Change daily turnover.
+HARVEST_AGRI = y
+# default = y
+
+# herbivores allowed?
+# With this variable, you can activate herbivores 
+HERBIVORES = n
+# default = n
+
+# treat expansion of PFTs across a grid cell?
+# With this variable, you can determine
+#  whether we treat expansion of PFTs across a
+#  grid cell.
+TREAT_EXPANSION = n
+# default = n
+
+#**************************************************************************
+
+# Time within the day simulated
+# This is the time spent simulating the current day. This variable is
+#  prognostic as it will trigger all the computations which are
+#  only done once a day.
+SECHIBA_DAY = 0.0
+# default = 0.0
+
+# Time step of STOMATE and other slow processes
+# Time step (s) of regular update of vegetation
+#  cover, LAI etc. This is also the time step
+#  of STOMATE.
+DT_SLOW = 86400.
+# default = un_jour = 86400.
+
+#**************************************************************************
+
+# Allows to switch on the multilayer hydrology of CWRR
+# This flag allows the user to decide if the vertical
+#  hydrology should be treated using the multi-layer 
+#  diffusion scheme adapted from CWRR by Patricia de Rosnay.
+#  by default the Choisnel hydrology is used.
+HYDROL_CWRR = n
+# default = n
+
+# do horizontal diffusion?
+# If TRUE, then water can diffuse horizontally between
+#  the PFTs' water reservoirs.
+HYDROL_OK_HDIFF = n
+# default = n
+ 
+
+# time scale (s) for horizontal diffusion of water
+# If HYDROL_OK_HDIFF
+# Defines how fast diffusion occurs horizontally between
+#  the individual PFTs' water reservoirs. If infinite, no
+#  diffusion.
+HYDROL_TAU_HDIFF = 1800.
+# default = 86400.
+
+# Percent by PFT of precip that is not intercepted by the canopy (since TAG 1.8).
+# During one rainfall event, PERCENT_THROUGHFALL_PFT% of the incident rainfall
+#  will get directly to the ground without being intercepted, for each PFT..
+PERCENT_THROUGHFALL_PFT = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30.
+# default = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30.
+
+# Decides if we route the water or not
+# This flag allows the user to decide if the runoff
+#  and drainage should be routed to the ocean
+#  and to downstream grid boxes.
+RIVER_ROUTING = y
+# default = n
+
+# Name of file which contains the routing information
+# The file provided here should allow the routing module to
+#  read the high resolution grid of basins and the flow direction 
+#  from one mesh to the other.
+ROUTING_FILE = routing.nc
+# default = routing.nc
+
+# Time step of the routing scheme
+# If RIVER_ROUTING
+# This values gives the time step in seconds of the routing scheme. 
+#   It should be multiple of the main time step of ORCHIDEE. One day
+#   is a good value.
+ROUTING_TIMESTEP = 86400
+# default = 86400
+
+# Number of rivers 
+# If RIVER_ROUTING
+# This parameter chooses the number of largest river basins
+#  which should be treated as independently as rivers and not
+#  flow into the oceans as diffusion coastal flow.
+ROUTING_RIVERS = 50
+# default = 50
+
+# Should we compute an irrigation flux 
+# This parameters allows the user to ask the model
+#  to compute an irigation flux. This performed for the
+#  on very simple hypothesis. The idea is to have a good
+#  map of irrigated areas and a simple function which estimates
+#  the need to irrigate.
+DO_IRRIGATION = n
+# default = n
+
+# Name of file which contains the map of irrigated areas
+# If IRRIGATE
+# The name of the file to be opened to read the field
+#  with the area in m^2 of the area irrigated within each
+#  0.5 0.5 deg grid box. The map currently used is the one
+#  developed by the Center for Environmental Systems Research 
+#  in Kassel (1995).
+IRRIGATION_FILE = irrigated.nc
+# default = irrigated.nc
+
+# Should we include floodplains 
+# This parameters allows the user to ask the model
+#  to take into account the flood plains and return 
+#  the water into the soil moisture. It then can go 
+#  back to the atmopshere. This tried to simulate 
+#  internal deltas of rivers.
+DO_FLOODPLAINS = n
+# default = n
+
+#**************************************************************************
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L19_AP
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L19_AP	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L19_AP	(revision 1738)
@@ -0,0 +1,215 @@
+## $Id$
+#
+# PARAMETRES ANCIENNEMENT DANS gcm.def
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=y
+##  Activation ou non de la parametrisation de Hines pour la strato
+ok_hines=n
+##  Orodr  ou  non   pour l orographie              
+ok_orodr=y
+##  Orolf  ou  non   pour l orographie              
+ok_orolf=y
+##   Si = .T. ,  lecture du fichier limit avec la bonne annee             
+ok_limitvrai=n
+## Nombre  d'appels des routines de rayonnements ( par jour)                 
+nbapp_rad=24
+##  Flag  pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect
+iflag_con=30
+## Facteur multiplication des precip convectives dans KE
+cvl_corr=1.0
+##  Facteur additif pour l'albedo
+# variable modife avec lmdz.card
+pmagic=999.9
+#
+#
+# Parametres fichiers de sortie
+#
+# Noms des fichiers
+phys_out_filenames=      histmth histday histhf  histhf3h histhf3hm histstn
+# Sortir ou non les fichiers
+phys_out_filekeys=       y       y       y       y       y          n
+# Niveaux de sorties
+phys_out_filelevels=     5       5       5       5       5          5
+### Type de fichier : global (n) ou stations (y)
+phys_out_filestations =  n       n       n       n       n          y
+### Frequences des sorties
+phys_out_filetimesteps = 1.day,  1.day,  0.25day, 0.125day, 0.125day, 1800.s
+
+# more variables for level 5 : vitw and pres
+# Vertical wind
+flag_vitw         =  2, 3, 5, 5, 1
+name_vitw         =  vitw
+# Air pressure
+flag_pres         =  2, 3, 5, 5, 1
+name_pres         =  pres
+#
+### parametres pour hist*NMC.nc
+#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
+#                histdayNMC: 8 ou 17
+lev_histdayNMC=8
+### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
+ok_histNMC = y, n, n
+
+### freqin_isccp = frequence input en secondes du simulateur ISCCP
+freq_ISCCP=10800.
+### freqout_isccp = frequence output en jours du simulateur ISCCP
+ecrit_ISCCP=_ecrit-ISCCP_
+### niveau du diagnostique de conservation d energie
+if_ebil=0
+#
+# parametres KE
+#
+### epmax = Efficacite precipitation maximale
+epmax = .99
+### ok_adj_ema = ?? pas utilise
+ok_adj_ema = n
+### iflag_clw Flag calcul eau liquide
+# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
+# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
+# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
+iflag_clw = 1
+# 
+# parametres nuages
+#
+### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp) 
+cld_lc_lsc = 2.6e-4
+### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp) 
+cld_lc_con = 2.6e-4
+### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) 
+cld_tau_lsc = 3600.
+### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) 
+cld_tau_con = 3600.
+### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) 
+ffallv_lsc = 1.
+### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) 
+ffallv_con = 1.
+### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
+coef_eva = 2.e-5
+### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) 
+reevap_ice = y
+### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) 
+# - iflag_cldcon<=-1 diag. rain_Tiedtke 
+# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
+# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
+# - iflag_cldcon=3,   ratqs=ratqss
+iflag_cldcon = 3
+### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages 
+# - iflag_pdf=0, version ratqs, 
+# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
+# -              des PDFS de Sandrine Bony 
+iflag_pdf = 1
+### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) 
+fact_cldcon = 1.
+### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
+facttemps = 0.
+## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) 
+ok_newmicro = y
+### iflag_ratqs=0 correspond a la version IPCC AR4 
+iflag_ratqs=0
+### ratqsbas     ratqs en bas si iflag_cldcon=1 
+ratqsbas = 0.005
+### ratqshaut    ratqs en haut pour ratqss "stables"  
+ratqshaut = 0.33
+### rad_froid    rayon cristaux des nuages de glace "froids" 
+rad_froid = 35
+### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev 
+rad_chau1=12
+### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3  
+rad_chau2=11
+# 
+# Coefficient et parametres sur les drags
+#
+f_cdrag_ter=1.
+f_cdrag_stable=1.
+f_cdrag_oce=0.8
+f_rugoro=0.
+### ksta_ter      coef.diffusion minimale sur terre/sic/lic 
+ksta_ter=1.e-7
+### cdmmax =     cdrag maximum pour le moment
+cdmmax = 2.5E-3
+### cdhmax =     cdrag maximum pour l energie
+cdhmax = 2.0E-3
+#
+# Parametres "orbitaux/ ere geologique"
+#
+### R_ecc =      Excentricite
+R_ecc = 0.016715
+### R_peri =     Equinoxe
+R_peri = 102.7
+### R_incl =     Inclinaison
+R_incl = 23.441
+### solaire =    Constante solaire
+solaire = 1365.
+#
+# Taux gaz a effet de serre
+#
+### co2_ppm =    taux CO2 en ppm
+co2_ppm = 348.
+### CH4_ppb =    taux CH4 en ppb
+CH4_ppb = 1650.
+### N2O_ppb =    taux N2O en ppb
+N2O_ppb = 306.
+### CFC11_ppt =  taux CFC11 en ppt
+CFC11_ppt = 280.
+### CFC12_ppt =  taux CFC12 en ppt
+CFC12_ppt = 484
+#
+# Parametres effets directs/indirects des "aerosols"
+#
+### ok_ade=y/n   flag Aerosol direct effect
+ok_ade=n
+### ok_aie=y/n   flag Aerosol indirect effect
+ok_aie=n
+### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
+# variable modife avec lmdz.card
+aer_type=xxx
+###  type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
+flag_aerosol=1 
+### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b0=1.7
+### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b1=0.2
+#
+# Parametre de lecture de l'ozone
+#
+# Allowed values are 0, 1 and 2
+# 0: do not read an ozone climatology
+# 1: read a single ozone climatology that will be used day and night
+# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
+read_climoz=0
+#
+# Parametres simulateur COSP (CFMIP Observational Simulator Package)
+#
+### ok_cosp=y/n flag simulateur COSP
+ok_cosp=n
+### freq_COSP = frequence d'appel de COSP en secondes
+freq_COSP=10800.
+### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
+ok_mensuelCOSP=y
+### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
+ok_journeCOSP=y
+### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
+ok_hfCOSP=n
+#
+# Parametres simulateur ISCCP
+#
+### ok_isccp=y/n flag simulateur ISCCP
+ok_isccp=n
+### top_height = flag choix calcul nuages par le simulateur en utilisant 
+# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
+# - top_height = 1 -> algo IR-VIS
+# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
+# - top_height = 3 -> algo IR
+top_height = 3
+### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
+# - overlap=1    Max overlap
+# - overlap=2    Random overlap
+# - overlap=3    Max/Random overlap
+overlap = 3
+config_inca=xxxx
+aerosol_couple=n
+type_trac=inca
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L39_AP
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L39_AP	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L39_AP	(revision 1738)
@@ -0,0 +1,215 @@
+## $Id$
+#
+# PARAMETRES ANCIENNEMENT DANS gcm.def
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+##  Soil Model  ou non               
+soil_model=y
+##  Choix ou non  de  New oliq               
+new_oliq=y
+##  Activation ou non de la parametrisation de Hines pour la strato
+ok_hines=y
+##  Orodr  ou  non   pour l orographie              
+ok_orodr=y
+##  Orolf  ou  non   pour l orographie              
+ok_orolf=y
+##   Si = .T. ,  lecture du fichier limit avec la bonne annee             
+ok_limitvrai=n
+## Nombre  d'appels des routines de rayonnements ( par jour)                 
+nbapp_rad=24
+##  Flag  pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect
+iflag_con=30
+## Facteur multiplication des precip convectives dans KE
+cvl_corr=1.0
+##  Facteur additif pour l'albedo
+# variable modife avec lmdz.card
+pmagic=999.9
+#
+#
+# Parametres fichiers de sortie
+#
+# Noms des fichiers
+phys_out_filenames=      histmth histday histhf  histhf3h histhf3hm histstn
+# Sortir ou non les fichiers
+phys_out_filekeys=       y       y       y       y       y          n
+# Niveaux de sorties
+phys_out_filelevels=     5       5       5       5       5          5
+### Type de fichier : global (n) ou stations (y)
+phys_out_filestations =  n       n       n       n       n          y
+### Frequences des sorties
+phys_out_filetimesteps = 1.day,  1.day,  0.25day, 0.125day, 0.125day, 1800.s
+
+# more variables for level 5 : vitw and pres
+# Vertical wind
+flag_vitw         =  2, 3, 5, 5, 1
+name_vitw         =  vitw
+# Air pressure
+flag_pres         =  2, 3, 5, 5, 1
+name_pres         =  pres
+#
+### parametres pour hist*NMC.nc
+#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
+#                histdayNMC: 8 ou 17
+lev_histdayNMC=8
+### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
+ok_histNMC = y, n, n
+
+### freqin_isccp = frequence input en secondes du simulateur ISCCP
+freq_ISCCP=10800.
+### freqout_isccp = frequence output en jours du simulateur ISCCP
+ecrit_ISCCP=_ecrit-ISCCP_
+### niveau du diagnostique de conservation d energie
+if_ebil=0
+#
+# parametres KE
+#
+### epmax = Efficacite precipitation maximale
+epmax = .999
+### ok_adj_ema = ?? pas utilise
+ok_adj_ema = n
+### iflag_clw Flag calcul eau liquide
+# - iflag_clw=0 : qcond_incld(i,l) = em_qcondc(l)
+# - iflag_clw=1 : qcond_incld(i,l) = em_qcond(l)
+# - iflag_clw=2 : eau liquide diagnostique en fonction de la Precip
+iflag_clw = 1
+# 
+# parametres nuages
+#
+### cld_lc_lsc  contenu en eau liquide des nuages large-scale (fisrtilp) 
+cld_lc_lsc = 4.16e-4
+### cld_lc_con  contenu en eau liquide des nuages convectifs (fisrtilp) 
+cld_lc_con = 4.16e-4
+### cld_tau_lsc cte de temps utilisee pour eliminer l eau large-scale (fisrtilp) 
+cld_tau_lsc = 1800.
+### cld_tau_con cte de temps utilisee pour eliminer l eau convective (fisrtilp) 
+cld_tau_con = 1800.
+### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace large-scale (fisrtilp) 
+ffallv_lsc = 0.5
+### ffallv_lsc  cte utilisee dans calcul vitesse de chute cristaux de glace convectifs (fisrtilp) 
+ffallv_con = 0.5
+### coef_eva coef evaporation precips eau/glace (fisrtilp/fisrtilp_tr?/conlmd?)
+coef_eva = 2.e-5
+### reevap_ice  reevaporation de toute la precip dans la couche du dessous pour la glace (fisrtilp) 
+reevap_ice = y
+### iflag_cldcon  flag pour calculer ratqsc=F(ratqsbas,fact_cldcon,q_seri) (physiq) 
+# - iflag_cldcon<=-1 diag. rain_Tiedtke 
+# - iflag_cldcon=1, ratqsc=ratqsbas+fact_cldcon*(q_seri(1)-q_seri(k))/q_seri(k)
+# - iflag_cldcon=1/2, ratqs=max(ratqs,ratqsc)
+# - iflag_cldcon=3,   ratqs=ratqss
+iflag_cldcon = 3
+### iflag_pdf :  flag calcul distribution sous-maille de l eau et des nuages 
+# - iflag_pdf=0, version ratqs, 
+# - iflag_pdf=1, calcul eau condensee, fraction nuageuse, eau nuageuse a partir
+# -              des PDFS de Sandrine Bony 
+iflag_pdf = 1
+### fact_cldcon  constante calcul ratqsc (voir iflag_cldcon) et proprietes nuages convectifs, clwcon0 (physiq.F) 
+fact_cldcon = 1.
+### facttemps=   facteur de relaxation de ratqs (iflag_cldcon=1/2) et rnebcon (iflag_cldcon=3)
+facttemps = 0.
+## ok_newmicro   =y appel newmicro , =n appel nuage (calcul epaisseur optique et emmissivite des nuages) 
+ok_newmicro = y
+### iflag_ratqs=0 correspond a la version IPCC AR4 
+iflag_ratqs=0
+### ratqsbas     ratqs en bas si iflag_cldcon=1 
+ratqsbas = 0.005
+### ratqshaut    ratqs en haut pour ratqss "stables"  
+ratqshaut = 0.33
+### rad_froid    rayon cristaux des nuages de glace "froids" 
+rad_froid = 35
+### rad_chau1    rayon goutelettes d eau chauds", en haut: k=4-klev 
+rad_chau1=12
+### rad_chau2    rayon goutelettes d eau chauds", en bas: k=1-3  
+rad_chau2=11
+# 
+# Coefficient et parametres sur les drags
+#
+f_cdrag_ter=1.
+f_cdrag_stable=1.
+f_cdrag_oce=0.8
+f_rugoro=0.
+### ksta_ter      coef.diffusion minimale sur terre/sic/lic 
+ksta_ter=1.e-7
+### cdmmax =     cdrag maximum pour le moment
+cdmmax = 2.5E-3
+### cdhmax =     cdrag maximum pour l energie
+cdhmax = 2.0E-3
+#
+# Parametres "orbitaux/ ere geologique"
+#
+### R_ecc =      Excentricite
+R_ecc = 0.016715
+### R_peri =     Equinoxe
+R_peri = 102.7
+### R_incl =     Inclinaison
+R_incl = 23.441
+### solaire =    Constante solaire
+solaire = 1365.
+#
+# Taux gaz a effet de serre
+#
+### co2_ppm =    taux CO2 en ppm
+co2_ppm = 348.
+### CH4_ppb =    taux CH4 en ppb
+CH4_ppb = 1650.
+### N2O_ppb =    taux N2O en ppb
+N2O_ppb = 306.
+### CFC11_ppt =  taux CFC11 en ppt
+CFC11_ppt = 280.
+### CFC12_ppt =  taux CFC12 en ppt
+CFC12_ppt = 484
+#
+# Parametres effets directs/indirects des "aerosols"
+#
+### ok_ade=y/n   flag Aerosol direct effect
+ok_ade=n
+### ok_aie=y/n   flag Aerosol indirect effect
+ok_aie=n
+### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
+# variable modife avec lmdz.card
+aer_type=xxx
+###  type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
+flag_aerosol=1 
+### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b0=1.7
+### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b1=0.2
+#
+# Parametre de lecture de l'ozone
+#
+# Allowed values are 0, 1 and 2
+# 0: do not read an ozone climatology
+# 1: read a single ozone climatology that will be used day and night
+# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
+read_climoz=0
+#
+# Parametres simulateur COSP (CFMIP Observational Simulator Package)
+#
+### ok_cosp=y/n flag simulateur COSP
+ok_cosp=n
+### freq_COSP = frequence d'appel de COSP en secondes
+freq_COSP=10800.
+### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
+ok_mensuelCOSP=y
+### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
+ok_journeCOSP=y
+### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
+ok_hfCOSP=n
+#
+# Parametres simulateur ISCCP
+#
+### ok_isccp=y/n flag simulateur ISCCP
+ok_isccp=n
+### top_height = flag choix calcul nuages par le simulateur en utilisant 
+# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
+# - top_height = 1 -> algo IR-VIS
+# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
+# - top_height = 3 -> algo IR
+top_height = 3
+### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
+# - overlap=1    Max overlap
+# - overlap=2    Random overlap
+# - overlap=3    Max/Random overlap
+overlap = 3
+config_inca=xxxx
+aerosol_couple=n
+type_trac=inca
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L39_NPv3.1
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L39_NPv3.1	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/physiq.def_L39_NPv3.1	(revision 1738)
@@ -0,0 +1,254 @@
+##############################################################################
+# PARAMETRES ET CLES DE CONFIGURATION LMDZ 
+# Cette parie sera bientot extrait dans un fichier config.def
+##############################################################################
+# D: Valeur par default
+#
+# Controle des sorties
+#*********************
+#ok_regdyn : y/n calcul/non des regymes dynamiques sur regions pre-definies
+ok_regdyn=y
+# A transformer en la ligne ci-dessous quand ce sera compatible avec libIGCM
+phys_out_filekeys=       y       y       y       y       y           n
+phys_out_filenames=      histmth histday histhf  histhf3h histhf3hm  histstn
+phys_out_filelevels=     5       5       5       5       5            5
+phys_out_filetypes=      ave(X)  ave(X)  ave(X) inst(X) inst(X)
+phys_out_filetimesteps=  1.mth,  1.day,  0.25day, 0.125day, 0.125day, 1800.s
+### parametres pour hist*NMC.nc
+#lev_histdayNMC: nombre de niveaux sur la verticale dans le fichier
+#                histdayNMC: 8 ou 17
+lev_histdayNMC=8
+### - ok_histNMC = y/n, y/n, y/n pour avoir les sorties NMC mensuelles, journalieres et haute frequence
+ok_histNMC = y, n, n
+#
+#Couplage avec autres modules
+#******************************
+##  Soil Model  ou non               
+soil_model=y
+##  Cycle diurne  ou non                 
+cycle_diurne=y
+
+#Code de transfert radiatif
+#**************************
+# ajout de ce flag pour activer/desactiver le rayonnement (MPL)
+# 0 : pas de rayonnement. 1: on active le rayonnement (D=1)
+iflag_radia=1
+## Nombre appels des routines de rayonnements ( par jour)                 
+nbapp_rad=24
+
+# Parametres "orbitaux/ ere geologique"
+#######################################
+#Parametres orbitaux
+#excentricite
+R_ecc=0.016715
+#equinoxe
+R_peri=102.7
+#inclinaison
+R_incl=23.441
+#constante solaire
+solaire = 1366.0896 
+##  Facteur additif pour l albedo
+pmagic=999.9
+#
+# Taux gaz a effet de serre
+#######################################
+#concentration des gaz (co2 ch4   n2o  cfc11 cfc12)
+#                  (D: 348. 1650. 306.  280.  484.)
+co2_ppm = 0.36886500E+03
+#RCO2=co2_ppm * 1.0e-06  * 44.011/28.97= 5.286789092164308E-04
+CH4_ppb = 0.17510225E+04
+#RCH4=1.65E-06* 16.043/28.97= 9.137366240938903E-07
+N2O_ppb = 0.31585000E+03
+#RN2O=306.E-09* 44.013/28.97= 4.648939592682085E-07
+CFC11_ppt = 5.18015181E+01
+#RCFC11=280.E-12* 137.3686/28.97= 1.327690990680013E-09
+CFC12_ppt = 0.99862742E+03
+#RCFC12=484.E-12* 120.9140/28.97= 2.020102726958923E-09
+#
+# Parametres effets directs/indirects des "aerosols"
+#######################################
+### ok_ade=y/n   flag Aerosol direct effect
+ok_ade=n
+### ok_aie=y/n   flag Aerosol indirect effect
+ok_aie=n
+### aer_type =   Aerosol variation type : actuel / preind / scenario / annuel
+aer_type=xxx
+###  type of coupled aerosol =1 (default) =2 => bc  only =3 => pom only =4 => seasalt only =5 => dust only =6 => all aerosol
+flag_aerosol=6 
+### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b0=1.7
+### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b1=0.2
+#
+# Parametre de lecture de l'ozone
+#######################################
+# Allowed values are 0, 1 and 2
+# 0: do not read an ozone climatology
+# 1: read a single ozone climatology that will be used day and night
+# 2: read two ozone climatologies, the average day and night climatology and the daylight climatology
+read_climoz=0
+#
+# Parametres simulateur COSP (CFMIP Observational Simulator Package)
+#######################################
+### ok_cosp=y/n flag simulateur COSP
+ok_cosp=n
+### freq_COSP = frequence d'appel de COSP en secondes
+freq_COSP=10800.
+### ok_mensuelCOSP=y/n sortir fichier mensuel COSP histmthCOSP.nc, =n pas de fichier histmthCOSP.nc
+ok_mensuelCOSP=y
+### ok_journeCOSP=y/n sortir fichier journalier histdayCOSP.nc, =n pas de fichier histdayCOSP.nc
+ok_journeCOSP=y
+### ok_hfCOSP=y/n, ecrire sorties  haute frequence histhfCOSP.nc, =n pas de fichier histhfCOSP.nc
+ok_hfCOSP=n
+#
+# Parametres simulateur ISCCP
+#######################################
+### ok_isccp=y/n flag simulateur ISCCP
+ok_isccp=n
+### freqin_isccp = frequence input en secondes du simulateur ISCCP
+freq_ISCCP=10800.
+### top_height = flag choix calcul nuages par le simulateur en utilisant 
+# -              les donnees IR et/ou VIS et l algorithme ISCCP-D1
+# - top_height = 1 -> algo IR-VIS
+# - top_height = 2 -> identique a 1, mais "ptop(ibox)=pfull(ilev)"
+# - top_height = 3 -> algo IR
+top_height = 1
+### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
+# - overlap=1    Max overlap
+# - overlap=2    Random overlap
+# - overlap=3    Max/Random overlap
+overlap = 3
+#
+
+# Specifique aquaplanete
+########################
+#qsol0=100.
+##############################################################################
+# PARAMETRES ET CLES DE LA PHYSIQUE LMDZ
+# Cette partie constitue la partie : physiq.def_L39_NPv3.0
+##############################################################################
+# D: Valeur par default
+#
+#Parametres Orographiques et cdrags
+#**********************************
+##   Si=.T. ,  lecture du fichier limit avec la bonne annee 
+ok_limitvrai=n
+#Cdrags 
+f_cdrag_stable=1.
+f_cdrag_ter=1.
+f_cdrag_oce=0.7
+#cdmmax
+cdmmax=2.5E-3
+#cdhmax
+cdhmax=2.0E-3
+##  Orodr  ou  non   pour l orographie
+ok_orodr=y
+##  Orolf  ou  non   pour l orographie              
+ok_orolf=y
+#Rugoro
+f_rugoro=0.
+#
+# Rayonnement
+#*********************
+# activation du nouveau code de rayonnement RRTM
+# 0 : Ancien code et 1 : RRTM (D=0)
+iflag_rrtm=0
+#
+##  Activation ou non de la parametrisation de Hines pour la strato
+#******************************************************************
+ok_hines=y
+#
+# Parametres nuages
+#******************
+#seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
+cld_lc_lsc=6.e-4
+cld_lc_con=6.e-4
+#constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.)
+cld_tau_lsc=1800.
+cld_tau_con=1800.
+#facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
+ffallv_lsc=1.35
+ffallv_con=1.35
+#coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
+coef_eva=1e-4
+reevap_ice=y
+#calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4) 
+iflag_cldcon=6
+fact_cldcon=1.
+facttemps=0.
+#calcul eau condensee et fraction nuageuse a partir des PDF (D:0, 0:version avec ratqs sinon nouvelles PDFs)  
+iflag_pdf=1
+#calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
+ok_newmicro=y
+iflag_ratqs=2
+ratqsbas=0.002
+ratqshaut=0.25
+#rayon effectif pour les nuages de glace (D:35)
+rad_froid=35
+#rayons effectifs pour les nuages eau liq (tailles des gouttes d eau nuageuse, D: 13 9)
+rad_chau1=12
+rad_chau2=11
+##  Choix ou non  de  New oliq               
+new_oliq=y
+#
+rei_min=20.
+rei_max=61.29
+#
+t_glace_min=258.
+t_glace_max=273.13
+#
+# Flags Convection 
+#*****************
+#flag  pour la convection (D:2, 1:LMD, 2:Tiedtke, 3:KE nouvelle physique, 30:KE IPCC)
+iflag_con=3
+#niveau de sorties des diagnos sur la conservation energie
+if_ebil=0
+#efficacite de precepitation maxim (D:.993)
+epmax=0.997
+#ajustement convectif sec au debut de la convection (D: n 0)
+ok_adj_ema=n
+iflag_clw=0
+#fermeture convective (D:1, 0:AR4, 2:ALE and ALP)
+iflag_clos=2
+#loi de melange a l entrainement  (D:1, 0:AR4=PDF plate, 1=PDF)
+iflag_mix=1
+#poids des PDFs plate et en cloche (D: 1 0)
+qqa1=1.
+qqa2=0.
+# Facteur multiplication des precip convectives dans KE
+cvl_corr=1.0
+#
+Fmax=0.65
+#
+# Flags Couche limite
+#********************
+#shema de CL des thermiques (D:0, 0:ajustement sec,=>1:versions thermiques) 
+iflag_thermals=15
+#pas de temps du spliting pour les thermiques
+nsplit_thermals=1
+tau_thermals=0.
+iflag_thermals_ed=10
+iflag_thermals_optflux=0
+#shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada) 
+iflag_pbl=8
+#diffusion turbulente
+ksta_ter=1.e-7
+ksta=1.e-10
+#ok_kzmin : calcul Kzmin dans la CL de surface
+ok_kzmin=y
+#couplage ou non avec la convection
+iflag_coupl=5 (0:AR4 et 1:nouvelle physique)
+seuil_inversion=-0.08
+#rendement sur la puissance fournie par les thermiques a la convection
+alp_bl_k=0.5
+#
+# Flags wakes
+#************
+#activation ou non des wakes (D:0, 0:AR4, 1:nouvelle physique) 
+iflag_wake=1
+alp_offset=0.
+
+nsplit_phys=2
+config_inca=xxxx
+type_trac=inca
+
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/run.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/run.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/run.def	(revision 1738)
@@ -0,0 +1,41 @@
+## $Id: run.def 71 2007-06-21 14:26:03Z acosce $
+INCLUDEDEF=physiq.def
+INCLUDEDEF=gcm.def
+INCLUDEDEF=orchidee.def
+INCLUDEDEF=guide.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_
+## Nombre de jours d'integration
+nday=_nday_
+## Remise a zero de la date initiale
+raz_date=_raz_date_
+## periode de sortie des variables de controle (en pas)
+iconser=480
+## periode d'ecriture du fichier histoire (en jour)
+iecri=1
+## flag de sortie dynzon
+ok_dynzon=n
+## periode de stockage fichier dynzon (en jour)
+periodav=_periodav_
+## activation du calcul d equilibrage de charge
+adjust=_adjust_
+## activation du filtre fft
+use_filtre_fft=y
+## niveau d impression de controle
+prt_level=0
+###
+### Informations sur la configuration utilisee
+###
+### type_ocean = force / slab  /couple
+type_ocean=force
+### VEGET= y si ORCHIDEE, =n si bucket
+VEGET=y
+### Parameters read only by ce0l :
+### create file grilles_gcm.nc
+grilles_gcm_netcdf=y
+
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/traceur.def
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/traceur.def	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/PARAM/traceur.def	(revision 1738)
@@ -0,0 +1,5 @@
+4
+14 14 H2Ov
+10 10 H2Ol
+10 10 RN
+10 10 PB
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/POST/monitoring01_inca.cfg
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/POST/monitoring01_inca.cfg	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/POST/monitoring01_inca.cfg	(revision 1738)
@@ -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/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/config.card
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/config.card	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/EXP_LOI/config.card	(revision 1738)
@@ -0,0 +1,144 @@
+#- $Id: config.card 69 2007-06-21 14:08:22Z acosce $
+#========================================================================
+#D-- Carte pour run AR5WERF -- 
+#D-- Les fichiers sflx - aircraft - start - startphy créés par S. Szopa
+#D--
+#========================================================================
+#D-- Compatibility -
+[Compatibility]
+libIGCM=1.0
+#========================================================================
+#D-- UserChoices -
+[UserChoices]
+#============================
+#D-- (<8 chars MAX for JobName)
+JobName=LOI3_AR5WERF
+LongName='Ajout de Stomate (NEW) - Fichiers sflx / aircraft / start / startphy  crees par S.Szopa '
+TagName=LMDZORINCA/${CHEM}
+#============================
+#D-- leap, noleap, 360d
+CalendarType=leap
+#D-- Début et fin de Job
+#D-- "YYYY-MM-DD"
+DateBegin=2006-01-01
+DateEnd=2006-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=32
+#============================
+#-- Run Options : here for oasis coupling
+#-- For Vargas
+#JobRunOptions=JobRunOptions=''
+#-- For Titane 
+JobRunOptions='"-np "${config_UserChoices_JobNumProcTot}'
+#========================================================================
+#D-- ListOfComponents -
+[ListOfComponents]
+#D- For each component, Name of component, Tag of component
+ATM= (lmdz, LMDZ4)
+SRF= (orchidee, orchidee_1_9_5)
+SBG= (stomate, ORCHIDEE_1_9_5)
+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)
+SRF= ("", "")
+SBG= ("", "")
+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=2005-12-31
+#D- Define restart simulation name (=> JOB_OS)
+RestartJobName=LOI_AR5WERF_E0_2006
+#D- Path Server Group Login (=> PSGL)
+RestartPath=/dmnfs/cont003/p86cozic/IGCM_OUT/LMDZORINCA/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
+RebuildFrequency=1Y
+#D- Do we rebuild parallel output from archive
+RebuildFromArchive=NONE
+#D- If you want to produce time series, this flag determines
+#D- frequency of post-processing submission
+TimeSeriesFrequency=NONE
+#D- If you want to produce seasonal average, this flag determines
+#D- the period of this average
+SeasonalFrequency=NONE
+
+#========================================================================
+#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/LMDZORINCA/NMHC_AER
+# Old component name for restart (if empty, use new name)
+OldName=""
+
+
+#========================================================================
+#D-- SRF -
+[SRF]
+WriteFrequency="1M"
+Restart=n
+##-- Last day of the experience used as restart
+RestartDate=1999-12-30
+# Define restart simulation name
+RestartJobName=JobName
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZORINCA/NMHC_AER
+#-- Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
+#D-- SBG - STOMATE
+[SBG]
+WriteFrequency="1M"
+# If config_Restarts_OverRule == 'n' next 4 params are read
+Restart=n
+#-- Last day of the experience used as restart for this component if Restart=y
+RestartDate=2000-01-31
+# Define restart simulation name for this component
+RestartJobName=EXP00
+RestartPath=${ARCHIVE}/IGCM_OUT/LMDZOR/DEVT/amip
+#-- 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/LMDZORINCA/NMHC_AER
+#D-- Old component name for restart (if empty, use new name)
+OldName=""
+
+#========================================================================
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/scripts/BB_make
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/scripts/BB_make	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/scripts/BB_make	(revision 1738)
@@ -0,0 +1,447 @@
+# $Id$
+# some variables
+TMP = ../../../tmp
+L_X = $(L_O) $(L_C)
+#-Q- sxnec  # Compiler options for the ORCA model
+#-Q- sxnec  F_O = -dw -Wf\"-A dbl4\" -sx5 -C vsafe -Wf"-P nh" -Wf,-pvctl noassume loopcnt=10000 -L transform -I $(MODDIR) -I../../../prism/SX/build/lib/psmile.$(LIB_MPI)
+#-Q- sx6nec  # Compiler options for the ORCA model
+#-Q- sx6nec  F_O = -size_t64 -dw -Wf\"-A dbl4\" -C vsafe -Wf"-P nh" -Wf,-pvctl noassume loopcnt=10000 -L transform -I $(MODDIR) -I../../../prism/SX/build/lib/psmile.$(LIB_MPI)
+#-Q- sx8brodie  # Compiler options for the ORCA model
+#-Q- sx8brodie  F_O = -dw -Wf\"-A dbl4\" -C vsafe -Wf"-P nh" -Wf,-pvctl noassume loopcnt=10000 -L transform -I $(MODDIR) -I../../../prism/SX/build/lib/psmile.$(LIB_MPI)
+#-Q- linux  # Compiler options for the ORCA model
+#-Q- linux  F_O = -Mlist -O2 -byteswapio -r8  -module $(MODDIR)
+PREF=
+#-Q- sxnec	PREF=i.
+#-Q- sx6nec	PREF=i.
+#-Q- sx8brodie	PREF=i.
+
+KEY_LISTE=key_antarctic key_arctic key_bimg key_bulk key_convevd key_convnpc key_coupled key_coupled_ipslcm4 key_coupled_prism key_ddmixing key_debug key_diaeiv key_diagap key_diahdy key_diahth key_diainstant key_diaspr key_diatrddyn key_diatrdmld key_diatrdtra key_diaznl key_diaznl_brokenline key_dynhdfbilap key_dynhdfcoef1d key_dynhdfcoef2d key_dynhdfcoef3d key_dynhdfgeop key_dynhdfiso key_eastobc key_eel2 key_eel6 key_fdir key_feti key_flx key_flx_clio key_flx_lmd key_flxqsr key_freesurf_cstvol key_hpgimplicit key_ice key_ice_igloo key_ice_lln key_icedta key_islands key_lim_age key_lim_cp1 key_lim_cp2 key_lim_cp3 key_lim_fdd key_monotasking key_mpp key_mpp_mpi key_mpp_pvm key_mpp_shmem key_northobc key_noslip_accurate key_openmp key_orca_r025 key_orca_r05 key_orca_r2 key_orca_r4 key_passivetrc key_s_coord key_saldta key_saldta_monthly key_southobc key_sst key_tau key_tau_clio key_tau_ers_ncep key_tau_lmd key_temdta key_tradmp key_trahdfbilap key_trahdfcoef1d key_trahdfcoef2d key_trahdfcoef3d key_trahdfeiv key_trahdfgeop key_trahdfiso key_trddyn_new key_vairmer key_vorcombined key_vorenergy key_vpp key_westobc key_zdfconstant key_zdfexplicit key_zdfrichardson key_zdftke 
+
+# .o file list in the library
+LIB_OBJ = \
+		$(MODEL_LIB)(Fetilib.o)\
+		$(MODEL_LIB)(Fioopa2.o)\
+		$(MODEL_LIB)(allcray.o)\
+		$(MODEL_LIB)(alllib.o)\
+		$(MODEL_LIB)(ctlopn.o)\
+		$(MODEL_LIB)(prihin.o)\
+		$(MODEL_LIB)(prihre.o)\
+		$(MODEL_LIB)(prizre.o)\
+		$(MODEL_LIB)(Fioopa.o)\
+		$(MODEL_LIB)(Fmpplib.o)\
+		$(MODEL_LIB)(acrlbq.o)\
+		$(MODEL_LIB)(adv.o)\
+		$(MODEL_LIB)(advx.o)\
+		$(MODEL_LIB)(advy.o)\
+		$(MODEL_LIB)(angle.o)\
+		$(MODEL_LIB)(bn2.o)\
+		$(MODEL_LIB)(closea.o)\
+		$(MODEL_LIB)(cofdis.o)\
+		$(MODEL_LIB)(cur.o)\
+		$(MODEL_LIB)(day.o)\
+		$(MODEL_LIB)(defcst.o)\
+		$(MODEL_LIB)(defgrid.o)\
+		$(MODEL_LIB)(diabort.o)\
+		$(MODEL_LIB)(diadyn.o)\
+		$(MODEL_LIB)(diafwb.o)\
+		$(MODEL_LIB)(diagap.o)\
+		$(MODEL_LIB)(diahdy.o)\
+		$(MODEL_LIB)(diahf.o)\
+		$(MODEL_LIB)(diahth.o)\
+		$(MODEL_LIB)(diamld.o)\
+		$(MODEL_LIB)(diamxl.o)\
+		$(MODEL_LIB)(dianam.o)\
+		$(MODEL_LIB)(diaspr.o)\
+		$(MODEL_LIB)(diatra.o)\
+		$(MODEL_LIB)(diawri.o)\
+		$(MODEL_LIB)(diaznl.o)\
+		$(MODEL_LIB)(diffus.o)\
+		$(MODEL_LIB)(div.o)\
+		$(MODEL_LIB)(domclo.o)\
+		$(MODEL_LIB)(domhgr.o)\
+		$(MODEL_LIB)(dommba.o)\
+		$(MODEL_LIB)(dommsk.o)\
+		$(MODEL_LIB)(domstp.o)\
+		$(MODEL_LIB)(domwri.o)\
+		$(MODEL_LIB)(domzgr.o)\
+		$(MODEL_LIB)(dtacof.o)\
+		$(MODEL_LIB)(dtaobc.o)\
+		$(MODEL_LIB)(dtasal.o)\
+		$(MODEL_LIB)(dtasst.o)\
+		$(MODEL_LIB)(dtatem.o)\
+		$(MODEL_LIB)(dtrlec.o)\
+		$(MODEL_LIB)(dtrsal.o)\
+		$(MODEL_LIB)(dtrtem.o)\
+		$(MODEL_LIB)(dtruvg.o)\
+		$(MODEL_LIB)(dynami.o)\
+		$(MODEL_LIB)(dynhdf.o)\
+		$(MODEL_LIB)(dynhpg.o)\
+		$(MODEL_LIB)(dynkeg.o)\
+		$(MODEL_LIB)(dynnxt.o)\
+		$(MODEL_LIB)(dynspg.o)\
+		$(MODEL_LIB)(dynvor.o)\
+		$(MODEL_LIB)(dynzad.o)\
+		$(MODEL_LIB)(dynzdf.o)\
+		$(MODEL_LIB)(eos.o)\
+		$(MODEL_LIB)(fetmat.o)\
+		$(MODEL_LIB)(fetsch.o)\
+		$(MODEL_LIB)(fetstr.o)\
+		$(MODEL_LIB)(flx.o)\
+		$(MODEL_LIB)(flxrnf.o)\
+		$(MODEL_LIB)(fontbc.o)\
+		$(MODEL_LIB)(forcat.o)\
+		$(MODEL_LIB)(forcng.o)\
+		$(MODEL_LIB)(forkjoin.o)\
+		$(MODEL_LIB)(foroutp.o)\
+		$(MODEL_LIB)(gather.o)\
+		$(MODEL_LIB)(geo2oce.o)\
+		$(MODEL_LIB)(hdfeiv.o)\
+		$(MODEL_LIB)(hdfght.o)\
+		$(MODEL_LIB)(hdfguv.o)\
+		$(MODEL_LIB)(hdfslp.o)\
+		$(MODEL_LIB)(hgrcoo.o)\
+		$(MODEL_LIB)(icdadv.o)\
+		$(MODEL_LIB)(icdyna.o)\
+		$(MODEL_LIB)(icelln.o)\
+		$(MODEL_LIB)(icemodel.o)\
+		$(MODEL_LIB)(icetooce.o)\
+		$(MODEL_LIB)(informe.o)\
+		$(MODEL_LIB)(inicmo.o)\
+		$(MODEL_LIB)(inidom.o)\
+		$(MODEL_LIB)(inidta.o)\
+		$(MODEL_LIB)(inidtr.o)\
+		$(MODEL_LIB)(inihdf.o)\
+		$(MODEL_LIB)(inimpp.o)\
+		$(MODEL_LIB)(inimpp2.o)\
+		$(MODEL_LIB)(iniobc.o)\
+		$(MODEL_LIB)(iniout.o)\
+		$(MODEL_LIB)(inipar.o)\
+		$(MODEL_LIB)(inisol.o)\
+		$(MODEL_LIB)(initice.o)\
+		$(MODEL_LIB)(inizdf.o)\
+		$(MODEL_LIB)(islbsf.o)\
+		$(MODEL_LIB)(isldom.o)\
+		$(MODEL_LIB)(islmat.o)\
+		$(MODEL_LIB)(islpri.o)\
+		$(MODEL_LIB)(islpth.o)\
+		$(MODEL_LIB)(lbc.o)\
+		$(MODEL_LIB)(lbc2.o)\
+		$(MODEL_LIB)(lbcice.o)\
+		$(MODEL_LIB)(mlbxinit.o)\
+		$(MODEL_LIB)(mpplnkice.o)\
+		$(MODEL_LIB)(obc.o)\
+		$(MODEL_LIB)(obcdom.o)\
+		$(MODEL_LIB)(obcfld.o)\
+		$(MODEL_LIB)(obclec.o)\
+		$(MODEL_LIB)(obcspg.o)\
+		$(MODEL_LIB)(obcwri.o)\
+		$(MODEL_LIB)(ocesla.o)\
+		$(MODEL_LIB)(ocetoice.o)\
+		$(MODEL_LIB)(ocfzpt.o)\
+		$(MODEL_LIB)(opa.o)\
+		$(MODEL_LIB)(outave.o)\
+		$(MODEL_LIB)(parcst.o)\
+		$(MODEL_LIB)(parctl.o)\
+		$(MODEL_LIB)(parlec.o)\
+		$(MODEL_LIB)(pdecli.o)\
+		$(MODEL_LIB)(redrunb.o)\
+		$(MODEL_LIB)(redrunnc.o)\
+		$(MODEL_LIB)(redtab.o)\
+		$(MODEL_LIB)(repcmo.o)\
+		$(MODEL_LIB)(repere.o)\
+		$(MODEL_LIB)(rstwri.o)\
+		$(MODEL_LIB)(savrunb.o)\
+		$(MODEL_LIB)(savrunnc.o)\
+		$(MODEL_LIB)(savtab.o)\
+		$(MODEL_LIB)(scatter.o)\
+		$(MODEL_LIB)(shine.o)\
+		$(MODEL_LIB)(solfet.o)\
+		$(MODEL_LIB)(solmat.o)\
+		$(MODEL_LIB)(solpcg.o)\
+		$(MODEL_LIB)(solrbs.o)\
+		$(MODEL_LIB)(solsor.o)\
+		$(MODEL_LIB)(sprmat.o)\
+		$(MODEL_LIB)(staocb.o)\
+		$(MODEL_LIB)(start.o)\
+		$(MODEL_LIB)(step.o)\
+		$(MODEL_LIB)(stpcmo.o)\
+		$(MODEL_LIB)(stpctl.o)\
+		$(MODEL_LIB)(tau.o)\
+		$(MODEL_LIB)(thersf.o)\
+		$(MODEL_LIB)(trabbc.o)\
+		$(MODEL_LIB)(trabbl.o)\
+		$(MODEL_LIB)(tradmp.o)\
+		$(MODEL_LIB)(trahad.o)\
+		$(MODEL_LIB)(trahdf.o)\
+		$(MODEL_LIB)(tranpc.o)\
+		$(MODEL_LIB)(tranxt.o)\
+		$(MODEL_LIB)(traqsr.o)\
+		$(MODEL_LIB)(trasbc.o)\
+		$(MODEL_LIB)(trazad.o)\
+		$(MODEL_LIB)(trazdf.o)\
+		$(MODEL_LIB)(wzv.o)\
+		$(MODEL_LIB)(zdfbfr.o)\
+		$(MODEL_LIB)(zdfddm.o)\
+		$(MODEL_LIB)(zdfevd.o)\
+		$(MODEL_LIB)(zdfric.o)\
+		$(MODEL_LIB)(zdftke.o)
+#
+.SUFFIXES: .F
+#
+# target
+all: 		$(EXEC_BIN)
+		@cp SRC_FILE_LIST.temp SRC_FILE_LIST
+		@echo OPA model is OK
+# cleaning objects, libraries and executables
+clean:
+		$(RM) modele.o $(MODEL_LIB) $(SXMODEL_LIB) $(EXEC_BIN)
+		@echo OPA model \(.o libraries and executables are removed\)
+# dependancies
+$(EXEC_BIN) :	libioipsl key keyverif src_file_list modele.o $(MODEL_LIB) $(IOIPSL_LIB)
+		$(F_L) $(L_X) -o $(EXEC_BIN) modele.o $(SXMODEL_LIB) $(OASISMPI2_LIB) $(IOIPSL_LIB) $(NCDF_LIB) $(LIBMPI) $(LIBMP)
+
+libioipsl :
+	(cd ../../IOIPSL/src ; $(M_K) -f Makefile )
+
+# main program
+modele.o :	modele.f
+		$(F_C) $(F_O) -c modele.f
+
+# if the make script is stopped, one doesn't destroy library
+.PRECIOUS : $(MODEL_LIB)
+#-Q- sxnec .PRECIOUS : $(SXMODEL_LIB)
+#-Q- sx6nec .PRECIOUS : $(SXMODEL_LIB)
+#-Q- sx8brodie .PRECIOUS : $(SXMODEL_LIB)
+
+# library dependancies
+$(MODEL_LIB):	$(LIB_OBJ)
+		@echo
+		@echo La bibliotheque est a jour
+		@echo
+
+# Rules
+.DEFAULTS:
+# rules for the .F files
+.F.a :           
+		@echo default rule .F.a used 
+		$(F_C) $(P_P) $(F_O) \
+#-Q- sxnec		-R2 \
+#-Q- sx6nec		-R2 \
+#-Q- sx8brodie		-R2 \
+		$*.F || { mv $(PREF)$*.L $(TMP) ; false ; exit ; }
+#-Q- sxnec		$(A_X) $(SXMODEL_LIB) $*.o
+#-Q- sx6nec		$(A_X) $(SXMODEL_LIB) $*.o
+#-Q- sx8brodie		$(A_X) $(SXMODEL_LIB) $*.o
+		$(A_C) $(MODEL_LIB) $*.o
+		$(RM) $*.[of]
+#-Q- linux		mv $(PREF)$*.lst $(PREF)$*.L
+		mv $(PREF)$*.L $(TMP)
+
+# rules for the library from the .f files 
+.f.a : 		
+		@echo default rule .f.a used 
+		$(F_C) $(F_O) \
+#-Q- sxnec		-R2 \
+#-Q- sx6nec		-R2 \
+#-Q- sx8brodie		-R2 \
+		$*.f|| { mv $*.L $(TMP) ; false ; exit ; }
+#-Q- sxnec		$(A_X) $(SXMODEL_LIB) $*.o
+#-Q- sx6nec		$(A_X) $(SXMODEL_LIB) $*.o
+#-Q- sx8brodie		$(A_X) $(SXMODEL_LIB) $*.o
+		$(A_C) $(MODEL_LIB) $*.o
+		$(RM) -f $*.o
+#-Q- linux	mv $(PREF)$*.lst $(PREF)$*.L
+		mv $*.L $(TMP)
+#
+$(MODEL_LIB)(Fioopa.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP Fioopa.F
+$(MODEL_LIB)(Fmpplib.o) :	common.coupled.h common.h common.ice.h common.mpp.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.mpp.h parameter.obc.h parameter.passivetrc.h KEY_CPP Fmpplib.F
+$(MODEL_LIB)(acrlbq.o) :	bloc.h const.h ice.h para.h thermo.h type.h KEY_CPP acrlbq.F
+$(MODEL_LIB)(adv.o) :	bloc.h const.h dynami.h para.h type.h KEY_CPP adv.F
+$(MODEL_LIB)(advx.o) :	bloc.h const.h dynami.h para.h type.h KEY_CPP advx.F
+$(MODEL_LIB)(advy.o) :	bloc.h const.h dynami.h para.h type.h KEY_CPP advy.F
+$(MODEL_LIB)(angle.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP angle.F
+$(MODEL_LIB)(bn2.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP bn2.F
+$(MODEL_LIB)(closea.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP closea.F
+$(MODEL_LIB)(cofdis.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP cofdis.F
+$(MODEL_LIB)(cur.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP cur.F
+$(MODEL_LIB)(day.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP day.F
+$(MODEL_LIB)(defcst.o) :	bloc.h const.h datadc.h dynami.h ice.h para.h reper.h type.h varno.h KEY_CPP defcst.F
+$(MODEL_LIB)(defgrid.o) :	bloc.h const.h dynami.h ice.h para.h reper.h type.h KEY_CPP defgrid.F
+$(MODEL_LIB)(diabort.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP diabort.F
+$(MODEL_LIB)(diadyn.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP diadyn.F
+$(MODEL_LIB)(diafwb.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP diafwb.F
+$(MODEL_LIB)(diagap.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP diagap.F
+$(MODEL_LIB)(diahdy.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP diahdy.F
+$(MODEL_LIB)(diahf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP diahf.F
+$(MODEL_LIB)(diahth.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP diahth.F
+$(MODEL_LIB)(diamld.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h diamld.ncinit.h diamld.ncwrite.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP diamld.F
+$(MODEL_LIB)(diamxl.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP diamxl.F
+$(MODEL_LIB)(dianam.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP dianam.F
+$(MODEL_LIB)(diaspr.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP diaspr.F
+$(MODEL_LIB)(diatra.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP diatra.F
+$(MODEL_LIB)(diawri.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h diawri.fdir.h diawri.vairmer.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP diawri.F
+$(MODEL_LIB)(diaznl.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h diaznl.brokenline.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP diaznl.F
+$(MODEL_LIB)(diffus.o) :	bloc.h const.h dynami.h para.h type.h KEY_CPP diffus.F
+$(MODEL_LIB)(div.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP div.F
+$(MODEL_LIB)(domclo.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP domclo.F
+$(MODEL_LIB)(domhgr.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP domhgr.F
+$(MODEL_LIB)(dommba.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dommba.F
+$(MODEL_LIB)(dommsk.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dommsk.F
+$(MODEL_LIB)(domstp.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP domstp.F
+$(MODEL_LIB)(domwri.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h domwri.fdir.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP domwri.F
+$(MODEL_LIB)(domzgr.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h domzgr.s.h domzgr.z.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP domzgr.F
+$(MODEL_LIB)(dtacof.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP dtacof.F
+$(MODEL_LIB)(dtaobc.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dtaobc.F
+$(MODEL_LIB)(dtasal.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h dtasal.fdir.h dtasal.monthly.fdir.h dtasal.monthly.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP dtasal.F
+$(MODEL_LIB)(dtasst.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h dtasst.forced.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dtasst.F
+$(MODEL_LIB)(dtatem.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h dtatem.fdir.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP dtatem.F
+$(MODEL_LIB)(dtrlec.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h dtrlec.fdir.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP dtrlec.F
+$(MODEL_LIB)(dtrsal.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dtrsal.F
+$(MODEL_LIB)(dtrtem.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP dtrtem.F
+$(MODEL_LIB)(dtruvg.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dtruvg.F
+$(MODEL_LIB)(dynami.o) :	bloc.h const.h dynami.h ice.h para.h type.h KEY_CPP dynami.F
+$(MODEL_LIB)(dynhdf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h dynhdf.bilapgeopot.h dynhdf.bilaplacian.h dynhdf.isopycnal.h dynhdf.laplacian.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dynhdf.F
+$(MODEL_LIB)(dynhpg.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dynhpg.F
+$(MODEL_LIB)(dynkeg.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dynkeg.F
+$(MODEL_LIB)(dynnxt.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dynnxt.F
+$(MODEL_LIB)(dynspg.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h dynspg.freesurf.cstvol.h dynspg.rigidlid.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dynspg.F
+$(MODEL_LIB)(dynvor.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h dynvor.combined.h dynvor.energy.h dynvor.enstrophy.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dynvor.F
+$(MODEL_LIB)(dynzad.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP dynzad.F
+$(MODEL_LIB)(dynzdf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h dynzdf.explicit.h dynzdf.implicit.h dynzdf.isopycnal.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h zdf.matrixsolver.h KEY_CPP dynzdf.F
+$(MODEL_LIB)(eos.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP eos.F
+$(MODEL_LIB)(fetmat.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP fetmat.F
+$(MODEL_LIB)(fetsch.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP fetsch.F
+$(MODEL_LIB)(fetstr.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP fetstr.F
+$(MODEL_LIB)(flx.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h flx.bulk.h flx.coupled.clio.h flx.coupled.h flx.forced.clio.fdir.h flx.forced.clio.h flx.forced.h flx.forced.lmd.h flx.sensitivity.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP flx.F
+$(MODEL_LIB)(flxrnf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h flxrnf.orca.r05.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP flxrnf.F
+$(MODEL_LIB)(fontbc.o) :	bloc.h const.h ice.h para.h thermo.h type.h KEY_CPP fontbc.F
+$(MODEL_LIB)(forcat.o) :	bloc.h const.h ice.h para.h type.h KEY_CPP forcat.F
+$(MODEL_LIB)(forcng.o) :	bloc.h const.h dynami.h ice.h para.h type.h KEY_CPP forcng.F
+$(MODEL_LIB)(forkjoin.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP forkjoin.F
+$(MODEL_LIB)(foroutp.o) :	bloc.h const.h datadc.h para.h type.h varno.h KEY_CPP foroutp.F
+$(MODEL_LIB)(gather.o) :	type.h KEY_CPP gather.F
+$(MODEL_LIB)(geo2oce.o) :	parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP geo2oce.F
+$(MODEL_LIB)(hdfeiv.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP hdfeiv.F
+$(MODEL_LIB)(hdfght.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP hdfght.F
+$(MODEL_LIB)(hdfguv.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP hdfguv.F
+$(MODEL_LIB)(hdfslp.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP hdfslp.F
+$(MODEL_LIB)(hgrcoo.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h hgrcoo.fdir.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP hgrcoo.F
+$(MODEL_LIB)(icdadv.o) :	bloc.h const.h dynami.h ice.h moment.h para.h type.h KEY_CPP icdadv.F
+$(MODEL_LIB)(icdyna.o) :	bloc.h const.h dynami.h ice.h para.h type.h KEY_CPP icdyna.F
+$(MODEL_LIB)(icelln.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h icelln.coupled.h icelln.lmd.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP icelln.F
+$(MODEL_LIB)(icemodel.o) :	bloc.h const.h dynami.h ice.h para.h type.h KEY_CPP icemodel.F
+$(MODEL_LIB)(icetooce.o) :	bloc.h const.h dynami.h ice.h para.h type.h KEY_CPP icetooce.F
+$(MODEL_LIB)(informe.o) :	bloc.h const.h dynami.h ice.h para.h reper.h type.h KEY_CPP informe.F
+$(MODEL_LIB)(inicmo.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP inicmo.F
+$(MODEL_LIB)(inidom.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP inidom.F
+$(MODEL_LIB)(inidta.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP inidta.F
+$(MODEL_LIB)(inidtr.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP inidtr.F
+$(MODEL_LIB)(inihdf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h inihdf.antarctic.h inihdf.arctic.h inihdf.dyn.coef1d.h inihdf.dyn.coef2d.h inihdf.dyn.coef3d.h inihdf.tra.coef1d.h inihdf.tra.coef2d.h inihdf.tra.coef3d.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP inihdf.F
+$(MODEL_LIB)(inimpp.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h ../../../lib/libioipsl.a KEY_CPP inimpp.F
+$(MODEL_LIB)(inimpp2.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h ../../../lib/libioipsl.a KEY_CPP inimpp2.F
+$(MODEL_LIB)(iniobc.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP iniobc.F
+$(MODEL_LIB)(iniout.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP iniout.F
+$(MODEL_LIB)(inipar.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP inipar.F
+$(MODEL_LIB)(inisol.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP inisol.F
+$(MODEL_LIB)(initice.o) :	bloc.h const.h dynami.h ice.h para.h type.h KEY_CPP initice.F
+$(MODEL_LIB)(inizdf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP inizdf.F
+$(MODEL_LIB)(islbsf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h islbsf.fdir.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP islbsf.F
+$(MODEL_LIB)(isldom.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP isldom.F
+$(MODEL_LIB)(islmat.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h islmat.fdir.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP islmat.F
+$(MODEL_LIB)(islpri.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP islpri.F
+$(MODEL_LIB)(islpth.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP islpth.F
+$(MODEL_LIB)(lbc.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP lbc.F
+$(MODEL_LIB)(lbc2.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP lbc2.F
+$(MODEL_LIB)(lbcice.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP lbcice.F
+$(MODEL_LIB)(mlbxinit.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP mlbxinit.F
+$(MODEL_LIB)(mpplnkice.o) :	common.coupled.h common.h common.ice.h common.mpp.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.mpp.h parameter.obc.h parameter.passivetrc.h KEY_CPP mpplnkice.F
+$(MODEL_LIB)(obc.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h obc.east.h obc.north.h obc.south.h obc.west.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP obc.F
+$(MODEL_LIB)(obcdom.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP obcdom.F
+$(MODEL_LIB)(obcfld.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP obcfld.F
+$(MODEL_LIB)(obclec.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP obclec.F
+$(MODEL_LIB)(obcspg.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h obcspg.east.h obcspg.north.h obcspg.south.h obcspg.west.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP obcspg.F
+$(MODEL_LIB)(obcwri.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP obcwri.F
+$(MODEL_LIB)(ocesla.o) :	bloc.h const.h para.h type.h KEY_CPP ocesla.F
+$(MODEL_LIB)(ocetoice.o) :	bloc.h const.h dynami.h ice.h ocetoice.coupled.h para.h type.h KEY_CPP ocetoice.F
+$(MODEL_LIB)(ocfzpt.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP ocfzpt.F
+$(MODEL_LIB)(opa.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP opa.F
+$(MODEL_LIB)(outave.o) :	bloc.h const.h dynami.h ice.h para.h thermo.h type.h ../../../lib/libioipsl.a KEY_CPP outave.F
+$(MODEL_LIB)(parcst.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP parcst.F
+$(MODEL_LIB)(parctl.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP parctl.F
+$(MODEL_LIB)(parlec.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP parlec.F
+$(MODEL_LIB)(pdecli.o) :	const.h type.h KEY_CPP pdecli.F
+$(MODEL_LIB)(redrunb.o) :	bloc.h const.h datadc.h dynami.h ice.h moment.h para.h type.h vareq.h varno.h KEY_CPP redrunb.F
+$(MODEL_LIB)(redrunnc.o) :	bloc.h const.h datadc.h dynami.h ice.h moment.h para.h type.h vareq.h varno.h ../../../lib/libioipsl.a KEY_CPP redrunnc.F
+$(MODEL_LIB)(redtab.o) :	type.h KEY_CPP redtab.F
+$(MODEL_LIB)(repcmo.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h KEY_CPP repcmo.F
+$(MODEL_LIB)(repere.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP repere.F
+$(MODEL_LIB)(rstwri.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h rstwri.fdir.h stafun.h ../../../lib/libioipsl.a KEY_CPP rstwri.F
+$(MODEL_LIB)(savrunb.o) :	bloc.h const.h datadc.h dynami.h ice.h moment.h para.h type.h vareq.h varno.h KEY_CPP savrunb.F
+$(MODEL_LIB)(savrunnc.o) :	bloc.h const.h datadc.h dynami.h ice.h moment.h para.h type.h vareq.h varno.h ../../../lib/libioipsl.a KEY_CPP savrunnc.F
+$(MODEL_LIB)(savtab.o) :	type.h KEY_CPP savtab.F
+$(MODEL_LIB)(scatter.o) :	type.h KEY_CPP scatter.F
+$(MODEL_LIB)(shine.o) :	type.h KEY_CPP shine.F
+$(MODEL_LIB)(solfet.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP solfet.F
+$(MODEL_LIB)(solmat.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP solmat.F
+$(MODEL_LIB)(solpcg.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP solpcg.F
+$(MODEL_LIB)(solrbs.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP solrbs.F
+$(MODEL_LIB)(solsor.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP solsor.F
+$(MODEL_LIB)(sprmat.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP sprmat.F
+$(MODEL_LIB)(staocb.o) :	bloc.h const.h dynami.h ice.h moment.h para.h type.h KEY_CPP staocb.F
+$(MODEL_LIB)(start.o) :	bloc.h const.h ice.h para.h type.h KEY_CPP start.F
+$(MODEL_LIB)(step.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP step.F
+$(MODEL_LIB)(stpcmo.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h ../../../lib/libioipsl.a KEY_CPP stpcmo.F
+$(MODEL_LIB)(stpctl.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP stpctl.F
+$(MODEL_LIB)(tau.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h tau.coupled.clio.h tau.coupled.h tau.forced.clio.fdir.h tau.forced.clio.h tau.forced.ers.ncep.h tau.forced.h tau.forced.lmd.h KEY_CPP tau.F
+$(MODEL_LIB)(thersf.o) :	bloc.h const.h dynami.h ice.h para.h thermo.h type.h KEY_CPP thersf.F
+$(MODEL_LIB)(trabbc.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP trabbc.F
+$(MODEL_LIB)(trabbl.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP trabbl.F
+$(MODEL_LIB)(tradmp.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP tradmp.F
+$(MODEL_LIB)(trahad.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP trahad.F
+$(MODEL_LIB)(trahdf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h diaznl.isopycnal.brokenline.h diaznl.isopycnal.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h trahdf.bilapgeopot.h trahdf.bilaplacian.h trahdf.isopycnal.h trahdf.laplacian.h KEY_CPP trahdf.F
+$(MODEL_LIB)(tranpc.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP tranpc.F
+$(MODEL_LIB)(tranxt.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP tranxt.F
+$(MODEL_LIB)(traqsr.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP traqsr.F
+$(MODEL_LIB)(trasbc.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP trasbc.F
+$(MODEL_LIB)(trazad.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP trazad.F
+$(MODEL_LIB)(trazdf.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h trazdf.explicit.h trazdf.implicit.h trazdf.isopycnal.h zdf.matrixsolver.h KEY_CPP trazdf.F
+$(MODEL_LIB)(wzv.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP wzv.F
+$(MODEL_LIB)(zdfbfr.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP zdfbfr.F
+$(MODEL_LIB)(zdfddm.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP zdfddm.F
+$(MODEL_LIB)(zdfevd.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP zdfevd.F
+$(MODEL_LIB)(zdfric.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h KEY_CPP zdfric.F
+$(MODEL_LIB)(zdftke.o) :	common.coupled.h common.h common.ice.h common.obc.h common.passivetrc.h parameter.coupled.h parameter.h parameter.ice.h parameter.obc.h parameter.passivetrc.h stafun.h zdf.matrixsolver.h KEY_CPP zdftke.F
+# key changing control
+key: 
+	@echo CHANGEMENT de CLES CPP oui/non ?
+	-@if [ "`cat KEY_CPP`" != "$(P_P)" ] ; then \
+		echo CPP options changed ; echo "$(P_P)" > KEY_CPP ; fi 
+
+keyverif:
+	@echo VERIFICATION DES CLES
+	@echo CLES UTILISEES  : 
+	@for v in ` echo $(P_P) | sed  -e s/-Wp,-D//g` ; \
+	do \
+	echo $$v ;\
+        w=`echo $$v | cut -c1-3` ; \
+        if [ $$w = `echo "key"` ] ; then \
+        r=n ; \
+        echo " $(KEY_LISTE) " | (grep $$v >/dev/null 2>&1) || r=o ; \
+           if [ $$r != "n" ] ; then \
+              echo "UNKNOWN KEY" ;\
+              false ;  \
+              exit ;  \
+           fi \
+        fi \
+	done
+# Checking the number of routines
+src_file_list: 
+	@echo CHECKING THE NUMBER AND NAMES OF SOURCE FILES 
+	@\$(RM) `ls | fgrep -v -f .patron` 
+	@if [ -d ../SRC_ORCA ] ; then ln -sf ../SRC_ORCA/*.[Ffh] . ; fi 
+	@if [ -d ../SRC_UCL ] ; then  ln -sf ../SRC_UCL/*[Ffh] . ; fi 
+	@if [ -d ../SRC_ORCA ] && [ -d ../SRC_TRC ] && [ ! -d ../SRC_OFFLINE_TRC ] ; then $(RM) *.passivetrc.h ; ln -sf ../SRC_ORCA/*.[Ffh] . ; ln -sf ../SRC_TRC/*.[Ffh] ../SRC_TRC/trc_sms/*.[Ffh] . ; fi 
+	@if [ -d ../SRC_OFFLINE_TRC ] ; then  ln -sf ../SRC_TRC/*.[Ffh] ../SRC_TRC/trc_sms/*.[Ffh] . ;ln -sf ../SRC_OFFLINE_TRC/*.[Ffh] . ; fi 
+	@$(RM) tmplist
+	@ls -1 | fgrep -v -f .patron  >tmplist
+	@cmp tmplist SRC_FILE_LIST.temp || { echo some routines have been added or suppressed ; sdiff -s tmplist SRC_FILE_LIST.temp ; cp tmplist SRC_FILE_LIST.temp ; $(RM) modele.o $(MODEL_LIB) $(SXMODEL_LIB) $(EXEC_BIN) ; echo Rebuilding the AA_make files ; cd ..  ; ../UTIL/fait_AA_make ; echo   ; echo THE ins_make COMMAND in modipsl/util MUST BE RERUN PLEASE DO SO BEFORE gmake ; echo  ==========================================  ; echo   ;  false ; exit ; }
Index: CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/scripts/BB_make.ldef
===================================================================
--- CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/scripts/BB_make.ldef	(revision 1738)
+++ CONFIG/LMDZORINCA/tags/LMDZORINCA_v4_2/scripts/BB_make.ldef	(revision 1738)
@@ -0,0 +1,21 @@
+# $Id$
+#---------------------------------------------------------------------
+#- OPA SYSTEM : Ocean model + Sea Ice Model
+#---------------------------------------------------------------------
+SHELL = /bin/sh
+#-Q- t3e     #--- Please use gmake instead of make
+#-Q- fjvpp   #--- Please use gmake instead of make
+#---------------------------------------------------------------------
+LIBDIR = ../../../lib
+MODDIR = $(LIBDIR)
+MODEL_LIB = $(LIBDIR)/libopa.a
+IOIPSL_LIB = $(LIBDIR)/libioipsl.a
+SXMODEL_LIB = $(MODEL_LIB)
+#-Q- sxnec  SXMODEL_LIB = $(LIBDIR)/libsxopa.a
+#-Q- sx6nec  SXMODEL_LIB = $(LIBDIR)/libsxopa.a
+#-Q- sx8brodie  SXMODEL_LIB = $(LIBDIR)/libsxopa.a
+OASISMPI2_LIB = 
+LIBMPI  =
+LIBMP   =
+EXEC_BIN = ../../../bin/opa
+#-
