source: trunk/NEMOGCM/NEMO/OPA_SRC/par_depth.h90 @ 32

Last change on this file since 32 was 32, checked in by cholod, 12 years ago

par_TROP12.h90 => par_TROP12.h90 + par_depth.h90 (also for BB12,PERU12)
File size: 6.3 KB
Line 
1   !!----------------------------------------------------------------------
2   !!                     ***  par_depth.h90  ***
3   !!                    Vertical domain Parameter
4   !!----------------------------------------------------------------------
5   !!----------------------------------------------------------------------
6   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
7   !! $Id: par_ORCA_R2.h90 2715 2011-03-30 15:58:35Z rblod $
8   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
9   !!----------------------------------------------------------------------
10
11   INTEGER     &
12#if !defined key_agrif
13      , PARAMETER  &
14#endif
15      :: &
16#if key_depth==75
17      jpkdta  = 75              !: number of levels           > or = to jpk
18#elif key_depth==50
19      jpkdta  = 50              !: number of levels           > or = to jpk
20#elif key_depth==46
21      jpkdta  = 46              !: number of levels           > or = to jpk
22#else
23      par_depth.h90 :  You should not be there !!! Vertical Level not supported (key_depth)
24#endif
25
26   !!  Coefficients associated with the vertical coordinate system
27   !!---------------------------------------------------------------------
28#if key_depth==75
29   !!---------------------------------------------------------------------
30   REAL(wp), PARAMETER  ::       &
31      &     ppsur =  -3958.951371276829_wp  ,  &  !: ORCA r1 coefficients
32      &     ppa0  =   103.9530096000000_wp  ,  &  !: (75 levels case)
33      &     ppa1  =   2.415951269000000_wp  ,  &  !:
34      &     ppkth =   15.35101370000000_wp  ,  &  !: (non dimensional): gives the approximate
35      !                                           !: layer number above which  stretching will
36      !                                           !: be maximum. Usually of order jpk/2.
37      &     ppacr =       7.00000000000_wp        !: (non dimensional): stretching factor
38      !                                           !: for the grid. The higher zacr, the smaller
39      !                                           !: the stretching.
40   !!
41   !!  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then
42   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
43   !!
44   REAL(wp), PARAMETER ::                      &
45      &     ppdzmin = pp_not_used           ,  &  !: (meters) vertical thickness of the top layer
46      &     pphmax  = pp_not_used                 !: (meters) Maximum depth of the ocean gdepw(jpk)
47   !!
48   LOGICAL,  PARAMETER ::                      &
49      &     ldbletanh = .TRUE.                    !: Use/do not use double tanf function for vertical coordinates
50   REAL(wp), PARAMETER ::                      &
51      &     ppa2  =   100.7609285000000_wp  ,  &  !: Double tanh function parameters
52      &     ppkth2=   48.02989372000000_wp  ,  &  !:
53      &     ppacr2=    13.00000000000_wp          !:
54
55   !!---------------------------------------------------------------------
56#elif key_depth==50
57   !!---------------------------------------------------------------------
58   REAL(wp), PARAMETER  ::       &
59      &     ppsur = -8494.48_wp ,   &  !: vertical grid Mercator
60      &     ppa0  =  257.609_wp ,   &  !: (default coefficients refine in surface)
61      &     ppa1  =  256.819_wp ,   &  !:
62      !
63      &     ppkth =    40_wp ,      &  !: (non dimensional): gives the approximate
64      !                                !: layer number above which  stretching will
65      !                                !: be maximum. Usually of order jpk/2.
66      &     ppacr =   10.00000000000_wp           !: (non dimensional): stretching factor
67      !                                           !: for the grid. The highest zacr, the smallest
68      !                                           !: the stretching.
69   !!
70   !!  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then
71   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
72   !!
73   REAL(wp), PARAMETER ::                      &
74      &     ppdzmin = 6._wp           ,        &  !: (meters) vertical thickness of the top layer
75      &     pphmax  = 5750._wp                    !: (meters) Maximum depth of the ocean gdepw(jpk)
76   !!
77   LOGICAL,  PARAMETER ::                      &
78      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates
79   REAL(wp), PARAMETER ::                      &
80      &     ppa2  =   100.760928500000_wp ,  &  ! IDC, 26/6/09
81      &     ppacr2 =   13.000000000000_wp ,  &  ! IDC, 26/6/09
82      &     ppkth2 =   48.029893720000_wp       ! IDC, 26/6/09
83
84   !!---------------------------------------------------------------------
85#elif key_depth==46
86   !!---------------------------------------------------------------------
87   REAL(wp), PARAMETER  ::       &
88      &     ppsur = pp_to_be_computed ,        &  !: Computed in domzgr, set ppdzmin and pphmax below
89      &     ppa0  = pp_to_be_computed ,        &  !:    "           "
90      &     ppa1  = pp_to_be_computed ,        &  !:    "           "
91      !
92      &     ppkth =  23.563_wp        ,        &  !: (non dimensional): gives the approximate
93      !                                           !: layer number above which  stretching will
94      !                                           !: be maximum. Usually of order jpk/2.
95      &     ppacr =   9.00000000000_wp            !: (non dimensional): stretching factor
96      !                                           !: for the grid. The highest zacr, the smallest
97      !                                           !: the stretching.
98   !!
99   !!  If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then
100   !!  they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr
101   !!
102   REAL(wp), PARAMETER ::                      &
103      &     ppdzmin = 6._wp           ,        &  !: (meters) vertical thickness of the top layer
104      &     pphmax  = 5750._wp                    !: (meters) Maximum depth of the ocean gdepw(jpk)
105   !!
106   LOGICAL,  PARAMETER ::                      &
107      &     ldbletanh = .FALSE.                   !: Use/do not use double tanf function for vertical coordinates
108   REAL(wp), PARAMETER ::                      &
109      &     ppa2    = pp_not_used           ,  &  !: Double tanh function parameters
110      &     ppkth2  = pp_not_used           ,  &  !:
111      &     ppacr2  = pp_not_used                 !:
112#endif
113   !!---------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.