Newer
Older
sparonuz
committed
!!=====================================================================
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
!! *** Include file tide.h90 ***
!!======================================================================
!! History : 3.2 ! 2007 (O. Le Galloudec) Original code
!! ! 2019 (S. Mueller, N. Bruneau) Alternative parameter set
!!----------------------------------------------------------------------
!!
!! ** Purpose : Inclusion of alternative variants of tidal-constituent
!! parameter definitions during code preprocessing: the default
!! variant includes the 34 constituents available in the FES2014
!! version of the Finite Element Solution - Global Tide data
!! product
!! (https://www.aviso.altimetry.fr/en/data/products/auxiliary-products/global-tide-fes.html);
!! also available is the default parameter set available in
!! previous NEMO versions
!!
!! ** References :
!! S58) Schureman, P. (1958): Manual of Harmonic Analysis and
!! Prediction of Tides (Revised (1940) Edition (Reprinted 1958
!! with corrections). Reprinted June 2001). U.S. Department of
!! Commerce, Coast and Geodetic Survey Special Publication
!! No. 98. Washington DC, United States Government Printing
!! Office. 317 pp. DOI: 10.25607/OBP-155.
!! CT71) Cartwright, D. E. and Tayler, R. J. (1971): New computations of
!! the Tide-generating Potential. Geophys. J. R. astr. Soc. 23,
!! pp. 45-74. DOI: 10.1111/j.1365-246X.1971.tb01803.x
!! CE73) Cartwright, D. E. and Edden, A. C. (1973): Corrected Tables of
!! Tidal Harmonics. Geophys. J. R. astr. Soc. 33,
!! pp. 253-264. DOI: 10.1111/j.1365-246X.1973.tb03420.x
!! FES2014) FES (Finite element Solution) - Global
!! tide. https://www.aviso.altimetry.fr/en/data/products/auxiliary-products/global-tide-fes.html
!!
!!----------------------------------------------------------------------
!!----------------------------------------------------------------------
!! NEMO/OCE 4.0 , NEMO Consortium (2019)
!! $Id: tide.h90 14502 2021-02-18 18:48:54Z smueller $
!! Software governed by the CeCILL license (see ./LICENSE)
!!----------------------------------------------------------------------
#ifndef TIDE_VAR_0
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
! | Name | Equilibrium | nt | ns | nh | np | np1 | Phase | nxi | nnu0 | nnu1 | nnu2 | R | Nodal | Equilibrium | Parameters source | Notes |
! | | tide | | | | | | shift | | | | | | correction | tide | | |
! | | | | | | | | | | | | | | formula | source/comment | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
! | Long-period tidal constituents |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 1) = tide( 'Mf' , 0.042054_wp , 0 , 2 , 1 , 0 , 0 , 0 , -2 , 0 , 0 , 0 , 0 , 74 ) ! CE73 | S54 (Table 2, A6) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 2) = tide( 'Mm' , 0.022187_wp , 0 , 1 , 0 , -1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 73 ) ! CE73 | S54 (Table 2, A2) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 3) = tide( 'Ssa' , 0.019572_wp , 0 , 0 , 2 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ) ! CE73 | S54 (Table 2, B6) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 4) = tide( 'Mtm' , 0.008052_wp , 0 , 3 , 0 , -1 , 0 , 0 , -2 , 0 , 0 , 0 , 0 , 74 ) ! CE73 | FES2014 (prediction | |
! | | | | | | | | | | | | | | | | algorithm); S54 | |
! | | | | | | | | | | | | | | | | (Table 2, A7) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 5) = tide( 'Msf' , 0.003677_wp , 0 , 2 , -2 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 73 ) ! CE73 | S54 (Table 2, A5) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 6) = tide( 'Msqm' , 0.001287_wp , 0 , 4 , -2 , 0 , 0 , 0 , -2 , 0 , 0 , 0 , 0 , 74 ) ! CE73 | FES2014 (prediction | |
! | | | | | | | | | | | | | | | | algorithm); S54 | |
! | | | | | | | | | | | | | | | | (Table 2, A12) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 7) = tide( 'Sa' , 0.000000_wp , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ) ! Meteorological | S54 (Table 2, B64) | |
! | | | | | | | | | | | | | | | tide only | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
! | Diurnal tidal constituents |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 8) = tide( 'K1' , 0.142486_wp , 1 , 0 , 1 , 0 , 0 , -90 , 0 , 0 , -1 , 0 , 0 , 227 ) ! CE73, sign | S54 (Table 2) | Note 1 |
! | | | | | | | | | | | | | | | change | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components( 9) = tide( 'O1' , 0.101316_wp , 1 , -2 , 1 , 0 , 0 , 90 , 2 , -1 , 0 , 0 , 0 , 75 ) ! CE73 | S54 (Table 2, A14) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(10) = tide( 'P1' , 0.047152_wp , 1 , 0 , -1 , 0 , 0 , 90 , 0 , 0 , 0 , 0 , 0 , 0 ) ! CE73 | S54 (Table 2, B14) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(11) = tide( 'Q1' , 0.019396_wp , 1 , -3 , 1 , 1 , 0 , 90 , 2 , -1 , 0 , 0 , 0 , 75 ) ! CE73 | S54 (Table 2, A15) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(12) = tide( 'J1' , 0.007967_wp , 1 , 1 , 1 , -1 , 0 , -90 , 0 , -1 , 0 , 0 , 0 , 76 ) ! CE73, sign | S54 (Table 2, A24) | Note 1 |
! | | | | | | | | | | | | | | | change | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(13) = tide( 'S1' , 0.000000_wp , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ) ! Meteorological | S54 (Table 2, B71) | |
! | | | | | | | | | | | | | | | tide only | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
! | Semidiurnal tidal constituents |
! +--------+-------------+-----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(14) = tide( 'M2' , 0.244081_wp , 2 , -2 , 2 , 0 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 ) ! CE73 | S54 (Table 2, A39) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(15) = tide( 'S2' , 0.110242_wp , 2 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ) ! CE73 | S54 (Table 2, B39) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(16) = tide( 'N2' , 0.046732_wp , 2 , -3 , 2 , 1 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 ) ! CE73 | S54 (Table 2, A40) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(17) = tide( 'K2' , 0.030905_wp , 2 , 0 , 2 , 0 , 0 , 0 , 0 , 0 , 0 , -2 , 0 , 235 ) ! CE73 | S54 (Table 2) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(18) = tide( 'nu2' , 0.008877_wp , 2 , -3 , 4 , -1 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 ) ! CE73 | S54 (Table 2, A43) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(19) = tide( 'mu2' , 0.007463_wp , 2 , -4 , 4 , 0 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 ) ! CE73 | S54 (Table 2, A45) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(20) = tide( '2N2' , 0.006184_wp , 2 , -4 , 2 , 2 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 ) ! CE73 | S54 (Table 2, A42) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(21) = tide( 'L2' , 0.006899_wp , 2 , -1 , 2 , -1 , 0 , 180 , 2 , -2 , 0 , 0 , -1 , 215 ) ! CE73, sign | S54 (Table 2) | Note 1 |
! | | | | | | | | | | | | | | | change | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(22) = tide( 'T2' , 0.006655_wp , 2 , 0 , -1 , 0 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ) ! CE73 | S54 (Table 2, B40) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(23) = tide( 'eps2' , 0.001804_wp , 2 , -5 , 4 , 1 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 ) ! CE73 | FES2014 (prediction | |
! | | | | | | | | | | | | | | | | algorithm) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(24) = tide( 'lam2' , 0.001800_wp , 2 , -1 , 0 , 1 , 0 , 180 , 2 , -2 , 0 , 0 , 0 , 78 ) ! CE73, sign | S54 (Table 2, A44) | Note 1 |
! | | | | | | | | | | | | | | | change | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(25) = tide( 'R2' , 0.000952_wp , 2 , 0 , 1 , 0 , -1 , 180 , 0 , 0 , 0 , 0 , 0 , 0 ) ! CE73, sign | S54 (Table 2, B41) | Note 1 |
! | | | | | | | | | | | | | | | change | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
! | Terdiurnal tidal constituents |
! +--------+-------------+-----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(26) = tide( 'M3' , 0.003192_wp , 3 , -3 , 3 , 0 , 0 , 0 , 3 , -3 , 0 , 0 , 0 , 149 ) ! CT71, sign | S54 (Table 2, A82) | Note 2 |
! | | | | | | | | | | | | | | | change | | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
! | Compound tides |
! +--------+-------------+-----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(27) = tide( 'MKS2' , 0.000000_wp , 2 , -2 , 4 , 0 , 0 , 0 , 2 , -2 , 0 , -2 , 0 , 4 ) ! Compound tide | FES2014 (prediction | |
! | | | | | | | | | | | | | | | | algorithm) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(28) = tide( 'MN4' , 0.000000_wp , 4 , -5 , 4 , 1 , 0 , 0 , 4 , -4 , 0 , 0 , 0 , 1 ) ! Compound tide | S54 (Table 2a) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(29) = tide( 'MS4' , 0.000000_wp , 4 , -2 , 2 , 0 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 ) ! Compound tide | FES2014 (prediction | Note 3 |
! | | | | | | | | | | | | | | | | algorithm); S54 | |
! | | | | | | | | | | | | | | | | (Table 2a) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
! | Overtides |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(30) = tide( 'M4' , 0.000000_wp , 4 , -4 , 4 , 0 , 0 , 0 , 4 , -4 , 0 , 0 , 0 , 1 ) ! Overtide | S54 | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(31) = tide( 'N4' , 0.000000_wp , 4 , -6 , 4 , 2 , 0 , 0 , 4 , -4 , 0 , 0 , 0 , 1 ) ! Overtide | FES2014 (prediction | |
! | | | | | | | | | | | | | | | | algorithm) | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(32) = tide( 'S4' , 0.000000_wp , 4 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ) ! Overtide | S54 | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(33) = tide( 'M6' , 0.000000_wp , 6 , -6 , 6 , 0 , 0 , 0 , 6 , -6 , 0 , 0 , 0 , 18 ) ! Overtide | S54 | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
tide_components(34) = tide( 'M8' , 0.000000_wp , 8 , -8 , 8 , 0 , 0 , 0 , 8 , -8 , 0 , 0 , 0 , 20 ) ! Overtide | S54 | |
! +--------+--------------+----+----+----+----+-----+-------+-----+------+------+------+----+------------+----------------+-----------------------+--------+
! Note 1: the negative sign of the equilibrium-tide value derived from CE73 has been changed to accomodate the phase shift from Table 2 of S54.
! Note 2: the negative sign of the equilibrium-tide value derived from CT71 has been changed to accomodate the phase shift from Table 2 of S54.
! Note 3: the nodal correction factor formulas from FES2014 and S54 differ; here, the version from FES2014 has been selected.
#else
! !! name_tide , equitide , nt , ns , nh , np , np1 , shift , nksi , nnu0 , nnu1 , nnu2 , R , formula !!
! !! ! ! ! ! ! ! ! ! ! ! ! ! ! !!
tide_components( 1) = tide( 'M2' , 0.242297 , 2 , -2 , 2 , 0 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 )
tide_components( 2) = tide( 'N2' , 0.046313 , 2 , -3 , 2 , 1 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 )
tide_components( 3) = tide( '2N2' , 0.006184 , 2 , -4 , 2 , 2 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 )
tide_components( 4) = tide( 'S2' , 0.113572 , 2 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 )
tide_components( 5) = tide( 'K2' , 0.030875 , 2 , 0 , 2 , 0 , 0 , 0 , 0 , 0 , 0 , -2 , 0 , 235 )
! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !
tide_components( 6) = tide( 'K1' , 0.142408 , 1 , 0 , 1 , 0 , 0 , -90 , 0 , 0 , -1 , 0 , 0 , 227 )
tide_components( 7) = tide( 'O1' , 0.101266 , 1 , -2 , 1 , 0 , 0 , +90 , 2 , -1 , 0 , 0 , 0 , 75 )
tide_components( 8) = tide( 'Q1' , 0.019387 , 1 , -3 , 1 , 1 , 0 , +90 , 2 , -1 , 0 , 0 , 0 , 75 )
tide_components( 9) = tide( 'P1' , 0.047129 , 1 , 0 , -1 , 0 , 0 , +90 , 0 , 0 , 0 , 0 , 0 , 0 )
! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !
tide_components(10) = tide( 'M4' , 0.000000 , 4 , -4 , 4 , 0 , 0 , 0 , 4 , -4 , 0 , 0 , 0 , 1 )
! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !
tide_components(11) = tide( 'Mf' , 0.042017 , 0 , 2 , 0 , 0 , 0 , 0 , -2 , 0 , 0 , 0 , 0 , 74 )
tide_components(12) = tide( 'Mm' , 0.022191 , 0 , 1 , 0 , -1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 73 )
tide_components(13) = tide( 'Msqm' , 0.000667 , 0 , 4 , -2 , 0 , 0 , 0 , -2 , 0 , 0 , 0 , 0 , 74 )
tide_components(14) = tide( 'Mtm' , 0.008049 , 0 , 3 , 0 , -1 , 0 , 0 , -2 , 0 , 0 , 0 , 0 , 74 )
! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !
tide_components(15) = tide( 'S1' , 0.000000 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 )
tide_components(16) = tide( 'MU2' , 0.005841 , 2 , -4 , 4 , 0 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 )
tide_components(17) = tide( 'NU2' , 0.009094 , 2 , -3 , 4 , -1 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 78 )
tide_components(18) = tide( 'L2' , 0.006694 , 2 , -1 , 2 , -1 , 0 , +180 , 2 , -2 , 0 , 0 , 0 , 215 )
tide_components(19) = tide( 'T2' , 0.006614 , 2 , 0 , -1 , 0 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 )
#endif