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_R2.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_R2.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 = 2,            &  !: resolution of the configuration (degrees) 
     21      ::                      & 
     22      jp_cfg  =  2     ,      &  !: Resolution of the configuration (degrees) 
    2323 
    24       ! data size              !!! * size of all input files * 
    25       jpidta  = 182,         &  !: 1st lateral dimension ( >= jpiglo ) 
    26       jpjdta  = 149,         &  !: 2nd    "       "      ( >= jpjglo ) 
    27       jpkdta  = 31              !: number of levels      ( >= jpk    )  
     24      ! Data domain size         !!! *  Size of all input files * 
     25      jpidta  =  182   ,      &  !: 1st lateral dimension ( >= jpiglo ) 
     26      jpjdta  =  149   ,      &  !: 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  =  50    ,      &  !: 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  =  142   ,      &  !: 1st dimension of global domain --> i 
     53      jpjglo  =  jpjdta-97+1, &  !: 2nd dimension of global domain --> j 
     54      ! Starting position of the zoom 
     55      jpizoom =  21    ,      &  !: Left bottom (i,j) indices of the zoom 
     56      jpjzoom =  97    ,      &  !: in data domain indices 
     57      ! Domain characteristics 
     58      jperio  =  3               !: 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  =  4               !: Lateral cond. type (between 0 and 6) 
    3468#endif 
    35       :: & 
    36       jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i 
    37       jpjglo  = 50,          &  !: 2nd    "                  "    --> j 
    38       ! zoom starting position 
    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       jpiglo  = 142   ,      &  !: 1st dimension of global domain --> i 
    52       jpjglo  = jpjdta-97+1, &  !: 2nd    "                  "    --> j 
    53       ! zoom starting position  
    54       jpizoom =  21   ,      &  !: left bottom (i,j) indices of the zoom 
    55       jpjzoom =  97   ,      &  !: in data domain indices 
    56       ! Domain characteristics 
    57       jperio  =   3             !: lateral cond. type (between 0 and 6) 
    58  
    59 #elif defined key_c1d 
    60       ! global domain size     !!! *  global domain  * 
    61    INTEGER    & 
    62 #if !defined key_agrif 
    63       , PARAMETER  & 
    64 #endif 
    65       :: & 
    66       jpiglo  = 3     ,      &  !: 1st dimension of global domain --> i 
    67       jpjglo  = 3     ,      &  !: 2nd    "                  "    --> j 
    68       ! starting position of the zoom  
    69       ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    70       ! jpjzoom =   133  ,    &  !: in data domain indices (160W,75N) 
    71       ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    72       ! jpjzoom =   110  ,    &  !: in data domain indices (160W,50N) 
    73       ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    74       ! jpjzoom =   97   ,    &  !: in data domain indices (160W,30N) 
    75       ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    76       ! jpjzoom =   86   ,    &  !: in data domain indices (160W,10N) 
    77       ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    78       ! jpjzoom =   49   ,    &  !: in data domain indices (160W,30S) 
    79       ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    80       ! jpjzoom =   27   ,    &  !: in data domain indices (160W,60S) 
    81       ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    82       ! jpjzoom =    7   ,    &  !: in data domain indices (160W,75S) 
    83       jpizoom =   110   ,    &  !: left bottom (i,j) indices of the zoom 
    84       jpjzoom =   97   ,    &  !: in data domain indices (64W,31.5N) BATS site 
    85       ! Domain characteristics 
    86       jperio  =   0            !: lateral cond. type (between 0 and 6) 
    87 #else 
    88       ! global domain size     !!! *  global domain  * 
    89    INTEGER    & 
    90 #if !defined key_agrif 
    91       , PARAMETER  & 
    92 #endif 
    93       :: & 
    94       jpiglo  = jpidta,      &  !: 1st dimension of global domain --> i 
    95       jpjglo  = jpjdta,      &  !: 2nd    "                  "    --> j 
    96       ! starting position of the zoom  
    97       jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    98       jpjzoom =   1   ,      &  !: in data domain indices 
    99       ! Domain characteristics 
    100       jperio  =   4             !: lateral cond. type (between 0 and 6) 
    101  
    102 #endif 
    10370 
    10471   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    10572   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    106    REAL(wp), PARAMETER ::   & 
    107       pp_not_used       = 999999_wp , &  !: 
    108       pp_to_be_computed = 0._wp          !: 
    10973 
    110    !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
     74   REAL(wp), PARAMETER ::              & 
     75      pp_not_used       = 999999._wp , &  !: 
     76      pp_to_be_computed = 0._wp           !: 
    11177 
    112    INTEGER,PARAMETER   ::    & ! 
    113       jphgr_msh = 0            !: type of horizontal mesh 
    114       !                        !  = 0 curvilinear coordinate on the sphere 
    115       !                        !      read in coordinate.nc file 
    116       !                        !  = 1 geographical mesh on the sphere 
    117       !                        !      with regular grid-spacing 
    118       !                        !  = 2 f-plane with regular grid-spacing 
    119       !                        !  = 3 beta-plane with regular grid-spacing 
    120       !                        !  = 4 Mercator grid with T/U point at the equator  with 
    121       !                        !      isotropic resolution (e1_deg) 
     78 
     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) 
    12291 
    12392      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    12493      !   The mercator grid starts only approximately at gphi0 because 
    12594      !   of the constraint that the equator be a T point. 
    126    REAL(wp) ,PARAMETER ::       &  ! 
    127       ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    128       ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    129       !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    130       ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
    131       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) 
    132102      ! 
    133       ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
    134       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 ) 
    135105 
    136    !! 
    137    !! Vertical grid parameter for domzgr 
    138    !! ================================== 
    139    !! 
    140    REAL(wp), PARAMETER  ::       & 
    141       &     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 
    142111      &     ppa0  =   255.58049070440_wp    ,  &  !: (default coefficients) 
    143112      &     ppa1  =   245.58132232490_wp    ,  &  !: 
    144       &     ppkth =    21.43336197938_wp    ,  &  !: (non dimensional): gives the approximate 
     113      &     ppkth =   21.43336197938_wp     ,  &  !: (non dimensional): gives the approximate 
    145114      !                                           !: layer number above which  stretching will 
    146115      !                                           !: be maximum. Usually of order jpk/2. 
    147       &     ppacr =     3.00000000000_wp          !: (non dimensional): stretching factor 
    148       !                                           !: 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 
    149118      !                                           !: the stretching. 
    150119 
    151    !! 
    152    !!  If both ppa0 ppa1 and ppsur are specified to 0, then 
    153    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    154    !! 
    155    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 ::                      & 
    156124      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer 
    157125      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk) 
    158    LOGICAL,  PARAMETER ::        & 
     126   LOGICAL,  PARAMETER ::                      & 
    159127      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates 
    160    REAL(wp), PARAMETER ::        & 
     128   REAL(wp), PARAMETER ::                      & 
    161129      &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters 
    162130      &     ppkth2  = pp_not_used           ,  &  !: 
Note: See TracChangeset for help on using the changeset viewer.