Changeset 32


Ignore:
Timestamp:
03/14/12 11:41:39 (12 years ago)
Author:
cholod
Message:

par_TROP12.h90 => par_TROP12.h90 + par_depth.h90 (also for BB12,PERU12)

Location:
trunk/NEMOGCM/NEMO/OPA_SRC
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/par_BB12.h90

    r17 r32  
    1  !!--------------------------------------------------------------------- 
    2   !!                     ***  par_BB12.h90  ***   
    3   !!   Ocean Domain : 0.12 degrees resolution bay of bengal 
    4   !!                  (BB12 configuration) 
    5   !!--------------------------------------------------------------------- 
    6   !!---------------------------------------------------------------------- 
    7   !!  OPA 9.0 , LOCEAN-IPSL (2005)  
    8   !! $Id$ 
    9   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    10   !!---------------------------------------------------------------------- 
    11   CHARACTER (len=16)      & 
     1   !!--------------------------------------------------------------------- 
     2   !!                     ***  par_BB12.h90  *** 
     3   !!   Ocean Domain : 1/12 degrees resolution bay of bengal 
     4   !!                  (BB12 configuration) 
     5   !!--------------------------------------------------------------------- 
     6   !!---------------------------------------------------------------------- 
     7   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     8   !! $Id$ 
     9   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     10   !!---------------------------------------------------------------------- 
     11 
     12   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
     13   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
     14   REAL,PARAMETER      ::  pp_not_used = 999999_wp , & 
     15      &                    pp_to_be_computed = 0._wp 
     16 
     17   !!--------------------------------------------------------------------- 
     18   !!   Vertical Domain 
     19   !!--------------------------------------------------------------------- 
     20# include "par_depth.h90" 
     21 
     22   !!--------------------------------------------------------------------- 
     23   !!   Horizontal Domain 
     24   !!--------------------------------------------------------------------- 
     25   CHARACTER (len=16)      & 
    1226#if !defined key_agrif 
    13      , PARAMETER  & 
     27      , PARAMETER  & 
    1428#endif 
    15      ::    &   
    16      cp_cfg = "bengal"           !: name of the configuration 
    17   INTEGER     & 
     29      ::    &   
     30      cp_cfg = "bengal"         !: name of the configuration 
     31   INTEGER     & 
    1832#if !defined key_agrif 
    19      , PARAMETER  & 
     33      , PARAMETER  & 
    2034#endif 
    21      :: & 
    22      jp_cfg = 12  ,        &  !: resolution of the configuration (degrees) 
    23      ! Original data size 
    24      jpidta  = 281,        &  !: first horizontal dimension > or = to jpi 
    25      jpjdta  = 256,        &  !: second                     > or = to jpj 
    26      jpkdta  = 50 ,         &  !: number of levels           > or = to jpk 
    27      ! total domain matrix size 
    28      jpiglo  = jpidta,      &  !: first  dimension of global domain --> i 
    29      jpjglo  = jpjdta,      &  !: second dimension of global domain --> j 
    30      ! starting position of the zoom 
    31      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    32      jpjzoom =   1   ,      &  !: in data indices 
    33      ! Domain characteristics 
    34      jperio  =    0  ,      &  !: lateral cond. type (between 0 and 6) 
    35      jpisl   =  150  ,      &  !: number of islands 
    36      jpnisl  = 3000            !: maximum number of points per island 
     35      :: & 
     36      jp_cfg = 12  ,        &  !: resolution of the configuration (degrees) 
     37      ! 
     38      jpidta    = 281,  & 
     39      jpjdta    = 256   
    3740 
    38   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    39   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    40   REAL,PARAMETER      ::  pp_not_used = 999999_wp , & 
    41      &                    pp_to_be_computed = 0._wp 
    42   !! 
    43   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
    44   !! 
    45   INTEGER, PARAMETER ::     & ! 
    46      jphgr_msh = 0            !: type of horizontal mesh 
    47      !                        !  = 0 curvilinear coordinate on the sphere 
    48      !                        !      read in coordinate.nc file 
    49      !                        !  = 1 geographical mesh on the sphere 
    50      !                        !      with regular grid-spacing 
    51      !                        !  = 2 f-plane with regular grid-spacing 
    52      !                        !  = 3 beta-plane with regular grid-spacing 
    53      !                        !  = 4 Mercator grid with T/U point at the equator  with 
    54      !                        !      isotropic resolution (e1_deg) 
     41   INTEGER     & 
     42#if !defined key_agrif 
     43      , PARAMETER  & 
     44#endif 
     45      :: & 
     46      ! total domain matrix size 
     47      jpiglo  = jpidta,      &  !: first  dimension of global domain --> i 
     48      jpjglo  = jpjdta,      &  !: second dimension of global domain --> j 
     49      ! starting position of the zoom 
     50      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
     51      jpjzoom =   1   ,      &  !: in data indices 
     52      ! Domain characteristics 
     53      jperio  =  0       !: lateral cond. type (between 0 and 6) 
     54      !                  !  = 0 closed                 ;   = 1 cyclic East-West 
     55      !                  !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
     56      !                  !  = 4 cyclic East-West AND North fold T-point pivot 
     57      !                  !  = 5 North fold F-point pivot 
     58      !                  !  = 6 cyclic East-West AND North fold F-point pivot 
    5559 
    56      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    57      !   The mercator grid starts only approximately at gphi0 because 
    58      !   of the constraint that the equator be a T point. 
    59   REAL(wp), PARAMETER ::       &  ! 
    60      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    61      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    62      !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    63      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
    64      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
    65      ! 
    66      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
    67      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     60   !! 
     61   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
     62   !! 
     63   INTEGER, PARAMETER ::     & ! 
     64      jphgr_msh = 0            !: type of horizontal mesh 
     65      !                        !  = 0 curvilinear coordinate on the sphere 
     66      !                        !      read in coordinate.nc file 
     67      !                        !  = 1 geographical mesh on the sphere 
     68      !                        !      with regular grid-spacing 
     69      !                        !  = 2 f-plane with regular grid-spacing 
     70      !                        !  = 3 beta-plane with regular grid-spacing 
     71      !                        !  = 4 Mercator grid with T/U point at the equator  with 
     72      !                        !      isotropic resolution (e1_deg) 
    6873 
    69   !!  Coefficients associated with the vertical coordinate system 
    70   REAL(wp), PARAMETER  ::       & 
    71 #ifdef key_zrefsurf 
    72      &     ppsur = -8494.48_wp    ,  &   !: vertical grid Mercator 
    73      &     ppa0  =   257.609_wp    ,  &  !: (default coefficients refine in surface) 
    74      &     ppa1  =   256.819_wp    ,  &  !: 
    75      &     ppkth =    40_wp    ,  &  !: (non dimensional): gives the approximate 
    76      !                                           !: layer number above which  stretching will 
    77      !                                           !: be maximum. Usually of order jpk/2. 
    78      &     ppacr =     10.00000000000_wp , &     !: (non dimensional): stretching factor 
    79      !                                           !: for the grid. The highest zacr, the smallest 
    80      !                                           !: the stretching. 
    81      &     ppa2  =   100.760928500000_wp ,  &  ! IDC, 26/6/09 
    82      &     ppacr2 =   13.000000000000_wp ,  &  ! IDC, 26/6/09 
    83      &     ppkth2 =   48.029893720000_wp       ! IDC, 26/6/09 
    84 #else 
    85      &     ppsur = pp_to_be_computed ,  &  !: Computed in domzgr, set ppdzmin and pphmax below 
    86      &     ppa0  = pp_to_be_computed ,  &  !:    "           " 
    87      &     ppa1  = pp_to_be_computed ,  &  !:    "           " 
    88      ! 
    89      &     ppkth =  23.563_wp        ,  &  !: (non dimensional): gives the approximate 
    90      !                                     !: layer number above which  stretching will 
    91      !                                     !: be maximum. Usually of order jpk/2. 
    92      &     ppacr =    7.0000000000000_wp ,  &  !: (non dimensional): stretching factor 
    93      !                                         !: for the grid. The highest zacr, the smallest 
    94      !                                         !: the stretching. 
    95      &     ppa2  =   100.760928500000_wp ,  &  ! IDC, 26/6/09 
    96      &     ppacr2 =   13.000000000000_wp ,  &  ! IDC, 26/6/09 
    97      &     ppkth2 =   48.029893720000_wp       ! IDC, 26/6/09 
    98 #endif 
    99   !! 
    100   !!  If both ppa0 ppa1 and ppsur are specified to 0, then 
    101   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    102   !! 
    103   REAL(wp), PARAMETER ::        & 
    104      &     ppdzmin = 6._wp           ,  &  !: (meters) vertical thickness of the top layer 
    105      &     pphmax  = 5750._wp              !: (meters) Maximum depth of the ocean gdepw(jpk) 
    106    LOGICAL , PUBLIC, PARAMETER ::   ldbletanh = .false.   !: Use/do not use double tanf function for vertical coordinates 
    107   !!--------------------------------------------------------------------- 
     74      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
     75      !   The mercator grid starts only approximately at gphi0 because 
     76      !   of the constraint that the equator be a T point. 
     77   REAL(wp), PARAMETER ::       &  ! 
     78      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
     79      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
     80      !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
     81      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
     82      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
     83      ! 
     84      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
     85      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     86 
     87   !!--------------------------------------------------------------------- 
  • trunk/NEMOGCM/NEMO/OPA_SRC/par_PERU12.h90

    r24 r32  
    1  !!--------------------------------------------------------------------- 
    2   !!                     ***  par_BB12.h90  ***   
    3   !!   Ocean Domain : 0.12 degrees resolution bay of bengal 
    4   !!                  (BB12 configuration) 
    5   !!--------------------------------------------------------------------- 
    6   !!---------------------------------------------------------------------- 
    7   !!  OPA 9.0 , LOCEAN-IPSL (2005)  
    8   !! $Id:$ 
    9   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    10   !!---------------------------------------------------------------------- 
    11   CHARACTER (len=16)      & 
     1   !!--------------------------------------------------------------------- 
     2   !!                     ***  par_peru12.h90  *** 
     3   !!   Ocean Domain : 1/12 degrees resolution Peru/Chile Configuration 
     4   !!                  (PERU12 configuration) 
     5   !!--------------------------------------------------------------------- 
     6   !!---------------------------------------------------------------------- 
     7   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     8   !! $Id: par_ORCA_R2.h90 2715 2011-03-30 15:58:35Z rblod $ 
     9   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     10   !!---------------------------------------------------------------------- 
     11 
     12   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
     13   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
     14   REAL,PARAMETER      ::  pp_not_used = 999999_wp , & 
     15      &                    pp_to_be_computed = 0._wp 
     16 
     17   !!--------------------------------------------------------------------- 
     18   !!   Vertical Domain 
     19   !!--------------------------------------------------------------------- 
     20# include "par_depth.h90" 
     21 
     22   !!--------------------------------------------------------------------- 
     23   !!   Horizontal Domain 
     24   !!--------------------------------------------------------------------- 
     25   CHARACTER (len=16)      & 
    1226#if !defined key_agrif 
    13      , PARAMETER  & 
     27      , PARAMETER  & 
    1428#endif 
    15      ::    &   
    16      cp_cfg = "peru"           !: name of the configuration 
    17   INTEGER     & 
     29      ::    &   
     30      cp_cfg = "peru"         !: name of the configuration 
     31   INTEGER     & 
    1832#if !defined key_agrif 
    19      , PARAMETER  & 
     33      , PARAMETER  & 
    2034#endif 
    21      :: & 
    22      jp_cfg = 12  ,        &  !: resolution of the configuration (degrees) 
    23      ! Original data size 
    24      jpidta  = 361,        &  !: first horizontal dimension > or = to jpi 
    25      jpjdta  = 505,        &  !: second                     > or = to jpj 
    26      jpkdta  = 46 ,         &  !: number of levels           > or = to jpk 
    27      ! total domain matrix size 
    28      jpiglo  = jpidta,      &  !: first  dimension of global domain --> i 
    29      jpjglo  = jpjdta,      &  !: second dimension of global domain --> j 
    30      jpk     = jpkdta,      &  !: number of vertical levels 
    31      ! starting position of the zoom 
    32      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    33      jpjzoom =   1   ,      &  !: in data indices 
    34      ! Domain characteristics 
    35      jperio  =    0  ,      &  !: lateral cond. type (between 0 and 6) 
    36      jpisl   =  150  ,      &  !: number of islands 
    37      jpnisl  = 3000            !: maximum number of points per island 
     35      :: & 
     36      jp_cfg = 12  ,        &  !: resolution of the configuration (degrees) 
     37      ! 
     38      jpidta    = 361,  & 
     39      jpjdta    = 505   
    3840 
    39   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    40   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    41   REAL,PARAMETER      ::  pp_not_used = 999999_wp , & 
    42      &                    pp_to_be_computed = 0._wp 
    43   !! 
    44   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
    45   !! 
    46   INTEGER, PARAMETER ::     & ! 
    47      jphgr_msh = 0            !: type of horizontal mesh 
    48      !                        !  = 0 curvilinear coordinate on the sphere 
    49      !                        !      read in coordinate.nc file 
    50      !                        !  = 1 geographical mesh on the sphere 
    51      !                        !      with regular grid-spacing 
    52      !                        !  = 2 f-plane with regular grid-spacing 
    53      !                        !  = 3 beta-plane with regular grid-spacing 
    54      !                        !  = 4 Mercator grid with T/U point at the equator  with 
    55      !                        !      isotropic resolution (e1_deg) 
     41   INTEGER     & 
     42#if !defined key_agrif 
     43      , PARAMETER  & 
     44#endif 
     45      :: & 
     46      ! total domain matrix size 
     47      jpiglo  = jpidta,      &  !: first  dimension of global domain --> i 
     48      jpjglo  = jpjdta,      &  !: second dimension of global domain --> j 
     49      ! starting position of the zoom 
     50      jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
     51      jpjzoom =   1   ,      &  !: in data indices 
     52      ! Domain characteristics 
     53      jperio  =  0       !: lateral cond. type (between 0 and 6) 
     54      !                  !  = 0 closed                 ;   = 1 cyclic East-West 
     55      !                  !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
     56      !                  !  = 4 cyclic East-West AND North fold T-point pivot 
     57      !                  !  = 5 North fold F-point pivot 
     58      !                  !  = 6 cyclic East-West AND North fold F-point pivot 
    5659 
    57      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    58      !   The mercator grid starts only approximately at gphi0 because 
    59      !   of the constraint that the equator be a T point. 
    60   REAL(wp), PARAMETER ::       &  ! 
    61      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    62      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    63      !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    64      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
    65      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
    66      ! 
    67      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
    68      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     60   !! 
     61   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
     62   !! 
     63   INTEGER, PARAMETER ::     & ! 
     64      jphgr_msh = 0            !: type of horizontal mesh 
     65      !                        !  = 0 curvilinear coordinate on the sphere 
     66      !                        !      read in coordinate.nc file 
     67      !                        !  = 1 geographical mesh on the sphere 
     68      !                        !      with regular grid-spacing 
     69      !                        !  = 2 f-plane with regular grid-spacing 
     70      !                        !  = 3 beta-plane with regular grid-spacing 
     71      !                        !  = 4 Mercator grid with T/U point at the equator  with 
     72      !                        !      isotropic resolution (e1_deg) 
    6973 
    70   !!  Coefficients associated with the vertical coordinate system 
    71   REAL(wp), PARAMETER  ::       & 
    72 #ifdef key_zrefsurf 
    73      &     ppsur = -8494.48_wp    ,  &   !: vertical grid Mercator 
    74      &     ppa0  =   257.609_wp    ,  &  !: (default coefficients refine in surface) 
    75      &     ppa1  =   256.819_wp    ,  &  !: 
    76      &     ppkth =    40_wp    ,  &  !: (non dimensional): gives the approximate 
    77      !                                           !: layer number above which  stretching will 
    78      !                                           !: be maximum. Usually of order jpk/2. 
    79      &     ppacr =     10.00000000000_wp , &     !: (non dimensional): stretching factor 
    80      !                                           !: for the grid. The highest zacr, the smallest 
    81      !                                           !: the stretching. 
    82      &     ppa2  =   100.760928500000_wp ,  &  ! IDC, 26/6/09 
    83      &     ppacr2 =   13.000000000000_wp ,  &  ! IDC, 26/6/09 
    84      &     ppkth2 =   48.029893720000_wp       ! IDC, 26/6/09 
    85 #else 
    86      &     ppsur = pp_to_be_computed ,  &  !: Computed in domzgr, set ppdzmin and pphmax below 
    87      &     ppa0  = pp_to_be_computed ,  &  !:    "           " 
    88      &     ppa1  = pp_to_be_computed ,  &  !:    "           " 
    89      ! 
    90      &     ppkth =  23.563_wp        ,  &  !: (non dimensional): gives the approximate 
    91      !                                     !: layer number above which  stretching will 
    92      !                                     !: be maximum. Usually of order jpk/2. 
    93      &     ppacr =    7.0000000000000_wp ,  &  !: (non dimensional): stretching factor 
    94      !                                         !: for the grid. The highest zacr, the smallest 
    95      !                                         !: the stretching. 
    96      &     ppa2  =   100.760928500000_wp ,  &  ! IDC, 26/6/09 
    97      &     ppacr2 =   13.000000000000_wp ,  &  ! IDC, 26/6/09 
    98      &     ppkth2 =   48.029893720000_wp       ! IDC, 26/6/09 
    99 #endif 
    100   !! 
    101   !!  If both ppa0 ppa1 and ppsur are specified to 0, then 
    102   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    103   !! 
    104   REAL(wp), PARAMETER ::        & 
    105      &     ppdzmin = 6._wp           ,  &  !: (meters) vertical thickness of the top layer 
    106      &     pphmax  = 5750._wp              !: (meters) Maximum depth of the ocean gdepw(jpk) 
    107    LOGICAL , PUBLIC, PARAMETER ::   ldbletanh = .false.   !: Use/do not use double tanf function for vertical coordinates 
    108   !!--------------------------------------------------------------------- 
     74      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
     75      !   The mercator grid starts only approximately at gphi0 because 
     76      !   of the constraint that the equator be a T point. 
     77   REAL(wp), PARAMETER ::       &  ! 
     78      ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
     79      ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
     80      !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
     81      ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
     82      ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
     83      ! 
     84      ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
     85      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     86 
     87   !!--------------------------------------------------------------------- 
  • trunk/NEMOGCM/NEMO/OPA_SRC/par_TROP12.h90

    r12 r32  
    11   !!--------------------------------------------------------------------- 
    2    !!                     ***  par_trop12.h90  ***   
    3    !!   Ocean Domain : 0.25 degrees resolution Tropical Channel 
    4    !!                  (0RCA_R025 configuration) 
     2   !!                     ***  par_trop12.h90  *** 
     3   !!   Ocean Domain : 1/12 degrees resolution Tropical Channel 
     4   !!                  (TROP12 configuration) 
    55   !!--------------------------------------------------------------------- 
    66   !!---------------------------------------------------------------------- 
    77   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    8    !! $Id: par_ORCA_R2.h90 2715 2011-03-30 15:58:35Z rblod $  
     8   !! $Id: par_ORCA_R2.h90 2715 2011-03-30 15:58:35Z rblod $ 
    99   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1010   !!---------------------------------------------------------------------- 
     11 
     12   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
     13   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
     14   REAL,PARAMETER      ::  pp_not_used = 999999_wp , & 
     15      &                    pp_to_be_computed = 0._wp 
     16 
     17   !!--------------------------------------------------------------------- 
     18   !!   Vertical Domain 
     19   !!--------------------------------------------------------------------- 
     20# include "par_depth.h90" 
     21 
     22   !!--------------------------------------------------------------------- 
     23   !!   Horizontal Domain 
     24   !!--------------------------------------------------------------------- 
    1125   CHARACTER (len=16)      & 
    1226#if !defined key_agrif 
     
    2337      ! 
    2438      jpidta    = 4322,  & 
    25       jpjdta    = 1249,  &   
    26 #if key_trop12==75 
    27       jpkdta  = 75 ,         &  !: number of levels           > or = to jpk 
    28 #else 
    29       jpkdta  = 46 ,         &  !: number of levels           > or = to jpk 
     39      jpjdta    = 1249   
     40 
     41   INTEGER     & 
     42#if !defined key_agrif 
     43      , PARAMETER  & 
    3044#endif 
     45      :: & 
    3146      ! total domain matrix size 
    3247      jpiglo  = jpidta,      &  !: first  dimension of global domain --> i 
     
    3651      jpjzoom =   1   ,      &  !: in data indices 
    3752      ! Domain characteristics 
    38       jperio  =   1             !: lateral cond. type (between 0 and 6) 
     53      jperio  =  1       !: lateral cond. type (between 0 and 6) 
     54      !                  !  = 0 closed                 ;   = 1 cyclic East-West 
     55      !                  !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
     56      !                  !  = 4 cyclic East-West AND North fold T-point pivot 
     57      !                  !  = 5 North fold F-point pivot 
     58      !                  !  = 6 cyclic East-West AND North fold F-point pivot 
    3959 
    40    !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    41    !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    42    REAL,PARAMETER      ::  pp_not_used = 999999_wp , & 
    43       &                    pp_to_be_computed = 0._wp 
    4460   !! 
    4561   !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
     
    6985      ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
    7086 
    71    !!  Coefficients associated with the vertical coordinate system 
    72  
    73 #if key_trop12==75 
    74    REAL(wp), PARAMETER  ::       & 
    75       &     ppsur =  -3958.951371276829_wp  ,  &  !: ORCA r1 coefficients 
    76       &     ppa0  =   103.9530096000000_wp  ,  &  !: (75 levels case) 
    77       &     ppa1  =   2.415951269000000_wp  ,  &  !: 
    78       &     ppkth =   15.35101370000000_wp  ,  &  !: (non dimensional): gives the approximate 
    79       !                                           !: layer number above which  stretching will 
    80       !                                           !: be maximum. Usually of order jpk/2. 
    81       &     ppacr =       7.00000000000_wp        !: (non dimensional): stretching factor 
    82       !                                           !: for the grid. The higher zacr, the smaller 
    83       !                                           !: the stretching. 
    84    !! 
    85    !!  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
    86    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    87    !! 
    88    REAL(wp), PARAMETER ::                      & 
    89       &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer 
    90       &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk) 
    91    !! 
    92    LOGICAL,  PARAMETER ::                      & 
    93       &     ldbletanh = .TRUE.                    !: Use/do not use double tanf function for vertical coordinates 
    94    REAL(wp), PARAMETER ::                      & 
    95       &     ppa2  =   100.7609285000000_wp  ,  &  !: Double tanh function parameters 
    96       &     ppkth2=   48.02989372000000_wp  ,  &  !: 
    97       &     ppacr2=    13.00000000000_wp          !: 
    98       ! 
    99 #else 
    100    REAL(wp), PARAMETER  ::       & 
    101       &     ppsur = pp_to_be_computed ,        &  !: Computed in domzgr, set ppdzmin and pphmax below 
    102       &     ppa0  = pp_to_be_computed ,        &  !:    "           " 
    103       &     ppa1  = pp_to_be_computed ,        &  !:    "           " 
    104       ! 
    105       &     ppkth =  23.563_wp        ,        &  !: (non dimensional): gives the approximate 
    106       !                                           !: layer number above which  stretching will 
    107       !                                           !: be maximum. Usually of order jpk/2. 
    108       &     ppacr =   9.00000000000_wp            !: (non dimensional): stretching factor 
    109       !                                           !: for the grid. The highest zacr, the smallest 
    110       !                                           !: the stretching. 
    111    !! 
    112    !!  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
    113    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    114    !! 
    115    REAL(wp), PARAMETER ::                      & 
    116       &     ppdzmin = 6._wp           ,        &  !: (meters) vertical thickness of the top layer 
    117       &     pphmax  = 5750._wp                    !: (meters) Maximum depth of the ocean gdepw(jpk) 
    118    !! 
    119    LOGICAL,  PARAMETER ::                      & 
    120       &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates 
    121    REAL(wp), PARAMETER ::                      & 
    122       &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters 
    123       &     ppkth2  = pp_not_used           ,  &  !: 
    124       &     ppacr2  = pp_not_used                 !: 
    125 #endif 
    12687   !!--------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.