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_R1.h90 – NEMO

Ignore:
Timestamp:
2013-10-28T14:50:08+01:00 (11 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_R1.h90

    r2715 r4144  
    1212   !! Use: key_orca_r1=75 to set 75 levels 
    1313   !!---------------------------------------------------------------------- 
    14    CHARACTER (len=16)      & 
     14   CHARACTER (len=16)         & 
    1515#if !defined key_agrif 
    16       , PARAMETER  & 
     16      , PARAMETER             & 
    1717#endif 
    18       ::     
    19       cp_cfg = "orca"           !: name of the configuration 
    20    INTEGER     & 
     18      ::                       
     19      cp_cfg  =  "orca"          !: Name of the configuration 
     20   INTEGER                    & 
    2121#if !defined key_agrif 
    22       , PARAMETER  & 
     22      , PARAMETER             & 
    2323#endif 
    24       :: & 
    25       jp_cfg = 1    ,        &  !: resolution of the configuration (degrees) 
    26       ! Original data size 
    27       jpidta  =  362,        &  !: first horizontal dimension > or = to jpi 
    28       jpjdta  =  292,        &  !: second                     > or = to jpj 
     24      ::                      & 
     25      jp_cfg  =  1     ,      &  !: Resolution of the configuration (degrees) 
     26 
     27      ! Data domain size         !!! *  Size of all input files  * 
     28      jpidta  =  362   ,      &  !: 1st lateral dimension ( >= jpiglo ) 
     29      jpjdta  =  292   ,      &  !: 2nd lateral dimension ( >= jpjglo ) 
    2930#if key_orca_r1==75 
    30       jpkdta  = 75 ,         &  !: number of levels           > or = to jpk 
     31      jpkdta  =  75    ,      &  !: Number of levels      ( >= jpk    ) 
    3132#else 
    32       jpkdta  = 46 ,         &  !: number of levels           > or = to jpk 
     33      jpkdta  =  46    ,      &  !: Number of levels      ( >= jpk    ) 
    3334#endif 
    34       ! total domain matrix size 
    35       jpiglo  = jpidta,      &  !: first  dimension of global domain --> i 
    36       jpjglo  = jpjdta,      &  !: second dimension of global domain --> j 
    37       ! starting position of the zoom 
    38       jpizoom =   1   ,      &  !: left bottom (i,j) indices of the zoom 
    39       jpjzoom =   1   ,      &  !: in data indices 
     35 
     36#if defined key_c1d 
     37      ! Zoom domain size         !!! *  C1D zoom  * 
     38      jpiglo  =  3     ,      &  !: 1st dimension of global domain --> i 
     39      jpjglo  =  3     ,      &  !: 2nd dimension of global domain --> j 
    4040      ! Domain characteristics 
    41       jperio  =   6             !: lateral cond. type (between 0 and 6) 
     41      jperio  =  0               !: Lateral cond. type (between 0 and 6) 
     42 
     43   INTEGER                    & 
     44      ::                      & 
     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#else 
     49      ! Global domain size       !!! *  Global domain  * 
     50      jpiglo  =  jpidta,      &  !: 1st dimension of global domain --> i 
     51      jpjglo  =  jpjdta,      &  !: 2nd dimension of global domain --> j 
     52      ! Starting position of the zoom 
     53      jpizoom =  1     ,      &  !: Left bottom (i,j) indices of the zoom 
     54      jpjzoom =  1     ,      &  !: in data domain indices 
     55      ! Domain characteristics 
     56      jperio  =  6               !: Lateral cond. type (between 0 and 6) 
     57#endif 
     58 
    4259 
    4360   !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    4461   !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    45    REAL,PARAMETER      ::  pp_not_used = 999999_wp , & 
    46       &                    pp_to_be_computed = 0._wp 
    47    !! 
    48    !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 
    49    !! 
    50    INTEGER, PARAMETER ::     & ! 
    51       jphgr_msh = 0            !: type of horizontal mesh 
    52       !                        !  = 0 curvilinear coordinate on the sphere 
    53       !                        !      read in coordinate.nc file 
    54       !                        !  = 1 geographical mesh on the sphere 
    55       !                        !      with regular grid-spacing 
    56       !                        !  = 2 f-plane with regular grid-spacing 
    57       !                        !  = 3 beta-plane with regular grid-spacing 
    58       !                        !  = 4 Mercator grid with T/U point at the equator  with 
    59       !                        !      isotropic resolution (e1_deg) 
     62 
     63   REAL(wp), PARAMETER ::              & 
     64      pp_not_used       = 999999._wp , &  !: 
     65      pp_to_be_computed = 0._wp           !: 
     66 
     67 
     68   !! Coefficients associated with the horizontal coordinate system 
     69 
     70   INTEGER, PARAMETER  ::     & 
     71      jphgr_msh = 0              !: type of horizontal mesh 
     72      !                          !: = 0 curvilinear coordinate on the sphere 
     73      !                          !:     read in coordinate.nc file 
     74      !                          !: = 1 geographical mesh on the sphere 
     75      !                          !:     with regular grid-spacing 
     76      !                          !: = 2 f-plane with regular grid-spacing 
     77      !                          !: = 3 beta-plane with regular grid-spacing 
     78      !                          !: = 4 Mercator grid with T/U point at the equator  with 
     79      !                          !:     isotropic resolution (e1_deg) 
    6080 
    6181      !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 
    6282      !   The mercator grid starts only approximately at gphi0 because 
    6383      !   of the constraint that the equator be a T point. 
    64    REAL(wp), PARAMETER ::       &  ! 
    65       ppglam0  = pp_not_used,   &  !: longitude of first raw and column T-point (jphgr_msh = 1) 
    66       ppgphi0  = pp_not_used,   &  !: latitude  of first raw and column T-point (jphgr_msh = 1) 
    67       !                            !  latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
    68       ppe1_deg = pp_not_used,   &  !: zonal      grid-spacing (degrees) 
    69       ppe2_deg = pp_not_used,   &  !: meridional grid-spacing (degrees) 
     84 
     85   REAL(wp), PARAMETER ::     & 
     86      ppglam0  = pp_not_used, &  !: longitude of first raw and column T-point   (jphgr_msh = 1) 
     87      ppgphi0  = pp_not_used, &  !: latitude  of first raw and column T-point   (jphgr_msh = 1) 
     88      !                          !: latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 
     89      ppe1_deg = pp_not_used, &  !: zonal      grid-spacing (degrees) 
     90      ppe2_deg = pp_not_used, &  !: meridional grid-spacing (degrees) 
    7091      ! 
    71       ppe1_m   = pp_not_used,   &  !: zonal      grid-spacing (meters ) 
    72       ppe2_m   = pp_not_used       !: meridional grid-spacing (meters ) 
     92      ppe1_m   = pp_not_used, &  !: zonal      grid-spacing (meters ) 
     93      ppe2_m   = pp_not_used     !: meridional grid-spacing (meters ) 
     94 
    7395 
    7496   !!  Coefficients associated with the vertical coordinate system 
    7597 
    7698#if key_orca_r1==75 
    77    REAL(wp), PARAMETER  ::       & 
    78       &     ppsur =  -3958.951371276829_wp ,  &  !: ORCA r1 coefficients 
     99   REAL(wp), PARAMETER  ::                     & 
     100      &     ppsur =   -3958.951371276829_wp ,  &  !: ORCA r1 coefficients 
    79101      &     ppa0  =   103.9530096000000_wp  ,  &  !: (75 levels case) 
    80102      &     ppa1  =   2.415951269000000_wp  ,  &  !: 
     
    82104      !                                           !: layer number above which  stretching will 
    83105      !                                           !: be maximum. Usually of order jpk/2. 
    84       &     ppacr =       7.00000000000_wp        !: (non dimensional): stretching factor 
     106      &     ppacr =   7.00000000000_wp            !: (non dimensional): stretching factor 
    85107      !                                           !: for the grid. The higher zacr, the smaller 
    86108      !                                           !: the stretching. 
    87    !! 
    88    !!  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
    89    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    90    !! 
     109 
     110      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
     111      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
     112 
    91113   REAL(wp), PARAMETER ::                      & 
    92114      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer 
    93115      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk) 
    94    !! 
    95116   LOGICAL,  PARAMETER ::                      & 
    96117      &     ldbletanh = .TRUE.                    !: Use/do not use double tanf function for vertical coordinates 
    97118   REAL(wp), PARAMETER ::                      & 
    98       &     ppa2  =   100.7609285000000_wp  ,  &  !: Double tanh function parameters 
    99       &     ppkth2=   48.02989372000000_wp  ,  &  !: 
    100       &     ppacr2=    13.00000000000_wp          !: 
    101       ! 
     119      &     ppa2    = 100.7609285000000_wp  ,  &  !: Double tanh function parameters 
     120      &     ppkth2  = 48.02989372000000_wp  ,  &  !: 
     121      &     ppacr2  = 13.00000000000_wp           !: 
    102122#else 
    103    REAL(wp), PARAMETER  ::       & 
    104       &     ppsur = pp_to_be_computed ,        &  !: Computed in domzgr, set ppdzmin and pphmax below 
    105       &     ppa0  = pp_to_be_computed ,        &  !:    "           " 
    106       &     ppa1  = pp_to_be_computed ,        &  !:    "           " 
    107       ! 
    108       &     ppkth =  23.563_wp        ,        &  !: (non dimensional): gives the approximate 
     123   REAL(wp), PARAMETER  ::                     & 
     124      &     ppsur =   pp_to_be_computed     ,  &  !: Computed in domzgr, set ppdzmin and pphmax below 
     125      &     ppa0  =   pp_to_be_computed     ,  &  !:    "           " 
     126      &     ppa1  =   pp_to_be_computed     ,  &  !:    "           " 
     127      &     ppkth =   23.563_wp             ,  &  !: (non dimensional): gives the approximate 
    109128      !                                           !: layer number above which  stretching will 
    110129      !                                           !: be maximum. Usually of order jpk/2. 
    111130      &     ppacr =   9.00000000000_wp            !: (non dimensional): stretching factor 
    112       !                                           !: for the grid. The highest zacr, the smallest 
     131      !                                           !: for the grid. The higher zacr, the smaller 
    113132      !                                           !: the stretching. 
    114    !! 
    115    !!  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
    116    !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
    117    !! 
     133 
     134      !  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 
     135      !  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 
     136 
    118137   REAL(wp), PARAMETER ::                      & 
    119       &     ppdzmin = 6._wp           ,        &  !: (meters) vertical thickness of the top layer 
     138      &     ppdzmin = 6._wp                 ,  &  !: (meters) vertical thickness of the top layer 
    120139      &     pphmax  = 5750._wp                    !: (meters) Maximum depth of the ocean gdepw(jpk) 
    121    !! 
    122140   LOGICAL,  PARAMETER ::                      & 
    123141      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates 
Note: See TracChangeset for help on using the changeset viewer.