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_POMME_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_POMME_R025.h90

    r3294 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 = "pomme"          !: POMME025 name of the configuration 
    17    INTEGER     & 
     15      ::                       
     16      cp_cfg  =  "pomme"         !: 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 
     21      ::                      & 
     22      jp_cfg  =  025   ,      &  !: Resolution of the configuration (degrees) 
     23 
    2424      ! ORCA025 global grid size 
    25       jpiglo_ORCA025 = 1442, & 
    26       jpjglo_ORCA025 = 1021, &  ! not used currently 
     25      jpiglo_ORCA025 = 1442,  & 
     26      jpjglo_ORCA025 = 1021,  &  ! not used currently 
    2727      ! POMME "global" domain localisation in the ORCA025 global grid 
    28       jpi_iw    = 1059,      &  
    29       jpi_ie    = 1088,      &   
    30       jpj_js    = 661,       & 
    31       jpj_jn    = 700,       & 
    32       jpidta  = ( jpi_ie - jpi_iw + 1 ), &   !: =30 first horizontal dimension > or = to jpi 
    33       jpjdta  = ( jpj_jn - jpj_js + 1 ), &   !: =40 second                     > or = to jpj 
    34       jpkdta  = 46 ,         &  !: number of levels           > or = to jpk 
    35       ! total domain matrix size 
    36       jpiglo  = jpidta,      &  !: first  dimension of global domain --> i 
    37       jpjglo  = jpjdta,      &  !: second dimension of global domain --> j 
    38       ! starting position of the zoom 
    39       jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    40       jpjzoom =   1   ,      &  !: in data indices 
     28      jpi_iw  =  1059  ,      &  
     29      jpi_ie  =  1088  ,      &   
     30      jpj_js  =  661   ,      & 
     31      jpj_jn  =  700   ,      & 
     32 
     33      ! Data domain size                      !!! *  Size of all input files  * 
     34      jpidta  =  ( jpi_ie - jpi_iw + 1 ), &   !: =30 1st lateral dimension ( >= jpiglo ) 
     35      jpjdta  =  ( jpj_jn - jpj_js + 1 ), &   !: =40 2nd lateral dimension ( >= jpjglo ) 
     36      jpkdta  =  46    ,      &               !: Number of levels          ( >= jpk    ) 
     37 
     38      ! Global domain size       !!! *  Global domain  * 
     39      jpiglo  =  jpidta,      &  !: 1st dimension of global domain --> i 
     40      jpjglo  =  jpjdta,      &  !: 2nd dimension of global domain --> j 
     41      ! Starting position of the zoom 
     42      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom 
     43      jpjzoom =  1     ,      &  !: in data domain indices 
    4144      ! Domain characteristics 
    42       jperio  =    0            !: lateral cond. type (between 0 and 6) 
     45      jperio  =  0               !: Lateral cond. type (between 0 and 6) 
     46 
    4347 
    4448   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    4549   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    46    REAL,PARAMETER      ::  pp_not_used = 999999_wp , & 
    47       &                    pp_to_be_computed = 0._wp 
    48    !! 
    49    !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
    50    !! 
    51    INTEGER, PARAMETER ::     & ! 
    52       jphgr_msh = 0            !: type of horizontal mesh 
    53       !                        !  = 0 curvilinear coordinate on the sphere 
    54       !                        !      read in coordinate.nc file 
    55       !                        !  = 1 geographical mesh on the sphere 
    56       !                        !      with regular grid-spacing 
    57       !                        !  = 2 f-plane with regular grid-spacing 
    58       !                        !  = 3 beta-plane with regular grid-spacing 
    59       !                        !  = 4 Mercator grid with T/U point at the equator  with 
    60       !                        !      isotropic resolution (e1_deg) 
     50 
     51   REAL(wp), PARAMETER ::              & 
     52      pp_not_used       = 999999._wp , &  !: 
     53      pp_to_be_computed = 0._wp           !: 
     54 
     55 
     56   !! Coefficients associated with the horizontal coordinate system 
     57 
     58   INTEGER, PARAMETER  ::     & 
     59      jphgr_msh = 0              !: type of horizontal mesh 
     60      !                          !: = 0 curvilinear coordinate on the sphere 
     61      !                          !:     read in coordinate.nc file 
     62      !                          !: = 1 geographical mesh on the sphere 
     63      !                          !:     with regular grid-spacing 
     64      !                          !: = 2 f-plane with regular grid-spacing 
     65      !                          !: = 3 beta-plane with regular grid-spacing 
     66      !                          !: = 4 Mercator grid with T/U point at the equator  with 
     67      !                          !:     isotropic resolution (e1_deg) 
    6168 
    6269      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    6370      !   The mercator grid starts only approximately at gphi0 because 
    6471      !   of the constraint that the equator be a T point. 
    65    REAL(wp), PARAMETER ::       &  ! 
    66       ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    67       ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    68       !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    69       ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
    70       ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
     72 
     73   REAL(wp), PARAMETER ::     & 
     74      ppglam0  = pp_not_used, &  !: longitude of first raw and column T-point   (jphgr_msh = 1) 
     75      ppgphi0  = pp_not_used, &  !: latitude  of first raw and column T-point   (jphgr_msh = 1) 
     76      !                          !: latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
     77      ppe1_deg = pp_not_used, &  !: zonal      grid-spacing (degrees) 
     78      ppe2_deg = pp_not_used, &  !: meridional grid-spacing (degrees) 
    7179      ! 
    72       ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
    73       ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     80      ppe1_m   = pp_not_used, &  !: zonal      grid-spacing (meters ) 
     81      ppe2_m   = pp_not_used     !: meridional grid-spacing (meters ) 
     82 
    7483 
    7584   !!  Coefficients associated with the vertical coordinate system 
    7685 
    77    REAL(wp), PARAMETER  ::       & 
    78       &     ppsur = pp_to_be_computed ,  &  !: Computed in domzgr, set ppdzmin and pphmax below 
    79       &     ppa0  = pp_to_be_computed ,  &  !:    "           " 
    80       &     ppa1  = pp_to_be_computed ,  &  !:    "           " 
    81       ! 
    82       &     ppkth =  23.563_wp        ,  &  !: (non dimensional): gives the approximate 
    83       !                                     !: layer number above which  stretching will 
    84       !                                     !: be maximum. Usually of order jpk/2. 
    85       &     ppacr =    9.00000000000_wp     !: (non dimensional): stretching factor 
    86       !                                     !: for the grid. The highest zacr, the smallest 
    87       !                                     !: the stretching. 
     86   REAL(wp), PARAMETER ::                      & 
     87      &     ppsur =   pp_to_be_computed     ,  &  !: Computed in domzgr, set ppdzmin and pphmax below 
     88      &     ppa0  =   pp_to_be_computed     ,  &  !:    "           " 
     89      &     ppa1  =   pp_to_be_computed     ,  &  !:    "           " 
     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 =   9.00000000000_wp            !: (non dimensional): stretching factor 
     94      !                                           !: for the grid. The higher zacr, the smaller 
     95      !                                           !: the stretching. 
    8896 
    89    !! 
    90    !!  If both ppa0 ppa1 and ppsur are specified to 0, then 
    91    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    92    !! 
    93    REAL(wp), PARAMETER ::        & 
    94       &     ppdzmin = 6._wp           ,  &  !: (meters) vertical thickness of the top layer 
    95       &     pphmax  = 5750._wp              !: (meters) Maximum depth of the ocean gdepw(jpk) 
    96    LOGICAL,  PARAMETER ::        & 
    97       &     ldbletanh = .FALSE.             !: Use/do not use double tanf function for vertical coordinates 
    98    REAL(wp), PARAMETER ::        & 
    99       &     ppa2    = pp_not_used     ,  &  !: Double tanh function parameters 
    100       &     ppkth2  = pp_not_used     ,  &  !: 
    101       &     ppacr2  = pp_not_used           !: 
     97      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
     98      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
     99 
     100   REAL(wp), PARAMETER ::                      & 
     101      &     ppdzmin = 6._wp                 ,  &  !: (meters) vertical thickness of the top layer 
     102      &     pphmax  = 5750._wp                    !: (meters) Maximum depth of the ocean gdepw(jpk) 
     103   LOGICAL,  PARAMETER ::                      & 
     104      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates 
     105   REAL(wp), PARAMETER ::                      & 
     106      &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters 
     107      &     ppkth2  = pp_not_used           ,  &  !: 
     108      &     ppacr2  = pp_not_used                 !: 
    102109   !!--------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.