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_R025.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_R025.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 = 025  ,        &  !: resolution of the configuration (degrees) 
    23       ! Original data size 
    24       jpidta  = 1442,        &  !: first horizontal dimension > or = to jpi 
    25       jpjdta  = 1021,        &  !: second                     > or = to jpj 
     21      ::                      & 
     22      jp_cfg  =  025   ,      &  !: Resolution of the configuration (degrees) 
     23 
     24      ! Data domain size         !!! *  Size of all input files  * 
     25      jpidta  =  1442  ,      &  !: 1st lateral dimension ( >= jpiglo ) 
     26      jpjdta  =  1021  ,      &  !: 2nd lateral dimension ( >= jpjglo ) 
    2627#if key_orca_r025==75 
    27       jpkdta  = 75 ,         &  !: number of levels           > or = to jpk 
     28      jpkdta  =  75    ,      &  !: Number of levels      ( >= jpk    ) 
    2829#else 
    29       jpkdta  = 46 ,         &  !: number of levels           > or = to jpk 
     30      jpkdta  =  46    ,      &  !: Number of levels      ( >= jpk    ) 
    3031#endif 
    31       ! total domain matrix size 
    32       jpiglo  = jpidta,      &  !: first  dimension of global domain --> i 
    33       jpjglo  = jpjdta,      &  !: second dimension of global domain --> j 
    34       ! starting position of the zoom 
    35       jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    36       jpjzoom =   1   ,      &  !: in data indices 
     32 
     33#if defined key_c1d 
     34      ! Zoom domain size         !!! *  C1D zoom  * 
     35      jpiglo  =  3     ,      &  !: 1st dimension of global domain --> i 
     36      jpjglo  =  3     ,      &  !: 2nd dimension of global domain --> j 
    3737      ! Domain characteristics 
    38       jperio  =    4            !: lateral cond. type (between 0 and 6) 
     38      jperio  =  0               !: Lateral cond. type (between 0 and 6) 
     39 
     40   INTEGER                    & 
     41      ::                      & 
     42      ! Starting position of the zoom 
     43      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom 
     44      jpjzoom =  1               !: in data domain indices 
     45#else 
     46      ! Global domain size       !!! *  Global domain  * 
     47      jpiglo  =  jpidta,      &  !: 1st dimension of global domain --> i 
     48      jpjglo  =  jpjdta,      &  !: 2nd 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 domain indices 
     52      ! Domain characteristics 
     53      jperio  =  4               !: Lateral cond. type (between 0 and 6) 
     54#endif 
     55 
    3956 
    4057   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    4158   !!  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 
    44    !! 
    45    !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
    46    !! 
    47    INTEGER, PARAMETER ::     & ! 
    48       jphgr_msh = 0            !: type of horizontal mesh 
    49       !                        !  = 0 curvilinear coordinate on the sphere 
    50       !                        !      read in coordinate.nc file 
    51       !                        !  = 1 geographical mesh on the sphere 
    52       !                        !      with regular grid-spacing 
    53       !                        !  = 2 f-plane with regular grid-spacing 
    54       !                        !  = 3 beta-plane with regular grid-spacing 
    55       !                        !  = 4 Mercator grid with T/U point at the equator  with 
    56       !                        !      isotropic resolution (e1_deg) 
     59 
     60   REAL(wp), PARAMETER ::              & 
     61      pp_not_used       = 999999._wp , &  !: 
     62      pp_to_be_computed = 0._wp           !: 
     63 
     64 
     65   !! Coefficients associated with the horizontal coordinate system 
     66 
     67   INTEGER, PARAMETER  ::     & 
     68      jphgr_msh = 0              !: type of horizontal mesh 
     69      !                          !: = 0 curvilinear coordinate on the sphere 
     70      !                          !:     read in coordinate.nc file 
     71      !                          !: = 1 geographical mesh on the sphere 
     72      !                          !:     with regular grid-spacing 
     73      !                          !: = 2 f-plane with regular grid-spacing 
     74      !                          !: = 3 beta-plane with regular grid-spacing 
     75      !                          !: = 4 Mercator grid with T/U point at the equator  with 
     76      !                          !:     isotropic resolution (e1_deg) 
    5777 
    5878      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    5979      !   The mercator grid starts only approximately at gphi0 because 
    6080      !   of the constraint that the equator be a T point. 
    61    REAL(wp), PARAMETER ::       &  ! 
    62       ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    63       ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    64       !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    65       ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
    66       ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
     81 
     82   REAL(wp), PARAMETER ::     & 
     83      ppglam0  = pp_not_used, &  !: longitude of first raw and column T-point   (jphgr_msh = 1) 
     84      ppgphi0  = pp_not_used, &  !: latitude  of first raw and column T-point   (jphgr_msh = 1) 
     85      !                          !: latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
     86      ppe1_deg = pp_not_used, &  !: zonal      grid-spacing (degrees) 
     87      ppe2_deg = pp_not_used, &  !: meridional grid-spacing (degrees) 
    6788      ! 
    68       ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
    69       ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     89      ppe1_m   = pp_not_used, &  !: zonal      grid-spacing (meters ) 
     90      ppe2_m   = pp_not_used     !: meridional grid-spacing (meters ) 
     91 
    7092 
    7193   !!  Coefficients associated with the vertical coordinate system 
    7294 
    7395#if key_orca_r025==75 
    74    REAL(wp), PARAMETER  ::       & 
    75       &     ppsur =  -3958.951371276829_wp ,  &  !: ORCA r025 coefficients 
     96   REAL(wp), PARAMETER  ::                     & 
     97      &     ppsur =   -3958.951371276829_wp ,  &  !: ORCA r025 coefficients 
    7698      &     ppa0  =   103.9530096000000_wp  ,  &  !: (75 levels case) 
    7799      &     ppa1  =   2.415951269000000_wp  ,  &  !: 
     
    79101      !                                           !: layer number above which  stretching will 
    80102      !                                           !: be maximum. Usually of order jpk/2. 
    81       &     ppacr =       7.00000000000_wp        !: (non dimensional): stretching factor 
     103      &     ppacr =   7.00000000000_wp            !: (non dimensional): stretching factor 
    82104      !                                           !: for the grid. The higher zacr, the smaller 
    83105      !                                           !: 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    !! 
     106 
     107      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
     108      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
     109 
    88110   REAL(wp), PARAMETER ::                      & 
    89111      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer 
    90112      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk) 
    91    !! 
    92113   LOGICAL,  PARAMETER ::                      & 
    93114      &     ldbletanh = .TRUE.                    !: Use/do not use double tanf function for vertical coordinates 
    94115   REAL(wp), PARAMETER ::                      & 
    95       &     ppa2  =   100.7609285000000_wp  ,  &  !: Double tanh function parameters 
    96       &     ppkth2=   48.02989372000000_wp  ,  &  !: 
    97       &     ppacr2=    13.00000000000_wp          !: 
    98       ! 
     116      &     ppa2    = 100.7609285000000_wp  ,  &  !: Double tanh function parameters 
     117      &     ppkth2  = 48.02989372000000_wp  ,  &  !: 
     118      &     ppacr2  = 13.00000000000_wp           !: 
    99119#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 
     120   REAL(wp), PARAMETER  ::                     & 
     121      &     ppsur =   pp_to_be_computed     ,  &  !: Computed in domzgr, set ppdzmin and pphmax below 
     122      &     ppa0  =   pp_to_be_computed     ,  &  !:    "           " 
     123      &     ppa1  =   pp_to_be_computed     ,  &  !:    "           " 
     124      &     ppkth =   23.563_wp             ,  &  !: (non dimensional): gives the approximate 
    106125      !                                           !: layer number above which  stretching will 
    107126      !                                           !: be maximum. Usually of order jpk/2. 
    108127      &     ppacr =   9.00000000000_wp            !: (non dimensional): stretching factor 
    109       !                                           !: for the grid. The highest zacr, the smallest 
     128      !                                           !: for the grid. The higher zacr, the smaller 
    110129      !                                           !: 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    !! 
     130 
     131      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
     132      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
     133 
    115134   REAL(wp), PARAMETER ::                      & 
    116       &     ppdzmin = 6._wp           ,        &  !: (meters) vertical thickness of the top layer 
     135      &     ppdzmin = 6._wp                 ,  &  !: (meters) vertical thickness of the top layer 
    117136      &     pphmax  = 5750._wp                    !: (meters) Maximum depth of the ocean gdepw(jpk) 
    118    !! 
    119137   LOGICAL,  PARAMETER ::                      & 
    120138      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates 
Note: See TracChangeset for help on using the changeset viewer.