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 | !!--------------------------------------------------------------------- |
---|