Index: LMDZ4/trunk/000-README
===================================================================
--- LMDZ4/trunk/000-README	(revision 5799)
+++ 	(revision )
@@ -1,20 +1,0 @@
-
-Logiciel LMDZ
--------------
-
-La documentation relative à LMDZ est accessible sur :
-http://lmdz.lmd.jussieu.fr/documentation
-
-Les quides d'installation et utilisation de LMDZ sont accessibles sur :
-http://lmdz.lmd.jussieu.fr/documentation/guides
-
-==========================================================================
-
-LMDZ software
--------------
-
-Documentation about the LMDZ software is available on the web at this address:
-http://lmdz.lmd.jussieu.fr/documentation
-
-Practical installation and user guides are available here:
-http://lmdz.lmd.jussieu.fr/documentation/guides
Index: LMDZ4/trunk/bld.cfg
===================================================================
--- LMDZ4/trunk/bld.cfg	(revision 5799)
+++ 	(revision )
@@ -1,105 +1,0 @@
-# ----------------------- FCM extract configuration file -----------------------
-cfg::type                           bld
-cfg::version                        1.0
-
-
-# ------------------------------------------------------------------------------
-# Build information
-# ------------------------------------------------------------------------------
-
-#Default value of FPP fortran preprocessor
-%FPP cpp
-
-inc arch.fcm
-inc config.fcm
-
-%CONFIG_NAME       %{ARCH}%SUFF_NAME
-%BASE_CONFIG_PATH  %LIBO/%CONFIG_NAME
-%CONFIG_PATH       %BASE_CONFIG_PATH/.config
-%SRC_PATH          %LIBF
-
-%FFLAGS            %BASE_FFLAGS %COMPIL_FFLAGS %PARA_FFLAGS
-%LD_FLAGS          %BASE_LD %PARA_LD
-
-src::dyn     %SRC_PATH/%DYN
-src::phys    %SRC_PATH/%PHYS
-src::grid    %SRC_PATH/grid
-src::filtrez %SRC_PATH/filtrez
-src::bibio   %SRC_PATH/bibio
-src::cosp    %COSP
-
-bld::lib::dyn      %DYN
-bld::lib::phys     %PHYS
-bld::lib::grid      grid
-bld::lib::filtrez   filtrez
-bld::lib::bibio     bibio
-bld::lib::cosp      cosp
-
-
-bld::outfile_ext::exe    %SUFF_NAME.e
-bld::target              lib%{DYN}.a lib%{PHYS}.a libgrid.a libfiltrez.a libbibio.a libcosp.a
-bld::target              %EXEC%SUFF_NAME.e
-bld::exe_dep             %{DYN} %{PHYS} grid filtrez bibio cosp
-
-
-dir::root            %CONFIG_PATH
-dir::lib             %BASE_CONFIG_PATH
-dir::bin             %ROOT_PATH/bin
-
-#search_src           1
-
-bld::tool::fpp       %FPP
-bld::tool::fc        %COMPILER 
-bld::tool::ld        %LINK
-bld::tool::ar        %AR
-bld::tool::make      %MAKE
-bld::tool::fflags    %FFLAGS %INCDIR 
-bld::tool::ldflags   %LD_FLAGS %LIB  
-
-bld::tool::cppflags  %FPP_FLAGS %INCDIR
-bld::tool::fppflags  %FPP_FLAGS %INCDIR
-bld::tool::fppkeys   %CPP_KEY %FPP_DEF
-
-
-#bld::tool::fflags::phys::readaerosol         %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt -pi auto
-#bld::tool::fflags::phys::aeropt_2bands       %BASE_FFLAGS %PROD_FFLAGS  %INCDIR
-#bld::tool::fflags::phys::radiation_AR4       %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -Wf,-O,extendreorder
-#bld::tool::fflags::phys::radiation_AR4_param %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -f3
-#bld::tool::fflags::phys::fisrtilp            %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
-#bld::tool::fflags::phys::cv30_routines       %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -Wf,-O,extendreorder
-#bld::tool::fflags::phys::cvltr               %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
-#bld::tool::fflags::phys::clouds_gno          %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
-#bld::tool::fflags::dyn::vlsplt_p             %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
-#bld::tool::fflags::dyn::groupeun_p           %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
-
-
-inc arch.opt
-
-# Pre-process code before analysing dependencies
-bld::pp              1
-
-
-# Ignore the following dependencies
-bld::excl_dep        inc::netcdf.inc
-bld::excl_dep        use::netcdf
-bld::excl_dep        use::typesizes
-bld::excl_dep        h::netcdf.inc
-bld::excl_dep        h::mpif.h
-bld::excl_dep        inc::mpif.h
-bld::excl_dep        use::ioipsl
-bld::excl_dep        use::intersurf
-bld::excl_dep        use::mod_prism_proto
-bld::excl_dep        use::mod_prism_def_partition_proto
-bld::excl_dep        use::mod_prism_get_proto
-bld::excl_dep        use::mod_prism_put_proto
-bld::excl_dep        use::mkl_dfti
-
-# Don't generate interface files
-bld::tool::geninterface none
-
-# Allow ".inc" as an extension for CPP include files
-bld::infile_ext::inc  CPP::INCLUDE
-
-# extension for module output
-bld::outfile_ext::mod .mod
-bld::tool::SHELL   /bin/bash
Index: LMDZ4/trunk/build_gcm
===================================================================
--- LMDZ4/trunk/build_gcm	(revision 5799)
+++ 	(revision )
@@ -1,24 +1,0 @@
-#!/bin/csh
-
-if ( -f '.lock' ) then
-  echo 'ATTENTION: vous etes sans doute en train de compiler le modele par ailleurs'
-  echo "Attendez que la premiere compilation soit terminee pour relancer la suivante."
-  echo "Si vous etes sur que vous ne compilez pas le modele par ailleurs,"
-  echo  vous pouvez continuer en repondant oui.
-  echo "Voulez-vous vraiment continuer?"
-
-  if ( $< == "oui" ) then
-    
-  else
-    exit
-  endif
-else
-  echo "compilation en cours..." > '.lock'
-endif
-
-#set arch=$1
-
-
-fcm build
-
-\rm -f '.lock' 
Index: LMDZ4/trunk/cosp_input_nl.txt
===================================================================
--- LMDZ4/trunk/cosp_input_nl.txt	(revision 5799)
+++ 	(revision )
@@ -1,104 +1,0 @@
-! (c) British Crown Copyright 2008, the Met Office.
-! All rights reserved.
-! 
-! Redistribution and use in source and binary forms, with or without modification, are permitted 
-! provided that the following conditions are met:
-! 
-!     * Redistributions of source code must retain the above copyright notice, this list 
-!       of conditions and the following disclaimer.
-!     * Redistributions in binary form must reproduce the above copyright notice, this list
-!       of conditions and the following disclaimer in the documentation and/or other materials 
-!       provided with the distribution.
-!     * Neither the name of the Met Office nor the names of its contributors may be used 
-!       to endorse or promote products derived from this software without specific prior written 
-!       permission.
-! 
-! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
-! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
-! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR 
-! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 
-! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, 
-! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER 
-! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT 
-! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
-
-
-! Namelist that sets up the main COSP options
-&COSP_INPUT
-  CMOR_NL='./cmor/cosp_cmor_nl.txt', ! CMOR namelist
-  NPOINTS=9026,! Number of gridpoints (klon dans LMDZi : ici correspond a klon de 96x95)
-  NPOINTS_IT=10000,! Max number of gridpoints to be processed in one iteration
-  NCOLUMNS=20,  ! Number of subcolumns
-  NLEVELS=39,   ! Number of model levels 
-  USE_VGRID=.true., ! Use fixed vertical grid for outputs? (if .true. then you need to define number of levels with Nlr)
-  NLR=40,       ! Number of levels in statistical outputs (only used if USE_VGRID=.true.)
-  CSAT_VGRID=.true., ! CloudSat vertical grid? (if .true. then the CloudSat standard grid is used for the outputs.
-                     !  USE_VGRID needs also be .true.)
-  FINPUT='histday.nc', ! NetCDF file with 1D inputs
-!  FINPUT='cosp_input_um_2d.nc', ! NetCDF file with 2D inputs
-  !----------------------------------------------------------------------------------
-  !--------------- Inputs related to radar simulations
-  !----------------------------------------------------------------------------------
-  RADAR_FREQ=94.0, ! CloudSat radar frequency (GHz)
-  SURFACE_RADAR=0, ! surface=1, spaceborne=0
-  use_mie_tables=0,! use a precomputed lookup table? yes=1,no=0
-  use_gas_abs=1,   ! include gaseous absorption? yes=1,no=0
-  do_ray=0,        ! calculate/output Rayleigh refl=1, not=0
-  melt_lay=0,      ! melting layer model off=0, on=1
-  k2=-1,           ! |K|^2, -1=use frequency dependent default
-  use_reff=.true., ! True if you want effective radius to be used by radar simulator (always used by lidar)
-  use_precipitation_fluxes=.true.,  ! True if precipitation fluxes are input to the algorithm 
-  !----------------------------------------------------------------------------------
-  !---------------- Inputs related to lidar simulations
-  !----------------------------------------------------------------------------------
-  Nprmts_max_hydro=12, ! Max number of parameters for hydrometeor size distributions
-  Naero=1,             ! Number of aerosol species (Not used)
-  Nprmts_max_aero=1,   ! Max number of parameters for aerosol size distributions (Not used)
-  lidar_ice_type=0,    ! Ice particle shape in lidar calculations (0=ice-spheres ; 1=ice-non-spherical)
-  OVERLAP=3,   !  overlap type: 1=max, 2=rand, 3=max/rand
-  !----------------------------------------------------------------------------------
-  !---------------- Inputs related to ISCCP simulator
-  !----------------------------------------------------------------------------------
-  ISCCP_TOPHEIGHT=1,  !  1 = adjust top height using both a computed
-                       !  infrared brightness temperature and the visible
-                       !  optical depth to adjust cloud top pressure. Note
-                       !  that this calculation is most appropriate to compare
-                       !  to ISCCP data during sunlit hours.
-                      !  2 = do not adjust top height, that is cloud top
-                       !  pressure is the actual cloud top pressure
-                       !  in the model
-                      !  3 = adjust top height using only the computed
-                       !  infrared brightness temperature. Note that this
-                       !  calculation is most appropriate to compare to ISCCP
-                       !  IR only algortihm (i.e. you can compare to nighttime
-                       !  ISCCP data with this option)
-  ISCCP_TOPHEIGHT_DIRECTION=1,   ! direction for finding atmosphere pressure level
-                                 ! with interpolated temperature equal to the radiance
-                                 ! determined cloud-top temperature
-                                 ! 1 = find the *lowest* altitude (highest pressure) level
-                                 ! with interpolated temperature equal to the radiance
-                                 ! determined cloud-top temperature
-                                 ! 2 = find the *highest* altitude (lowest pressure) level
-                                 ! with interpolated temperature equal to the radiance 
-                                 ! determined cloud-top temperature
-                                 ! ONLY APPLICABLE IF top_height EQUALS 1 or 3
-                                 ! 1 = default setting, and matches all versions of 
-                                 ! ISCCP simulator with versions numbers 3.5.1 and lower
-                                 ! 2 = experimental setting  
-  !----------------------------------------------------------------------------------
-  !-------------- RTTOV inputs
-  !----------------------------------------------------------------------------------
-  Platform=1,    ! satellite platform
-  Satellite=15,  ! satellite
-  Instrument=0,  ! instrument
-  Nchannels=8,   ! Number of channels to be computed
-  Channels=1,3,5,6,8,10,11,13,        ! Channel numbers (please be sure that you supply Nchannels)
-  Surfem=0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,  ! Surface emissivity (please be sure that you supply Nchannels)
-  ZenAng=50.0, ! Satellite Zenith Angle
-  CO2=5.241e-04, ! Mixing ratios of trace gases
-  CH4=9.139e-07,
-  N2O=4.665e-07,
-  CO=2.098e-07
-/
-
-
Index: LMDZ4/trunk/cosp_output_nl.txt
===================================================================
--- LMDZ4/trunk/cosp_output_nl.txt	(revision 5799)
+++ 	(revision )
@@ -1,63 +1,0 @@
-! (c) British Crown Copyright 2008, the Met Office.
-! All rights reserved.
-! 
-! Redistribution and use in source and binary forms, with or without modification, are permitted 
-! provided that the following conditions are met:
-! 
-!     * Redistributions of source code must retain the above copyright notice, this list 
-!       of conditions and the following disclaimer.
-!     * Redistributions in binary form must reproduce the above copyright notice, this list
-!       of conditions and the following disclaimer in the documentation and/or other materials 
-!       provided with the distribution.
-!     * Neither the name of the Met Office nor the names of its contributors may be used 
-!       to endorse or promote products derived from this software without specific prior written 
-!       permission.
-! 
-! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR 
-! IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND 
-! FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR 
-! CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 
-! DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, 
-! DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER 
-! IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT 
-! OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
-
-! Namelist that sets up output-related variables. It controls 
-! the instrument simulators to run and the list of variables 
-! to be written to file
-&COSP_OUTPUT
-  ! Simulator flags
-  Lradar_sim=.false.,
-  Llidar_sim=.true.,
-  Lisccp_sim=.true.,
-  Lmisr_sim=.false.,
-  ! Output variables
-  Lalbisccp=.true.,
-  Latb532=.true.,
-  Lboxptopisccp=.true.,
-  Lboxtauisccp=.true.,
-  Lcfad_dbze94=.false.,
-  Lcfad_lidarsr532=.true.,
-  Lclcalipso=.true.,
-  Lclhcalipso=.true.,
-  Lclisccp2=.true.,
-  Lcllcalipso=.true.,
-  Lclmcalipso=.true.,
-  Lcltcalipso=.true.,
-  Lctpisccp=.true.,
-  Ldbze94=.false.,
-  Ltauisccp=.true.,
-  Ltclisccp=.true.,
-  Llongitude=.false.,
-  Llatitude=.false.,
-  Lparasol_refl=.true.,
-  LclMISR=.false.,
-  Lmeantbisccp=.true.,
-  Lmeantbclrisccp=.true.,
-  ! Use lidar and radar
-  Lclcalipso2=.false.,
-  Lcltlidarradar=.false.,
-  ! These are provided for debugging or special purposes
-  Lfrac_out=.false.,
-  Lbeta_mol532=.true.,  
-/
Index: LMDZ4/trunk/create_make_gcm
===================================================================
--- LMDZ4/trunk/create_make_gcm	(revision 5799)
+++ 	(revision )
@@ -1,253 +1,0 @@
-#!/bin/sh
-#
-# $Header$
-#
-#set -xv
-machine=`hostname`
-os=`uname`
-gcm=`pwd`
-libf=$gcm/libf
-libo=$gcm/libo
-CRAY=0
-if [ "$machine" = "atlas" -o "$machine" = "etoile" -o "$machine" = "axis" ] ; then
-  CRAY=1
-fi
-XNEC=0
-if [ "$machine" = "rhodes" ] ; then
-  XNEC=1
-fi
-X6NEC=0
-if [ "$machine" = "mercure" ] ; then
-  X6NEC=1
-fi
-X8BRODIE=0
-if [ "$machine" = "brodie" ] ; then
-  X8BRODIE=1
-fi
-VPP=0
-if [ "$machine" = "nymphea0" ] ; then
-  VPP=1
-fi
-#
-echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
-echo "# Definitions de Macros pour Make"
-echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
-echo
-echo "# Repertoires :"
-echo
-echo "GCM     = "$gcm
-echo 'LIBF    = $(GCM)/libf'
-if [ "$CRAY" = '0' ] ; then
-#   echo 'LIBO    = $(GCM)/libo/$(MACHINE)'
-   echo 'LIBO    = $(LIBOGCM)/$(MACHINE)'
-else
-   echo 'LIBO    = $(GCM)/libo'
-fi
-#echo 'LOCAL_DIR=$(GCM)'
-#echo $localdir
-echo "LOCAL_DIR=`echo $localdir`"
-echo 'BIBIO    = $(LIBF)/bibio'
-echo "FILTRE   = filtre"
-echo "PHYS  = "
-echo "DYN  = dyn "
-echo 'LIBPHY = $(LIBO)/libphy$(PHYS).a'
-echo 'DIRMAIN=dyn$(DIM)d$(FLAG_PARA)'
-echo 'RM=rm'
-echo
-echo "OPLINK = "
-echo
-echo '# Les differentes librairies pour l"edition des liens:'
-echo
-if ( [ "$XNEC" = '1' ] || [ "$X6NEC" = '1' ] || [ "$X8BRODIE" = '1' ] ) ; then
-  echo 'dyn3d      = $(LIBO)/libsxdyn3d.a $(LIBO)/libsx$(FILTRE).a'
-  echo 'dyn3dpar     = $(LIBO)/libsxdyn3dpar.a $(LIBO)/libsx$(FILTRE).a'
-  echo 'dyn2d      = $(LIBO)/libsxdyn2d.a'
-  echo 'dyn1d      = $(LIBO)/libsxdyn1d.a'
-  echo 'L_DYN      = -lsxdyn$(DIM)d$(FLAG_PARA)'
-  echo 'L_FILTRE   = -lsx$(FILTRE)'
-  echo 'L_PHY = -lsxphy$(PHYS) '
-  echo 'L_BIBIO    = -lsxbibio'
-  echo 'L_ADJNT    ='
-  echo 'L_COSP     = -lsxcosp'
-else
-  echo 'dyn3d            = $(LIBO)/libdyn3d.a $(LIBO)/lib$(FILTRE).a'
-  echo 'dyn3dpar      = $(LIBO)/libdyn3dpar.a $(LIBO)/lib$(FILTRE).a'
-  echo 'dyn2d            = $(LIBO)/libdyn2d.a'
-  echo 'dyn1d            = $(LIBO)/libdyn1d.a'
-  echo 'L_DYN      = -ldyn$(DIM)d$(FLAG_PARA)'
-  echo 'L_FILTRE   = -l$(FILTRE)'
-  echo 'L_PHY = -lphy$(PHYS) '
-  echo 'L_BIBIO    = -lbibio'
-  echo 'L_ADJNT    ='
-  echo 'L_COSP     = -lcosp'
-fi
-
-echo
-echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
-echo "# Option de compilation FORTRAN"
-echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
-echo
-   echo 'COMPILE = $(F77) $(OPTIM) $(INCLUDE) -c'
-   echo 'COMPILE90 = $(F90) $(OPTIM90) $(INCLUDE) -c'
-   echo 'COMPTRU90 = $(F90) $(OPTIMTRU90) $(INCLUDE) -c'
-   echo "LINK    = $LINK"
-   echo "AR      = $AR"
-echo
-echo
-echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
-echo "# Creation des differents executables"
-echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
-echo
-echo "# Executables:"
-echo "# ------------"
-echo
-echo "PROG = code"
-echo
-#echo 'main : chimie $(DYN) bibio phys $(OPTION_DEP) '
-echo 'main : $(DYN) bibio phys $(OPTION_DEP) '
-echo '	cd $(LIBO) ; $(RANLIB) lib*.a ; cd $(GCM) ;\'
-echo '	cd $(LOCAL_DIR); \'
-echo '	$(COMPILE90) $(LIBF)/$(DIRMAIN)/$(SOURCE) -o $(PROG).o ; \'
-echo '	$(LINK) $(PROG).o -L$(LIBO) $(L_DYN) $(L_ADJNT) $(L_COSP) $(L_FILTRE) $(L_PHY) $(L_DYN) $(L_BIBIO) $(L_DYN) $(OPLINK) $(OPTION_LINK) -o $(LOCAL_DIR)/$(PROG).e ; $(RM) $(PROG).o '
-echo
-echo 'dyn : $(LIBO)/libdyn$(DIM)d$(FLAG_PARA).a $(FILTRE)$(DIM)d'
-echo
-echo 'phys : $(LIBPHY)'
-echo
-#echo 'chimie : $(LIBO)/libchimie.a'
-echo
-echo 'bibio : $(LIBO)/libbibio.a'
-echo
-echo 'adjnt : $(LIBO)/libadjnt.a'
-echo
-echo 'cosp : $(LIBO)/libcosp.a'
-echo
-echo '$(FILTRE)3d : $(LIBO)/lib$(FILTRE).a'
-echo
-echo '$(FILTRE)2d :'
-echo
-echo '$(FILTRE)1d :'
-echo
-echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
-echo "# Contenu des differentes bibliotheques"
-echo "#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"
-echo
-echo
-cd $libf >/dev/null 2>&1
-for diri in ` ls `
-do
-   if [ -d $diri ] ; then
-   if [ "`ls $diri/*.F`" != "" ] || [ "`ls $diri/*.F90`" != "" ]  ; then 
-      cd $diri >/dev/null 2>&1
-      echo
-      listlib=""
-      for i in `ls *.F`
-      do
-         fili=`basename $i .F`
-         test=` (  head $i | grep '      PROGRAM' ) `
-         if [ "$test" = "" ] ; then 
-            listlib=$listlib" "$fili
-         fi
-      done
-      for i in `ls *.F90`
-      do
-         fili=`basename $i .F90`
-         test=` (  head $i | grep '      PROGRAM' ) `
-         if [ "$test" = "" ] ; then
-            listlib=$listlib" "$fili
-         fi 
-      done
-#
-      echo
-      echo
-      echo "#======================================================================="
-      echo "# Contenu de la bibliotheque correspondant au Directory "$diri
-      echo "#======================================================================="
-      echo
-      for fili in $listlib
-      do
-         echo '$(LIBO)/lib'$diri".a : " '$(LIBO)/lib'$diri".a("$fili".o)"
-         echo
-      done
-      echo '.PRECIOUS	: $(LIBO)/lib'$diri'.a'
-      echo
-      echo
-      echo "# Compilation des membres de la bibliotheque lib"$diri".a"
-      echo
-      for fili in $listlib
-      do
-         if [ -f $fili.F90 ] ; then
-           trufile=$fili.F90
-         else
-           trufile=$fili.F
-         fi
-         F90=0 ; egrep -i '^ *use ' $trufile > /dev/null 2>&1 && F90=1
-                 egrep -i '^ *module ' $trufile > /dev/null 2>&1 && F90=1
-                 egrep -i '#include*.inc ' $trufile > /dev/null 2>&1 && F90=1
-         str1='$(LIBO)/lib'$diri'.a('$fili'.o) : $(LIBF)/'$diri/$trufile
-         [ "$fili" = "chem.subs" ] && str1=$str1' $(LIBF)/'$diri/chem.mods.F
-         for stri in ` ( sed -n "/\#include/s/\#include//p" $trufile | sed 's/\"//g' ; egrep -i '^ *use ' $trufile | sed -e 's/,/ /' | awk ' { print $2 } ' ) `
-         do
-
-
-# Differents cas de dependance correspondant a des include ou des
-# use module.
-# soit dans le repertoire local soit dans un autre.
-
-            stri=`echo $stri | tr [A-Z] [a-z]`
-            if [ -f $stri ] ; then
-               echo $str1 \\
-               str1='$(LIBF)/'$diri'/'$stri
-            else
-               if [ -f $stri.F ] || [ -f $stri.F90 ] ; then
-                  echo $str1 \\
-                  str1='$(LIBO)/lib'$diri'.a('$stri'.o)'
-               else
-                  for dirinc in dyn3d grid bibio filtrez
-                  do
-                     if [ -f ../$dirinc/$stri ] ; then
-                        echo $str1 \\
-                        str1='$(LIBF)/'`cd .. ; ls */$stri | head -1`
-                     fi
-                     if [ -f ../$dirinc/$stri.F90 ] ; then
-                        echo $str1 \\
-                        str1='$(LIBO)/lib'$dirinc'.a('$stri'.o)'
-                     fi
-                  done
-               fi
-            fi
-         done
-         echo $str1
-	 if [ "$F90" -eq '0' ] ; then
-           echo '	cd $(LOCAL_DIR); \'
-	   echo '	$(COMPILE) $(LIBF)/'$diri'/'$trufile' ; \'
-	 else
-           echo '	cd $(LOCAL_DIR); \'
-           if [ -f $fili.F90 ] ; then
-	      echo '	$(COMPTRU90) $(LIBF)/'$diri'/'$trufile' ; \'
-           else
-	      echo '	$(COMPILE90) $(LIBF)/'$diri'/'$trufile' ; \'
-           fi
-           MODU=0; egrep -i '^ *module ' $trufile> /dev/null 2>&1 && MODU=1
-            if [ "$MODU" -eq '1' -a "$CRAY" != '1' ] ; then
-              if [ "$os" = 'UNIX_System_V' ] ; then
-                echo '	cp $(MOD_LOC_DIR)/*.$(MOD_SUFFIX) $(LIBO)/ ; \'
-              else
-                echo '	mv $(MOD_LOC_DIR)/'$fili'.$(MOD_SUFFIX) $(LIBO)/'$fili'.$(MOD_SUFFIX) ; \'
-#                echo '	if [ "$(MOD_LOC_DIR)" ne "$(LIBO)" ] ; then mv $(MOD_LOC_DIR)/'*'.$(MOD_SUFFIX) $(LIBO) ; fi ; \'
-              fi
-            fi
-	 fi
-	 if ( [ "$XNEC" -eq '1' ] || [ "$X6NEC" = '1' ] || [ "$X8BRODIE" = '1' ] ) ; then
-           echo '	sxar r $(LIBO)/libsx'$diri'.a '$fili'.o ; \'
-         fi
-         echo '	$(AR) r $(LIBO)/lib'$diri'.a '$fili'.o ; $(RM) '$fili'.o ; \'
-         echo '	cd $(GCM)'
-	 echo
-      done
-#	 
-      echo
-   cd $libf
-   fi
-   fi
-done
Index: LMDZ4/trunk/gcm.def
===================================================================
--- LMDZ4/trunk/gcm.def	(revision 5799)
+++ 	(revision )
@@ -1,62 +1,0 @@
-## $Id$
-## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
-day_step=480
-## periode pour le pas Matsuno (en pas)
-iperiod=5
-## periode de la dissipation (en pas)
-idissip=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:y ou Matsuno-leapfrog:n) 
-purmats=n
-## avec ou sans physique
-## 0: pas de physique (e.g. en mode Shallow Water)
-## 1: avec physique (e.g. physique phylmd)
-## 2: avec rappel newtonien dans la dynamique                                         
-iflag_phys=1
-## avec ou sans fichiers de demarrage (start.nc, startphy.nc) ?
-## (sans fichiers de demarrage, initialisation des champs par iniacademic
-##  dans la dynamique)
-read_start=y
-## periode de la physique (en pas dynamiques, n'a de sens que si iflag_phys=1)                                       
-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=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: LMDZ4/trunk/guide.def
===================================================================
--- LMDZ4/trunk/guide.def	(revision 5799)
+++ 	(revision )
@@ -1,37 +1,0 @@
-ok_guide=y
-# guidage sur niveaux modèle (y) ou standard
-guide_modele=y
-# inversion de l'ordre des niveaux verticaux
-ok_invertp=y
-ncep=y
- ######################################
- #### guidage de u #####
- guide_u=y
- ######################################
- #### guidage de v #####
- guide_v=y
- ######################################
- #### guidage de T #####
- guide_T=y
- ######################################
- #### guidage de P #####
- guide_P=n
- ######################################
- #### guidage de Q (hr=y:hum.rel,n:hum.spec) #####
- guide_Q=n
- guide_hr=n
- ######################################
- ## guidage dans la couche limite
- guide_BL=n
- ######################################
-ini_anal=n
-tau_min_u=0.04166667
-tau_max_u=0.125
-tau_min_v=0.04166667
-tau_max_v=0.125
-tau_min_T=0.04166667
-tau_max_T=10.
-tau_min_Q=0.2
-tau_max_Q=10.
-# gamma limité
-gamma4=n
Index: LMDZ4/trunk/makegcm
===================================================================
--- LMDZ4/trunk/makegcm	(revision 5799)
+++ 	(revision )
@@ -1,1102 +1,0 @@
-#!/bin/csh
-#
-# $Id$
-#
-#set verbose echo
-########################################################################
-# options par defaut pour la commande make
-########################################################################
-set dim="96x71x19"
-set physique=lmd
-set phys="PHYS=$physique"
-set include='-I$(LIBF)/grid -I$(LIBF)/bibio -I$(LIBF)/filtrez -I. '
-set filtre=filtrez
-set grille=reg
-set couple=false
-set veget=false
-set chimie=false
-set psmile=true
-set parallel=false
-set vampir=false
-set OPT_STACK='-Wf,-init stack=nan'
-set OPT_STACK=' '
-set OPTIMI='-C debug -eC '
-set OPTIMI=' -ftrace '
-set OPT_LINUX='-O3'
-set OPT_LINUX="-i4 -r8 -O3"
-set io=ioipsl
-set cosp=false
-
-set FC_LINUX=g95
-set FC_LINUX=gfortran
-#set FC_LINUX=pgf90
-if ( $FC_LINUX == g95 ) then
-   set OPT_LINUX="-O3"
-   set OPT_LINUX="-O3"
-else if ( $FC_LINUX == gfortran ) then
-   set OPT_LINUX="-fdefault-real-8 -O3"
-#   set OPT_LINUX="-O3 -fno-second-underscore"
-   set OPT_LINUX="-O3 "
-else
-   # pgf90 options
-   set OPT_LINUX="-i4 -r8 -O2 -Munroll -Mnoframe -Mautoinline -Mcache_align"
-endif
-
-########################################################################
-# path a changer contenant les sources et les objets du modele
-########################################################################
-
-###### VERSION LMDZ.4
-set INCALIB=../INCA3/config/lib
-set LMDGCM="`pwd`"
-setenv LIBOGCM $LMDGCM/libo
-#
-#
-setenv IOIPSLDIR /d4/fairhead/gfortran/ioipsl_v2_1_9
-setenv MODIPSLDIR /d4/fairhead/gfortran/ioipsl_v2_1_9
-setenv NCDFINC /d4/fairhead/gfortran_4.4/netcdf-4.1.1/include
-setenv NCDFLIB /d4/fairhead/gfortran_4.4/netcdf-4.1.1/lib
-
-
-
-
-
-setenv localdir "`pwd`"
-set MODIPSL=0
-echo $localdir | grep modipsl >& /dev/null
-if ( ! $status ) then
-  set MODIPSL=1
-  setenv LMDGCM $localdir
-  cd ../..
-  setenv LIBOGCM "`pwd`/lib"
-  setenv IOIPSLDIR $LIBOGCM
-  setenv MODIPSLDIR $LIBOGCM
-  cd $localdir
-  if ( `hostname` == rhodes ) then
-    set NCDFINC=`grep sxnec ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"`
-    set NCDFLIB=`grep sxnec ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'`
-  else
-    if ( `hostname` == nymphea0 ) then
-      set NCDFINC=`grep fjvpp ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"`
-      set NCDFLIB=`grep fjvpp ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'`
-    else if ( `hostname` == mercure ) then
-      set NCDFINC=`grep sx6nec ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"`
-      set NCDFLIB=`grep sx6nec ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'`
-    else  if ( `hostname` == brodie ) then
-      set NCDFINC=`grep sx8brodie ../../util/AA_make.gdef| grep NCDF_INC|sed -e "s/^.* =//"`
-      set NCDFLIB=`grep sx8brodie ../../util/AA_make.gdef| grep NCDF_LIB|sed -e 's/^.* =//'`
-    else
-      echo 'Probleme de definition des variables NCDFINC et NCDFLIB'
-    endif 
-  endif 
-else
-  if ( ! $?LMDGCM ) then
-    echo You must initialize the variable LMDGCM in your environnement
-    echo for instance: "setenv LMDGCM /usr/myself/supergcm" in .cshrc
-    exit
-  endif
-  if ( ! $?LIBOGCM ) then
-    set LIBOGCM=$LMDGCM/libo
-  endif
-  if ( ! $?IOIPSLDIR ) then
-    echo You must initialize the variable IOIPSLDIR in your environnement
-    echo for instance: "setenv IOIPSLDIR /usr/myself/ioipsl" in .cshrc
-    exit
-  else
-      setenv MODIPSLDIR $IOIPSLDIR
-  endif
-  if ( ! $?NCDFLIB ) then
-    echo You must initialize the variable NCDFLIB in your environnement
-    echo for instance: "setenv NCDFLIB /usr/myself/netcdf" in .cshrc
-    exit
-  endif
-  if ( ! $?NCDFINC ) then
-    echo You must initialize the variable NCDFINC in your environnement
-    echo for instance: "setenv NCDFINC /usr/myself/netcdf" in .cshrc
-    exit
-  endif
-endif
-set model=$LMDGCM
-set libo=$LIBOGCM
-
-########################################################################
-#  Les differentes platformes reconnues
-########################################################################
-
-set HP=0
-set IBM=0
-set SUN=0
-set VPP=0
-set CRAY=0
-set DEC=0
-set LINUX=0
-set NEC=0
-set XNEC=0
-set X6NEC=0
-set X8BRODIE=0
-if ( `uname` == HP-UX ) then
-   set machine=HP
-   set HP=1
-else if (`uname` == UNIX_System_V ) then
-   set machine=VPP
-   set VPP=1
-else if (`uname` == SunOS ) then
-   set machine=SUN
-   set SUN=1
-else if ( `uname` == AIX ) then
-   set machine=IBM
-   set IBM=1
-else if ( `uname` == OSF1 ) then
-   set machine=ALPHA
-   set DEC=1
-else if ( `uname` == Linux && `hostname` != mercure  && `hostname` != brodie ) then
-   set machine=LINUX
-   set LINUX=1
-else if ( `hostname` == atlas || `hostname` == axis  || `hostname` == etoile ) then
-   set machine=CRAY
-   set CRAY=1
-else if ( `uname` == SUPER-UX ) then
-   set machine=NEC
-   set NEC=1
-else if ( `hostname` == rhodes) then
-   set machine=XNEC
-   set XNEC=1
-else if ( `hostname` == mercure) then
-   set machine=X6NEC
-   set X6NEC=1
-else if ( `hostname` == brodie) then
-   set machine=X8BRODIE
-   set X8BRODIE=1
-else
-   echo Vous travaillez sur une machine non prevue par le reglement
-   exit
-endif
-
-if ( ! -d $libo )  then
-   mkdir $libo
-endif
-
-
-if $VPP then
-set netcdf=netcdf_v
-else 
-set netcdf=netcdf
-endif
-########################################################################
-#  Quelques initialisations de variables du shell.
-########################################################################
-
-set dyn=
-set opt_link=""
-set adjnt=""
-set lcosp=""
-set opt_dep=""
-set libchimie=""
-
-set optim=""
-set optimbis=""
-set optim90=""
-set oplink=""
-
-########################################################################
-#  Optimisations par defaut suivant les machines
-########################################################################
-
-echo "Optimisations par defaut suivant les machines"
-set libf=$model/libf
-#setenv localdir "LOCAL_DIR=`pwd`"
-#setenv localdir "`pwd`"
-cd $model
-#############
-if $CRAY then
-#############
-#   set optim="-Wf'-ei' -dp -Wf'-a static'"
-   set optimbis=" -DCRAY "
-   set optim90="-Wp'-P' -DCRAY -p$IOIPSLDIR "'-p$(LIBO) -eiv '
-   set optim="$optim90"
-   if ( $io == "ioipsl" ) then
-     set oplink="-Wl'-DSTACK=128 -f indef' -L$IOIPSLDIR -lioipsl  -L$NCDFLIB -lnetcdf "
-   else
-     set oplink="-Wl'-DSTACK=128 -f indef' -L$IOIPSLDIR -L$NCDFLIB -lnetcdf "
-   endif
-   set mod_loc_dir=" "
-   set mod_suffix=" "
-#################
-else if $SUN then
-#################
-   set optim=" -fast "
-   set optimbis=" "
-   set optim90=" -fast -fixed "
-   set optimtru90=" -fast -free "
-   if ( $io == "ioipsl" ) then
-     set opt_link="-lf77compat -L$MODIPSLDIR -lsechiba -lparameters -lstomate -lioipsl -L$NCDFLIB -lnetcdf "
-   else
-     set opt_link="-lf77compat -L$MODIPSLDIR -lsechiba -lparameters -lstomate -L$NCDFLIB -lnetcdf "
-   endif
-   set mod_loc_dir=$localdir
-   set mod_suffix=mod
-#################
-else if $HP then
-#################
-   set optim=" +U77 -O +E1 "
-   set optimbis=" "
-#################
-else if $IBM then
-#################
-   set optim=" -O3 -qtune=pwr2 -qarch=pwr2"
-   set optimbis=" "
-#################
-else if $VPP then
-#################
-#   set optim="-Dasuxm  -On, -g -Ad -Potilax -Eciplume -Si"
-#   set optimbis="  -Wv,-m3 -Wp,-DVPP -Z $LMDGCM/listage"
-   set optimbis=" -Wp,-DNC_DOUBLE -Ad -Z $LMDGCM/listage -X9"
-   set optim90="$optim $optimbis -X9 -w"
-   set mod_loc_dir=$MODIPSLDIR
-   set mod_suffix=mod
-#################
-else if $DEC then
-#################
-   set optim=" "
-   set optimbis=" "
-#################
-else if $LINUX then
-#################
-   if ( $FC_LINUX == pgf90 || $FC_LINUX == g95 || $FC_LINUX == gfortran ) then
-     set optim=" $OPT_LINUX "
-     set optim90=" $OPT_LINUX "
-     set optimtru90=" $OPT_LINUX "
-   else
-     echo 'compilateur linux non reconnu'
-     exit
-   endif
-   set mod_loc_dir=$MODIPSLDIR
-   set mod_suffix=mod
-#################
-else if $NEC then
-#################
-   set optim90=' -clear -C hopt -float0 -ew -P static -Wf,"-pvctl fullmsg noassume "'
-   set optimtru90=' -clear -f4 -C hopt -float0 -ew -P static -Wf,"-pvctl fullmsg noassume "'
-   set optim="$optim90"
-   set optimbis=" "
-   if ( $io == "ioipsl" ) then
-     set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR -lioipsl  $NCDFLIB -lnetcdf_i8r8_v "
-   else
-     set opt_link=" -C hopt -float0 -ew -P static -L$MODIPSLDIR $NCDFLIB -lnetcdf_i8r8_v "
-   endif
-   set mod_loc_dir="."
-   set mod_suffix="mod"
-#################
-else if $XNEC then
-#################
-   set optdbl='-dw -Wf\"-A dbl4\"'
-   set optim90=' -clear -float0 -f3 -Ep -DNC_DOUBLE -dw -Wf\"-A dbl4\" -R5 -Wf,"-pvctl loopcnt=40000 fullmsg noassume "'
-   set optimtru90=' -clear -f4 -float0 -Ep -DNC_DOUBLE -dw -Wf\"-A dbl4\" -R2 -R3 -R4 -R5 -Wf,"-pvctl loopcnt=40000 fullmsg noassume"'
-   set optim="$optim90"
-   set optimbis=" "
-   set mod_suffix="mod"
-   set mod_loc_dir="./"
-#################
-else if $X6NEC then
-#################
-   set optdbl='-dw -Wf\"-A dbl4\"'  
-   set optim90=' -clear -float0 -size_t64 -P stack -Wf "-ptr byte -init stack=nan -init heap=nan" -Ep -DNC_DOUBLE -dw -Wf\"-A dbl4\" -R5 -Wf,"-pvctl loopcnt=40000 fullmsg noassume "'
-   set optimtru90=' -clear -f4 -float0 -size_t64 -P stack -Wf "-ptr byte -init stack=nan -init heap=nan" -Ep -DNC_DOUBLE -dw -Wf\"-A dbl4\" -R2 -R3 -R4 -R5 -Wf,"-pvctl loopcnt=40000 fullmsg noassume"'
-   set optim="$optim90"
-   set optimbis=" "
-   set mod_suffix="mod"
-   set mod_loc_dir="./"
-#################
-else if $X8BRODIE then
-##################
-   set optdbl='-dw -Wf\"-A dbl4\"'  
-#   set optim90='-P stack -Wf,-pvctl res=whole,-A dbl4,-init stack=nan,-init heap=nan,-ptr byte -EP -R5 -float0 -dw -Wf,"-pvctl loopcnt=999999 fullmsg noassume" -I/SX/usr/include'
-   set optim90='-C vopt -Wf,-pvctl res=whole,-A dbl4,-init stack=nan,-init heap=nan,-ptr byte -EP -DNC_DOUBLE -R5 -float0 -dw -Wf,"-pvctl loopcnt=999999 noassume" -I/SX/usr/include'
-#   set optim90='-C vsafe -P stack -Wf,-pvctl res=whole,-A dbl4,-ptr byte -EP -R5 -float0 -dw -Wf,"-pvctl loopcnt=999999 fullmsg noassume" -I/SX/usr/include'
-   set optimtru90="$optim90"
-   set optim90="$optim90"
-   set optim="$optim90"
-   set optimbis=" "
-   set mod_suffix="mod"
-   set mod_loc_dir="./"
-else
-   set optim=""
-   set optimbis=" "
-endif
-
-set nomlib=$machine
-
-########################################################################
-# lecture des options de mymake
-########################################################################
-
-top:
-if ($#argv > 0) then
-    switch ($1:q)
-
-    case -h:
-
-########################################################################
-# Manuel en ligne
-########################################################################
-more <<eod
-
-
-makegcm [Options] prog
-
-
-
-
-Par default, la commande makegcm:
----------------------------------
-
-1. compile une serie de sous programmes se trouvant dans des sous-repertoires
-de $LMDGCM/libf.
-Les sous programmes sont ensuite stokes sur dans des librairies FORTRAN
-sur $LIBOGCM.
-
-2. Ensuite, makegcm compile le programme prog.f se trouvant par default sur
-$LMDGCM/libf/dyn3d et effectue le lien avec l ensemble des librairies.
-
-La variable '$LMDGCM' doit etre initialisee dans votre .cshrc ou en dur
-dans la comande makegcm.
-
-La commande makegcm est faite pour permettre de gerer en parallele des
-versions differentes du modele, compilees avec des options de compilation
-et des dimensions differentes sans avoir a chaque fois a recompiler tout
-le modele.
-
-Les librairies FORTRAN sont stoquees sur le directory $LIBOGCM.
-
-
-OPTIONS:
---------
-
-Les options suivantes peuvent etre definies soit par defaut en editant le
-"script" makegcm, soit en interactif:
-
--d imxjmxlm  ou im, jm, et lm sont resp. le nombre de longitudes, latitudes
-             et couches verticales. 
-             L'effet des options -d est d'ecraser le fichier 
-             $LMDGCM/libf/grid/dimensions.h
-             qui contient sous forme de 3 PARAMETER FORTRAN les 3 dimensions
-             de la grille horizontale im, jm et verticale lm, par un nouveu fichier
-             $LMDGCM/libf/grid/dimension/dimensions.im.jm.lm
-             Si ce fichier n'existe pas encore, il est cree par le script
-             $LMDGCM/libf/grid/dimension/makdim
-
--p PHYS    pour selectionner le jeu de parametrisations physiques avec
-           lequel on veut compiler le modele.
-           Le modele sera alors compile en prenant les sources des
-           parametrisations physiques dans le repertoire:
-            $LMDGCM/libf/phyPHYS
-
--c false|true
-           pour selectionner le mode force (par defaut) ou couple
-
--io ioipsl|noioipsl
-           pour selectionner le logiciel IO : IOIPSL par defaut
-
--psmile false|true
-           pour selectionner le mode psmile ou non (par defaut)
-
--parallel  false|true
-           pour selectionner le mode parallele ou non (false par defaut)
-
--v true|false
-           pour selectionner la vegetation (par defaut) ou non
-
--chimie INCA|false
-	   pour selectionner ou non la chimie (par defaut sans)
-
--g grille  selectionne le type de grille qu'on veut utiliser.
-           L'effet de cette option est d'ecraser le fichier
-           $LMDGCM/libf/grid/fxyprim.h
-           avec le fichier
-           $LMDGCM/libf/grid/fxy_grille.h
-           grille peut prendre les valeurs:
-           1. reg pour la grille reguliere
-           2. sin pour avoir des points equidistants en sinus de la latitude
-           3. new pour pouvoir zoomer sur une partie du globe
-
--O "optimisation fortran" ou les optimisations fortran sont les options de la
-            commande f77
-
--include path
-           Dans le cas ou on a dans des sous programmes des fichiers 
-           #include (cpp) qui se trouve sur des repertoires non references
-           par defaut
-
--adjnt     Pour compiler la l'adjoint du code dynamique
-
--cosp true|false      
-           Pour compiler avec cosp
-
--filtre  filtre
-           Pour choisir le filtre en longitude dans les regions polaires.
-           "filtre" correspond au nom d'un repertoire se trouvant sur
-           $LMDGCM/libf. Le filtre standard du modele est "filtrez" qui peut
-           etre utilise aussi bien pour une grille reguliere que pour une 
-           grille zoomee en longitude.
-
--link "-Ldir1 -lfile1 -Ldir2 -lfile2 ..."
-          Pour rajouter un lien avec les librairies FORTRAN
-          libfile1.a, libfile2.a ... se trouvant respectivement sur les
-          repertoires dir1, dir2 ...
-          Si dirn est un repertoire dont le chemin est automatique (comme
-          par exemple /usr/lib ...) il n'est pas besoin de specifier -Ldirn.
-
-Auteur: Frederic Hourdin  (hourdin@lmd.jussieu.fr)
-eod
-exit
-
-########################################################################
-# Lecture des differentes options
-########################################################################
-
-    case -d:
-        set dim=$2 ; shift ; shift ; goto top
-                        
-    case -O:
-        set optim="$2" ; shift ; shift ; goto top
-
-     case -p
-        set physique="$2" ; set phys="PHYS=$physique" ; shift ; shift ; goto top
-
-     case -g
-        set grille="$2" ; shift ; shift ; goto top
-
-     case -c
-        set couple="$2" ; shift ; shift ; goto top
-
-     case -io
-        set io="$2" ; shift ; shift ; goto top
-
-     case -v
-        set veget="$2" ; shift ; shift ; goto top
-
-     case -chimie
-	set chimie="$2" ; shift ; shift ; goto top
-
-     case -parallel
-        set parallel="$2" ; shift ; shift ; goto top
-  
-     case -include
-        set include="$include -I$2" ; shift ; shift ; goto top
-
-     case -adjnt
-        set opt_dep="$opt_dep adjnt" ; set adjnt="-ladjnt -ldyn3d "
-        set optim="$optim -Dadj" ; shift ; goto top
-
-     case -cosp
-        set cosp="$2"; shift ; shift ; goto top
-
-     case -filtre
-        set filtre=$2 ; shift ; shift ; goto top
-
-     case -link
-        set opt_link="$opt_link $2" ; shift ; shift ; goto top
-
-     case -debug
-        if $HP then
-           set optim=" -g "
-        else if $SUN then
-           setenv PARALLEL 2
-# Il faut rajouter l'option -dalign a -g pour pouvoir editer les liens
-# avec des programmes compiles avec -fast
-           set optim=" -g -dalign "
-           set optim90=" -fixed -g "
-           set optimtru90=" -free -g -C -dalign "
-        else if $CRAY then
-           set optim="$optim"" -g "
-           set optim90="$optim90"" -G1 "
-        else if $LINUX then
-           if ( $FC_LINUX == "pgf90" ) then
-             set optim="$optim"" -g -C -Mbounds "
-             set optim90="$optim90"" -g -C -Mbounds "
-             set optimtru90="$optimtru90"" -g -C -Mbounds "
-           else if ( $FC_LINUX == 'g95' ) then
-             set optim="$optim"" -g -fbounds-check "
-             set optim90="$optim90"" -g -fbounds-check "
-             set optimtru90="$optimtru90"" -g -fbounds-check "
-           else
-             echo 'compilateur linux non reconnu'
-             exit
-           endif
-        else 
-           echo pas d option debug predefinie pour cette machine
-           exit
-        endif
-        shift ; goto top
-
-     default
-        set code="$1" ; shift ; goto top
-
-   endsw
-endif
-
-########################################################################
-# Definition des clefs CPP
-########################################################################
-
-set cppflags=''
-
-if $X8BRODIE then
-  set cppflags="$cppflags -DNC_DOUBLE -DBLAS -DSGEMV=DGEMV -DSGEMM=DGEMM"
-endif
-
-if ( $io == ioipsl ) then
-   set cppflags="$cppflags -DCPP_IOIPSL"
-endif
-
-if ( "$cosp" == 'true' ) then
-    set cppflags="$cppflags -DCPP_COSP"
-    set include="$include"' -I$(LIBF)/cosp '
-    set opt_dep="$opt_dep cosp"
-#    set lcosp="-lcosp -lphy$physique "
-     set lcosp="-lcosp "
-   if ( $XNEC || $X8BRODIE || $X6NEC) then
-#    set lcosp="-lsxcosp -lsxphy$physique "
-     set lcosp="-lsxcosp "
-   endif
-endif
-
-if ( "$physique" == 'nophys' ) then
-   set phys="L_PHY= LIBPHY="
-else
-   #Default planet type is Earth
-   set cppflags="$cppflags -DCPP_EARTH"
-endif
-
-set link_veget=" "
-if ( "$veget" == 'true' ) then
-   set cppflags="$cppflags -DCPP_VEGET"
-#   set link_veget=" -lsechiba -lparameters -lstomate -lorglob -lparallel"
-   set link_veget=" -lsechiba -lparameters -lstomate"
-   if ( $XNEC || $X8BRODIE || $X6NEC) then
-#      set link_veget=" -lsxsechiba -lsxparameters -lsxstomate -lsxorglob -lsxparallel"
-      set link_veget=" -lsxsechiba -lsxparameters -lsxstomate "
-   endif
-endif
-
-if ( "$chimie" == 'INCA' ) then
-    set cppflags="$cppflags -DINCA" 
-    set libchimie=" -L$INCALIB -lchimie"
-    set opt_link="$opt_link  -L$INCALIB -lchimie"
-endif
-    
-if ( "$couple" == 'true' ) then
-   set cppflags="$cppflags -DCPP_COUPLE"
-endif
-
-set FLAG_PARA=''
-if ( "$parallel" == 'true' ) then
-   set cppflags="$cppflags -DCPP_PARA"
-   set FLAG_PARA='par'
-endif
-
-set optim="$optim $cppflags"
-set optim90="$optim90 $cppflags"
-set optimtru90="$optimtru90 $cppflags"
-
-
-########################################################################
-# cas special sans physique
-########################################################################
-if ( "$physique" == 'nophys' ) then
-   set phys="L_PHY= LIBPHY="
-endif
-
-########################################################################
-#subtilites sur le nom de la librairie
-########################################################################
-
-\rm -f tmp ; touch tmp
-\rm -f tmp90 ; touch tmp90
-foreach i ( $optim )
-   echo $i | sed -e 's/\"//g' -e "s/\'//g" -e 's/-//g'  >> tmp
-end
-set suf=
-foreach i ( `sort tmp | uniq ` )
-   set suf=$suf$i
-end
-if ( ! $IBM ) then
-   set nomlib="$nomlib$suf"
-endif
-if ( $DEC ) then
-   set nomlib=DEC
-endif
-if ( $IBM ) then
-   set dim=`echo $dim | sed -en 's/[^0-9]/ /g'`
-   set dim_=`echo $dim | sed -en 's/[^0-9]/_/g'`
-else if ( $SUN ) then
-   set dim=`echo $dim | sed -e 's/[^0-9]/ /g'` 
-   set dim_=`echo $dim | sed -e 's/[^0-9]/_/g'`
-else
-   set dim_=`echo $dim | sed -e 's/[^0-9]/_/g'`
-   set dim=`echo $dim | sed -e 's/[^0-9]/ /g'`
-endif
-set nomlib=${nomlib}${physique}_${dim_}_$grille
-## M-A-F nomlib trop long sur CRAY pour ar
-if ( $CRAY ) then
-    set nomlib=F90_${dim_}
-endif
-if ( $NEC || $XNEC || $X6NEC || $X8BRODIE ) then
-    set nomlib=F90_${dim_}_'phy'${physique}${FLAG_PARA}
-endif
-echo calcul de la dimension
-set dimc=`echo $dim | wc -w`
-
-if ( "$dimc" == "2" ) then
-set include="$include "'-I$(LIBF)/dyn2d '
-set dimh=$dim
-else
-set include="$include "'-I$(LIBF)/dyn3d${FLAG_PARA} '
-set dimh=`echo $dim | awk ' { print $1 "." $2 } '`
-endif
-echo $dimc
-
-########################################################################
-# path pour les #include
-########################################################################
-
-if ( $XNEC ) then
-  set include="$include -I$NCDFINC -I$IOIPSLDIR"
-else
-  set include="$include -I$NCDFINC -I$IOIPSLDIR"
-endif
-echo $include
-
-########################################################################
-# Gestion des dimensions du modele.
-# on cree ou remplace le fichier des dimensions
-########################################################################
-
-cd $libf/grid
-if ( -f dimensions.h ) then
-echo 'ATTENTION: vous etes sans doute en train de compiler le modele par ailleurs'
-echo "Attendez que la premiere compilation soit terminee pour relancer la suivante."
-echo "Si vous etes sur que vous ne compilez pas le modele par ailleurs,"
-echo  vous pouvez continuer en repondant oui.
-echo "Voulez-vous vraiment continuer?"
-if ( $< == "oui" ) then
-\rm -f $libf/grid/dimensions.h
-else
-exit
-endif
-endif
-
-cd dimension
-./makdim $dim
-cat $libf/grid/dimensions.h
-
-cd $LMDGCM
-set libo=$libo/$nomlib
-if ( ! -d $libo )  then
-   mkdir $libo
-   cd $model
-endif
-
-########################################################################
-# Differentes dynamiques (3d, 2d, 1d)
-########################################################################
-
-set dimension=`echo $dim | wc -w`
-echo dimension $dimension
-if ( $dimension == 1 ) then
-echo pas de dynamique
-set dyn="L_DYN= DYN= L_FILTRE= DIRMAIN=phy$physique "
-endif
-endif
-cd $model
-if ( $dimension == 3 ) then
-cd libf/grid
-\rm fxyprim.h
-cp -p fxy_${grille}.h fxyprim.h
-endif
-
-######################################################################
-#   Traitement special pour le nouveau rayonnement de Laurent Li.
-######################################################################
-
-#if ( -f $libf/phy$physique/raddim.h ) then
-# if ( -f $libf/phy$physique/raddim.$dimh.h ) then
-#  \rm -f $libf/phy$physique/raddim.h
-#  cp -p $libf/phy$physique/raddim.$dimh.h $libf/phy$physique/raddim.h
-#  echo $libf/phy$physique/raddim.$dimh.h 
-#  cat $libf/phy$physique/raddim.$dimh.h 
-#  cat $libf/phy$physique/raddim.h
-# else
-#  echo On peut diminuer la taille de l executable en creant
-#  echo le fichier $libf/phy$physique/raddim.$dimh.h
-#  \cp -p $libf/phy$physique/raddim.defaut.h $libf/phy$physique/raddim.h
-# endif
-#endif
-
-######################################################################
-# Gestion du filtre qui n'existe qu'en 3d.
-######################################################################
-
-if ( `expr $dimc \> 2` == 1 ) then
-   set filtre="FILTRE=$filtre"
-else
-   set filtre="FILTRE= L_FILTRE= "
-endif
-echo MACRO FILTRE $filtre
-
-echo $dimc
-
-########################################################################
-#  Avant de lancer le make, on recree le makefile si necessaire
-########################################################################
-########################################################################
-# c'est a dire dans 3 cas:
-# 1. si la liste des fichiers .F et .h a ete modifiee depuis la
-#    derniere creation du makefile
-# 2. si le fichier contenant cette liste "liste_des_sources"
-#    n'existe pas.
-# 3. Si le makefile n'existe pas.
-########################################################################
-##########################################
-# On adapte d'abord certains include à F90
-##########################################
-##########################################
-cd $model
-find libf -name '*.[Fh]' -print >! tmp77
-find libf -name '*.[Fh]' -exec egrep -i " *use *ioipsl" {} \; -print >! tmp90
-find libf -name '*.[Fh]90' -print >> tmp90
-
-if (    `diff tmp77 liste_des_sources_f77 | wc -w` \
-     || `diff tmp90 liste_des_sources_f90 | wc -w` \
-     || ! -f makefile \
-     || ! -f liste_des_sources_f90 \
-     || ! -f liste_des_sources_f77 ) then
-        echo les fichiers suivants ont ete crees ou detruits
-        echo ou les fichiers suivants sont passes ou ne sont plus en Fortran 90
-        diff liste_des_sources_f77 tmp77
-        diff liste_des_sources_f90 tmp90
-        \cp tmp77 liste_des_sources_f77
-        \cp tmp90 liste_des_sources_f90
-        echo On recree le makefile
-        ./create_make_gcm >! tmp
-        \mv tmp makefile
-        echo Nouveau makefile cree.
-endif
-
-########################################################################
-#  Execution de la comande make
-########################################################################
-
-echo PHYSIQUE $phys
-echo dynamique $dyn $dimension
-echo OPTIM="$optim" $filtre LIBO=$libo $dyn PHYS=$phys DIM=$dimc PROG=$code
-echo PATH pour les fichiers INCLUDE $include
-echo OPLINK="$oplink"
-
-#################
-if $HP then
-#################
-   set f77='fort77 +OP'
-   set f90='jensaisrien'
-   set opt_link="$opt_link -lm"
-#################
-else  if $VPP then
-#################
-   set f77=frt
-   set f90=$f77
-   if ($couple == true) then
-     set opt_link="-Wg,-c $MODIPSLDIR/liboasis2.4_mpi2.a /usr/lang/mpi2/lib64/libmpi.a /usr/lang/mpi2/lib64/libmp.a -L$MODIPSLDIR -lioipsl /usr/local/lib/lib64/libnetcdf_cc_31.a"
-     set oplink="-Wl,-t,-P,-dy "
-   else
-     set opt_link="-Wg,-c -L$MODIPSLDIR -lioipsl /usr/local/lib/lib64/libnetcdf_cc_31.a"
-     set oplink="-Wl,-t,-dy "
-   endif
-   if ($veget == true) then
-     set opt_link="$opt_link $link_veget -lioipsl /usr/local/lib/lib64/libnetcdf_cc_31.a"
-   endif
-#################
-else if $CRAY then
-#################
-   set f77=f90
-   set f90=f90
-#################
-else if $LINUX then
-#################
-#   set f77=pgf90
-#   set f90=pgf90
-   set f77=$FC_LINUX
-   set f90=$FC_LINUX
-   if ( $FC_LINUX == 'pgf90' ) then
-     if ( $io == "ioipsl" ) then
-       set opt_link=" -L$MODIPSLDIR $link_veget -L$NCDFLIB -lioipsl -lnetcdf "
-     else
-       set opt_link=" -L$MODIPSLDIR $link_veget -L$NCDFLIB -lnetcdf "
-     endif
-   else if ($FC_LINUX == 'g95' || $FC_LINUX == 'gfortran' ) then
-     if ( $io == "ioipsl" ) then
-       set opt_link="-L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf -lioipsl -lnetcdf "
-     else
-       set opt_link="-L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf -lnetcdf "
-     endif
-   else
-     set opt_link=" "
-   endif
-#################
-else if $SUN then
-#################
-   set f77=f90
-   set f90=f90
-   if ( $io == "ioipsl" ) then
-     set opt_link="-lf77compat -L$MODIPSLDIR $link_veget -lioipsl -L$NCDFLIB -lnetcdf "
-   else
-     set opt_link="-lf77compat -L$MODIPSLDIR $link_veget -L$NCDFLIB -lnetcdf "
-   endif
-#################
-else if $NEC then
-#################
-   set f77=f90 -ftrace
-   set f90=f90 -ftrace
-   set opt_link="-L$MODIPSLDIR"
-   if ($veget == true) then
-     set opt_link="$opt_link $link_veget"
-   endif
-   if ($couple == true) then
-     set opt_link="$opt_link -lioipsl -loasis2.4_mpi2 -float0 -ew -P static $NCDFLIB "
-   else
-     set opt_link="$opt_link -L$MODIPSLDIR -lioipsl -float0 -ew -P static $NCDFLIB "
-   endif
-   set mod_loc_dir="./"
-#################
-else if $XNEC then
-#################
-   set f77="sxmpif90 -ftrace"
-   set f90="sxmpif90 -ftrace"
-   if $MODIPSL then
-     set opt_link="-L$MODIPSLDIR"
-     if ($veget == true) then
-       set opt_link="$opt_link $link_veget"
-     endif
-     if ($couple == true) then
-       if ($psmile == true) then
-       set opt_link="$opt_link -lsxioipsl -float0 $optdbl -P static $NCDFLIB "
-       else
-       set opt_link="$opt_link -lsxioipsl -loasis2.4_mpi2 -float0 $optdbl -P static $NCDFLIB "
-       endif
-     else
-       set opt_link="$opt_link -lsxioipsl -float0 $optdbl -P static $NCDFLIB "
-     endif
-   else
-     if ($couple == true) then
-       set opt_link="-L$MODIPSLDIR"
-       set opt_link="$opt_link $link_veget -lsxioipsl -loasis2.4_mpi2 -float0 $optdbl -P static $NCDFLIB "
-     else
-       set opt_link=" -C hopt -float0 $optdbl -P static -L$MODIPSLDIR $link_veget -lsxioipsl $NCDFLIB "
-     endif
-   endif
-   set mod_loc_dir="./"
-##################
-else if $X6NEC then
-##################
-   set f77=sxmpif90
-   set f90=sxmpif90
-   if $MODIPSL then
-     set opt_link="$opt_link -L$MODIPSLDIR"
-     if ($veget == true) then
-       set opt_link="$opt_link $link_veget"
-     endif
-     if ($couple == true) then
-	if ($psmile == true) then
-	set opt_link="$opt_link -lsxioipsl -float0 -size_t64 $optdbl -P static $NCDFLIB "
-	else
-	set opt_link="$opt_link -lsxioipsl -loasis2.4_mpi2 -float0 -size_t64 $optdbl -P static $NCDFLIB "
-	endif
-     else
-       set opt_link="$opt_link -lsxioipsl -float0 -size_t64 $optdbl -P static $NCDFLIB "
-     endif
-   else
-#     set opt_link=" -float0 -size_t64 $optdbl -P static -L$MODIPSLDIR -lsxsechiba -lsxparameters -lsxstomate -lsxioipsl $NCDFLIB "
-     set opt_link=" $opt_link -float0 -size_t64 $optdbl -P static -L$MODIPSLDIR -lsxioipsl $NCDFLIB "
-
-   endif
-   set mod_loc_dir="./"
-##################
-else if $X8BRODIE then
-##################
-   set f77=sxmpif90
-   set f90=sxmpif90 
-   if $MODIPSL then
-     set opt_link="$opt_link -float0 -Wf,-A dbl4 -L$MODIPSLDIR -lblas"
-     if ($veget == true) then
-       set opt_link="$opt_link $link_veget"
-     endif
-     if ($couple == true) then
-       set opt_link="$opt_link -lsxioipsl -float0 $optdbl -P static $NCDFLIB "
-     else
-       set opt_link="$opt_link -lsxioipsl -float0 $optdbl -P static $NCDFLIB "
-     endif
-   else
-#     set opt_link=" -float0 $optdbl -P static -L$MODIPSLDIR -lsxsechiba -lsxparameters -lsxstomate -lsxioipsl $NCDFLIB "
-     set opt_link=" -float0 $optdbl -P static -L$MODIPSLDIR -lsxioipsl $NCDFLIB -lblas"
-
-   endif
-   set mod_loc_dir="./"
-#################
-else
-#################
-   set f77=f77
-   set f90=f90
-endif
-
-cd $model
-
-if $VPP then
-set make="gmake RANLIB=ls"
-else if $CRAY then
-set make="make RANLIB=ls"
-else if $NEC then
-set make="make RANLIB=ls"
-else if $LINUX then
-set make="make -k RANLIB=ranlib"
-else if $XNEC then
-set make="gmake RANLIB=ls"
-else if $X6NEC then
-set make="gmake RANLIB=ls"
-else if $X8BRODIE then
-set make="gmake RANLIB=ls"
-else
-set make="make RANLIB=ranlib"
-endif
-
-
-
-
-#
-# etat0_netcdf a besoin d'info de la physique
-# A revoir
-set include="$include"' -I$(LIBF)/phy'"$physique"
-#
-# le programme principal create_limit a besoin de l'info du module
-# startvar: on met donc libo dans les include pour Nec
-set include="$include"' -I$(LIBO)'
-
-
-#################################################################
-# Execution de la comande make... ENFIN!
-#################################################################
-
-if $VPP then
-  set optim90=" $optim90 -Am -M$libo"
-  set optimtru90="$optim90"
- \cp $IOIPSLDIR/*.mod $libo
-else if $SUN then
- set optim90=" $optim90 -M$libo -M$MODIPSLDIR "
- set optimtru90=" $optimtru90 -M$libo -M$MODIPSLDIR "
- set optim="$optim90"
- \cp $IOIPSLDIR/*.mod $libo
-else if $NEC then
- set optim90=" $optim90 -I$libo "
-else if $XNEC then
- set optim90=" $optim90 -I$libo "
- set optimtru90=" $optimtru90 -I$libo "
-else if $X6NEC then
- set optim90=" $optim90 -I$libo "
- set optimtru90=" $optimtru90 -I$libo "
-else if $X8BRODIE then
- set optim90=" $optim90 -I$libo "
- set optimtru90=" $optimtru90 -I$libo "
-else if $LINUX then
- if ( $FC_LINUX == "pgf90" ) then
-   set optimtru90=" $optimtru90 -module $libo "
-   set optim90=" $optim90 -module $libo "
- else if ( $FC_LINUX == 'g95' ) then
-   set optimtru90=" $optimtru90 -fmod=$libo  "
-   set optim90=" $optim90 -fmod=$libo  "
- else if ( $FC_LINUX == 'gfortran' ) then
-   set optimtru90=" $optimtru90 -M $libo  "
-   set optim90=" $optim90 -M $libo  "
- endif
- set optim="$optim90"
- set mod_loc_dir=$libo
-# \cp /d3/fairhead/sechiba/ioipsl/*.mod $libo
-# \cp $IOIPSLDIR/*.mod $libo
-endif
-
-set link="$f90 $optim90"
-
-set ar=ar
-
-if $XNEC then
-  set link="sxld $opt_link"
-  set link="$f90 "
-#  set ar=sxar
-else if $X6NEC then
-  set link="sxld $opt_link"
-  set link="$f90 -Wl,-hlib_cyclic "
-else if $X8BRODIE then
-  set link="sxld $opt_link"
-  set link="$f90 -Wl,-hlib_cyclic "
-endif
-
-
-cd $localdir
-
-set source_code=${code}.F
-if ( -f $LMDGCM/libf/dyn${dimc}d${FLAG_PARA}/${code}.F90 ) then
-  set source_code=${code}.F90
-endif
-
-echo $make -f $LMDGCM/makefile \
-OPTION_DEP="$opt_dep" OPTION_LINK="$opt_link" \
-OPTIM90="$optim90" \
-OPTIMTRU90="$optimtru90" \
-OPTIM="$optim$optimbis" \
-INCLUDE="$include" \
-$filtre \
-LIBO=$libo \
-$dyn \
-$phys \
-DIM=$dimc \
-FLAG_PARA="$FLAG_PARA"\
-L_ADJNT="$adjnt" \
-L_COSP="$lcosp" \
-L_CHIMIE="$libchimie" \
-LOCAL_DIR="$localdir"  \
-F77="$f77" \
-F90="$f90" \
-OPLINK="$oplink" \
-LINK="$link" \
-GCM="$LMDGCM" \
-MOD_LOC_DIR=$mod_loc_dir \
-MOD_SUFFIX=$mod_suffix \
-AR=$ar \
-SOURCE=$source_code \
-PROG=$code
-
-$make -f $LMDGCM/makefile \
-OPTION_DEP="$opt_dep" OPTION_LINK="$opt_link" \
-OPTIM90="$optim90" \
-OPTIMTRU90="$optimtru90" \
-OPTIM="$optim$optimbis" \
-INCLUDE="$include" \
-$filtre \
-LIBO=$libo \
-$dyn \
-$phys \
-DIM=$dimc \
-FLAG_PARA="$FLAG_PARA"\
-L_ADJNT="$adjnt" \
-L_COSP="$lcosp" \
-L_CHIMIE="$libchimie" \
-LOCAL_DIR="$localdir"  \
-F77="$f77" \
-F90="$f90" \
-OPLINK="$oplink" \
-LINK="$link" \
-GCM="$LMDGCM" \
-MOD_LOC_DIR=$mod_loc_dir \
-MOD_SUFFIX=$mod_suffix \
-AR=$ar \
-SOURCE=$source_code \
-PROG=$code
-
-\rm -f $libf/grid/dimensions.h
Index: LMDZ4/trunk/makelmdz_fcm
===================================================================
--- LMDZ4/trunk/makelmdz_fcm	(revision 5799)
+++ 	(revision )
@@ -1,411 +1,0 @@
-#!/bin/bash
-# $Id$
-# This is a script in Bash.
-
-# FH : on ne crée plus le fichier arch.mk qui est supposé exister par
-# FH : ailleurs.
-# FH : ulterieurement, ce fichier sera pré-existant pour une série
-# FH : de configurations en versions optimisées et debug qui seront
-# FH : liés (ln -s) avec arch.mk en fonction de l'architecture.
-# FH : Pour le moment, cette version est en test et on peut créer les
-# FH : arch.mk en lançant une première fois makegcm.
-#
-##set -x
-########################################################################
-# options par defaut pour la commande make
-########################################################################
-
-dim="96x72x19"
-physique=lmd
-filtre=filtrez
-grille=reg
-couple=false
-veget=false
-chimie=false
-parallel=none
-compil_mod=prod
-io=ioipsl
-LIBPREFIX=""
-fcm_path=none
-cosp=false
-
-LMDGCM=`/bin/pwd`
-LIBOGCM=$LMDGCM/libo
-LIBFGCM=$LMDGCM/libf
-COSP_PATH=$LMDGCM/.void_dir
-
-########################################################################
-#  Quelques initialisations de variables du shell.
-########################################################################
-
-CPP_KEY="" 
-INCLUDE=""
-LIB=""
-adjnt=""
-COMPIL_FFLAGS="%PROD_FFLAGS"
-PARA_FFLAGS=""
-PARA_LD=""
-EXT_SRC=""
-
-########################################################################
-# lecture des options de mymake
-########################################################################
-
-while (($# > 0))
-  do
-  case $1 in
-      "-h") cat <<fin
-manuel complet sur http://...
-Usage :
-makegcm [options] -m arch exec
-[-h]                       : manuel abrÃ©gÃ©
-[-d [[IMx]JMx]LM]          : IM, JM, LM sont les dims en x, y, z (def: $dim)
-[-p PHYS]                  : compilation avec la physique libf/phyPHYS, (def: lmd)
-[-prod / -dev / -debug]    : compilation en mode production (default) / developpement / debug .
-[-c false/MPI1/MPI2]       : couplÃ© ocÃ©an : MPI1/MPI2/false (def: false)
-[-v false/true]            : avec ou sans vÃ©gÃ©tation (def: false)
-[-chimie INCA/false]       : avec ou sans model de chimie INCA (def: false)
-[-parallel none/mpi/omp/mpi_omp] : parallelisation (default: none) : mpi, openmp ou mixte mpi_openmp
-[-g GRI]                   : conf. grille dans dyn3d/GRI_xy.h  (def: reg inclue un zoom)
-[-io IO]                   : choix d'une librairie I/O, experts (def: ioipsl)
-[-include INCLUDES]        : variables supplementaires pour include
-[-cpp CPP_KEY]             : cle cpp supplementaires
-[-adjnt]                   : adjoint, a remettre en route ...
-[-filtre NOMFILTRE]        : prend le filtre dans libf/NOMFILTRE (def: filtrez)
-[-link LINKS]              : liens optionels avec d'autres librairies
-[-fcm_path path]           : chemin pour fcm (def: le chemin est suppose deja exister dans le PATH)
-[-ext_src path]            : chemin d'un repertoire source avec des sources externe a compiler avec le modele
- -arch nom_arch            : nom de l'architecture cible
- exec                      : exÃ©cutable gÃ©nÃ©rÃ©
-fin
-	  exit;;
-
-      "-d")
-	  dim=$2 ; shift ; shift ;;
-      
-      "-O")
-	  echo "option obsolete dans cette version intermediaire de makegcm"
-	  exit;;
-
-      "-p")
-	  physique="$2" ;  shift ; shift ;;
-
-      "-g")
-	  grille="$2" ; shift ; shift ;;
-
-      "-c")
-	  couple="$2" ; shift ; shift ;;
-
-      "-prod")
-	  compil_mod="prod" ; shift ;;
-
-      "-dev")
-	  compil_mod="dev" ; shift ;;
-
-      "-debug")
-	  compil_mod="debug" ; shift ;;
-
-      "-io")
-	  io="$2" ; shift ; shift ;;
-
-      "-v")
-	  veget="$2" ; shift ; shift ;;
-
-      "-chimie")
-	  chimie="$2" ; shift ; shift ;;
-
-      "-parallel")
-	  parallel="$2" ; shift ; shift ;;
-      
-      "-include")
-	  INCLUDE="$INCLUDE -I$2" ; shift ; shift ;;
-
-      "-cpp")
-	  CPP_KEY="$CPP_KEY $2" ; shift ; shift ;;
-
-      "-adjnt")
-	  echo "otpion a reactiver ";exit
-	  opt_dep="$opt_dep adjnt" ; adjnt="-ladjnt -ldyn3d "
-	  optim="$optim -Dadj" ; shift ;;
-
-      "-cosp")
-          cosp="$2" ; shift ; shift ;;
-
-
-      "-filtre")
-	  filtre=$2 ; shift ; shift ;;
-
-      "-link")
-	  LIB="$LIB $2" ; shift ; shift ;;
-
-      "-fcm_path")
-	  fcm_path=$2 ; shift ; shift ;;
-
-      "-ext_src")
-	  EXT_SRC=$2 ; shift ; shift ;;
-
-      "-arch")
-	  arch=$2 ; shift ; shift ;;
-
-      *)
-	  code="$1" ; shift ;;
-  esac
-done
-
-###############################################################
-# mettre le chemin du fcm dans le path
-###############################################################
-if  [[ "$fcm_path" != "none" ]]
-then
-    export PATH=${fcm_path}:${PATH}
-fi
-
-echo "Chemin du fcm utlise :" 
-which fcm
-
-###############################################################
-# lecture des chemins propres à l'architecture de la machine #
-###############################################################
-rm -f .void_file
-echo > .void_file
-rm -rf .void_dir
-mkdir .void_dir
-rm -f arch.path
-ln -s arch/arch-${arch}.path ./arch.path
-source arch.path
-
-########################################################################
-# Definition des clefs CPP, des chemins des includes et modules
-#  et des libraries
-########################################################################
-
-if [[ "$compil_mod" == "prod" ]]
-then
-  COMPIL_FFLAGS="%PROD_FFLAGS"
-elif [[ "$compil_mod" == "dev" ]]
-then
-  COMPIL_FFLAGS="%DEV_FFLAGS"
-elif [[ "$compil_mod" == "debug" ]]
-then
-  COMPIL_FFLAGS="%DEBUG_FFLAGS"
-fi
-
-if [[ "$physique" != "nophys" ]]
-then
-   #Default planet type is Earth
-   CPP_KEY="$CPP_KEY CPP_EARTH"
-fi
-
-if [[ "$chimie" == "INCA" ]]
-then
-   CPP_KEY="$CPP_KEY INCA"
-   INCLUDE="$INCLUDE -I${INCA_INCDIR}"
-   LIB="$LIB -L${INCA_LIBDIR} -lchimie"
-fi
-
-if [[ "$couple" != "false" ]]
-then
-   CPP_KEY="$CPP_KEY CPP_COUPLE"
-   INCLUDE="$INCLUDE -I${OASIS_INCDIR}"
-   LIB="$LIB -L${OASIS_LIBDIR} -lpsmile.${couple} -lmpp_io"
-fi
-
-if [[ "$parallel" == "mpi" ]]
-then
-   CPP_KEY="$CPP_KEY CPP_PARA CPP_MPI"
-   PARA_FFLAGS="%MPI_FFLAGS"
-   PARA_LD="%MPI_LD"
-elif [[ "$parallel" == "omp" ]]
-then
-   CPP_KEY="$CPP_KEY CPP_PARA CPP_OMP"
-   PARA_FFLAGS="%OMP_FFLAGS"
-   PARA_LD="%OMP_LD"
-elif [[ "$parallel" == "mpi_omp" ]]
-then
-   CPP_KEY="$CPP_KEY CPP_PARA CPP_MPI CPP_OMP"
-   PARA_FFLAGS="%MPI_FFLAGS %OMP_FFLAGS"
-   PARA_LD="%MPI_LD %OMP_LD"
-fi
-
-if [[ ( "$parallel" == "omp" || "$parallel" == "mpi_omp" ) \
-   && "$compil_mod" == "debug" ]]
-then
-    echo "Usually, parallelization with OpenMP requires some optimization."
-    echo "We suggest switching to \"-dev\"."
-fi
-
-if [[ "$veget" == "true" ]]
-then
-   CPP_KEY="$CPP_KEY CPP_VEGET"
-   INCLUDE="${INCLUDE} -I${ORCH_INCDIR}"
-   LIB="${LIB} -L${ORCH_LIBDIR} -l${LIBPREFIX}sechiba -l${LIBPREFIX}parameters -l${LIBPREFIX}stomate -l${LIBPREFIX}parallel -l${LIBPREFIX}orglob"
-fi
-
-if [[ $io == ioipsl ]]
-then
-   CPP_KEY="$CPP_KEY CPP_IOIPSL"
-   INCLUDE="$INCLUDE -I${IOIPSL_INCDIR}"
-   LIB="$LIB -L${IOIPSL_LIBDIR} -l${LIBPREFIX}ioipsl"
-fi
-if [[ "$cosp" == "true" ]]
-then
-   CPP_KEY="$CPP_KEY CPP_COSP"
-   COSP_PATH="$LIBFGCM/cosp"
-#   LIB="${LIB} -l${LIBPREFIX}cosp"
-fi
-INCLUDE="$INCLUDE -I${NETCDF_INCDIR}"
-LIB="$LIB -L${NETCDF_LIBDIR} -lnetcdf"
-
-########################################################################
-# calcul du nombre de dimensions
-########################################################################
-
-
-dim_full=$dim
-dim=`echo $dim | sed -e 's/[^0-9]/ /g'` 
-set $dim
-dimc=$#
-echo calcul de la dimension
-echo dim $dim
-echo dimc $dimc
-
-
-########################################################################
-# Gestion des dimensions du modele.
-# on cree ou remplace le fichier des dimensions
-########################################################################
-
-cd $LIBFGCM/grid/dimension
-./makdim $dim
-cat $LIBFGCM/grid/dimensions.h
-cd $LMDGCM
-
-
-########################################################################
-# Differentes dynamiques (3d, 2d, 1d)
-########################################################################
-
-dimension=`echo $dim | wc -w`
-echo dimension $dimension
-
-if (( $dimension == 3 ))
-then
-  cd $LIBFGCM/grid
-  \rm fxyprim.h
-  cp -p fxy_${grille}.h fxyprim.h
-else
-  echo "Probleme dans les dimensions de la dynamique !!"
-  echo "Non reactive pour l'instant !!!"
-fi
-
-######################################################################
-#   Traitement special pour le nouveau rayonnement de Laurent Li.
-#   ---> YM desactive pour le traitemement en parallele
-######################################################################
-
-#if [[ -f $libf/phy$physique/raddim.h ]]
-#then
-# if [[ -f $libf/phy$physique/raddim.$dimh.h ]]
-#then
-#  \rm -f $libf/phy$physique/raddim.h
-#  cp -p $libf/phy$physique/raddim.$dimh.h $libf/phy$physique/raddim.h
-#  echo $libf/phy$physique/raddim.$dimh.h 
-#  cat $libf/phy$physique/raddim.h
-# else
-#  echo On peut diminuer la taille de l executable en creant
-#  echo le fichier $libf/phy$physique/raddim.$dimh.h
-#  \cp -p $libf/phy$physique/raddim.defaut.h $libf/phy$physique/raddim.h
-# fi
-#fi
-
-######################################################################
-# Gestion du filtre qui n'existe qu'en 3d.
-######################################################################
-
-if (( `expr $dimc \> 2` == 1 ))
-then
-   filtre="FILTRE=$filtre"
-else
-   filtre="FILTRE= L_FILTRE= "
-fi
-echo MACRO FILTRE $filtre
-
-echo $dimc
-
-
-
-######################################################################
-# Creation du suffixe de la configuration
-######################################################################
-
-
-SUFF_NAME=_${dim_full}
-SUFF_NAME=${SUFF_NAME}_phy${physique}
-
-if [[ "$parallel" != "none" ]]
-then
-  SUFF_NAME=${SUFF_NAME}_para
-  DYN=dyn${dimc}dpar
-else
-  SUFF_NAME=${SUFF_NAME}_seq
-  DYN=dyn${dimc}d
-fi
-
-if [[ $veget == "true" ]]
-then
-  SUFF_NAME=${SUFF_NAME}_orch
-fi
-
-if [[ $couple != "false" ]]
-then
-  SUFF_NAME=${SUFF_NAME}_couple
-fi
-
-if [[ $chimie == "INCA" ]]
-then
-  SUFF_NAME=${SUFF_NAME}_inca
-fi
-
-cd $LMDGCM
-config_fcm="config.fcm"
-rm -f $config_fcm
-touch $config_fcm
-rm -f bin/${code}${SUFF_NAME}.e
-rm -f arch.fcm
-rm -f arch.opt
-
-echo "%ARCH          $arch"          >> $config_fcm
-echo "%INCDIR        $INCLUDE"       >> $config_fcm 
-echo "%LIB           $LIB"           >> $config_fcm
-echo "%ROOT_PATH     $PWD"           >> $config_fcm
-echo "%LIBF          $LIBFGCM"       >> $config_fcm
-echo "%LIBO          $LIBOGCM"       >> $config_fcm
-echo "%DYN           $DYN"           >> $config_fcm
-echo "%PHYS          phy${physique}" >> $config_fcm
-echo "%COSP          $COSP_PATH"     >> $config_fcm
-echo "%CPP_KEY       $CPP_KEY"       >> $config_fcm
-echo "%EXEC          $code"          >> $config_fcm
-echo "%SUFF_NAME     $SUFF_NAME"     >> $config_fcm
-echo "%COMPIL_FFLAGS $COMPIL_FFLAGS" >> $config_fcm
-echo "%PARA_FFLAGS   $PARA_FFLAGS"   >> $config_fcm
-echo "%PARA_LD       $PARA_LD"       >> $config_fcm
-echo "%EXT_SRC       $EXT_SRC"       >> $config_fcm
-
-
-
-ln -s arch/arch-${arch}.fcm arch.fcm
-if test -f arch/arch-${arch}.opt &&  [ $compil_mod = "prod" ]
-  then
-  ln -s arch/arch-${arch}.opt arch.opt
-else
-  ln -s .void_file arch.opt
-fi
-
-
-rm -f $LIBOGCM/${arch}${SUFF_NAME}/.config/fcm.bld.lock
-./build_gcm
-
-rm -rf tmp_src
-rm -rf config
-ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config config
-ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config/tmp tmp_src
Index: LMDZ4/trunk/offline.def
===================================================================
--- LMDZ4/trunk/offline.def	(revision 5799)
+++ 	(revision )
@@ -1,12 +1,0 @@
-#
-# $Header$
-#
-T
-4
-T
--2.
-48.1
-1
-T
-6
-2
Index: LMDZ4/trunk/orchidee.def
===================================================================
--- LMDZ4/trunk/orchidee.def	(revision 5799)
+++ 	(revision )
@@ -1,63 +1,0 @@
-#
-#
-# Parameter file for LMDZ4OR_v2 configuration
-# See comments : http://forge.ipsl.jussieu.fr/orchidee/
-#
-STOMATE_OK_CO2=TRUE
-# STOMATE_OK_STOMATE is not set
-# STOMATE_OK_DGVM is not set
-# STOMATE_WATCHOUT is not set
-SECHIBA_restart_in=default
-SECHIBA_rest_out=sechiba_rest.nc
-SECHIBA_reset_time=y
-#
-OUTPUT_FILE=sechiba_out.nc
-WRITE_STEP=2592000
-SECHIBA_HISTLEVEL=5
-#
-SECHIBA_HISTFILE2 = FALSE
-SECHIBA_OUTPUT_FILE2 = sechiba_out_2.nc
-WRITE_STEP2 = 86400.0
-SECHIBA_HISTLEVEL2 = 1
-#
-STOMATE_OUTPUT_FILE=stomate_history.nc
-STOMATE_HIST_DT=10.
-STOMATE_HISTLEVEL=0
-SECHIBA_DAY=0.0
-SECHIBA_ZCANOP=0.5
-DT_SLOW=86400.
-# IMPOSE_VEG is not set
-VEGETATION_FILE=carteveg5km.nc
-# VEGETATION_FILE=pft_new.nc
-DIFFUCO_LEAFCI=233.
-CONDVEG_SNOWA=default
-# IMPOSE_AZE is not set
-SOILALB_FILE=soils_param.nc
-SOILTYPE_FILE=soils_param.nc 
-ENERBIL_TSURF=280.
-HYDROL_SNOW=0.0
-HYDROL_SNOWAGE=0.0
-HYDROL_SNOWICE=0.0
-HYDROL_SNOWICEAGE=0.0
-HYDROL_HDRY=1.0
-HYDROL_HUMR=1.0
-HYDROL_BQSB=default
-HYDROL_GQSB=0.0
-HYDROL_DSG=0.0
-HYDROL_DSP=default
-HYDROL_QSV=0.0
-HYDROL_OK_HDIFF=n
-HYDROL_TAU_HDIFF=1800.
-THERMOSOIL_TPRO=280.
-RIVER_ROUTING=y
-ROUTING_FILE=routing.nc
-LAI_MAP=y
-LAI_FILE=lai2D.nc
-SECHIBA_QSINT=0.02
-ALB_BARE_MODEL = FALSE
-PERCENT_THROUGHFALL_PFT = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30.
-RVEG_PFT = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5
-CDRAG_FROM_GCM = .TRUE.
-#LAND_USE=y
-#VEGET_YEAR=0
-#VEGET_UPDATE=1Y
Index: LMDZ4/trunk/output.def
===================================================================
--- LMDZ4/trunk/output.def	(revision 5799)
+++ 	(revision )
@@ -1,794 +1,0 @@
-######## Abderrahmane le 20 11 08 ########################
-# Niveaux de sorties et nom pour chaque variable dans les#
-# fichiers histmth histday histhf histins histLES        #
-##########################################################
-# Surface geop.height
-flag_phis         =  1, 1, 10, 1, 1    
-name_phis         =  phis
-# Grid area
-flag_aire         =  1, 1, 10,  1, 1    
-name_aire         =  aire
-# Surfac ter+lic
-flag_contfracATM  =  10, 1,  1, 10, 10    
-name_contfracATM  =  contfracATM
-# Surface terre OR  
-flag_contfracOR   =  10, 1,  1, 10, 10    
-name_contfracOR   =  contfracOR
-# Grid area CONT
-flag_aireTER      =  10, 10, 1, 10, 10    
-name_aireTER      =  aireTER
-# Latent heat flux
-flag_flat         =  10, 1, 10, 10, 1    
-name_flat         =  flat
-# Sea Level Pressure
-flag_slp          =  1, 1, 1, 10, 1    
-name_slp          =  slp
-# Surface Temperature
-flag_tsol         =  1, 1, 1, 1, 1    
-name_tsol         =  tsol
-# Temperature 2m
-flag_t2m          =  1, 1, 1, 1, 1    
-name_t2m          =  t2m
-# Temperature min 2m
-flag_t2m_min      =  1, 1, 10, 10, 10    
-name_t2m_min      =  t2m_min
-# Temperature max 2m
-flag_t2m_max      =  1, 1, 10, 10, 10    
-name_t2m_max      =  t2m_max
-# Temperature ter lic oce sic at 2m
-flag_t2m_ter      =  10, 4, 10, 10, 10    
-name_t2m_ter      =  t2m_ter 
-flag_t2m_lic      =  10, 4, 10, 10, 10    
-name_t2m_lic      =  t2m_lic
-flag_t2m_oce      =  10, 4, 10, 10, 10    
-name_t2m_oce      =  t2m_oce
-flag_t2m_sic      =  10, 4, 10, 10, 10    
-name_t2m_sic      =  t2m_sic 
-# 10-m wind speed
-flag_wind10m      =  1, 1, 1, 10, 10    
-name_wind10m      =  wind10m
-# 10-m wind speed max
-flag_wind10max    =  10, 1, 10, 10, 10    
-name_wind10max    =  wind10max
-# Sea-ice fraction
-flag_sicf         =  1, 1, 10, 10, 10    
-name_sicf         =  sicf
-# Specific humidity 2m
-flag_q2m          =  1, 1, 1, 1, 1    
-name_q2m          =  q2m
-# 10m zonal wind
-flag_u10m         =  1, 1, 1, 1, 1    
-name_u10m         =  u10m
-# 10m meridional wind
-flag_v10m         =  1, 1, 1, 1, 1    
-name_v10m         =  v10m
-# Surface Pressure
-flag_psol         =  1, 1, 1, 1, 1    
-name_psol         =  psol
-# Surface Air humidity
-flag_qsurf        =  1, 10, 10, 10, 10    
-name_qsurf        =  qsurf
-# 10m zonal wind (ter lic oce sic)
-flag_u10m_ter     =  10, 4, 10, 10, 10    
-name_u10m_ter     =  u10m_ter
-flag_u10m_lic     =  10, 4, 10, 10, 10    
-name_u10m_lic     =  u10m_lic
-flag_u10m_oce     =  10, 4, 10, 10, 10    
-name_u10m_oce     =  u10m_oce
-flag_u10m_sic     =  10, 4, 10, 10, 10    
-name_u10m_sic     =  u10m_sic
-# 10m meridien  wind (ter oce ice lic)
-flag_v10m_ter     =  10, 4, 10, 10, 10    
-name_v10m_ter     =  v10m_ter
-flag_v10m_lic     =  10, 4, 10, 10, 10    
-name_v10m_lic     =  v10m_lic
-flag_v10m_oce     =  10, 4, 10, 10, 10    
-name_v10m_oce     =  v10m_oce
-flag_v10m_sic     =  10, 4, 10, 10, 10    
-name_v10m_sic     =  v10m_sic
-# Soil watter content
-flag_qsol         =  1, 10, 10, 1, 1    
-name_qsol         =  qsol
-# Number of dayrain(liq+sol)
-flag_ndayrain     =  1, 10, 10, 10, 10    
-name_ndayrain     =  ndayrain
-# Precip Totale liq+sol
-flag_precip       =  1, 1, 1, 1, 1    
-name_precip       =  precip
-# Large-scale Precip
-flag_plul         =  1, 1, 1, 1, 10    
-name_plul         =  plul 
-# Convective Precip
-flag_pluc         =  1, 1, 1, 1, 10    
-name_pluc         =  pluc
-# Snow fall
-flag_snow         =  1, 1, 10, 1, 10    
-name_snow         =  snow
-# Evaporation
-flag_evap         =  1, 1, 10, 1, 10    
-name_evap         =  evap
-# Solar rad. at TOA
-flag_tops         =  1, 1, 10, 10, 10    
-name_tops         =  tops
-# CS Solar rad. at TOA
-flag_tops0        =  1, 5, 10, 10, 10    
-name_tops0        =  tops0
-# IR rad. at TOA
-flag_topl         =  1, 1, 10, 1, 10    
-name_topl         =  topl
-# CR IR rad. at TOA
-flag_topl0        =  1, 5, 10, 10, 10    
-name_topl0        =  topl0
-# SWup at TOA
-flag_SWupTOA      =  1, 4, 10, 10, 10    
-name_SWupTOA      =  SWupTOA
-# CR SWup at TOA
-flag_SWupTOAclr   =  1, 4, 10, 10, 10    
-name_SWupTOAclr   =  SWupTOAclr
-# SWdn at TOA
-flag_SWdnTOA      =  1, 4, 10, 10, 10    
-name_SWdnTOA      =  SWdnTOA
-# CR SWdn at TOA
-flag_SWdnTOAclr   =  1, 4, 10, 10, 10    
-name_SWdnTOAclr   =  SWdnTOAclr
-# SWup at 200hPa
-flag_SWup200      =  1, 10, 10, 10, 10    
-name_SWup200      =  SWup200
-# CR SWup at 200hPa
-flag_SWup200clr   =  10, 1, 10, 10, 10    
-name_SWup200clr   =  SWup200clr
-# SWdn at 200hPa
-flag_SWdn200      =  1, 10, 10, 10, 10    
-name_SWdn200      =  SWdn200
-# CR SWdn at 200hPa
-flag_SWdn200clr   =  10, 1, 10, 10, 10    
-name_SWdn200clr   =  SWdn200clr
-# LWup at 200mb
-flag_LWup200      =  1, 10, 10, 10, 10    
-name_LWup200      =  LWup200
-# CR LWup at 200mb
-flag_LWup200clr   =  1, 10, 10, 10, 10    
-name_LWup200clr   =  LWup200clr
-# LWdn at 200mb
-flag_LWdn200      =  1, 10, 10, 10, 10    
-name_LWdn200      =  LWdn200
-# CR LWdn at 200mb
-flag_LWdn200clr   =  1, 10, 10, 10, 10    
-name_LWdn200clr   =  LWdn200clr
-# Solar rad. at surf
-flag_sols         =  1, 1, 10, 1, 10    
-name_sols         =  sols
-# CR Solar rad. at surf
-flag_sols0        =  1, 5, 10, 10, 10    
-name_sols0        =  sols0
-# IR rad. at surface
-flag_soll         =  1, 1, 10, 1, 10    
-name_soll         =  soll
-# CR IR rad. at surface
-flag_soll0        =  1, 5, 10, 10, 10    
-name_soll0        =  soll0
-# Rayonnement au sol
-flag_radsol       =  1, 1, 10, 10, 10    
-name_radsol       =  radsol
-# SWup at surface
-flag_SWupSFC      =  1, 4, 10, 10, 10    
-name_SWupSFC      =  SWupSFC
-# CR SWup at surface
-flag_SWupSFCclr   =  1, 4, 10, 10, 10    
-name_SWupSFCclr   =  SWupSFCclr
-# SWdn at surface
-flag_SWdnSFC      =  1, 1, 10, 10, 10    
-name_SWdnSFC      =  SWdnSFC
-# CR at surface
-flag_SWdnSFCclr   =  1, 4, 10, 10, 10    
-name_SWdnSFCclr   =  SWdnSFCclr
-# LWup at surface
-flag_LWupSFC      =  1, 4, 10, 10, 10    
-name_LWupSFC      =  LWupSFC
-# CR LWup at surface
-flag_LWupSFCclr   =  1, 4, 10, 10, 10    
-name_LWupSFCclr   =  LWupSFCclr
-# LWdn  at surface
-flag_LWdnSFC      =  1, 4, 10, 10, 10    
-name_LWdnSFC      =  LWdnSFC
-# CR LWdn  at surface
-flag_LWdnSFCclr   =  1, 4, 10, 10, 10    
-name_LWdnSFCclr   =  LWdnSFCclr
-# Surf. total heat flux
-flag_bils         =  1, 2, 10, 1, 10    
-name_bils         =  bils
-# Sensible heat flux
-flag_sens         =  1, 1, 10, 1, 1    
-name_sens         =  sens
-# Heat flux derivation
-flag_fder         =  1, 2, 10, 1, 10    
-name_fder         =  fder
-# Thermal flux for snow melting
-flag_ffonte       =  1, 10, 10, 10, 10    
-name_ffonte       =  ffonte
-# Ice Calving
-flag_fqcalving    =  1, 10, 10, 10, 10    
-name_fqcalving    =  fqcalving
-# Land ice melt
-flag_fqfonte      =  1, 10, 10, 10, 10    
-name_fqfonte      =  fqfonte
-# Zonal wind stress (ter ice liq oce)
-flag_taux_ter     =  1, 4, 10, 1, 10    
-name_taux_ter     =  taux_ter
-flag_taux_lic     =  1, 4, 10, 1, 10    
-name_taux_lic     =  taux_lic                                          
-flag_taux_oce     =  1, 4, 10, 1, 10    
-name_taux_oce     =  taux_oce
-flag_taux_sic     =  1, 4, 10, 1, 10    
-name_taux_sic     =  taux_sic
-# Meridien wind stress (ter ice liq oce)
-flag_tauy_ter     =  1, 4, 10, 1, 10    
-name_tauy_ter     =  tauy_ter
-flag_tauy_lic     =  1, 4, 10, 1, 10    
-name_tauy_lic     =  tauy_lic
-flag_tauy_oce     =  1, 4, 10, 1, 10    
-name_tauy_oce     =  tauy_oce
-flag_tauy_sic     =  1, 4, 10, 1, 10    
-name_tauy_sic     =  tauy_sic
-# % surface (ter ice liq oce)
-flag_pourc_ter    =  1, 4, 10, 1, 10    
-name_pourc_ter    =  pourc_ter  
-flag_pourc_lic    =  1, 4, 10, 1, 10    
-name_pourc_lic    =  pourc_lic
-flag_pourc_oce    =  1, 4, 10, 1, 10    
-name_pourc_oce    =  pourc_oce
-flag_pourc_sic    =  1, 4, 10, 1, 10    
-name_pourc_sic    =  pourc_sic
-# Fraction (ter ice liq oce)
-flag_fract_ter    =  1, 4, 10, 1, 10    
-name_fract_ter    =  fract_ter
-flag_fract_lic    =  1, 4, 10, 1, 10    
-name_fract_lic    =  fract_lic
-flag_fract_oce     =  1, 4, 10, 1, 10    
-name_fract_oce    =  fract_oce
-flag_fract_sic    =  1, 4, 10, 1, 10    
-name_fract_sic    =  fract_sic
-# Surface temperature (ter ice liq oce)
-flag_tsol_ter     =  1, 4, 10, 1, 10    
-name_tsol_ter     =  tsol_ter
-flag_tsol_lic     =  1, 4, 10, 1, 10    
-name_tsol_lic     =  tsol_lic
-flag_tsol_oce     =  1, 4, 10, 1, 10    
-name_tsol_oce     =  tsol_oce
-flag_tsol_sic     =  1, 4, 10, 1, 10    
-name_tsol_sic     =  tsol_sic
-# Sensible heat flux (ter ice liq oce)
-flag_sens_ter     =  1, 4, 10, 1, 10    
-name_sens_ter     =  sens_ter
-flag_sens_lic     =  1, 4, 10, 1, 10    
-name_sens_lic     =  sens_lic
-flag_sens_oce     =  1, 4, 10, 1, 10    
-name_sens_oce     =  sens_oce
-flag_sens_sic     =  1, 4, 10, 1, 10    
-name_sens_sic     =  sens_sic
-# Latent heat flux (ter ice liq oce)
-flag_lat_ter      =  1, 4, 10, 1, 10    
-name_lat_ter      =  lat_ter 
-flag_lat_lic      =  1, 4, 10, 1, 10    
-name_lat_lic      =  lat_lic
-flag_lat_oce      =  1, 4, 10, 1, 10    
-name_lat_oce      =  lat_oce
-flag_lat_sic      =  1, 4, 10, 1, 10    
-name_lat_sic      =  lat_sic 
-# LW (ter ice liq oce)
-flag_flw_ter      =  1, 10, 10, 10, 10    
-name_flw_ter      =  flw_ter
-flag_flw_lic      =  1, 10, 10, 10, 10    
-name_flw_lic      =  flw_lic
-flag_flw_oce      =  1, 10, 10, 10, 10    
-name_flw_oce      =  flw_oce
-flag_flw_sic      =  1, 10, 10, 10, 10    
-name_flw_sic      =  flw_sic
-# SW (ter ice liq oce)
-flag_fsw_ter      =  1, 10, 10, 10, 10    
-name_fsw_ter      =  fsw_ter
-flag_fsw_lic      =  1, 10, 10, 10, 10    
-name_fsw_lic      =  fsw_lic
-flag_fsw_oce      =  1, 10, 10, 10, 10    
-name_fsw_oce      =  fsw_oce
-flag_fsw_sic      =  1, 10, 10, 10, 10    
-name_fsw_sic      =  fsw_sic
-# Bilan sol (ter ice liq oce)
-flag_wbils_ter    =  1, 10, 10, 10, 10    
-name_wbils_ter    =  wbils_ter
-flag_wbils_lic    =  1, 10, 10, 10, 10    
-name_wbils_lic    =  wbils_lic
-flag_wbils_oce    =  1, 10, 10, 10, 10    
-name_wbils_oce    =  wbils_oce
-flag_wbils_sic   =  1, 10, 10, 10, 10    
-name_wbils_sic    =  wbils_sic 
-# Bilan eau (ter ice liq oce)
-flag_wbilo_ter    =  1, 10, 10, 10, 10    
-name_wbilo_ter    =  wbilo_ter
-flag_wbilo_lic    =  1, 10, 10, 10, 10    
-name_wbilo_lic    =  wbilo_lic
-flag_wbilo_oce    =  1, 10, 10, 10, 10    
-name_wbilo_oce    =  wbilo_oce
-flag_wbilo_sic    =  1, 10, 10, 10, 10    
-name_wbilo_sic    =  wbilo_sic
-# Momentum drag coef
-flag_cdrm         =  1, 10, 10, 1, 10    
-name_cdrm         =  cdrm
-# Heat drag coef
-flag_cdrh         =  1, 10, 10, 1, 10    
-name_cdrh         =  cdrh 
-# Low-level cloudiness
-flag_cldl         =  1, 1, 10, 10, 10    
-name_cldl         =  cldl
-# Mid-level cloudiness
-flag_cldm         =  1, 1, 10, 10, 10    
-name_cldm         =  cldm
-# High-level cloudiness
-flag_cldh         =  1, 1, 10, 10, 10    
-name_cldh         =  cldh
-# Total cloudiness
-flag_cldt         =  1, 1, 2, 10, 10    
-name_cldt         =  cldt
-# Cloud liquid water path
-flag_cldq         =  1, 1, 10, 10, 10    
-name_cldq         =  cldq
-# Cloud water path
-flag_lwp          =  1, 5, 10, 10, 10    
-name_lwp          =  lwp
-# Cloud ice water path
-flag_iwp          =  1, 5, 10, 10, 10    
-name_iwp          =  iwp
-# Zonal energy transport
-flag_ue           =  1, 10, 10, 10, 10    
-name_ue           =  ue
-# Merid energy transport
-flag_ve           =  1, 10, 10, 10, 10    
-name_ve           =  ve
-# Zonal humidity transport
-flag_uq           =  1, 10, 10, 10, 10    
-name_uq           =  uq
-# Merid humidity transport
-flag_vq           =  1, 10, 10, 10, 10    
-name_vq           =  vq
-# Conv avlbl pot ener
-flag_cape         =  1, 10, 10, 10, 10    
-name_cape         =  cape
-# Cld base pressure
-flag_pbase        =  1, 10, 10, 10, 10    
-name_pbase        =  pbase
-# Cld top pressure
-flag_ptop         =  1, 4, 10, 10, 10    
-name_ptop         =  ptop
-# Cld base mass flux
-flag_fbase        =  1, 10, 10, 10, 10    
-name_fbase        =  fbase
-# Precipitable water
-flag_prw          =  1, 1, 10, 10, 10    
-name_prw          =  prw
-# Boundary Layer Height
-flag_s_pblh       =  1, 10, 10, 1, 1    
-name_s_pblh       =  pblh
-# t at Boundary Layer Height
-flag_s_pblt       =  1, 10, 10, 1, 1    
-name_s_pblt       =  pblt
-# Condensation level
-flag_s_lcl        =  1, 10, 10, 1, 10    
-name_s_lcl        =  lcl
-# Conv avlbl pot enerfor ABL
-flag_s_capCL      =  1, 10, 10, 1, 10    
-name_s_capCL      =  capCL
-# Liq Water in BL
-flag_s_oliqCL     =  1, 10, 10, 1, 10    
-name_s_oliqCL     =  oliqCL
-# Instability criteria(ABL)
-flag_s_cteiCL     =  1, 10, 10, 1, 1    
-name_s_cteiCL     =  cteiCL
-# Exces du thermique
-flag_s_therm      =  1, 10, 10, 1, 1    
-name_s_therm      =  therm
-# deep_cape(HBTM2)
-flag_s_trmb1      =  1, 10, 10, 1, 10    
-name_s_trmb1      =  trmb1
-# inhibition (HBTM2)
-flag_s_trmb2      =  1, 10, 10, 1, 10    
-name_s_trmb2      =  trmb2
-# Point Omega (HBTM2)
-flag_s_trmb3      =  1, 10, 10, 1, 10    
-name_s_trmb3      =  trmb3
-# Bilan au sol sur ocean slab
-flag_slab_bils    =  1, 1, 10, 10, 10    
-name_slab_bils    =  slab_bils 
-# ALE BL
-flag_ale_bl       =  1, 1, 1, 1, 10    
-name_ale_bl       =  ale_bl 
-# alp_bl
-flag_alp_bl       =  1, 1, 1, 1, 10    
-name_alp_bl       =  alp_bl
-# ale_wk
-flag_ale_wk       =  1, 1, 1, 1, 10    
-name_ale_wk       =  ale_wk
-# alp_wk
-flag_alp_wk       =  1, 1, 1, 1, 10    
-name_alp_wk       =  alp_wk
-# ale
-flag_ale          =  1, 1, 1, 1, 10    
-name_ale          =  ale
-# alp
-flag_alp          =  1, 1, 1, 1, 10    
-name_alp          =  alp
-# Convective INhibition
-flag_cin          =  1, 1, 1, 1, 10    
-name_cin          =  cin
-# WAPE
-flag_wape         =  1, 1, 1, 1, 10    
-name_wape         =  wape
-# u, v w t q et phi aux niveaux 200, 500, 700 et 850 hPa
-flag_u850         =  1, 1, 3, 10, 10    
-name_u850         =  u850 
-flag_u700         =  1, 1, 3, 10, 10    
-name_u700         =  u700
-flag_u500         =  1, 1, 3, 10, 10    
-name_u500         =  u500
-flag_u200         =  1, 1, 3, 10, 10    
-name_u200         =  u200
-flag_v850         =  1, 1, 3, 10, 10    
-name_v850         =  v850 
-flag_v700         =  1, 1, 3, 10, 10    
-name_v700         =  v700
-flag_v500         =  1, 1, 3, 10, 10    
-name_v500         =  v500
-flag_v200         =  1, 1, 3, 10, 10    
-name_v200         =  v200
-flag_w850         =  1, 1, 3, 10, 10    
-name_w850         =  w850
-flag_w700         =  1, 1, 3, 10, 10    
-name_w700         =  w700
-flag_w500         =  1, 1, 3, 10, 10    
-name_w500         =  w500
-flag_w200         =  1, 1, 3, 10, 10    
-name_w200         =  w200
-flag_t850         =  1, 1, 3, 10, 10    
-name_t850         =  t850
-flag_t700         =  1, 1, 3, 10, 10    
-name_t700         =  t700
-flag_t500         =  1, 1, 3, 10, 10    
-name_t500         =  t500
-flag_t200         =  1, 1, 3, 10, 10    
-name_t200         =  t200
-flag_q850         =  1, 1, 3, 10, 10    
-name_q850         =  q850
-flag_q700         =  1, 1, 3, 10, 10    
-name_q700         =  q700
-flag_q500         =  1, 1, 3, 10, 10    
-name_q500         =  q500
-flag_q200         =  1, 1, 3, 10, 10    
-name_q200         =  q200
-flag_phi850       =  1, 1, 3, 10, 10    
-name_phi850       =  phi850 
-flag_phi700       =  1, 1, 3, 10, 10    
-name_phi700       =  phi700
-flag_phi500       =  1, 1, 3, 10, 10    
-name_phi500       =  phi500
-flag_phi200       =  1, 1, 3, 10, 10    
-name_phi200       =  phi200
-# Temp mixte oce-sic
-flag_t_oce_sic    =  1, 10, 10, 10, 10    
-name_t_oce_sic    =  t_oce_sic
-# Weak inversion
-flag_weakinv      =  10, 1, 10, 10, 10    
-name_weakinv      =  weakinv
-# dTheta mini
-flag_dthmin       =  10, 1, 10, 10, 10    
-name_dthmin       =  dthmin
-# 10m zonal and meriden wind (ter sic lic oce)
-flag_u10_ter      =  10, 4, 10, 10, 10    
-name_u10_ter      =  u10_ter 
-flag_u10_lic      =  10, 4, 10, 10, 10    
-name_u10_lic      =  u10_lic
-flag_u10_oce      =  10, 4, 10, 10, 10    
-name_u10_oce      =  u10_oce
-flag_u10_sic      =  10, 4, 10, 10, 10    
-name_u10_sic      =  u10_sic 
-flag_v10_ter      =  10, 4, 10, 10, 10    
-name_v10_ter      =  v10_ter
-flag_v10_lic      =  10, 4, 10, 10, 10    
-name_v10_lic      =  v10_lic
-flag_v10_oce      =  10, 4, 10, 10, 10    
-name_v10_oce      =  v10_oce
-flag_v10_sic      =  10, 4, 10, 10, 10    
-name_v10_sic      =  v10_sic
-# Cloud optical thickness
-flag_cldtau       =  10, 5, 10, 10, 10    
-name_cldtau       =  cldtau
-# Cloud optical emissivity
-flag_cldemi       =  10, 5, 10, 10, 10    
-name_cldemi       =  cldemi
-# Relative humidity at 2m
-flag_rh2m         =  10, 5, 10, 10, 10    
-name_rh2m         =  rh2m
-# Saturant humidity at 2m
-flag_qsat2m       =  10, 5, 10, 10, 10    
-name_qsat2m       =  qsat2m
-# Surface air potential temperature
-flag_tpot         =  10, 5, 10, 10, 10    
-name_tpot         =  tpot
-# Surface air equivalent potential temperature
-flag_tpote        =  10, 5, 10, 10, 10    
-name_tpote        =  tpote
-# TKE, tke max and tke (ter sic lic oce)
-flag_tke          =  4, 10, 10, 10, 10    
-name_tke          =  tke
-flag_tke_max      =  4, 10, 10, 10, 10    
-name_tke_max      =  tke_max
-flag_tke_ter      =  10, 4, 10, 10, 10    
-name_tke_ter      =  tke_ter 
-flag_tke_lic      =  10, 4, 10, 10, 10    
-name_tke_lic      =  tke_lic
-flag_tke_oce      =  10, 4, 10, 10, 10    
-name_tke_oce      =  tke_oce
-flag_tke_sic      =  10, 4, 10, 10, 10    
-name_tke_sic      =  tke_sic
-flag_tke_max_ter  =  10, 4, 10, 10, 10    
-name_tke_max_ter  =  tke_max_ter
-flag_tke_max_lic  =  10, 4, 10, 10, 10    
-name_tke_max_lic  =  tke_max_lic
-flag_tke_max_oce  =  10, 4, 10, 10, 10    
-name_tke_max_oce  =  tke_max_oce
-flag_tke_max_sic  =  10, 4, 10, 10, 10    
-name_tke_max_sic  =  tke_max_sic
-# Kz melange
-flag_kz           =  4, 10, 10, 10, 10    
-name_kz           =  kz
-# Kz max melange
-flag_kz_max       =  4, 10, 10, 10, 10    
-name_kz_max       =  kz_max
-# Sfce net SW radiation OR
-flag_SWnetOR      =  10, 10, 2, 10, 10    
-name_SWnetOR      =  SWnetOR
-# Sfce incident SW radiation OR
-flag_SWdownOR     =  10, 10, 2, 10, 10    
-name_SWdownOR     =  SWdownOR
-# Sfce incident LW radiation OR
-flag_LWdownOR     =  10, 10, 2, 10, 10    
-name_LWdownOR     =  LWdownOR
-# Solid Large-scale Precip
-flag_snowl        =  10, 1, 10, 10, 10    
-name_snowl        =  snowl
-# cape max
-flag_cape_max     =  10, 1, 10, 10, 10    
-name_cape_max     =  cape_max
-# Down. IR rad. at surface
-flag_solldown     =  10, 1, 10, 1, 10    
-name_solldown     =  solldown
-# Boundary-layer dTs(o)
-flag_dtsvdfo      =  10, 10, 10, 1, 10    
-name_dtsvdfo      =  dtsvdfo
-# Boundary-layer dTs(t)
-flag_dtsvdft      =  10, 10, 10, 1, 10    
-name_dtsvdft      =  dtsvdft
-# Boundary-layer dTs(g)
-flag_dtsvdfg      =  10, 10, 10, 1, 10    
-name_dtsvdfg      =  dtsvdfg
-# Boundary-layer dTs(g)
-flag_dtsvdfi      =  10, 10, 10, 1, 10    
-name_dtsvdfi      =  dtsvdfi
-# rugosity
-flag_rugs         =  10, 10, 10, 1, 1    
-name_rugs         =  rugs
-# Cloud liquid water content
-flag_lwcon        =  2, 5, 10, 10, 1    
-name_lwcon        =  lwcon
-# Cloud ice water content
-flag_iwcon        =  2, 5, 10, 10, 10    
-name_iwcon        =  iwcon
-# Air temperature
-flag_temp         =  2, 3, 4, 1, 1    
-name_temp         =  temp
-# Potential air temperature
-flag_theta        =  2, 3, 4, 1, 1    
-name_theta        =  theta
-# Specific humidity
-flag_ovap         =  2, 3, 4, 1, 1    
-name_ovap         =  ovap
-# ?
-flag_ovapinit     =  2, 3, 4, 1, 1
-name_ovapinit     =  ovapinit
-# ?
-flag_wvapp        =  2, 10, 10, 10, 10    
-name_wvapp        =  wvapp
-# Geopotential height
-flag_geop         =  2, 3, 10, 1, 1    
-name_geop         =  geop
-# Zonal wind
-flag_vitu         =  2, 3, 4, 1, 1    
-name_vitu         =  vitu
-# Meridional wind
-flag_vitv         =  2, 3, 4, 1, 1    
-name_vitv         =  vitv
-# Vertical wind
-flag_vitw         =  2, 3, 10, 10, 1    
-name_vitw         =  vitw
-# Air pressure
-flag_pres         =  2, 3, 10, 1, 1    
-name_pres         =  pres
-# Cloud Fraction
-flag_rneb         =  2, 5, 10, 10, 1    
-name_rneb         =  rneb
-# Convective Cloud Fraction
-flag_rnebcon      =  2, 5, 10, 10, 1    
-name_rnebcon      =  rnebcon
-# Relative humidity
-flag_rhum         =  2, 10, 10, 10, 10    
-name_rhum         =  rhum
-# Ozone concentration
-flag_ozone        =  2, 10, 10, 10, 10    
-name_ozone        =  ozone
-# saturated updraft
-flag_upwd         =  2, 10, 10, 10, 10    
-name_upwd         =  upwd
-# Physics dT
-flag_dtphy        =  2, 10, 10, 10, 1    
-name_dtphy        =  dtphy
-# Physics dq
-flag_dqphy        =  2, 10, 10, 10, 1    
-name_dqphy        =  dqphy
-# Convective precipitation lic and ice
-flag_pr_con_l     =  2, 10, 10, 10, 10    
-name_pr_con_l     =  pr_con_l
-flag_pr_con_i     =  2, 10, 10, 10, 10    
-name_pr_con_i     =  pr_con_i
-# Large scale precipitation lic and ice
-flag_pr_lsc_l     =  2, 10, 10, 10, 10    
-name_pr_lsc_l     =  pr_lsc_l
-flag_pr_lsc_i     =  2, 10, 10, 10, 10    
-name_pr_lsc_i     =  pr_lsc_i
-# Albedo surf, Snow age and rugosity (ter sic lic oce)
-flag_albe_ter     =  3, 4, 10, 1, 10    
-name_albe_ter     =  albe_ter 
-flag_albe_lic     =  3, 4, 10, 1, 10    
-name_albe_lic     =  albe_lic
-flag_albe_oce     =  3, 4, 10, 1, 10    
-name_albe_oce     =  albe_oce
-flag_albe_sic     =  3, 4, 10, 1, 10    
-name_albe_sic     =  albe_sic 
-flag_ages_ter     =  3, 10, 10, 10, 10    
-name_ages_ter     =  ages_ter
-flag_ages_lic     =  3, 10, 10, 10, 10    
-name_ages_lic     =  ages_lic
-flag_ages_oce     =  3, 10, 10, 10, 10    
-name_ages_oce     =  ages_oce
-flag_ages_sic     =  3, 10, 10, 10, 10    
-name_ages_sic     =  ages_sic
-flag_rugs_ter     =  3, 4, 10, 1, 10    
-name_rugs_ter     =  rugs_ter 
-flag_rugs_lic     =  3, 4, 10, 1, 10    
-name_rugs_lic     =  rugs_lic
-flag_rugs_oce     =  3, 4, 10, 1, 10    
-name_rugs_oce     =  rugs_oce
-flag_rugs_sic     =  3, 4, 10, 1, 10    
-name_rugs_sic     =  rugs_sic
-# Surface albedo
-flag_albs         =  3, 10, 10, 1, 10    
-name_albs         =  albs
-# Surface albedo LW
-flag_albslw       =  3, 10, 10, 1, 10    
-name_albslw       =  albslw
-# Convective Cloud Liquid water content
-flag_clwcon       =  4, 10, 10, 10, 10    
-name_clwcon       =  clwcon
-# undilute adiab updraft
-flag_Ma           =  4, 10, 10, 10, 10    
-name_Ma           =  Ma
-# saturated downdraft
-flag_dnwd         =  4, 10, 10, 10, 10    
-name_dnwd         =  dnwd
-# unsat. downdraft
-flag_dnwd0        =  4, 10, 10, 10, 10    
-name_dnwd0        =  dnwd0
-# Dynamics dT dQ dU dV, .....
-flag_dtdyn        =  4, 10, 10, 10, 1
-name_dtdyn        =  dtdyn
-flag_dqdyn        =  4, 10, 10, 10, 1    
-name_dqdyn        =  dqdyn
-flag_dudyn        =  4, 10, 10, 10, 1    
-name_dudyn        =  dudyn
-flag_dvdyn        =  4, 10, 10, 10, 1    
-name_dvdyn        =  dvdyn 
-flag_dtcon        =  4, 5, 10, 10, 10    
-name_dtcon        =  dtcon
-flag_ducon        =  4, 10, 10, 10, 10    
-name_ducon        =  ducon
-flag_dqcon        =  4, 5, 10, 10, 10    
-name_dqcon        =  dqcon
-flag_dtwak        =  4, 5, 10, 10, 10    
-name_dtwak        =  dtwak
-flag_dqwak        =  4, 5, 10, 10, 10    
-name_dqwak        =  dqwak
-flag_wake_h       =  4, 5, 10, 10, 10    
-name_wake_h       =  wake_h
-flag_wake_s       =  4, 5, 10, 10, 10    
-name_wake_s       =  wake_s
-flag_wake_deltat  =  4, 5, 10, 10, 10    
-name_wake_deltat  =  wake_deltat
-flag_wake_deltaq  =  4, 5, 10, 10, 10    
-name_wake_deltaq  =  wake_deltaq
-flag_wake_omg     =  4, 5, 10, 10, 10    
-name_wake_omg     =  wake_omg 
-flag_Vprecip      =  10, 10, 10, 10, 10    
-name_Vprecip      =  Vprecip
-flag_ftd          =  4, 5, 10, 10, 10    
-name_ftd          =  ftd
-flag_fqd          =  4, 5, 10, 10, 10    
-name_fqd          =  fqd
-flag_dtlsc        =  4, 10, 10, 10, 10    
-name_dtlsc        =  dtlsc
-flag_dtlschr      =  4, 10, 10, 10, 10    
-name_dtlschr      =  dtlschr
-flag_dqlsc        =  4, 10, 10, 10, 10    
-name_dqlsc        =  dqlsc
-flag_dtvdf        =  4, 10, 10, 1, 10    
-name_dtvdf        =  dtvdf
-flag_dqvdf        =  4, 10, 10, 1, 10    
-name_dqvdf        =  dqvdf
-flag_dteva        =  4, 10, 10, 10, 10    
-name_dteva        =  dteva
-flag_dqeva        =  4, 10, 10, 10, 10    
-name_dqeva        =  dqeva
-flag_ptconv       =  4, 10, 10, 10, 10    
-name_ptconv       =  ptconv
-flag_ratqs        =  4, 10, 10, 10, 10    
-name_ratqs        =  ratqs
-flag_dtthe        =  4, 10, 10, 10, 10    
-name_dtthe        =  dtthe
-flag_f_th         =  4, 10, 10, 10, 10    
-name_f_th         =  f_th
-flag_e_th         =  4, 10, 10, 10, 10    
-name_e_th         =  e_th
-flag_w_th         =  4, 10, 10, 10, 10    
-name_w_th         =  w_th
-flag_lambda_th    =  4, 10, 10, 10, 10    
-name_lambda_th    =  lambda_th
-flag_q_th         =  4, 10, 10, 10, 10    
-name_q_th         =  q_th
-flag_a_th         =  4, 10, 10, 10, 10    
-name_a_th         =  a_th
-flag_d_th         =  4, 10, 10, 10, 10    
-name_d_th         =  d_th
-flag_f0_th        =  4, 10, 10, 10, 10    
-name_f0_th        =  f0_th
-flag_zmax_th      =  4, 10, 10, 10, 10    
-name_zmax_th      =  zmax_th
-flag_dqthe        =  4, 10, 10, 10, 1    
-name_dqthe        =  dqthe
-flag_dtajs        =  4, 10, 10, 10, 10    
-name_dtajs        =  dtajs
-flag_dqajs        =  4, 10, 10, 10, 10    
-name_dqajs        =  dqajs
-flag_dtswr        =  4, 10, 10, 10, 1    
-name_dtswr        =  dtswr
-flag_dtsw0        =  4, 10, 10, 10, 10    
-name_dtsw0        =  dtsw0 
-flag_dtlwr        =  4, 10, 10, 10, 1    
-name_dtlwr        =  dtlwr
-flag_dtlw0        =  4, 10, 10, 10, 10    
-name_dtlw0        =  dtlw0
-flag_dtec         =  4, 10, 10, 10, 10    
-name_dtec         =  dtec
-flag_duvdf        =  4, 10, 10, 10, 10    
-name_duvdf        =  duvdf
-flag_dvvdf        =  4, 10, 10, 10, 10    
-name_dvvdf        =  dvvdf
-flag_duoro        =  4, 10, 10, 10, 10    
-name_duoro        =  duoro
-flag_dvoro        =  4, 10, 10, 10, 10    
-name_dvoro        =  dvoro
-flag_dulif        =  4, 10, 10, 10, 10    
-name_dulif        =  dulif
-flag_dvlif        =  4, 10, 10, 10, 10    
-name_dvlif        =  dvlif
-###! Attention a refaire correctement
-flag_trac01       =  4, 10, 10, 10, 10    
-name_trac01       =  trac01 
-flag_trac02       =  4, 10, 10, 10, 10    
-name_trac02       =  trac02
-
Index: LMDZ4/trunk/physiq.def
===================================================================
--- LMDZ4/trunk/physiq.def	(revision 5799)
+++ 	(revision )
@@ -1,219 +1,0 @@
-## $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
-pmagic=0.008
-#
-#
-#
-# Parametres fichiers de sortie
-#
-### type_run = type run par rapport aux fichiers et variables de sortie
-# - type_run = CLIM/ENSP (=1)
-# - type_run = AMIP/CFMI (=2)
-type_run=AMIP
-### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
-OK_journe=n
-### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc
-OK_mensuel=y
-### OK_instan=y, ecrire sorties "instantannees" (chaque pas de temps de la  physique)
-OK_instan=n
-### OK_hf=y, ecrire sorties hautes frequence histhf.nc, =n pas de fichier histhf.nc
-ok_hf=n
-#
-# Parametres niveau de sorties differents fichiers 
-#
-### lev_histhf=0-4, niveau de sortie fichier "histhf.nc" 
-# - lev_histhf=0 => pas de sorties histhf.nc
-# - lev_histhf=2 => defaut
-# - lev_histhf=3 => variables sur niveaux standards
-# - lev_histhf=4 => histhf3d.nc champs 3d niveaux modele => fichier. histhf3d.nc
-lev_histhf=2
-### lev_histday=0-5, niveau de sortie fichier "histday.nc"
-# - lev_histday=0 => pas de sorties lev_histday.nc
-# - lev_histday=2 => defaut
-# - lev_histday=3 => + champs 3D => F. Lott
-# - lev_histday=4 => + champs sous-surfaces
-# - lev_histday=5 => + champs F. Aires
-lev_histday=2
-### lev_histmth=0-4, niveau de sortie fichier "histmth.nc"
-# - lev_histmth=0 => pas de sorties lev_histmth.nc
-# - lev_histmth=2 => defaut
-# - lev_histmth=3 => albedo, rugosite sous-surfaces
-# - lev_histmth=4 => champs tendances 3d
-lev_histmth=2
-### ecrit_hf = frequence ecriture fichier histhf.nc en jours
-ecrit_hf=0.250
-### ecrit_day = frequence ecriture fichier histday.nc en jours
-ecrit_day=1.
-### ecrit_mth = frequence ecriture fichier histmth.nc en jours
-ecrit_mth=30
-### freqin_isccp = frequence input en secondes du simulateur ISCCP
-freq_ISCCP=10800.
-### freqout_isccp = frequence output en jours du simulateur ISCCP
-ecrit_ISCCP=30
-### 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 = 1361.
-#
-# Taux gaz a effet de serre
-#
-### co2_ppm =    taux CO2 en ppm
-co2_ppm = 367.
-### CH4_ppb =    taux CH4 en ppb
-CH4_ppb = 1760.
-### N2O_ppb =    taux N2O en ppb
-N2O_ppb = 316.
-### CFC11_ppt =  taux CFC11 en ppt
-CFC11_ppt = 741.2
-### CFC12_ppt =  taux CFC12 en ppt
-CFC12_ppt = 535.
-#
-# 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=actuel
-###  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.
-#
-# 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 = 1
-### overlap =    Hypothese de Recouvrement (HR) utilisee pour le simulateur ISCCP
-# - overlap=1    Max overlap
-# - overlap=2    Random overlap
-# - overlap=3    Max/Random overlap
-overlap = 3
Index: LMDZ4/trunk/run.def
===================================================================
--- LMDZ4/trunk/run.def	(revision 5799)
+++ 	(revision )
@@ -1,47 +1,0 @@
-# $Id$
-#
-## Fichier de configuration general
-## 
-INCLUDEDEF=physiq.def
-INCLUDEDEF=gcm.def
-INCLUDEDEF=orchidee.def
-INCLUDEDEF=output.def
-## Type de calendrier utilise
-## valeur possible: earth_360d (defaut), earth_365d, earth_366d
-calend=earth_360d
-## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
-dayref=1
-##  Annee de l'etat  initial (   avec  4  chiffres   )
-anneeref=1980
-## Nombre de jours d'integration
-nday=5
-## Remise a zero de la date initiale
-raz_date=0
-## periode de sortie des variables de controle (en pas)
-iconser=240
-## sorties instantanees dans la dynamique (fichiers dyn_hist.nc and co.)
-ok_dyn_ins=n
-## periode d'ecriture des sorties instantanees dans la dynamique
-## (en pas dynamiques)
-iecri=960
-## sorties de valeurs moyennes dans la dynamique (fichiers dyn_hist_ave.nc and co.)
-ok_dyn_ave=n
-## periode de stockage des moyennes dans la dynamique et dans dynzon
-periodav=30
-## flag de sortie dynzon
-ok_dynzon=n
-## activation du calcul d equilibrage de charge
-adjust=n
-## activation du filtre fft
-use_filtre_fft=n
-## niveau d'impression de controle
-prt_level=1
-##
-## Informations sur la configuration utilisee
-## 
-### type_ocean = force / slab  /couple
-type_ocean=force
-### version_ocean = nemo / opa8 
-version_ocean=nemo
-### VEGET= y si ORCHIDEE, =n si bucket
-VEGET=n
Index: LMDZ4/trunk/traceur.def
===================================================================
--- LMDZ4/trunk/traceur.def	(revision 5799)
+++ 	(revision )
@@ -1,5 +1,0 @@
-4
-14 14 H2Ov
-10 10 H2Ol
-10 10 RN
-10 10 PB
