New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4144 for branches/2013/dev_r3987_UKMO6_C1D/NEMOGCM/NEMO/OPA_SRC/par_ORCA_R05.h90 – NEMO

Ignore:
Timestamp:
2013-10-28T14:50:08+01:00 (10 years ago)
Author:
rfurner
Message:

Commit for 2013 changes; see #1085

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3987_UKMO6_C1D/NEMOGCM/NEMO/OPA_SRC/par_ORCA_R05.h90

    r2715 r4144  
    99   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1010   !!---------------------------------------------------------------------- 
    11    CHARACTER (len=16)      & 
     11   CHARACTER (len=16)         & 
    1212#if !defined key_agrif 
    13       , PARAMETER  & 
     13      , PARAMETER             & 
    1414#endif 
    15       ::     
    16       cp_cfg = "orca"           !: name of the configuration 
    17    INTEGER     & 
     15      ::                       
     16      cp_cfg  =  "orca"          !: Name of the configuration 
     17   INTEGER                    & 
    1818#if !defined key_agrif 
    19       , PARAMETER  & 
     19      , PARAMETER             & 
    2020#endif 
    21       :: & 
    22       jp_cfg = 05  ,         &  !: resolution of the configuration (degrees) 
     21      ::                      & 
     22      jp_cfg  =  05    ,      &  !: Resolution of the configuration (degrees) 
    2323 
    24       ! data size              !!! * size of all the input files * 
    25       jpidta  = 722,         &  !: 1st lateral dimension > or = to jpiglo 
    26       jpjdta  = 511,         &  !: 2nd   "         "     > or = to jpjglo 
    27       jpkdta  =  31             !: number of levels      > or = to jpkglo 
     24      ! Data domain size         !!! *  Size of all input files * 
     25      jpidta  =  722   ,      &  !: 1st lateral dimension ( >= jpiglo ) 
     26      jpjdta  =  511   ,      &  !: 2nd lateral dimension ( >= jpjglo ) 
     27      jpkdta  =  31    ,      &  !: Number of levels      ( >= jpk    ) 
    2828 
    29 #if defined key_antarctic 
    30       ! zoom domain size       !!! *  antarctic zoom  *  
    31    INTEGER     & 
    32 #if !defined key_agrif 
    33       , PARAMETER  & 
     29#if defined key_c1d 
     30      ! Zoom domain size         !!! *  C1D zoom  * 
     31      jpiglo  =  3     ,      &  !: 1st dimension of global domain --> i 
     32      jpjglo  =  3     ,      &  !: 2nd dimension of global domain --> j 
     33      ! Domain characteristics 
     34      jperio  =  0               !: Lateral cond. type (between 0 and 6) 
     35 
     36   INTEGER                    & 
     37      ::                      & 
     38      ! Starting position of the zoom 
     39      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom 
     40      jpjzoom =  1               !: in data domain indices 
     41#elif defined key_antarctic 
     42      ! Zoom domain size         !!! *  Antarctic zoom  *  
     43      jpiglo  =  jpidta,      &  !: 1st dimension of global domain --> i 
     44      jpjglo  =  187   ,      &  !: 2nd dimension of global domain --> j 
     45      ! Starting position of the zoom 
     46      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom 
     47      jpjzoom =  1     ,      &  !: in data domain indices 
     48      ! Domain characteristics 
     49      jperio  =  1               !: Lateral cond. type (between 0 and 6) 
     50#elif defined key_arctic 
     51      ! Zoom domain size         !!! *  Arctic zoom  * 
     52      jpiglo  =  562   ,      &  !: 1st dimension of global domain --> i 
     53      jpjglo  =  jpjdta-301+1,&  !: 2nd dimension of global domain --> j 
     54      ! Starting position of the zoom 
     55      jpizoom =  81    ,      &  !: Left bottom (i,j) indices of the zoom 
     56      jpjzoom =  301   ,      &  !: in data domain indices 
     57      ! Domain characteristics 
     58      jperio  =  5               !: Lateral cond. type (between 0 and 6) 
     59#else 
     60      ! Global domain size       !!! *  Global domain  * 
     61      jpiglo  =  jpidta,      &  !: 1st dimension of global domain --> i 
     62      jpjglo  =  jpjdta,      &  !: 2nd dimension of global domain --> j 
     63      ! Starting position of the zoom 
     64      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom 
     65      jpjzoom =  1     ,      &  !: in data domain indices 
     66      ! Domain characteristics 
     67      jperio  =  6               !: Lateral cond. type (between 0 and 6) 
    3468#endif 
    35       :: & 
    36       jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i 
    37       jpjglo  = 187   ,      &  !: 2nd     "                 "    --> j  
    38       ! starting position of the zoom 
    39       jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    40       jpjzoom =   1   ,      &  !: in data domain indices 
    41       ! Domain characteristics 
    42       jperio  =   1             !: lateral cond. type (between 0 and 6) 
    4369 
    44 #elif defined key_arctic 
    45       ! zoom domain size       !!! *  arctic zoom  * 
    46    INTEGER    & 
    47 #if !defined key_agrif 
    48       , PARAMETER  & 
    49 #endif 
    50       :: & 
    51       ! zoom domain size       !!! *  arctic zoom  * 
    52       jpiglo  = 562,         &  !: 1st dimension of global domain --> i 
    53       jpjglo  = jpjdta-301+1,&  !: 2nd     "                 "    --> j 
    54       ! zoom starting position 
    55       jpizoom =  81   ,      &  !: left bottom (i,j) indices of the zoom 
    56       jpjzoom = 301   ,      &  !: in data domain indices 
    57       ! Domain characteristics 
    58       jperio  =   5             !: lateral cond. type (between 0 and 6) 
    59  
    60 #else 
    61       ! global domain size     !!! *  global domain  * 
    62    INTEGER    & 
    63 #if !defined key_agrif 
    64       , PARAMETER  & 
    65 #endif 
    66       :: & 
    67       jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i 
    68       jpjglo  = jpjdta,      &  !: 2nd     "                 "    --> j 
    69       ! zoom starting position     
    70       jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    71       jpjzoom =   1   ,      &  !: in data domain indices 
    72       ! Domain characteristics 
    73       jperio  =    6            !: lateral cond. type (between 0 and 6) 
    74 #endif 
    7570 
    7671   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    7772   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    78    REAL(wp), PARAMETER ::   & 
     73 
     74   REAL(wp), PARAMETER ::              & 
    7975      pp_not_used       = 999999._wp , &  !: 
    80       pp_to_be_computed = 0._wp          !: 
     76      pp_to_be_computed = 0._wp           !: 
    8177 
    82    !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
    8378 
    84    INTEGER, PARAMETER   ::   & ! 
    85       jphgr_msh = 0            !: type of horizontal mesh 
    86       !                        !  = 0 curvilinear coordinate on the sphere 
    87       !                        !      read in coordinate.nc file 
    88       !                        !  = 1 geographical mesh on the sphere 
    89       !                        !      with regular grid-spacing 
    90       !                        !  = 2 f-plane with regular grid-spacing 
    91       !                        !  = 3 beta-plane with regular grid-spacing 
    92       !                        !  = 4 Mercator grid with T/U point at the equator  with 
    93       !                        !      isotropic resolution (e1_deg) 
     79   !! Coefficients associated with the horizontal coordinate system 
     80 
     81   INTEGER, PARAMETER  ::     & 
     82      jphgr_msh = 0              !: type of horizontal mesh 
     83      !                          !: = 0 curvilinear coordinate on the sphere 
     84      !                          !:     read in coordinate.nc file 
     85      !                          !: = 1 geographical mesh on the sphere 
     86      !                          !:     with regular grid-spacing 
     87      !                          !: = 2 f-plane with regular grid-spacing 
     88      !                          !: = 3 beta-plane with regular grid-spacing 
     89      !                          !: = 4 Mercator grid with T/U point at the equator  with 
     90      !                          !:     isotropic resolution (e1_deg) 
    9491 
    9592      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    9693      !   The mercator grid starts only approximately at gphi0 because 
    9794      !   of the constraint that the equator be a T point. 
    98    REAL(wp) , PARAMETER ::      &  ! 
    99       ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    100       ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    101       !                            ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    102       ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
    103       ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
     95 
     96   REAL(wp), PARAMETER ::     & 
     97      ppglam0  = pp_not_used, &  !: longitude of first raw and column T-point   (jphgr_msh = 1) 
     98      ppgphi0  = pp_not_used, &  !: latitude  of first raw and column T-point   (jphgr_msh = 1) 
     99      !                          !: latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
     100      ppe1_deg = pp_not_used, &  !: zonal      grid-spacing (degrees) 
     101      ppe2_deg = pp_not_used, &  !: meridional grid-spacing (degrees) 
    104102      ! 
    105       ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
    106       ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     103      ppe1_m   = pp_not_used, &  !: zonal      grid-spacing (meters ) 
     104      ppe2_m   = pp_not_used     !: meridional grid-spacing (meters ) 
    107105 
    108    !! 
    109    !! Vertical grid parameter for domzgr 
    110    !! ===================================== 
    111    !! 
    112    REAL(wp), PARAMETER  ::       & 
    113       &     ppsur = -4762.96143546300_wp    ,  &  !: ORCA r4, r2 and r05 coefficients 
     106 
     107   !!  Coefficients associated with the vertical coordinate system 
     108 
     109   REAL(wp), PARAMETER  ::                     & 
     110      &     ppsur =   -4762.96143546300_wp  ,  &  !: ORCA r4, r2 and r05 coefficients 
    114111      &     ppa0  =   255.58049070440_wp    ,  &  !: (default coefficients) 
    115112      &     ppa1  =   245.58132232490_wp    ,  &  !: 
    116       &     ppkth =    21.43336197938_wp    ,  &  !: (non dimensional): gives the approximate 
     113      &     ppkth =   21.43336197938_wp     ,  &  !: (non dimensional): gives the approximate 
    117114      !                                           !: layer number above which  stretching will 
    118115      !                                           !: be maximum. Usually of order jpk/2. 
    119       &     ppacr =     3.00000000000_wp          !: (non dimensional): stretching factor 
    120       !                                           !: for the grid. The highest zacr, the smallest 
     116      &     ppacr =   3.00000000000_wp            !: (non dimensional): stretching factor 
     117      !                                           !: for the grid. The higher zacr, the smaller 
    121118      !                                           !: the stretching. 
    122119 
    123    !! 
    124    !!  If both ppa0 ppa1 and ppsur are specified to 0, then 
    125    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    126    !! 
    127    REAL(wp), PARAMETER ::        & 
     120      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
     121      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
     122 
     123   REAL(wp), PARAMETER ::                      & 
    128124      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer 
    129125      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk) 
    130    LOGICAL,  PARAMETER ::        & 
     126   LOGICAL,  PARAMETER ::                      & 
    131127      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates 
    132    REAL(wp), PARAMETER ::        & 
     128   REAL(wp), PARAMETER ::                      & 
    133129      &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters 
    134130      &     ppkth2  = pp_not_used           ,  &  !: 
Note: See TracChangeset for help on using the changeset viewer.