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_R4.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_R4.h90

    r2715 r4144  
    11   !!--------------------------------------------------------------------- 
    2    !!                     ***  par_ORCA_R4.h90  ***    
     2   !!                     ***  par_ORCA_R4.h90  ***   
    33   !!   Ocean Domain : 4 degrees resolution global ocean 
    44   !!                  (0RCA_R4 configuration) 
     
    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 = 4      ,      &  !: resolution of the configuration (degrees) 
    23       ! Original data size 
    24       jpidta  =  92   ,      &  !: first horizontal dimension > or = to jpi 
    25       jpjdta  =  76   ,      &  !: second                     > or = to jpj 
    26       jpkdta  =  31   ,      &  !: number of levels           > or = to jpk 
    27       ! global 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 
     21      ::                      & 
     22      jp_cfg  =  4     ,      &  !: Resolution of the configuration (degrees) 
     23 
     24      ! Data domain size         !!! *  Size of all input files  * 
     25      jpidta  =  92    ,      &  !: 1st lateral dimension ( >= jpiglo ) 
     26      jpjdta  =  76    ,      &  !: 2nd lateral dimension ( >= jpjglo ) 
     27      jpkdta  =  31    ,      &  !: Number of levels      ( >= jpk    ) 
     28 
     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 
    3333      ! Domain characteristics 
    34       jperio  =   4             !: lateral cond. type (between 0 and 6) 
     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#else 
     42      ! Global domain size       !!! *  Global domain  * 
     43      jpiglo  =  jpidta,      &  !: 1st dimension of global domain --> i 
     44      jpjglo  =  jpjdta,      &  !: 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  =  4               !: Lateral cond. type (between 0 and 6) 
     50#endif 
     51 
    3552 
    3653   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    3754   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    38    REAL(wp), PARAMETER ::   & 
    39       pp_not_used       = 999999_wp , &  !: 
    40       pp_to_be_computed = 0._wp          !: 
    4155 
    42    !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
    43    INTEGER, PARAMETER ::     & ! 
    44       jphgr_msh = 0            !: type of horizontal mesh 
    45       !                        ! = 0 curvilinear coordinate on the sphere 
    46       !                        !     read in coordinate.nc file 
    47       !                        ! = 1 geographical mesh on the sphere 
    48       !                        !     with regular grid-spacing 
    49       !                        ! = 2 f-plane with regular grid-spacing 
    50       !                        ! = 3 beta-plane with regular grid-spacing 
    51       !                        ! = 4 Mercator grid with T/U point at the equator  with 
    52       !                        !     isotropic resolution (e1_deg) 
     56   REAL(wp), PARAMETER ::              & 
     57      pp_not_used       = 999999._wp , &  !: 
     58      pp_to_be_computed = 0._wp           !: 
     59 
     60 
     61   !! Coefficients associated with the horizontal coordinate system 
     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) 
    5373 
    5474      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    5575      !   The mercator grid starts only approximately at gphi0 because 
    5676      !   of the constraint that the equator be a T point. 
    57    REAL(wp) , PARAMETER ::      &  ! 
    58       ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    59       ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    60       !                            ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    61       ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
    62       ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
     77 
     78   REAL(wp), PARAMETER ::     & 
     79      ppglam0  = pp_not_used, &  !: longitude of first raw and column T-point   (jphgr_msh = 1) 
     80      ppgphi0  = pp_not_used, &  !: latitude  of first raw and column T-point   (jphgr_msh = 1) 
     81      !                          !: latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
     82      ppe1_deg = pp_not_used, &  !: zonal      grid-spacing (degrees) 
     83      ppe2_deg = pp_not_used, &  !: meridional grid-spacing (degrees) 
    6384      ! 
    64       ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
    65       ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     85      ppe1_m   = pp_not_used, &  !: zonal      grid-spacing (meters ) 
     86      ppe2_m   = pp_not_used     !: meridional grid-spacing (meters ) 
    6687 
    6788 
    68    !! Vertical grid parameter for domzgr 
    69    !! ================================== 
     89   !!  Coefficients associated with the vertical coordinate system 
    7090 
    71    REAL(wp), PARAMETER  ::       & 
    72       &     ppsur = -4762.96143546300_wp    ,  &  !: ORCA r4, r2 and r05 coefficients 
     91   REAL(wp), PARAMETER  ::                     & 
     92      &     ppsur =   -4762.96143546300_wp  ,  &  !: ORCA r4, r2 and r05 coefficients 
    7393      &     ppa0  =   255.58049070440_wp    ,  &  !: (default coefficients) 
    7494      &     ppa1  =   245.58132232490_wp    ,  &  !: 
    75       &     ppkth =    21.43336197938_wp    ,  &  !: (non dimensional): gives the approximate 
     95      &     ppkth =   21.43336197938_wp     ,  &  !: (non dimensional): gives the approximate 
    7696      !                                           !: layer number above which  stretching will 
    7797      !                                           !: be maximum. Usually of order jpk/2. 
    78       &     ppacr =     3.00000000000_wp          !: (non dimensional): stretching factor 
    79       !                                           !: for the grid. The highest zacr, the smallest 
     98      &     ppacr =   3.00000000000_wp            !: (non dimensional): stretching factor 
     99      !                                           !: for the grid. The higher zacr, the smaller 
    80100      !                                           !: the stretching. 
    81101 
    82    !!  If both ppa0 ppa1 and ppsur are specified to 0, then 
    83    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
     102      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
     103      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    84104 
    85    REAL(wp), PARAMETER ::        & 
     105   REAL(wp), PARAMETER ::                      & 
    86106      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer 
    87107      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk) 
    88    LOGICAL,  PARAMETER ::        & 
     108   LOGICAL,  PARAMETER ::                      & 
    89109      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates 
    90    REAL(wp), PARAMETER ::        & 
     110   REAL(wp), PARAMETER ::                      & 
    91111      &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters 
    92112      &     ppkth2  = pp_not_used           ,  &  !: 
Note: See TracChangeset for help on using the changeset viewer.