Index: anches/dev_001_GM/NEMO/OPA_SRC/par_EEL_R2.h90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_EEL_R2.h90	(revision 777)
+++ 	(revision )
@@ -1,98 +1,0 @@
-   !!---------------------------------------------------------------------
-   !!                     ***  par_EEL_R2.h90  ***
-   !!   Ocean Domain : 2 km resolution Channel (EEL_R2 configuration)
-   !!---------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-   !!----------------------------------------------------------------------
-   CHARACTER (len=16)      &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::    &  
-      cp_cfg = "eel"            !: name of the configuration
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jp_cfg = 2   ,         &  !: resolution of the configuration (km)
-
-      ! data size              !!! * size of all the input files *
-      jpidta  = 83,          &  !: 1st horizontal dimension ( >= jpi )
-      jpjdta  = 242,         &  !: 2nd    "            "    ( >= jpj )
-      jpkdta  = 30,          &  !: number of levels         ( >= jpk )
-
-      ! global domain size     !!! * full domain *
-      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
-      jpjglo  = jpjdta,      &  !: 2nd    "                  "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! zoom starting position
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data indices
-
-      ! Domain characteristics
-      jperio  =     1 ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =     1 ,      &  !: number of islands
-      jpnisl  = jpiglo          !: maximum number of points per island
-
-   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-   REAL(wp), PARAMETER ::   &  !:
-      pp_not_used       = 999999._wp  , & !: ???
-      pp_to_be_computed =      0._wp      !: ???
-   !!
-   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
-   !!
-   INTEGER,PARAMETER   ::    & !:
-      jphgr_msh = 3            !: type of horizontal mesh
-      !                        ! = 0 curvilinear coordinate on the sphere
-      !                        !     read in coordinate.nc file
-      !                        ! = 1 geographical mesh on the sphere
-      !                        !     with regular grid-spacing
-      !                        ! = 2 f-plane with regular grid-spacing
-      !                        ! = 3 beta-plane with regular grid-spacing
-      !                        ! = 4 Mercator grid with T/U point at the equator  with
-      !                        !     isotropic resolution (e1_deg)
-
-      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
-      !   The mercator grid starts only approximately at gphi0 because
-      !   of the constraint that the equator be a T point.
-   REAL(wp) ,PARAMETER ::     &  !:
-      ppglam0  =    0.0_wp,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  =   35.0_wp,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                          ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg = pp_not_used ,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg = pp_not_used ,   &  !: meridional grid-spacing (degrees)
-      !
-      ppe1_m   = 2000.0_wp,   &  !: zonal      grid-spacing (meters )
-      ppe2_m   = 2000.0_wp       !: meridional grid-spacing (meters )
-   !!
-   !!  Coefficients associated with the vertical coordinate system
-   !!
-   REAL(wp) &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::     &  !:
-      &     ppsur = -2033.194295283385_wp   ,  &  !: Computed in domzgr
-      &     ppa0  =  155.8325369664153_wp   ,  &  !:
-      &     ppa1  =  146.3615918601890_wp   ,  &  !:
-      !
-      &     ppkth =  17.28520372419791_wp   ,  &  !: (non dimensional): gives the approximate
-      !                                           !: layer number above which  stretching will
-      !                                           !: be maximum. Usually of order jpk/2.
-      &     ppacr =  5.000000000000000_wp         !: (non dimensional): stretching factor
-      !                                           !: for the grid. The highest zacr, the smallest
-      !                                           !: the stretching.
-
-   !!
-   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-   !!
-   REAL(wp), PARAMETER ::        &
-      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
-      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
-   !!---------------------------------------------------------------------
Index: anches/dev_001_GM/NEMO/OPA_SRC/par_EEL_R5.h90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_EEL_R5.h90	(revision 777)
+++ 	(revision )
@@ -1,98 +1,0 @@
-   !!---------------------------------------------------------------------
-   !!                     ***  par_EEL_R5.h90  ***   
-   !!   Ocean Domain : 5 km resolution Channel (EEL_R5 configuration)
-   !!---------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-   !!----------------------------------------------------------------------
-   CHARACTER (len=16)      &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::    &  
-      cp_cfg = "eel"            !: name of the configuration
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jp_cfg = 5      ,      &  !: resolution of the configuration (km)
-
-      ! data size              !!! * size of all the input files
-      jpidta  =  66   ,      &  !: first horizontal dimension > or = to jpi
-      jpjdta  =  66   ,      &  !: second                     > or = to jpj
-      jpkdta  =  31   ,      &  !: number of levels           > or = to jpk
-
-      ! total domain size      !!! * full domain *
-      jpiglo  = jpidta,      &  !: first  dimension of global domain --> i
-      jpjglo  = jpjdta,      &  !: second dimension of global domain --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! zoom starting position
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data indices
-
-      ! Domain characteristics
-      jperio  =   1   ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =   1   ,      &  !: number of islands
-      jpnisl  = jpiglo          !: maximum number of points per island
-
-   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-   REAL(wp), PARAMETER ::   &  !:
-      pp_not_used       = 999999._wp ,  &  !:
-      pp_to_be_computed =      0._wp       !:
-
-   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
-
-   INTEGER,PARAMETER   ::    & !:
-      jphgr_msh = 2            !: type of horizontal mesh
-      !                        ! = 0 curvilinear coordinate on the sphere
-      !                        !     read in coordinate.nc file
-      !                        ! = 1 geographical mesh on the sphere
-      !                        !     with regular grid-spacing
-      !                        ! = 2 f-plane with regular grid-spacing
-      !                        ! = 3 beta-plane with regular grid-spacing
-      !                        ! = 4 Mercator grid with T/U point at the equator  with
-      !                        !     isotropic resolution (e1_deg)
-
-      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
-      !   The mercator grid starts only approximately at gphi0 because
-      !   of the constraint that the equator be a T point.
-   REAL(wp) &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::     &  !:
-      ppglam0  =    0.0_wp,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  = 43.436430714_wp,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                          ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees)
-      !
-      ppe1_m   = 8000.0_wp,   &  !: zonal      grid-spacing (meters)
-      ppe2_m   = 8000.0_wp       !: meridional grid-spacing (meters)
-   !!
-   !!  Coefficients associated with the vertical coordinate system
-   !!
-   REAL(wp), PARAMETER  ::       &   !:
-      &     ppsur = -4762.96143546300_wp    ,  &  !: Computed in domzgr, set ppdzmin, pphmax below
-      &     ppa0  =   255.58049070440_wp    ,  &  !:
-      &     ppa1  =   245.58132232490_wp    ,  &  !:
-      !
-      &     ppkth =    21.43336197938_wp    ,  &  !: (non dimensional): gives the approximate
-      !                                           !: layer number above which  stretching will
-      !                                           !: be maximum. Usually of order jpk/2.
-      &     ppacr =     3.00000000000_wp          !: (non dimensional): stretching factor
-      !                                           !: for the grid. The highest zacr, the smallest
-      !                                           !: the stretching.
-
-   !!
-   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-   !!
-   REAL(wp), PARAMETER ::        &  !:
-      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
-      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
-   !!---------------------------------------------------------------------
Index: anches/dev_001_GM/NEMO/OPA_SRC/par_EEL_R6.h90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_EEL_R6.h90	(revision 777)
+++ 	(revision )
@@ -1,98 +1,0 @@
-   !!---------------------------------------------------------------------
-   !!                     ***  par_EEL_R6.h90  ***   
-   !!   Ocean Domain : 6 km resolution Channel (EEL_R6 configuration)
-   !!---------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-   !!----------------------------------------------------------------------
-   CHARACTER (len=16)      &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::    &  
-      cp_cfg = "eel"            !: name of the configuration
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::     &
-      jp_cfg = 6      ,      &  !: resolution of the configuration (km)
-
-      ! data size              !!! * size of all the input files *
-      jpidta  = 29,          &  !: 1st lateral dimension ( >= jpi )
-      jpjdta  = 83,          &  !: 2nd    "         "    ( >= jpj )
-      jpkdta  = 30,          &  !: number of levels      ( >= jpk )
-
-      ! global domain size     !!! * full domain *
-      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
-      jpjglo  = jpjdta,      &  !: 2nd    "                  "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! starting position of the zoom
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data domain indices
-
-      ! Domain characteristics
-      jperio  =      1,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =      1,      &  !: number of islands
-      jpnisl  = jpiglo          !: maximum number of points per island
-
-   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-   REAL(wp), PARAMETER ::   &  !:
-      pp_not_used       = 999999._wp ,  &  !: ???
-      pp_to_be_computed =      0._wp       !: ???
-
-   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
-
-   INTEGER,PARAMETER   ::    & !:
-      jphgr_msh = 3            !: type of horizontal mesh
-      !                        ! = 0 curvilinear coordinate on the sphere
-      !                        !     read in coordinate.nc file
-      !                        ! = 1 geographical mesh on the sphere
-      !                        !     with regular grid-spacing
-      !                        ! = 2 f-plane with regular grid-spacing
-      !                        ! = 3 beta-plane with regular grid-spacing
-      !                        ! = 4 Mercator grid with T/U point at the equator  with
-      !                        !     isotropic resolution (e1_deg)
-
-      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
-      !   The mercator grid starts only approximately at gphi0 because
-      !   of the constraint that the equator be a T point.
-   REAL(wp) &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::     &  !:
-      ppglam0  =    0.0_wp,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  =   35.0_wp,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                          !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees)
-      !
-      ppe1_m   = 6000.0_wp,   &  !: zonal      grid-spacing (meters )
-      ppe2_m   = 6000.0_wp       !: meridional grid-spacing (meters )
-
-   !!  Coefficients associated with the vertical coordinate system
-
-   REAL(wp), PARAMETER  ::       &  !:
-      &     ppsur = -2033.194295283385_wp   ,  &  !: Computed in domzgr
-      &     ppa0  =  155.8325369664153_wp   ,  &  !:
-      &     ppa1  =  146.3615918601890_wp   ,  &  !:
-      !
-      &     ppkth =  17.28520372419791_wp   ,  &  !: (non dimensional): gives the approximate
-      !                                           !: layer number above which  stretching will
-      !                                           !: be maximum. Usually of order jpk/2.
-      &     ppacr =  5.000000000000000_wp         !: (non dimensional): stretching factor
-      !                                           !: for the grid. The highest zacr, the smallest
-      !                                           !: the stretching.
-
-
-   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-   !!
-   REAL(wp), PARAMETER ::        &  !:
-      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
-      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
-   !!---------------------------------------------------------------------
Index: anches/dev_001_GM/NEMO/OPA_SRC/par_GYRE.h90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_GYRE.h90	(revision 777)
+++ 	(revision )
@@ -1,103 +1,0 @@
-   !!---------------------------------------------------------------------
-   !!                     ***  par_GYRE.h90  ***
-   !!   Ocean Domain : GYRE configuration at 1/jp_cfg degree resolution 
-   !!---------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-   !!----------------------------------------------------------------------
-   CHARACTER (len=16)      &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::    &  
-      cp_cfg = "gyre"           !: name of the configuration
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jp_cfg =  1   ,        &  !: 
-
-      ! data size              !!! * size of all the input files *
-      jpidta  = 30*jp_cfg+2, &  !: 1st horizontal dimension ( >= jpi )
-      jpjdta  = 20*jp_cfg+2, &  !: 2nd    "            "    ( >= jpj )
-      jpkdta  = 31,          &  !: number of levels         ( >= jpk )
-
-      ! global domain size     !!! * full domain *
-      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
-      jpjglo  = jpjdta,      &  !: 2nd    "                  "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! zoom starting position
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data indices
-
-      ! Domain characteristics
-      jperio  =     0 ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =     1 ,      &  !: number of islands
-      jpnisl  = jpiglo          !: maximum number of points per island
-
-   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-   REAL(wp), PARAMETER ::   &  !:
-      pp_not_used       = 999999._wp  , & !: ???
-      pp_to_be_computed =      0._wp      !: ???
-   !!
-   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
-   !!
-   INTEGER,PARAMETER   ::    & !:
-      jphgr_msh = 5            !: type of horizontal mesh
-      !                        ! = 0 curvilinear coordinate on the sphere
-      !                        !     read in coordinate.nc file
-      !                        ! = 1 geographical mesh on the sphere
-      !                        !     with regular grid-spacing
-      !                        ! = 2 f-plane with regular grid-spacing
-      !                        ! = 3 beta-plane with regular grid-spacing
-      !                        ! = 4 Mercator grid with T/U point at the equator  with
-      !                        !     isotropic resolution (e1_deg)
-      !                        ! =5  beta-plane with regular grid-spacing and rotated domain (GYRE configuration)
-
-      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
-      !   The mercator grid starts only approximately at gphi0 because
-      !   of the constraint that the equator be a T point.
-   REAL(wp) &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::     &  !:
-      ppglam0  =    0.0_wp,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  =   29.0_wp,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                          ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg = pp_not_used ,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg = pp_not_used ,   &  !: meridional grid-spacing (degrees)
-      !
-      ppe1_m   = pp_not_used,    &  !: zonal      grid-spacing (meters )
-      ppe2_m   = pp_not_used        !: meridional grid-spacing (meters )
-   !!
-   !!  Coefficients associated with the vertical coordinate system
-   !!
-
-   REAL(wp), PARAMETER  ::       &  !:
-      &     ppsur = -2033.194295283385_wp   ,  &  !: 
-      &     ppa0  =  155.8325369664153_wp   ,  &  !:
-      &     ppa1  =  146.3615918601890_wp   ,  &  !:
-      !
-      &     ppkth =  17.28520372419791_wp   ,  &  !: (non dimensional): gives the approximate
-      !                                           !    layer number above which  stretching will
-      !                                           !    be maximum. Usually of order jpk/2.
-      &     ppacr =  5.000000000000000_wp         !: (non dimensional): stretching factor
-      !                                           !    for the grid. The highest zacr, the smallest
-      !                                           !    the stretching.
-
-   !!
-   !!  If all ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-   !!
-   REAL(wp), PARAMETER ::        &  !:
-      &     ppdzmin = pp_not_used   ,  &  !: (meters): depth of the top (first) model layer
-      !                             !            depth of second "w" level
-      &     pphmax  = pp_not_used         !: (meters): maximum depth of the ocean
-      !                             !            depth of the last "w" level
-
-   !!---------------------------------------------------------------------
Index: anches/dev_001_GM/NEMO/OPA_SRC/par_ORCA_R025.h90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_ORCA_R025.h90	(revision 777)
+++ 	(revision )
@@ -1,91 +1,0 @@
-   !!---------------------------------------------------------------------
-   !!                     ***  par_ORCA_R025.h90  ***  
-   !!   Ocean Domain : 0.25 degrees resolution global ocean
-   !!                  (0RCA_R025 configuration)
-   !!---------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-   !!----------------------------------------------------------------------
-   CHARACTER (len=16)      &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::    &  
-      cp_cfg = "orca"           !: name of the configuration
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jp_cfg = 025  ,        &  !: resolution of the configuration (degrees)
-      ! Original data size
-      jpidta  = 1442,        &  !: first horizontal dimension > or = to jpi
-      jpjdta  = 1021,        &  !: second                     > or = to jpj
-      jpkdta  = 46 ,         &  !: number of levels           > or = to jpk
-      ! total domain matrix size
-      jpiglo  = jpidta,      &  !: first  dimension of global domain --> i
-      jpjglo  = jpjdta,      &  !: second dimension of global domain --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! starting position of the zoom
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data indices
-      ! Domain characteristics
-      jperio  =    4  ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =  150  ,      &  !: number of islands
-      jpnisl  = 3000            !: maximum number of points per island
-
-   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-   REAL,PARAMETER      ::  pp_not_used = 999999_wp , &
-      &                    pp_to_be_computed = 0._wp
-   !!
-   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
-   !!
-   INTEGER, PARAMETER ::     & !
-      jphgr_msh = 0            !: type of horizontal mesh
-      !                        !  = 0 curvilinear coordinate on the sphere
-      !                        !      read in coordinate.nc file
-      !                        !  = 1 geographical mesh on the sphere
-      !                        !      with regular grid-spacing
-      !                        !  = 2 f-plane with regular grid-spacing
-      !                        !  = 3 beta-plane with regular grid-spacing
-      !                        !  = 4 Mercator grid with T/U point at the equator  with
-      !                        !      isotropic resolution (e1_deg)
-
-      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
-      !   The mercator grid starts only approximately at gphi0 because
-      !   of the constraint that the equator be a T point.
-   REAL(wp), PARAMETER ::       &  !
-      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees)
-      !
-      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters )
-      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters )
-
-   !!  Coefficients associated with the vertical coordinate system
-
-   REAL(wp), PARAMETER  ::       &
-      &     ppsur = pp_to_be_computed ,  &  !: Computed in domzgr, set ppdzmin and pphmax below
-      &     ppa0  = pp_to_be_computed ,  &  !:    "           "
-      &     ppa1  = pp_to_be_computed ,  &  !:    "           "
-      !
-      &     ppkth =  23.563_wp        ,  &  !: (non dimensional): gives the approximate
-      !                                     !: layer number above which  stretching will
-      !                                     !: be maximum. Usually of order jpk/2.
-      &     ppacr =    9.00000000000_wp     !: (non dimensional): stretching factor
-      !                                     !: for the grid. The highest zacr, the smallest
-      !                                     !: the stretching.
-
-   !!
-   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-   !!
-   REAL(wp), PARAMETER ::        &
-      &     ppdzmin = 6._wp           ,  &  !: (meters) vertical thickness of the top layer
-      &     pphmax  = 5750._wp              !: (meters) Maximum depth of the ocean gdepw(jpk)
-   !!---------------------------------------------------------------------
Index: anches/dev_001_GM/NEMO/OPA_SRC/par_ORCA_R05.h90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_ORCA_R05.h90	(revision 777)
+++ 	(revision )
@@ -1,139 +1,0 @@
-   !!---------------------------------------------------------------------
-   !!                     ***  par_ORCA_R05.h90  ***  
-   !!   Ocean Domain : 0.5 degrees resolution global ocean
-   !!                  (0RCA_R05 configuration)
-   !!---------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-   !!----------------------------------------------------------------------
-   CHARACTER (len=16)      &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::    &  
-      cp_cfg = "orca"           !: name of the configuration
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jp_cfg = 05  ,         &  !: resolution of the configuration (degrees)
-
-      ! data size              !!! * size of all the input files *
-      jpidta  = 722,         &  !: 1st lateral dimension > or = to jpiglo
-      jpjdta  = 511,         &  !: 2nd   "         "     > or = to jpjglo
-      jpkdta  =  31             !: number of levels      > or = to jpkglo
-
-#if defined key_antarctic
-      ! zoom domain size       !!! *  antarctic zoom  * 
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
-      jpjglo  = 187   ,      &  !: 2nd     "                 "    --> j 
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! starting position of the zoom
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data domain indices
-      ! Domain characteristics
-      jperio  =   1   ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =   1   ,      &  !: number of islands
-      jpnisl  =   1             !: maximum number of points per island
-
-#elif defined key_arctic
-      ! zoom domain size       !!! *  arctic zoom  *
-   INTEGER    &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      ! zoom domain size       !!! *  arctic zoom  *
-      jpiglo  = 562,         &  !: 1st dimension of global domain --> i
-      jpjglo  = jpjdta-301+1,&  !: 2nd     "                 "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! zoom starting position
-      jpizoom =  81   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom = 301   ,      &  !: in data domain indices
-      ! Domain characteristics
-      jperio  =   5   ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =   1   ,      &  !: number of islands
-      jpnisl  =   1             !: maximum number of points per island
-
-#else
-      ! global domain size     !!! *  global domain  *
-   INTEGER    &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
-      jpjglo  = jpjdta,      &  !: 2nd     "                 "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! zoom starting position    
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data domain indices
-      ! Domain characteristics
-      jperio  =    6  ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =   79  ,      &  !: number of islands
-      jpnisl  = 2000            !: maximum number of points per island
-#endif
-
-   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-   REAL(wp), PARAMETER ::   &
-      pp_not_used       = 999999._wp , &  !:
-      pp_to_be_computed = 0._wp          !:
-
-   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
-
-   INTEGER, PARAMETER   ::   & !
-      jphgr_msh = 0            !: type of horizontal mesh
-      !                        !  = 0 curvilinear coordinate on the sphere
-      !                        !      read in coordinate.nc file
-      !                        !  = 1 geographical mesh on the sphere
-      !                        !      with regular grid-spacing
-      !                        !  = 2 f-plane with regular grid-spacing
-      !                        !  = 3 beta-plane with regular grid-spacing
-      !                        !  = 4 Mercator grid with T/U point at the equator  with
-      !                        !      isotropic resolution (e1_deg)
-
-      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
-      !   The mercator grid starts only approximately at gphi0 because
-      !   of the constraint that the equator be a T point.
-   REAL(wp) , PARAMETER ::      &  !
-      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                            ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees)
-      !
-      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters )
-      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters )
-
-   !!
-   !! Vertical grid parameter for domzgr
-   !! =====================================
-   !!
-   REAL(wp), PARAMETER  ::       &
-      &     ppsur = -4762.96143546300_wp    ,  &  !: ORCA r4, r2 and r05 coefficients
-      &     ppa0  =   255.58049070440_wp    ,  &  !: (default coefficients)
-      &     ppa1  =   245.58132232490_wp    ,  &  !:
-      &     ppkth =    21.43336197938_wp    ,  &  !: (non dimensional): gives the approximate
-      !                                           !: layer number above which  stretching will
-      !                                           !: be maximum. Usually of order jpk/2.
-      &     ppacr =     3.00000000000_wp          !: (non dimensional): stretching factor
-      !                                           !: for the grid. The highest zacr, the smallest
-      !                                           !: the stretching.
-
-   !!
-   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-   !!
-   REAL(wp), PARAMETER ::        &
-      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
-      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
-   !!---------------------------------------------------------------------
Index: anches/dev_001_GM/NEMO/OPA_SRC/par_ORCA_R2.h90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_ORCA_R2.h90	(revision 777)
+++ 	(revision )
@@ -1,168 +1,0 @@
-   !!---------------------------------------------------------------------
-   !!                     ***  par_ORCA_R2.h90  ***  
-   !!   Ocean Domain : 2 degrees resolution global ocean
-   !!                  (0RCA_R2 configuration)
-   !!---------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-   !!----------------------------------------------------------------------
-   CHARACTER (len=16)      &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      ::    &  
-      cp_cfg = "orca"           !: name of the configuration 
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jp_cfg = 2,            &  !: resolution of the configuration (degrees)
-
-      ! data size              !!! * size of all input files *
-      jpidta  = 182,         &  !: 1st lateral dimension ( >= jpiglo )
-      jpjdta  = 149,         &  !: 2nd    "       "      ( >= jpjglo )
-      jpkdta  = 31              !: number of levels      ( >= jpk    ) 
-
-#if defined key_antarctic
-      ! zoom domain size       !!! *  antarctic zoom  * 
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
-      jpjglo  = 50,          &  !: 2nd    "                  "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! zoom starting position
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data domain indices
-      ! Domain characteristics
-      jperio  =   1   ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =   3   ,      &  !: number of islands
-      jpnisl  = 400             !: maximum number of points per island
-
-#elif defined key_arctic
-      ! zoom domain size       !!! *  arctic zoom  *
-   INTEGER    &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jpiglo  = 142   ,      &  !: 1st dimension of global domain --> i
-      jpjglo  = jpjdta-97+1, &  !: 2nd    "                  "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! zoom starting position 
-      jpizoom =  21   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =  97   ,      &  !: in data domain indices
-      ! Domain characteristics
-      jperio  =   3   ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =   7   ,      &  !: number of islands
-      jpnisl  = 400             !: maximum number of points per island
-
-#elif defined key_cfg_1d
-      ! global domain size     !!! *  global domain  *
-   INTEGER    &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jpiglo  = 3     ,      &  !: 1st dimension of global domain --> i
-      jpjglo  = 3     ,      &  !: 2nd    "                  "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! starting position of the zoom 
-      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
-      ! jpjzoom =   133  ,    &  !: in data domain indices (160W,75N)
-      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
-      ! jpjzoom =   110  ,    &  !: in data domain indices (160W,50N)
-      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
-      ! jpjzoom =   97   ,    &  !: in data domain indices (160W,30N)
-      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
-      ! jpjzoom =   86   ,    &  !: in data domain indices (160W,10N)
-      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
-      ! jpjzoom =   49   ,    &  !: in data domain indices (160W,30S)
-      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
-      ! jpjzoom =   27   ,    &  !: in data domain indices (160W,60S)
-      jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =    7   ,    &  !: in data domain indices (160W,75S)
-      ! Domain characteristics
-      jperio  =   0   ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =  18   ,      &  !: number of islands
-      jpnisl  = 800             !: maximum number of points per island
-#else
-      ! global domain size     !!! *  global domain  *
-   INTEGER    &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i
-      jpjglo  = jpjdta,      &  !: 2nd    "                  "    --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! starting position of the zoom 
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data domain indices
-      ! Domain characteristics
-      jperio  =   4   ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =  18   ,      &  !: number of islands
-      jpnisl  = 800             !: maximum number of points per island
-
-#endif
-
-   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-   REAL(wp), PARAMETER ::   &
-      pp_not_used       = 999999_wp , &  !:
-      pp_to_be_computed = 0._wp          !:
-
-   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
-
-   INTEGER,PARAMETER   ::    & !
-      jphgr_msh = 0            !: type of horizontal mesh
-      !                        !  = 0 curvilinear coordinate on the sphere
-      !                        !      read in coordinate.nc file
-      !                        !  = 1 geographical mesh on the sphere
-      !                        !      with regular grid-spacing
-      !                        !  = 2 f-plane with regular grid-spacing
-      !                        !  = 3 beta-plane with regular grid-spacing
-      !                        !  = 4 Mercator grid with T/U point at the equator  with
-      !                        !      isotropic resolution (e1_deg)
-
-      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
-      !   The mercator grid starts only approximately at gphi0 because
-      !   of the constraint that the equator be a T point.
-   REAL(wp) ,PARAMETER ::       &  !
-      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees)
-      !
-      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters )
-      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters )
-
-   !!
-   !! Vertical grid parameter for domzgr
-   !! ==================================
-   !!
-   REAL(wp), PARAMETER  ::       &
-      &     ppsur = -4762.96143546300_wp    ,  &  !: ORCA r4, r2 and r05 coefficients
-      &     ppa0  =   255.58049070440_wp    ,  &  !: (default coefficients)
-      &     ppa1  =   245.58132232490_wp    ,  &  !:
-      &     ppkth =    21.43336197938_wp    ,  &  !: (non dimensional): gives the approximate
-      !                                           !: layer number above which  stretching will
-      !                                           !: be maximum. Usually of order jpk/2.
-      &     ppacr =     3.00000000000_wp          !: (non dimensional): stretching factor
-      !                                           !: for the grid. The highest zacr, the smallest
-      !                                           !: the stretching.
-
-   !!
-   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-   !!
-   REAL(wp), PARAMETER ::        &
-      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
-      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
-   !!---------------------------------------------------------------------
Index: anches/dev_001_GM/NEMO/OPA_SRC/par_ORCA_R4.h90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_ORCA_R4.h90	(revision 777)
+++ 	(revision )
@@ -1,91 +1,0 @@
-   !!---------------------------------------------------------------------
-   !!                     ***  par_ORCA_R4.h90  ***   
-   !!   Ocean Domain : 4 degrees resolution global ocean
-   !!                  (0RCA_R4 configuration)
-   !!---------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
-   !!----------------------------------------------------------------------
-   CHARACTER (len=16)      &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-       ::    &  
-     cp_cfg = "orca"           !: name of the configuration
-   INTEGER     &
-#if !defined key_agrif
-      , PARAMETER  &
-#endif
-      :: &
-      jp_cfg = 4      ,      &  !: resolution of the configuration (degrees)
-      ! Original data size
-      jpidta  =  92   ,      &  !: first horizontal dimension > or = to jpi
-      jpjdta  =  76   ,      &  !: second                     > or = to jpj
-      jpkdta  =  31   ,      &  !: number of levels           > or = to jpk
-      ! global domain matrix size
-      jpiglo  = jpidta,      &  !: first  dimension of global domain --> i
-      jpjglo  = jpjdta,      &  !: second dimension of global domain --> j
-      jpk     = jpkdta,      &  !: number of vertical levels
-      ! starting position of the zoom
-      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,      &  !: in data indices
-      ! Domain characteristics
-      jperio  =   4   ,      &  !: lateral cond. type (between 0 and 6)
-      jpisl   =  11   ,      &  !: number of islands
-      jpnisl  = 800             !: maximum number of points per island
-
-   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-   REAL(wp), PARAMETER ::   &
-      pp_not_used       = 999999_wp , &  !:
-      pp_to_be_computed = 0._wp          !:
-
-   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
-   INTEGER, PARAMETER ::     & !
-      jphgr_msh = 0            !: type of horizontal mesh
-      !                        ! = 0 curvilinear coordinate on the sphere
-      !                        !     read in coordinate.nc file
-      !                        ! = 1 geographical mesh on the sphere
-      !                        !     with regular grid-spacing
-      !                        ! = 2 f-plane with regular grid-spacing
-      !                        ! = 3 beta-plane with regular grid-spacing
-      !                        ! = 4 Mercator grid with T/U point at the equator  with
-      !                        !     isotropic resolution (e1_deg)
-
-      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
-      !   The mercator grid starts only approximately at gphi0 because
-      !   of the constraint that the equator be a T point.
-   REAL(wp) , PARAMETER ::      &  !
-      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                            ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees)
-      !
-      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters )
-      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters )
-
-
-   !! Vertical grid parameter for domzgr
-   !! ==================================
-
-   REAL(wp), PARAMETER  ::       &
-      &     ppsur = -4762.96143546300_wp    ,  &  !: ORCA r4, r2 and r05 coefficients
-      &     ppa0  =   255.58049070440_wp    ,  &  !: (default coefficients)
-      &     ppa1  =   245.58132232490_wp    ,  &  !:
-      &     ppkth =    21.43336197938_wp    ,  &  !: (non dimensional): gives the approximate
-      !                                           !: layer number above which  stretching will
-      !                                           !: be maximum. Usually of order jpk/2.
-      &     ppacr =     3.00000000000_wp          !: (non dimensional): stretching factor
-      !                                           !: for the grid. The highest zacr, the smallest
-      !                                           !: the stretching.
-
-   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-
-   REAL(wp), PARAMETER ::        &
-      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
-      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
-   !!---------------------------------------------------------------------
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r2.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r2.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r2.F90	(revision 778)
@@ -0,0 +1,131 @@
+MODULE par_eel_r2
+   !!======================================================================
+   !!                        ***  par_eel_r2  ***
+   !! Ocean Domain parameter : 2 km resolution cyclic Channel (EEL_R2 configuration)
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_EEL_R2.h90
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Header: $ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_eel_r2
+   !!---------------------------------------------------------------------
+   !!   'key_eel_r2'   :                            EEL configuration : R2
+   !!---------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                       AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                  no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution      !!! *  GYRE R2  *
+   CHARACTER (len=16) ::   cp_cfg = "eel"    !: name of the configuration 
+   INTEGER            ::   jp_cfg = 2        !: resolution of the configuration [km]
+
+   ! input data size                        !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  =  83                !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  = 242                !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  =  30                !: number of levels      ( >= jpk    ) 
+   
+
+#  if defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :            1D configuration of GYRE_R2
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                   !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                   !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   16                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   11                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   0                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   0                 !: maximum number of points per island
+
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :               Full global domain of GYRE_R2
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta              !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta              !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1                 !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                 !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   1                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   1                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = jpiglo              !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 3   !: type of horizontal mesh
+   !                                       ! = 3 beta-plane with regular grid-spacing
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator be a T point.
+   REAL(wp) ,PARAMETER ::   ppglam0  =    0.0_wp     !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  =   35.0_wp     !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
+   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = 2000.0_wp     !: zonal      grid-spacing (meters )
+   REAL(wp) ,PARAMETER ::   ppe2_m   = 2000.0_wp     !: meridional grid-spacing (meters )
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = -2033.194295283385_wp  !: ORCA r4, r2 and r05 coefficients
+   REAL(wp), PARAMETER ::   ppa0  =  155.8325369664153_wp  !: (default coefficients)
+   REAL(wp), PARAMETER ::   ppa1  =  146.3615918601890_wp  !:
+   REAL(wp), PARAMETER ::   ppkth =  17.28520372419791_wp  !: (non dimensional): gives the approximate
+   !                                                       !: layer number above which  stretching will
+   !                                                       !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     5.00000000000_wp   !: (non dimensional): stretching factor
+   !                                                       !: for the grid. The highest zacr, the smallest
+   !                                                       !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: (meters) vertical thickness of the top layer
+   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: (meters) Maximum depth of the ocean gdepw(jpk)
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_eel_r2
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r5.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r5.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r5.F90	(revision 778)
@@ -0,0 +1,131 @@
+MODULE par_eel_r5
+   !!======================================================================
+   !!                        ***  par_eel_r5  ***
+   !! Ocean Domain parameter : 5 km resolution cyclic Channel (EEL_R5 configuration)
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_EEL_R5.h90
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Header: $ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_eel_r5
+   !!---------------------------------------------------------------------
+   !!   'key_eel_r5'   :                            EEL configuration : R5
+   !!---------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                       AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                  no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution      !!! *  GYRE R5  *
+   CHARACTER (len=16) ::   cp_cfg = "eel"    !: name of the configuration 
+   INTEGER            ::   jp_cfg = 5        !: resolution of the configuration [km]
+
+   ! input data size                        !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  =  66                !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  =  66                !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  =  31                !: number of levels      ( >= jpk    ) 
+   
+
+#  if defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :            1D configuration of GYRE_R5
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                   !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                   !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   16                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   11                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   0                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   0                 !: maximum number of points per island
+
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :               Full global domain of GYRE_R5
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta              !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta              !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1                 !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                 !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   1                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   1                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = jpiglo              !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 3   !: type of horizontal mesh
+   !                                       ! = 3 beta-plane with regular grid-spacing
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator be a T point.
+   REAL(wp) ,PARAMETER ::   ppglam0  =    0.0_wp       !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  = 43.436430714_wp !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                   !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used     !: zonal      grid-spacing [degrees]
+   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used     !: meridional grid-spacing [degrees]
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = 8000.0_wp       !: zonal      grid-spacing [meters]
+   REAL(wp) ,PARAMETER ::   ppe2_m   = 8000.0_wp       !: meridional grid-spacing [meters]
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = -4762.96143546300_wp  !: ORCA r4, r2 and r05 coefficients
+   REAL(wp), PARAMETER ::   ppa0  =   255.58049070440_wp  !: (default coefficients)
+   REAL(wp), PARAMETER ::   ppa1  =   245.58132232490_wp  !:
+   REAL(wp), PARAMETER ::   ppkth =    21.43336197938_wp  !: (non dimensional): gives the approximate
+   !                                                      !: layer number above which  stretching will
+   !                                                      !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     3.00000000000_wp  !: (non dimensional): stretching factor
+   !                                                      !: for the grid. The highest zacr, the smallest
+   !                                                      !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: (meters) vertical thickness of the top layer
+   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: (meters) Maximum depth of the ocean gdepw(jpk)
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_eel_r5
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r6.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r6.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r6.F90	(revision 778)
@@ -0,0 +1,131 @@
+MODULE par_eel_r6
+   !!======================================================================
+   !!                        ***  par_eel_r6  ***
+   !! Ocean Domain parameter : 6 km resolution cyclic Channel (EEL_R6 configuration)
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_EEL_R6.h90
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Header: $ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_eel_r6
+   !!---------------------------------------------------------------------
+   !!   'key_eel_r6'   :                            EEL configuration : R6
+   !!---------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                       AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                  no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution      !!! *  GYRE R6  *
+   CHARACTER (len=16) ::   cp_cfg = "eel"    !: name of the configuration 
+   INTEGER            ::   jp_cfg = 6        !: resolution of the configuration [km]
+
+   ! input data size                        !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  =  29                !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  =  83                !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  =  30                !: number of levels      ( >= jpk    ) 
+   
+
+#  if defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :            1D configuration of GYRE_R6
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                   !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                   !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   16                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   11                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   0                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   0                 !: maximum number of points per island
+
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :               Full global domain of GYRE_R6
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta              !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta              !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1                 !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                 !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   1                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   1                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = jpiglo              !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 3   !: type of horizontal mesh
+   !                                       ! = 3 beta-plane with regular grid-spacing
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator be a T point.
+   REAL(wp) ,PARAMETER ::   ppglam0  =    0.0_wp     !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  =   35.0_wp     !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
+   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = 6000.0_wp     !: zonal      grid-spacing (meters )
+   REAL(wp) ,PARAMETER ::   ppe2_m   = 6000.0_wp     !: meridional grid-spacing (meters )
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = -2033.194295283385_wp  !: ORCA r4, r2 and r05 coefficients
+   REAL(wp), PARAMETER ::   ppa0  =  155.8325369664153_wp  !: (default coefficients)
+   REAL(wp), PARAMETER ::   ppa1  =  146.3615918601890_wp  !:
+   REAL(wp), PARAMETER ::   ppkth =  17.28520372419791_wp  !: (non dimensional): gives the approximate
+   !                                                       !: layer number above which  stretching will
+   !                                                       !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     5.00000000000_wp   !: (non dimensional): stretching factor
+   !                                                       !: for the grid. The highest zacr, the smallest
+   !                                                       !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: (meters) vertical thickness of the top layer
+   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: (meters) Maximum depth of the ocean gdepw(jpk)
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_eel_r6
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_gyre.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_gyre.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_gyre.F90	(revision 778)
@@ -0,0 +1,133 @@
+MODULE par_gyre
+   !!======================================================================
+   !!                        ***  par_gyre  ***
+   !! Ocean Domain parameter : rotated square bassin on a beta-plane at 1/jp_cfg 
+   !!    degree resolution (GYRE configuration) and a 1D configuration case
+   !!    NB: with AGRIF, values are not set as parameters
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_GYRE.h90
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Header: $ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_gyre
+   !!----------------------------------------------------------------------
+   !!   'key_gyre'   :                                   GYRE configuration
+   !!----------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                       AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                  no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution                 !!! *  GYRE  *
+   CHARACTER (len=16), PARAMETER ::   cp_cfg = "gyre"   !: name of the configuration 
+   INTEGER, PARAMETER            ::   jp_cfg = 1        !: resolution of the configuration (degrees)
+   !                                                    !  control the size of GYRE configurations
+
+   ! input data size                                   !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  = 30 * jp_cfg + 2    !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  = 20 * jp_cfg + 2    !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  = 31                 !: number of levels      ( >= jpk    ) 
+
+#  if defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :             1D configuration of MY_CFG
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                   !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                   !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   16                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   11                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   0                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   0                 !: maximum number of points per island
+
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :                  Full global domain of GYRE
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta             !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta             !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   1                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   1                !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 5   !: type of horizontal mesh
+   !                                       !  = 5  beta-plane with regular grid-spacing and rotated domain
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator be a T point.
+   REAL(wp) ,PARAMETER ::   ppglam0  =    0.0_wp     !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  =   29.0_wp     !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
+   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = pp_not_used   !: zonal      grid-spacing (meters )
+   REAL(wp) ,PARAMETER ::   ppe2_m   = pp_not_used   !: meridional grid-spacing (meters )
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = -2033.194295283385_wp     !: GYRE coefficients
+   REAL(wp), PARAMETER ::   ppa0  =   155.8325369664153_wp    !:
+   REAL(wp), PARAMETER ::   ppa1  =   146.3615918601890_wp    !:
+   REAL(wp), PARAMETER ::   ppkth =    17.28520372419791_wp   !: (non dimensional): gives the approximate
+   !                                                          !: layer number above which  stretching will
+   !                                                          !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     5.00000000000_wp      !: (non dimensional): stretching factor
+   !                                                          !: for the grid. The highest zacr, the smallest
+   !                                                          !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: vertical thickness of the top layer   [meters]
+   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: maximum depth of the ocean gdepw(jpk) [meters]
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_gyre
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_kind.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_kind.F90	(revision 777)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_kind.F90	(revision 778)
@@ -4,10 +4,10 @@
    !! Ocean :  define the kind of real for the whole model
    !!======================================================================
-   !! History :
-   !!   8.5   02/06  (G. Madec)  Original code
+   !! History :   8.5   2002/06  (G. Madec)  Original code
+   !!             9.0   2007/11  (G. Madec)  Coding style
    !!----------------------------------------------------------------------
-   !!   OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
+   !! NEMO 2.4, LOCEAN-IPSL (2007) 
+   !! $Id:$ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
@@ -15,8 +15,7 @@
    PRIVATE
 
-   INTEGER, PUBLIC, PARAMETER ::    &  !:
-      jpbyt   = 8       ,           &  !: real size for mpp communications
-      jpbytda = 4       ,           &  !: real size in input data files 4 or 8
-      jpbi3e  = 4                      !: real size for T3E
+   INTEGER, PUBLIC, PARAMETER ::   jpbyt   = 8        !: real size for mpp communications
+   INTEGER, PUBLIC, PARAMETER ::   jpbytda = 4        !: real size in input data files 4 or 8
+   INTEGER, PUBLIC, PARAMETER ::   jpbi3e  = 4        !: real size for T3E
 
    ! Number model from which the SELECTED_*_KIND are requested:
@@ -27,12 +26,12 @@
    !            exponent = 37     exponent = 307
 
-   INTEGER, PUBLIC, PARAMETER ::        &  !: Floating point section
-      sp = SELECTED_REAL_KIND( 6, 37),  &  !: single precision (real 4)
-      dp = SELECTED_REAL_KIND(12,307),  &  !: double precision (real 8)
-      wp = dp                              !: working precision
+   !                                       !: Floating point section
+   INTEGER, PUBLIC, PARAMETER ::   sp = SELECTED_REAL_KIND( 6, 37)    !: single precision (real 4)
+   INTEGER, PUBLIC, PARAMETER ::   dp = SELECTED_REAL_KIND(12,307)    !: double precision (real 8)
+   INTEGER, PUBLIC, PARAMETER ::   wp = dp                            !: working precision
 
-   INTEGER, PUBLIC, PARAMETER ::        &  !: Integer section
-      i4 = SELECTED_INT_KIND(9) ,       &  !: single precision (integer 4)
-      i8 = SELECTED_INT_KIND(14)           !: double precision (integer 8)
+   !                                       !: Integer section
+   INTEGER, PUBLIC, PARAMETER ::   i4 = SELECTED_INT_KIND( 9)         !: single precision (integer 4)
+   INTEGER, PUBLIC, PARAMETER ::   i8 = SELECTED_INT_KIND(14)         !: double precision (integer 8)
 
 !!----------------------------------------------------------------------
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_my_cfg.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_my_cfg.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_my_cfg.F90	(revision 778)
@@ -0,0 +1,137 @@
+MODULE par_my_cfg
+   !!======================================================================
+   !!                        ***  par_my_cfg  ***
+   !! Ocean Domain parameter : user defined configuration (MY_CFG configuration)
+   !!    including  a 1D configuration.
+   !!    NB: with AGRIF, values are not set as parameters
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code
+   !!----------------------------------------------------------------------
+   !!  OPA 9.0 , LOCEAN-IPSL (2007) 
+   !! $Id: $ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_my_cfg
+   !!---------------------------------------------------------------------
+   !!   'key_my_cfg'   :                       user configuration : MY_CFG
+   !!---------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                        AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                   no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution                  !!! *  MY_CFG  *
+   CHARACTER (len=16), PARAMETER ::   cp_cfg = "my_cfg"  !: name of the configuration 
+   INTEGER, PARAMETER            ::   jp_cfg = 1         !: resolution of the configuration (degrees)
+
+   ! input data size                          !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  =  10                 !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  =  12                 !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  =  31                 !: number of levels      ( >= jpk    ) 
+
+#  if defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :             1D configuration of MY_CFG
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                   !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                   !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =    5                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =    6                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   0                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   0                 !: maximum number of points per island
+   
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :                Full global domain of MY_CFG
+   !------------!!--------------------------------------------------------
+   ! global domain size                                 !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta              !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta              !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta              !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1                 !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                 !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                 !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   0                 !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   0                 !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 0   !: type of horizontal mesh
+   !                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc file
+   !                                       !  = 1 geographical mesh on the sphere with regular grid-spacing
+   !                                       !  = 2 f-plane with regular grid-spacing
+   !                                       !  = 3 beta-plane with regular grid-spacing
+   !                                       !  = 4 Mercator grid with T/U point at the equator
+   !                                       !  = 5  beta-plane with regular grid-spacing and rotated domain
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator be a T point.
+   REAL(wp) ,PARAMETER ::   ppglam0  =    0.0_wp     !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  =  -35.0_wp     !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg =    1.0_wp     !: zonal      grid-spacing (degrees)
+   REAL(wp) ,PARAMETER ::   ppe2_deg =    0.5_wp     !: meridional grid-spacing (degrees)
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = 5000.0_wp     !: zonal      grid-spacing (meters )
+   REAL(wp) ,PARAMETER ::   ppe2_m   = 5000.0_wp     !: meridional grid-spacing (meters )
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = -4762.96143546300_wp   !: ORCA r4, r2 and r05 coefficients
+   REAL(wp), PARAMETER ::   ppa0  =   255.58049070440_wp   !: (default coefficients)
+   REAL(wp), PARAMETER ::   ppa1  =   245.58132232490_wp   !:
+   REAL(wp), PARAMETER ::   ppkth =    21.43336197938_wp   !: (non dimensional): gives the approximate
+   !                                                       !: layer number above which  stretching will
+   !                                                       !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     3.00000000000_wp   !: (non dimensional): stretching factor
+   !                                                       !: for the grid. The highest zacr, the smallest
+   !                                                       !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin = 10._wp        !: (meters) vertical thickness of the top layer
+   REAL(wp), PARAMETER ::   pphmax  = 5000._wp      !: (meters) Maximum depth of the ocean gdepw(jpk)
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_my_cfg
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_oce.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_oce.F90	(revision 777)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_oce.F90	(revision 778)
@@ -4,198 +4,74 @@
    !! Ocean :   set the ocean parameters
    !!======================================================================
-   !! History :
-   !!   4.0  !  91     (Imbard, Levy, Madec)  Original code
-   !!   9.0  !  04-01  (G. Madec, J.-M. Molines)  Free form and module
-   !!    "   !  05-11  (V. Garnier) Surface pressure gradient organization
+   !! History :    -   !  1991     (Imbard, Levy, Madec)  Original code
+   !!             1.0  !  2004-01  (G. Madec, J.-M. Molines)  Free form and module
+   !!             2.4  !  2007-11  (G. Madec)  Suppress .h90 + style with AGRIF
    !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Id:$ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
-   !! * Modules used
-   USE par_kind          ! kind parameters
+   USE par_kind        ! kind parameters
+   USE par_orca_r4     ! ORCA R4   configuration
+   USE par_orca_r2     ! ORCA R2   configuration (including 1D or arctic/antarctic zoom)
+   USE par_orca_r05    ! ORCA R05  configuration (including 1D or arctic/antarctic zoom)
+   USE par_orca_r025   ! ORCA R025 configuration (including 1D or arctic/antarctic zoom)
+   USE par_gyre        ! GYRE      configuration
+   USE par_eel_r2      ! GYRE R2   configuration
+   USE par_eel_r5      ! GYRE R5   configuration
+   USE par_eel_r6      ! GYRE R6   configuration
+   USE par_my_cfg      ! MY_CFG    user defined configuration (including 1D)
 
    IMPLICIT NONE
    PUBLIC
 
-   !!----------------------------------------------------------------------
-   !!   Domain decomposition
-   !!----------------------------------------------------------------------
-   !! * if we dont use massively parallel computer (parameters jpni=jpnj=1)
-   !!      so jpiglo=jpi and jpjglo=jpj
-
-#if ! defined key_mpp_dyndist 
-   INTEGER, PUBLIC, PARAMETER ::    &  !: 
-      jpni   = 1,                   &  !: number of processors following i 
-      jpnj   = 1,                   &  !: number of processors following j
-      jpnij  = 1                       !: nb of local domain = nb of processors 
-      !                                !  ( <= jpni x jpnj )
-#else
-   INTEGER, PUBLIC ::               &  !
-      jpni      ,                   &  !: number of processors following i 
-      jpnj      ,                   &  !: number of processors following j
-      jpnij                            !: nb of local domain = nb of processors 
-      !                                !  ( <= jpni x jpnj )
-#endif
-
-   INTEGER, PUBLIC, PARAMETER ::    &  !:
-      jpr2di = 0,                   &  !: number of columns for extra outer halo 
-      jpr2dj = 0,                   &  !: number of rows    for extra outer halo 
-      jpreci = 1,                   &  !: number of columns for overlap 
-      jprecj = 1                       !: number of rows    for overlap 
-
-   !! Ocean Domain sizes
-   !! ------------------
-   !!   data           domain   (jpidta,jpjdta)
-   !!   global or zoom domain   (jpiglo,jpjglo)
-   !!   local          domain   ( jpi  , jpj  )
-   
-#if   defined key_orca_r4
    !!---------------------------------------------------------------------
-   !!   'key_orca_r4'   :                           global ocean : ORCA R4
+   !! Extra halos for mpp computing 
    !!---------------------------------------------------------------------
-#             include "par_ORCA_R4.h90"
-#elif defined key_orca_r2
-   !!---------------------------------------------------------------------
-   !!   'key_orca_r2'   :                           global ocean : ORCA R4
-   !!---------------------------------------------------------------------
-#             include "par_ORCA_R2.h90"
-#elif defined key_orca_r05
-   !!---------------------------------------------------------------------
-   !!   'key_orca_r05'  :                          global ocean : ORCA R05
-   !!---------------------------------------------------------------------
-#             include "par_ORCA_R05.h90"
-#elif defined key_orca_r025
-   !!---------------------------------------------------------------------
-   !!   'key_orca_r025' :                         global ocean : ORCA R025
-   !!---------------------------------------------------------------------
-#             include "par_ORCA_R025.h90"
-#elif defined key_eel_r2
-   !!---------------------------------------------------------------------
-   !!   'key_eel_r2'    :                                 channel : EEL R2
-   !!---------------------------------------------------------------------
-#             include "par_EEL_R2.h90"
-#elif defined key_eel_r5
-   !!---------------------------------------------------------------------
-   !!   'key_eel_r5'    :                                 channel : EEL R5
-   !!---------------------------------------------------------------------
-#             include "par_EEL_R5.h90"
-#elif defined key_eel_r6
-   !!---------------------------------------------------------------------
-   !!   'key_eel_r6'    :                                 channel : EEL R6
-   !!---------------------------------------------------------------------
-#             include "par_EEL_R6.h90"
-#elif defined key_gyre
-   !!---------------------------------------------------------------------
-   !!   'key_gyre'      :                        mid-latitude basin : GYRE
-   !!---------------------------------------------------------------------
-#             include "par_GYRE.h90"
-#else
-   !!---------------------------------------------------------------------
-   !!   default option  :                               small closed basin
-   !!---------------------------------------------------------------------
-   CHARACTER(len=16), PUBLIC, PARAMETER ::   &  !:
-      cp_cfg = "default"               !: name of the configuration
-   INTEGER, PARAMETER ::            &  !:
-      jp_cfg = 0  ,                 &  !: resolution of the configuration
-
-      ! data size                     !!! * size of all input files *
-      jpidta  = 10,                 &  !: 1st lateral dimension ( >= jpi )
-      jpjdta  = 12,                 &  !: 2nd    "         "    ( >= jpj )
-      jpkdta  = 31,                 &  !: number of levels      ( >= jpk )
-
-      ! global or zoom domain size    !!! * computational domain *
-      jpiglo  = jpidta,             &  !: 1st dimension of global domain --> i
-      jpjglo  = jpjdta,             &  !: 2nd    "                  "    --> j
-      jpk     = jpkdta,             &  !: number of vertical levels
-      ! zoom starting position 
-      jpizoom =   1   ,             &  !: left bottom (i,j) indices of the zoom
-      jpjzoom =   1   ,             &  !: in data domain indices
-
-      ! Domain characteristics
-      jperio  =  0,                 &  !: lateral cond. type (between 0 and 6)
-         !                             !  = 0 closed
-         !                             !  = 1 cyclic East-West
-         !                             !  = 2 equatorial symmetric
-         !                             !  = 3 North fold T-point pivot
-         !                             !  = 4 cyclic East-West AND North fold T-point pivot
-         !                             !  = 5 North fold F-point pivot
-         !                             !  = 6 cyclic East-West AND North fold F-point pivot
-      jpisl   =  0,                 &  !: number of islands (rigid-lid only)
-      jpnisl  =  0                     !: maximum number of points per island
-
-      !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
-      !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
-      REAL(wp), PARAMETER ::   &  !:
-         pp_not_used       = 999999._wp , &  !:
-         pp_to_be_computed = 999999._wp      !:
-
-
-   !! Horizontal grid parameters for domhgr
-   !! =====================================
-
-   INTEGER, PUBLIC, PARAMETER   ::   &  !:
-      jphgr_msh = 0            !: type of horizontal mesh
-      !                        !  = 0 curvilinear coordinate on the sphere
-      !                        !      read in coordinate.nc file
-      !                        !  = 1 geographical mesh on the sphere
-      !                        !      with regular grid-spacing
-      !                        !  = 2 f-plane with regular grid-spacing
-      !                        !  = 3 beta-plane with regular grid-spacing
-      !                        !  = 4 Mercator grid with T/U point at the equator  with
-      !                        !      isotropic resolution (e1_deg)
-
-   REAL(wp) , PUBLIC, PARAMETER ::   &   !:
-      ppglam0  =    0.0_wp,   &  !: longitude of first raw and column T-point (jphgr_msh = 1)
-      ppgphi0  =  -35.0_wp,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1)
-      !                          !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)
-      ppe1_deg =    1.0_wp,   &  !: zonal      grid-spacing (degrees)
-      ppe2_deg =    0.5_wp,   &  !: meridional grid-spacing (degrees)
-      ppe1_m   = 5000.0_wp,   &  !: zonal      grid-spacing (degrees)
-      ppe2_m   = 5000.0_wp       !: meridional grid-spacing (degrees)
-
-   !! Vertical grid parameter for domzgr
-   !! ==================================
-
-   REAL(wp), PUBLIC, PARAMETER  ::   &  !:
-      &     ppsur = -4762.96143546300_wp ,  &  !: ORCA r4, r2 and r05 coefficients
-      &     ppa0  =   255.58049070440_wp ,  &  !: (default coefficients)
-      &     ppa1  =   245.58132232490_wp ,  &  !:
-      &     ppkth =    21.43336197938_wp ,  &  !:
-      &     ppacr =     3.00000000000_wp       !:
-
-   !!  If both ppa0 ppa1 and ppsur are specified to 0, then
-   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
-
-   REAL(wp), PUBLIC, PARAMETER ::   &  !:
-      &     ppdzmin = 10._wp             ,  &  !: Minimum vertical spacing
-      &     pphmax  = 5000._wp                 !: Maximum depth
-
-   !!---------------------------------------------------------------------
-#endif
-
-   !!---------------------------------------------------------------------
-   !! Domain Matrix size
-   !!---------------------------------------------------------------------
-   INTEGER  &  !:
-#if !defined key_agrif
-      ,PARAMETER  &
-#endif
-    :: &
-      jpi = ( jpiglo-2*jpreci + (jpni-1) ) / jpni + 2*jpreci ,   &  !: first  dimension
-      jpj = ( jpjglo-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj ,   &  !: second dimension
-      jpim1 = jpi-1,                                             &  !: inner domain indices
-      jpjm1 = jpj-1,                                             &  !:   "            "
-      jpkm1 = jpk-1,                                             &  !:   "            "
-      jpij  = jpi*jpj                                               !:  jpi x jpj
+   INTEGER, PUBLIC, PARAMETER ::   jpr2di = 0   !: number of columns for extra outer halo (solsor only)
+   INTEGER, PUBLIC, PARAMETER ::   jpr2dj = 0   !: number of rows    for extra outer halo (solsor only)
+   INTEGER, PUBLIC, PARAMETER ::   jpreci = 1   !: number of columns for overlap 
+   INTEGER, PUBLIC, PARAMETER ::   jprecj = 1   !: number of rows    for overlap 
 
 #if defined key_agrif
    !!---------------------------------------------------------------------
-   !! Agrif variables
+   !!   'key_agrif'   :                                              AGRIF
    !!---------------------------------------------------------------------
-   INTEGER, PUBLIC, PARAMETER :: nbghostcells = 1
-   INTEGER, PUBLIC :: nbcellsx = jpiglo - 2 - 2*nbghostcells
-   INTEGER, PUBLIC :: nbcellsy = jpjglo - 2 - 2*nbghostcells
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+#if defined key_mpp_dyndist 
+   INTEGER, PUBLIC ::   jpni, jpnj, jpnij                      !: domain decomposition (AGRIF case)
 #endif
+   INTEGER, PUBLIC ::   jpi, jpj, jpim1, jpjm1, jpkm1, jpij    !: local domain size    (AGRIF case)
+
+   !! AGRIF ghost cells
+   INTEGER, PUBLIC, PARAMETER :: nbghostcells = 1                           !: ???
+   INTEGER, PUBLIC            :: nbcellsx = jpiglo - 2 - 2 * nbghostcells   !: ???
+   INTEGER, PUBLIC            :: nbcellsy = jpjglo - 2 - 2 * nbghostcells   !: ???
+
+#else
+   !!---------------------------------------------------------------------
+   !!   Default       :                         no AGRIF: model parameters
+   !!---------------------------------------------------------------------
+
+   ! Domain decomposition (NO massively parallel computer: jpni=jpnj=1 so jpi=jpiglo and jpj=jpjglo)
+   INTEGER, PUBLIC, PARAMETER ::   jpni   = 1   !: number of processors following i 
+   INTEGER, PUBLIC, PARAMETER ::   jpnj   = 1   !: number of processors following j
+   INTEGER, PUBLIC, PARAMETER ::   jpnij  = 1   !: nb of local domain = nb of processors ( <= jpni x jpnj )
+
+   !! Domain Matrix size on each processor
+   !!   data           domain   (jpidta,jpjdta)   defined in par_"config name" modules
+   !!   global or zoom domain   (jpiglo,jpjglo)   defined in par_"config name" modules
+   !!   local          domain   ( jpi  , jpj  )   defined below (= size of the domain on each processor)
+   INTEGER, PARAMETER ::   jpi = ( jpiglo-2*jpreci + (jpni-1) ) / jpni + 2*jpreci   !: first  dimension
+   INTEGER, PARAMETER ::   jpj = ( jpjglo-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj   !: second dimension
+   INTEGER, PARAMETER ::   jpim1 = jpi - 1                                          !: inner domain indices
+   INTEGER, PARAMETER ::   jpjm1 = jpj - 1                                          !: inner domain indices
+   INTEGER, PARAMETER ::   jpkm1 = jpk - 1                                          !: inner domain indices
+   INTEGER, PARAMETER ::   jpij  = jpi * jpj                                        !:  jpi x jpj
+#endif
+
+
    !!---------------------------------------------------------------------
    !! Optimization/control flags
@@ -206,5 +82,4 @@
    LOGICAL, PUBLIC, PARAMETER ::   lk_esopa     = .FALSE.  !: flag to activate the all options
 #endif
-
 #if defined key_vectopt_memory
    LOGICAL, PUBLIC, PARAMETER ::   lk_vopt_mem  = .TRUE.   !: vector optimization flag
@@ -212,5 +87,4 @@
    LOGICAL, PUBLIC, PARAMETER ::   lk_vopt_mem  = .FALSE.  !: vector optimization flag
 #endif
-
 #if defined key_vectopt_loop
    LOGICAL, PUBLIC, PARAMETER ::   lk_vopt_loop = .TRUE.   !: vector optimization flag
@@ -218,5 +92,4 @@
    LOGICAL, PUBLIC, PARAMETER ::   lk_vopt_loop = .FALSE.  !: vector optimization flag
 #endif
-
 #if defined key_mpp_omp
    LOGICAL, PUBLIC, PARAMETER ::   lk_jki = .TRUE.   !: j-k-i loop flag
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r025.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r025.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r025.F90	(revision 778)
@@ -0,0 +1,136 @@
+MODULE par_orca_r025
+   !!======================================================================
+   !!                        ***  par_ORCA_R025  ***
+   !! Ocean Domain parameter : 0.25 degrees ORCA global ocean (0RCA_R025 configuration)
+   !!    including a 1D configuration case.
+   !!    NB: with AGRIF, values are not set as parameters
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_ORCA_R025.h90
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Id$ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_orca_r025
+   !!---------------------------------------------------------------------
+   !!   'key_orca_r025'   :                       global ocean : ORCA R025
+   !!---------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                       AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                  no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution                 !!! *  ORCA R025  *
+   CHARACTER (len=16), PARAMETER ::   cp_cfg = "orca"   !: name of the configuration 
+   INTEGER           , PARAMETER ::   jp_cfg = 025      !: resolution of the configuration (degrees)
+
+   ! input data size                                   !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  = 1442               !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  = 1021               !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  =   46               !: number of levels      ( >= jpk    ) 
+   
+#  if defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :          1D configuration of ORCA R025
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                  !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                  !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   61               !: left bottom (i,j) indices of the zoom in data domain indices
+   INTEGER, PARAMETER ::   jpjzoom =  110               !: here : (jpizoom,jpjzoom) = (61,110)  for (??W,??N)
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =  18                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = 800                !: maximum number of points per island
+   
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :             Full global domain of ORCA R025
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta            !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta            !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta            !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1               !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1               !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =    4              !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =  150              !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = 3000              !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 0   !: type of horizontal mesh
+   !                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc file
+   !                                       !  = 1 geographical mesh on the sphere with regular grid-spacing
+   !                                       !  = 2 f-plane with regular grid-spacing
+   !                                       !  = 3 beta-plane with regular grid-spacing
+   !                                       !  = 4 Mercator grid with T/U point at the equator
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator be a T point.
+   REAL(wp) ,PARAMETER ::   ppglam0  = pp_not_used   !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  = pp_not_used   !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
+   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = pp_not_used   !: zonal      grid-spacing (meters )
+   REAL(wp) ,PARAMETER ::   ppe2_m   = pp_not_used   !: meridional grid-spacing (meters )
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = pp_to_be_computed      !: ORCA r4, r2 and r05 coefficients
+   REAL(wp), PARAMETER ::   ppa0  = pp_to_be_computed      !: (default coefficients)
+   REAL(wp), PARAMETER ::   ppa1  = pp_to_be_computed      !:
+   REAL(wp), PARAMETER ::   ppkth =    23.563_wp           !: (non dimensional): gives the approximate
+   !                                                       !: layer number above which  stretching will
+   !                                                       !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     9.00000000000_wp   !: (non dimensional): stretching factor
+   !                                                       !: for the grid. The highest zacr, the smallest
+   !                                                       !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin =    6._wp             !: (meters) vertical thickness of the top layer
+   REAL(wp), PARAMETER ::   pphmax  = 5750._wp             !: (meters) Maximum depth of the ocean gdepw(jpk)
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_orca_r025
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r05.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r05.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r05.F90	(revision 778)
@@ -0,0 +1,165 @@
+MODULE par_ORCA_R05
+   !!======================================================================
+   !!                        ***  par_ORCA_R05  ***
+   !! Ocean Domain parameter : 0.5 degrees orca global ocean (0RCA_R05 configuration)
+   !!    including antarctic, arctic or 1D configurations cases.
+   !!    NB: with AGRIF, values are not set as parameters
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_ORCA_R05.h90
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Id$ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_orca_r05
+   !!---------------------------------------------------------------------
+   !!   'key_orca_r05'   :                         global ocean : ORCA R05
+   !!---------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                       AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                  no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution                 !!! *  ORCA R05  *
+   CHARACTER (len=16), PARAMETER ::   cp_cfg = "orca"   !: name of the configuration 
+   INTEGER           , PARAMETER ::   jp_cfg = 05       !: resolution of the configuration (degrees)
+
+   ! input data size                                   !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  = 722                !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  = 511                !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  =  31                !: number of levels      ( >= jpk    ) 
+   
+
+#  if defined key_antarctic
+   !------------!!--------------------------------------------------------
+   !            !!   'key_antarctic'   :         antartic zoom of ORCA R05
+   !------------!!--------------------------------------------------------
+   ! zoom domain size                                  !!! *  antarctic zoom  * 
+   INTEGER, PARAMETER ::   jpiglo  = jpidta             !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 187                !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! zoom starting position
+   INTEGER, PARAMETER ::   jpizoom =   1                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   1                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   1                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   1                !: maximum number of points per island
+
+#  elif defined key_arctic
+   !------------!!--------------------------------------------------------
+   !            !!   'key_arctic' :                 artic zoom of ORCA R05
+   !------------!!--------------------------------------------------------
+   ! zoom domain size                                  !!! *  arctic zoom  *
+   INTEGER, PARAMETER ::   jpiglo  = 562                !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta - 301 + 1   !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! zoom starting position 
+   INTEGER, PARAMETER ::   jpizoom =  81                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom = 301                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   5                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   1                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   1                !: maximum number of points per island
+
+#  elif defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :           1D configuration of ORCA R05
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                  !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                  !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   61               !: left bottom (i,j) indices of the zoom in data domain indices
+   INTEGER, PARAMETER ::   jpjzoom =  110               !: here : (jpizoom,jpjzoom) = (61,110)  for (??W,??S)
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   1                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   1                !: maximum number of points per island
+   
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :              Full global domain of ORCA R05
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta            !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta            !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta            !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1               !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1               !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =    6              !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   79              !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = 2000              !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 0   !: type of horizontal mesh
+   !                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc file
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T-point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator is a T-point.
+   REAL(wp) ,PARAMETER ::   ppglam0  = pp_not_used   !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  = pp_not_used   !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
+   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = pp_not_used   !: zonal      grid-spacing (meters )
+   REAL(wp) ,PARAMETER ::   ppe2_m   = pp_not_used   !: meridional grid-spacing (meters )
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = -4762.96143546300_wp   !: ORCA r4, r2 and r05 coefficients
+   REAL(wp), PARAMETER ::   ppa0  =   255.58049070440_wp   !: (default coefficients)
+   REAL(wp), PARAMETER ::   ppa1  =   245.58132232490_wp   !:
+   REAL(wp), PARAMETER ::   ppkth =    21.43336197938_wp   !: (non dimensional): gives the approximate
+   !                                                       !: layer number above which  stretching will
+   !                                                       !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     3.00000000000_wp   !: (non dimensional): stretching factor
+   !                                                       !: for the grid. The highest zacr, the smallest
+   !                                                       !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: (meters) vertical thickness of the top layer
+   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: (meters) Maximum depth of the ocean gdepw(jpk)
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_ORCA_R05
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r2.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r2.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r2.F90	(revision 778)
@@ -0,0 +1,165 @@
+MODULE par_orca_r2
+   !!======================================================================
+   !!                        ***  par_orca_r2  ***
+   !! Ocean Domain parameter : 2 degrees orca global ocean (0RCA_R2 configuration)
+   !!    including antarctic, arctic or 1D configuration cases.
+   !!    NB: with AGRIF, values are not set as parameters
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_ORCA_R2.h90
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Id$ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_orca_r2
+   !!---------------------------------------------------------------------
+   !!   'key_orca_r2'   :                           global ocean : ORCA R2
+   !!---------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                        AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                   no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution                 !!! *  ORCA R2  *
+   CHARACTER (len=16), PARAMETER ::   cp_cfg = "orca"   !: name of the configuration 
+   INTEGER           , PARAMETER ::   jp_cfg = 2        !: resolution of the configuration (degrees)
+
+   ! input data size                                   !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  = 182                !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  = 149                !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  = 31                 !: number of levels      ( >= jpk    ) 
+   
+
+#  if defined key_antarctic
+   !------------!!--------------------------------------------------------
+   !            !!   'key_antarctic'   :          antartic zoom of ORCA R2
+   !------------!!--------------------------------------------------------
+   ! zoom domain size                                  !!! *  antarctic zoom  * 
+   INTEGER, PARAMETER ::   jpiglo  = jpidta             !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 50                 !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! zoom starting position
+   INTEGER, PARAMETER ::   jpizoom =   1                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   1                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   3                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = 400                !: maximum number of points per island
+
+#  elif defined key_arctic
+   !------------!!--------------------------------------------------------
+   !            !!   'key_arctic' :                  artic zoom of ORCA R2
+   !------------!!--------------------------------------------------------
+   ! zoom domain size                                  !!! *  arctic zoom  *
+   INTEGER, PARAMETER ::   jpiglo  = 142                !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta - 97 + 1    !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! zoom starting position 
+   INTEGER, PARAMETER ::   jpizoom =  21                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =  97                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   3                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   7                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = 400                !: maximum number of points per island
+
+#  elif defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :            1D configuration of ORCA R2
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                  !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                  !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   61               !: left bottom (i,j) indices of the zoom in data domain indices
+   INTEGER, PARAMETER ::   jpjzoom =  110               !: here : (jpizoom,jpjzoom) = (61,110)  for (160W,50N)
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =  18                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = 800                !: maximum number of points per island
+   
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :               Full global domain of ORCA R2
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta             !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta             !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   4                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =  18                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = 800                !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 0   !: type of horizontal mesh
+   !                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc file
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T-point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator is a T-point.
+   REAL(wp) ,PARAMETER ::   ppglam0  = pp_not_used   !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  = pp_not_used   !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
+   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = pp_not_used   !: zonal      grid-spacing (meters )
+   REAL(wp) ,PARAMETER ::   ppe2_m   = pp_not_used   !: meridional grid-spacing (meters )
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = -4762.96143546300_wp   !: ORCA r4, r2 and r05 coefficients
+   REAL(wp), PARAMETER ::   ppa0  =   255.58049070440_wp   !: (default coefficients)
+   REAL(wp), PARAMETER ::   ppa1  =   245.58132232490_wp   !:
+   REAL(wp), PARAMETER ::   ppkth =    21.43336197938_wp   !: (non dimensional): gives the approximate
+   !                                                       !: layer number above which  stretching will
+   !                                                       !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     3.00000000000_wp   !: (non dimensional): stretching factor
+   !                                                       !: for the grid. The highest zacr, the smallest
+   !                                                       !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: (meters) vertical thickness of the top layer
+   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: (meters) Maximum depth of the ocean gdepw(jpk)
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_orca_r2
Index: /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r4.F90
===================================================================
--- /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r4.F90	(revision 778)
+++ /branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r4.F90	(revision 778)
@@ -0,0 +1,132 @@
+MODULE par_orca_r4
+   !!======================================================================
+   !!                        ***  par_ORCA_R4  ***
+   !! Ocean Domain parameter : 4 degrees orca global ocean (0RCA_R4 configuration)
+   !!    including a 1D configuration case.
+   !!    NB: with AGRIF, values are not set as parameters
+   !!======================================================================
+   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_ORCA_R4.h90
+   !!----------------------------------------------------------------------
+   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 
+   !! $Id$ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
+#if defined key_orca_r4
+   !!---------------------------------------------------------------------
+   !!   'key_orca_r4'   :                           global ocean : ORCA R4
+   !!---------------------------------------------------------------------
+   USE par_kind          ! kind parameters
+
+   IMPLICIT NONE
+   PUBLIC
+
+# if defined key_agrif
+   !------!!--------------------------------------------------------------
+   !      !!   'key_agrif'   :                                       AGRIF 
+   !------!!--------------------------------------------------------------
+
+   ! Configuration, and domain size are not assigned to a parameter value
+   !    as they are different in mother and children models
+   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
+   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
+   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
+   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
+
+# else
+   !------!!--------------------------------------------------------------
+   !      !!   Default       :                  no AGRIF: model parameters
+   !------!!--------------------------------------------------------------
+
+   ! Configuration name and resolution                 !!! *  ORCA R4  *
+   CHARACTER (len=16), PARAMETER ::   cp_cfg = "orca"   !: name of the configuration 
+   INTEGER, PARAMETER            ::   jp_cfg = 4        !: resolution of the configuration (degrees)
+
+   ! input data size                                   !!! *  size of all input files  *
+   INTEGER, PARAMETER ::   jpidta  =  92                !: 1st lateral dimension ( >= jpiglo )
+   INTEGER, PARAMETER ::   jpjdta  =  76                !: 2nd    "       "      ( >= jpjglo )
+   INTEGER, PARAMETER ::   jpkdta  =  31                !: number of levels      ( >= jpk    ) 
+
+#  if defined key_cfg_1d
+   !------------!!--------------------------------------------------------
+   !            !!   'key_cfg_1d' :            1D configuration of ORCA R4
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  1D configuration  *
+   INTEGER, PARAMETER ::   jpiglo  = 3                  !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = 3                  !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   61               !: left bottom (i,j) indices of the zoom in data domain indices
+   INTEGER, PARAMETER ::   jpjzoom =   30               !: here : (jpizoom,jpjzoom) = (61,30)  for (??W,??N)
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   0                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =   1                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  =   1                !: maximum number of points per island
+
+#  else
+   !------------!!--------------------------------------------------------
+   !            !!   Default :               Full global domain of ORCA R4
+   !------------!!--------------------------------------------------------
+   ! global domain size                                !!! *  global domain  *
+   INTEGER, PARAMETER ::   jpiglo  = jpidta             !: 1st dimension of global domain --> i
+   INTEGER, PARAMETER ::   jpjglo  = jpjdta             !: 2nd    "                  "    --> j
+   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
+   ! starting position of the zoom 
+   INTEGER, PARAMETER ::   jpizoom =   1                !: left bottom (i,j) indices of the zoom
+   INTEGER, PARAMETER ::   jpjzoom =   1                !: in data domain indices
+   ! Domain characteristics
+   INTEGER, PARAMETER ::   jperio  =   4                !: lateral cond. type (between 0 and 6)
+   INTEGER, PARAMETER ::   jpisl   =  11                !: number of islands
+   INTEGER, PARAMETER ::   jpnisl  = 400                !: maximum number of points per island
+#  endif
+
+# endif
+
+   !!---------------------------------------------------------------------
+   !! grid parameters for domhgr and domzgr
+   !!---------------------------------------------------------------------
+   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
+   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
+   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
+   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
+
+   ! horizontal mesh (domhgr)
+   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
+   INTEGER, PARAMETER ::   jphgr_msh = 0   !: type of horizontal mesh
+   !                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc file
+
+   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T-point of the grid.
+   !   The mercator grid starts only approximately at gphi0 because
+   !   of the constraint that the equator is a T-point.
+   REAL(wp) ,PARAMETER ::   ppglam0  = pp_not_used   !: longitude of first raw and column T-point (jphgr_msh=1)
+   REAL(wp) ,PARAMETER ::   ppgphi0  = pp_not_used   !: latitude  of first raw and column T-point (jphgr_msh=1)
+   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
+   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
+   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
+   !
+   REAL(wp) ,PARAMETER ::   ppe1_m   = pp_not_used   !: zonal      grid-spacing (meters )
+   REAL(wp) ,PARAMETER ::   ppe2_m   = pp_not_used   !: meridional grid-spacing (meters )
+
+   ! Vertical mesh (domzgr)
+   !    parameter used in the analytical function defining the repartition of the level and their thickness
+   REAL(wp), PARAMETER ::   ppsur = -4762.96143546300_wp   !: ORCA r4, r2 and r05 coefficients
+   REAL(wp), PARAMETER ::   ppa0  =   255.58049070440_wp   !: (default coefficients)
+   REAL(wp), PARAMETER ::   ppa1  =   245.58132232490_wp   !:
+   REAL(wp), PARAMETER ::   ppkth =    21.43336197938_wp   !: (non dimensional): gives the approximate
+   !                                                       !: layer number above which  stretching will
+   !                                                       !: be maximum. Usually of order jpk/2.
+   REAL(wp), PARAMETER ::   ppacr =     3.00000000000_wp   !: (non dimensional): stretching factor
+   !                                                       !: for the grid. The highest zacr, the smallest
+   !                                                       !: the stretching.
+   !  If both ppa0 ppa1 and ppsur are specified to 0, then
+   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
+   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: (meters) vertical thickness of the top layer
+   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: (meters) Maximum depth of the ocean gdepw(jpk)
+
+#else
+   !!---------------------------------------------------------------------
+   !!   default option  :                                     Empty module
+   !!---------------------------------------------------------------------
+#endif
+
+   !!======================================================================
+END MODULE par_orca_r4
