1977 - manual for tidal heights analysis and prediction - foreman
TRANSCRIPT
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
1/66
MANUAL FOR TIDAL HEIGHTS
ANALYSIS AND PREDICTION
M .G . G . Fo re m a n
b y
Pacific Marine Science Report 77-10
IN STITUTE O F O C EA N SC IEN C ES, PA TRIC IA BA Y
Sid ne y, B.C .
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
2/66
F o r a d d i t i o n a l c o p i e s o r f u r t h e r i n f o r m a t i o n p l e a s e w r i t e t o :
D e p a r t m e n t o f F i s h e r i e s a n d O c e a n s
I n s t i t u t e o f O c e a n S c i e n c e s
P . O . B o x 6 0 0 0
S i d n e y , B . C . V 8 L 4 B 2
C a n a d a
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
3/66
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
4/66
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
5/66
i
C O N T E N T S
P r e f a c e . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . i i i
A c k n o w l e d g e m e n t s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . i v
1 U S E O F T H E T I D A L H E I G H T S A N A L Y S I S
C O M P U T E R P R O G R A M . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1 . 1 G e n e r a l D e s c r i p t i o n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1 . 2 R o u t i n e s R e q u i r e d . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1 . 3 D a t a I n p u t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1 . 4 O u t p u t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1 . 5 P r o g r a m C o n v e r s i o n , M o d i c a t i o n s , S t o r a g e a n d D i m e n s i o n G u i d e l i n e s . . . . . . 6
2 T I D A L H E I G H T S A N A L Y S I S P R O G R A M D E T A I L S . . . . . . . . . . . . 8
2 . 1 C o n s t i t u e n t D a t a P a c k a g e . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 . 1 . 1 A s t r o n o m i c a l v a r i a b l e s . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 . 1 . 2 C h o i c e o f c o n s t i t u e n t s a n d R a y l e i g h c o m p a r i s o n p a i r s . . . . . . . . . . . 9
2 . 1 . 3 S a t e l l i t e c o n s t i t u e n t s a n d n o d a l m o d u l a t i o n . . . . . . . . . . . . . . 1 5
2 . 1 . 4 S h a l l o w w a t e r c o n s t i t u e n t s . . . . . . . . . . . . . . . . . . . . . . . 1 6
2 . 2 T h e L e a s t S q u a r e s M e t h o d o f A n a l y s i s . . . . . . . . . . . . . . . . . . . . 1 6
2 . 2 . 1 F o r m u l a t i o n o f t h e p r o b l e m . . . . . . . . . . . . . . . . . . . . . . 1 6
2 . 2 . 2 S o l u t i o n o f t h e m a t r i x e q u a t i o n , t h e c o n d i t i o n n u m b e r a n d
s t a t i s t i c a l p r o p e r t i e s . . . . . . . . . . . . . . . . . . . . . . . . . 2 0
2 . 3 M o d i c a t i o n s t o t h e L e a s t S q u a r e s A n a l y s i s R e s u l t s . . . . . . . . . . . . . 2 3
2 . 3 . 1 A s t r o n o m i c a l a r g u m e n t a n d G r e e n w i c h p h a s e l a g . . . . . . . . . . . . 2 3
2 . 3 . 2 N o d a l c o r r e c t i o n s . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 4
2 . 3 . 3 F i n a l a m p l i t u d e a n d p h a s e r e s u l t s . . . . . . . . . . . . . . . . . . . 2 7
2 . 3 . 4 I n f e r r e d c o n s t i t u e n t s . . . . . . . . . . . . . . . . . . . . . . . . . 2 7
3 U S E O F T H E T I D A L H E I G H T S P R E D I C T I O N
C O M P U T E R P R O G R A M . . . . . . . . . . . . . . . . . . . . . . . . . 2 9
3 . 1 G e n e r a l D e s c r i p t i o n . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 9
3 . 2 R o u t i n e s R e q u i r e d . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 9
3 . 3 D a t a I n p u t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 9
3 . 4 O u t p u t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 1
3 . 5 P r o g r a m C o n v e r s i o n , M o d i c a t i o n s , S t o r a g e a n d D i m e n s i o n G u i d e l i n e s . . . . . 3 2
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
6/66
i i
4 T I D A L H E I G H T S P R E D I C T I O N P R O G R A M D E T A I L S . . . . . . . . . . 3 3
4 . 1 P r o b l e m F o r m u l a t i o n a n d t h e E q u a l l y S p a c e d P r e d i c t i o n s M e t h o d . . . . . . . 3 3
4 . 2 T h e H i g h a n d L o w T i d e P r e d i c t i o n M e t h o d . . . . . . . . . . . . . . . . . 3 4
5 C O N S I S T E N C Y O F T H E A N A L Y S I S A N D P R E D I C T I O N
P R O G R A M S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 8
6 R E F E R E N C E S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 9
7 A P P E N D I C E S
7 . 1 S t a n d a r d C o n s t i t u e n t I n p u t D a t a f o r t h e T i d a l H e i g h t s
A n a l y s i s C o m p u t e r P r o g r a m . . . . . . . . . . . . . . . . . . . . . . . . 4 0
7 . 2 S a m p l e T i d a l S t a t i o n I n p u t D a t a f o r t h e A n a l y s i s P r o g r a m . . . . . . . . . . 4 7
7 . 3 F i n a l A n a l y s i s R e s u l t s A r i s i n g f r o m t h e I n p u t D a t a o f
A p p e n d i c e s 7 . 1 a n d 7 . 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 0
7 . 4 S a m p l e I n p u t f o r t h e T i d a l H e i g h t s P r e d i c t i o n P r o g r a m . . . . . . . . . . . . 5 1
7 . 5 T i d a l H e i g h t s P r e d i c t i o n R e s u l t s A r i s i n g f r o m t h e I n p u t D a t a o f
A p p e n d i x 7 . 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 6
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
7/66
i i i
P R E F A C E
T h i s r e p o r t i s i n t e n d e d t o s e r v e a s a u s e r ' s m a n u a l t o G . G o d i n ' s t i d a l h e i g h t s a n a l y s i s
a n d p r e d i c t i o n s p r o g r a m s , r e v i s e d a l o n g l i n e s s u g g e s t e d b y G o d i n . I n a d d i t i o n t o d e s c r i b i n g
i n p u t a n d o u t p u t o f t h e s e p r o g r a m s , t h e r e p o r t g i v e s a n o u t l i n e o f t h e m e t h o d s u s e d a f u l l
p r e s e n t a t i o n o f w h i c h c a n b e f o u n d i n G o d i n ( 1 9 7 2 ) a n d G o d i n a n d T a y l o r ( 1 9 7 3 ) .
U s e r s w h o w i s h t o r e c e i v e u p d a t e s o f t h e s e p r o g r a m s a n d m a n u a l s h o u l d s e n d t h e i r n a m e s ,
a d d r e s s e s , a n d t y p e o f c o m p u t e r u s e d , t o t h e a u t h o r .
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
8/66
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
9/66
1
1 U S E O F T H E T I D A L H E I G H T S A N A L Y S I S
C O M P U T E R P R O G R A M
1 . 1 G e n e r a l D e s c r i p t i o n
T h i s p r o g r a m a n a l y s e s t h e h o u r l y h e i g h t t i d e g a u g e d a t a f o r a g i v e n p e r i o d o f t i m e . A m p l i -
t u d e s a n d G r e e n w i c h p h a s e l a g s a r e c a l c u l a t e d v i a a l e a s t s q u a r e s t m e t h o d c o u p l e d w i t h n o d a l
m o d u l a t i o n f o r o n l y t h o s e c o n s t i t u e n t s t h a t c a n b e r e s o l v e d o v e r t h e l e n g t h o f t h e r e c o r d . U n l e s s
s p e c i e d o t h e r w i s e , a s t a n d a r d d a t a p a c k a g e o f 6 9 c o n s t i t u e n t s w i l l b e c o n s i d e r e d f o r i n c l u s i o n
i n t h e a n a l y s i s . H o w e v e r , u p t o 7 7 a d d i t i o n a l s h a l l o w w a t e r c o n s t i t u e n t s c a n b e r e q u e s t e d . I f
t h e r e c o r d l e n g t h i s s u c h t h a t c e r t a i n i m p o r t a n t c o n s t i t u e n t s a r e n o t i n c l u d e d d i r e c t l y i n t h e
a n a l y s i s , p r o v i s i o n i s m a d e f o r t h e i n f e r e n c e o f t h e a m p l i t u d e a n d p h a s e o f t h e s e c o n s t i t u e n t s
f r o m o t h e r s . G a p s w i t h i n t h e t i d a l r e c o r d a r e p e r m i t t e d .
1 . 2 R o u t i n e s R e q u i r e d
( 1 ) M A I N : : : : : : r e a d s i n s o m e o f d a t a , c o n t r o l s m o s t o f t h e o u t p u t a n d c a l l s o t h e r
r o u t i n e s .
( 2 ) I N P U T : : : : : : r e a d s i n t h e h o u r l y h e i g h t d a t a f o r t h e d e s i r e d t i m e p e r i o d a n d
c h e c k s f o r e r r o r s .
( 3 ) U C O N : : : : : : c h o o s e s t h e c o n s t i t u e n t s t o b e i n c l u d e d i n t h e a n a l y s i s v i a t h e
R a y l e i g h c r i t e r i o n
( 4 ) S C F I T 2 : : : : : : n d s t h e l e a s t s q u a r e s t t o a n e q u a l l y s p a c e d t i m e s e r i e s u s i n g
s i n e s a n d c o s i n e s o f s p e c i e d f r e q u e n c i e s a s t t i n g f u n c t i o n s .
( 5 ) V U F : : : : : : r e a d s r e q u i r e d i n f o r m a t i o n a n d c a l c u l a t e s t h e n o d a l c o r r e c t i o n s f o r
a l l c o n s t i t u e n t s .
( 6 ) I N F E R : : : : : : r e a d s r e q u i r e d i n f o r m a t i o n a n d c a l c u l a t e s t h e a m p l i t u d e a n d p h a s e
o f i n f e r r e d c o n s t i t u e n t s , a s w e l l a s a d j u s t i n g t h e a m p l i t u d e a n d
p h a s e o f t h e c o n s t i t u e n t u s e d f o r t h e i n f e r e n c e .
( 7 ) C H L S K Y : : : : : : s o l v e s t h e s y m m e t r i c p o s i t i v e d e n i t e m a t r i x e q u a t i o n r e s u l t i n g
f r o m a l i n e a r l e a s t s q u a r e s t .
( 8 ) G D A Y : : : : : : r e t u r n s t h e c o n s e c u t i v e d a y n u m b e r f r o m a s p e c i c o r i g i n f o r a n y
g i v e n d a t e a n d v i c e v e r s a .
( 9 ) A S T R : : : : : : c a l c u l a t e s e p h e r m i d e s f o r t h e s u n a n d m o o n .
( 1 0 ) O U T P U T : : : : : : w r i t e s p r e d i c t e d h o u r l y h e i g h t s t o t h e o u t p u t l e .
( 1 1 ) S C U L P : : : : : : s c a l e s u p a m p l i t u d e s t o c o m p e n s a t e f o r m o v i n g a v e r a g e l t e r s .
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
10/66
2
1 . 3 D a t a I n p u t
F o r a c o m p u t e r r u n o f t h e t i d a l h e i g h t s a n a l y s i s p r o g r a m , t w o l o g i c a l u n i t s a r e u s e d f o r
d a t a i n p u t . L o g i c a l u n i t n u m b e r 8 c o n t a i n s t h e t i d a l c o n s t i t u e n t i n f o r m a t i o n w h i l e l o g i c a l u n i t 4
c o n t a i n s t h e h o u r l y h e i g h t s a n d i n f o r m a t i o n r e l a t i n g t o t h e t y p e o f a n a l y s i s a n d o u t p u t r e q u i r e d .
A l i s t i n g o f t h e s t a n d a r d c o n s t i t u e n t i n f o r m a t i o n f o r l o g i c a l u n i t 8 a n d a s a m p l e s e t o f i n p u t
f o r l o g i c a l u n i t 4 a r e g i v e n i n A p p e n d i c e s 7 . 1 a n d 7 . 2 r e s p e c t i v e l y .
L o g i c a l u n i t 8 e x p e c t s f o u r t y p e s o f d a t a :
( i ) O n e c a r d e a c h f o r a l l p o s s i b l e c o n s t i t u e n t s , K O N T A B , t o b e i n c l u d e d i n t h e a n a l y s i s a l o n g w i t h
t h e i r f r e q u e n c i e s , F R E Q , i n c y c l e s / h a n d t h e c o n s t i t u e n t w i t h w h i c h t h e y s h o u l d b e c o m p a r e d
u n d e r t h e R a y l e i g h c r i t e r i o n , K M P R . T h e f o r m a t u s e d i s ( 4 X , A 5 , 3 X , F 1 3 . 1 0 , 4 X , A 5 ) . U n l e s s
K O N T A B i s s p e c i c a l l y d e s i g n a t e d o n l o g i c a l u n i t 4 f o r i n c l u s i o n , a b l a n k d a t a e l d f o r K M P R
r e s u l t s i n t h e c o n s t i t u e n t n o t b e i n g i n c l u d e d i n t h e a n a l y s i s .
A b l a n k c a r d t e r m i n a t e s t h i s d a t a t y p e .
( i i ) T w o c a r d s s p e c i f y i n g v a l u e s f o r t h e a s t r o n o m i c a l a r g u m e n t s S O , H O , P O , E N P O , P P O , D S , D H , D P ,
D N P , D P P i n t h e f o r m a t ( 5 F 1 3 . 1 0 )
S O = m e a n l o n g i t u d e o f t h e m o o n ( c y c l e s ) a t t h e r e f e r e n c e t i m e o r i g i n
H O = m e a n l o n g i t u d e o f t h e s u n ( c y c l e s ) a t t h e r e f e r e n c e t i m e o r i g i n
P O = m e a n l o n g i t u d e o f t h e l u n a r p e r i g e e ( c y c l e s ) a t t h e r e f e r e n c e t i m e o r i g i n
E N P O = n e g a t i v e o f t h e m e a n l o n g i t u d e o f t h e a s c e n d i n g n o d e ( c y c l e s ) a t t h e
r e f e r e n c e t i m e o r i g i n
P P O = m e a n l o n g i t u d e o f t h e s o l a r p e r i g e e ( p e r i h e l i o n ) a t t h e r e f e r e n c e t i m e o r i g i n .
D S , D H , D P , D N P , D P P a r e t h e i r r e s p e c t i v e r a t e s o f c h a n g e o v e r a 3 6 5 - d a y p e r i o d a t t h e r e f e r -
e n c e t i m e o r i g i n .
A l t h o u g h t h e s e a r g u m e n t v a l u e s a r e n o t u s e d b y t h e p r o g r a m t h a t w a s r e v i s e d i n O c t o b e r
1 9 9 2 , i n o r d e r t o m a i n t a i n c o n s i s t e n c y w i t h e a r l i e r p r o g r a m s , t h e y a r e s t i l l r e q u i r e d a s i n p u t .
P o l y n o m i a l a p p r o x i m a t i o n s a r e n o w e m p l o y e d t o m o r e a c c u r a t e l y e v a l u a t e t h e a s t r o n o m i c a l
a r g u m e n t s a n d t h e i r r a t e s o f c h a n g e .
( i i i ) A t l e a s t o n e c a r d f o r a l l t h e m a i n t i d a l c o n s t i t u e n t s s p e c i f y i n g t h e i r D o o d s o n n u m b e r s a n d
p h a s e s h i f t s a l o n g w i t h a s m a n y c a r d s a s a r e n e c e s s a r y f o r t h e s a t e l l i t e c o n s t i t u e n t s . T h e
r s t c a r d f o r e a c h s u c h c o n s t i t u e n t i s i n t h e f o r m a t ( 6 X , A 5 , 1 X , 6 I 3 , F 5 . 2 , I 4 ) a n d c o n t a i n s
t h e f o l l o w i n g i n f o r m a t i o n :
K O N = c o n s t i t u e n t n a m e
I I , J J , K K , L L , M M , N N = t h e s i x D o o d s o n n u m b e r s f o r K O N
S E M I = t h e p h a s e c o r r e c t i o n f o r K O N
N J = t h e n u m b e r o f s a t e l l i t e c o n s t i t u e n t s .
A b l a n k c a r d t e r m i n a t e s t h i s d a t a t y p e .
I f N J > 0 , i n f o r m a t i o n o n t h e s a t e l l i t e c o n s t i t u e n t s f o l l o w s , t h r e e s a t e l l i t e s p e r c a r d , i n
t h e f o r m a t ( 1 1 X , 3 ( 3 I 3 , F 4 . 2 , F 7 . 4 , 1 X , I 1 , 1 X ) ) . F o r e a c h s a t e l l i t e t h e v a l u e s r e a d a r e :
L D E L , M D E L , N D E L = t h e l a s t t h r e e D o o d s o n n u m b e r s o f t h e m a i n c o n s t i t u e n t
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
11/66
3
s u b t r a c t e d f r o m t h e l a s t t h r e e D o o d s o n n u m b e r s o f t h e
s a t e l l i t e c o n s t i t u e n t
P H = p h a s e c o r r e c t i o n o f t h e s a t e l l i t e c o n s t i t u e n t r e l a t i v e t o t h e p h a s e
o f t h e m a i n c o n s t i t u e n t
E E = a m p l i t u d e r a t i o o f t h e s a t e l l i t e t i d a l p o t e n t i a l t o t h a t o f t h e
m a i n c o n s t i t u e n t
I R = 1 i f t h e a m p l i t u d e r a t i o h a s t o b e m u l t i p l i e d b y t h e l a t i t u d e
c o r r e c t i o n f a c t o r f o r d i u r n a l c o n s t i t u e n t s ,
= 2 i f t h e a m p l i t u d e r a t i o h a s t o b e m u l t i p l i e d b y t h e l a t i t u d e
c o r r e c t i o n f a c t o r f o r s e m i d i u r n a l c o n s t i t u e n t s ,
= o t h e r w i s e i f n o c o r r e c t i o n i s r e q u i r e d t o t h e a m p l i t u d e r a t i o .
( i v ) O n e c a r d s p e c i f y i n g e a c h o f t h e s h a l l o w w a t e r c o n s t i t u e n t s a n d t h e m a i n c o n s t i t u e n t s f r o m
w h i c h t h e y a r e d e r i v e d . T h e f o r m a t i s ( 6 X , A 5 , I 1 , 2 X , 4 ( F 5 . 2 , A 5 , 5 X ) ) a n d t h e r e s p e c t i v e
v a l u e s a r e :
K O N = n a m e o f t h e s h a l l o w w a t e r c o n s t i t u e n t
N J = n u m b e r o f m a i n c o n s t i t u e n t s f r o m w h i c h i t i s d e r i v e d
C O E F , K O N C O = c o m b i n a t i o n n u m b e r a n d n a m e o f t h e s e m a i n c o n s t i t u e n t s .
T h e e n d o f t h e s e s h a l l o w w a t e r c o n s t i t u e n t s i s d e n o t e d b y a b l a n k c a r d .
L o g i c a l u n i t 4 c o n t a i n s s i x t y p e s o f d a t a :
( i ) O n e c a r d f o r t h e v a r i a b l e s I O U T 1 , R A Y O P T , Z O F F , I C H K , O B S F A C , I N D P R , N S T R P i n t h e f o r m a t
( I 2 , 2 X , F 4 . 2 , 2 X , F 1 0 . 0 , I 2 , 3 X , F 1 0 . 7 , 2 1 5 )
I O U T 1 = 6 i f t h e o n l y o u t p u t d e s i r e d i s a l i n e p r i n t e r l i s t i n g o f r e s u l t s ,
= 2 i f b o t h a n a l y s i s o u t p u t a n d l i s t i n g a r e d e s i r e d
R A Y O P T = R a y l e i g h c r i t e r i o n c o n s t a n t v a l u e i f d i e r e n t f r o m 1 . 0
Z O F F = c o n s t a n t t o b e s u b t r a c t e d f r o m a l l t h e h o u r l y h e i g h t s
I C H K = 0 i f t h e h o u r l y h e i g h t i n p u t d a t a i s t o b e c h e c k e d f o r f o r m a t e r r o r s ,
= o t h e r w i s e i f t h i s c h e c k i n g t o b e w a i v e d
O B S F A C = s c a l i n g f a c t o r , i f d i e r e n t f r o m 0 . 0 1 , w h i c h w i l l m u l t i p l y t h e h o u r l y
o b s e r v a t i o n s , i n o r d e r t o p r o d u c e t h e d e s i r e d u n i t s f o r t h e n a l c o n s t i t u e n t
a m p l i t u d e s . ( e . g . i f t h e h o u r l y o b s e r v a t i o n s a r e i n m m / s a n d t h e n a l u n i t s
a r e t o b e f t / s e c , t h e n t h i s v a r i a b l e w o u l d b e s e t t o 0 . 0 0 3 2 8 0 8 . )
I N D P R = 1 i f h o u r l y h e i g h t p r e d i c t i o n s b a s e d o n t h e a n a l y s i s r e s u l t s a r e t o b e
c a l c u l a t e d a n d w r i t t e n o n t o d e v i c e n u m b e r 1 0 . I f t h e r e i s i n f e r e n c e , t h i s
p a r a m e t e r v a l u e w i l l a l s o g i v e t h e r m s r e s i d u a l e r r o r a f t e r i n f e r e n c e
a d j u s t m e n t s h a v e b e e n m a d e ,
= 0 i f n o s u c h p r e d i c t i o n s a r e d e s i r e d
N S T R P = n u m b e r o f s u c c e s s i v e m o v i n g a v e r a g e l t e r s t h a t h a v e b e e n a p p l i e d t o t h e
o r i g i n a l d a t a .
I f N S T R P > 0 , t h e n T I M I N T a n d ( L S T R P ( I ) , I = 1 , N S T R P ) w i l l b e r e a d o n a f o l l o w i n g c a r d , i n
t h e f o r m a t ( F 1 0 . 0 , 1 0 1 5 ) , a n d s u i t a b l e a m p l i t u d e c o r r e c t i o n s w i l l b e a p p l i e d t o c o m p e n s a t e f o r
t h e s m o o t h i n g e e c t o f t h e s e l t e r s .
T I M I N T = s a m p l i n g i n t e r v a l , i n m i n u t e s , o f t h e o r i g i n a l u n l t e r e d
r e c o r d
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
12/66
4
( L S T R P ( J ) , J = 1 , N S T R P ) = n u m b e r o f c o n s e c u t i v e o b s e r v a t i o n s u s e d i n c o m p u t i n g e a c h
o f t h e N S T R P m o v i n g a v e r a g e l t e r s .
( i i ) O n e c a r d f o r e a c h p o s s i b l e i n f e r e n c e p a i r . T h e f o r m a t i s ( 2 ( 4 X , A 5 , E 1 6 . 1 0 ) , 2 F 1 0 . 3 ) a n d
t h e r e s p e c t i v e v a l u e s r e a d a r e :
K O N A N & S I G A N = n a m e a n d f r e q u e n c y o f t h e a n a l y s e d c o n s t i t u e n t t o b e u s e d f o r t h e
i n f e r e n c e
K O N I N & S I G I N = n a m e a n d f r e q u e n c y o f t h e i n f e r r e d c o n s t i t u e n t
R = a m p l i t u d e r a t i o o f K O N I N t o K O N A N
Z E T A = G r e e n w i c h p h a s e l a g o f t h e i n f e r r e d c o n s t i t u e n t s u b t r a c t e d f r o m t h e
G r e e n w i c h p h a s e l a g o f t h e a n a l y s e d c o n s t i t u e n t .
T h e s e a r e t e r m i n a t e d b y o n e b l a n k c a r d .
( i i i ) O n e c a r d f o r e a c h s h a l l o w w a t e r c o n s t i t u e n t , o t h e r t h a n t h o s e i n t h e s t a n d a r d 6 9 c o n s t i t u e n t
d a t a p a c k a g e , t o b e c o n s i d e r e d f o r i n c l u s i o n i n t h e a n a l y s i s . T h e R a y l e i g h c o m p a r i s o n
c o n s t i t u e n t i s a l s o r e q u i r e d a n d t h e a d d i t i o n a l s h a l l o w w a t e r c o n s t i t u e n t m u s t b e f o u n d i n
d a t a t y p e ( i ) o f l o g i c a l u n i t 8 , b u t h a v e a b l a n k d a t a e l d w h e r e t h e R a y l e i g h c o m p a r i s o n
c o n s t i t u e n t i s e x p e c t e d . T h e f o r m a t i s ( 6 X , A 5 , 4 X , A 5 ) a n d a b l a n k c a r d i s r e q u i r e d a t t h e
e n d .
( i v ) O n e c a r d i n t h e f o r m a t ( I 1 , 1 X , 1 0 I 2 ) s p e c i f y i n g t h e f o l l o w i n g i n f o r m a t i o n o n t h e p e r i o d
o f t h e a n a l y s i s :
I N D Y = 8 i n d i c a t e s a n a n a l y s i s i s d e s i r e d f o r t h e u p c o m i n g p e r i o d
= 0 i n d i c a t e s n o f u r t h e r a n a l y s e s a r e r e q u i r e d
I H H 1 , I D D 1 , I M M 1 , I Y Y 1 , I C C 1 = h o u r , d a y , m o n t h , y e a r a n d c e n t u r y o f t h e b e g i n n i n g o f
t h e a n a l y s i s ( m e a s u r e d i n t i m e I T Z O N E o f i n p u t d a t a ( v ) )
I H H L , I D D L , I M M L , I Y Y L , I C C L = h o u r , d a y , m o n t h , y e a r a n d c e n t u r y o f t h e e n d o f t h e
a n a l y s i s .
I f I C C 1 o r I C C L a r e z e r o , t h e i r v a l u e i s r e s e t t o 1 9 .
( v ) O n e c a r d i n t h e f o r m a t ( I 1 , 4 X , A 5 , 3 A 6 , A 4 , A 3 , 1 X , 2 I 2 , I 3 , I 2 , 5 X , A 5 ) c o n t a i n i n g t h e f o l -
l o w i n g i n f o r m a t i o n o n t h e t i d a l s t a t i o n :
I N D I C = 1 i f J c a r d o u t p u t i s d e s i r e d ( n o l o n g e r u s e d ) ,
= o t h e r w i s e i f n o t
K S T N = t i d a l s t a t i o n n u m b e r
( N A ( J ) , J = 1 , 4 ) = t i d a l s t a t i o n n a m e ( 2 2 c h a r a c t e r s m a x i m u m l e n g t h )
I T Z O N E = t i m e z o n e o f t h e h o u r l y o b s e r v a t i o n s
L A D , L A M = s t a t i o n l a t i t u d e i n d e g r e e s a n d m i n u t e s
L O D , L O M = s t a t i o n l o n g i t u d e i n d e g r e e s a n d m i n u t e s
I R E F = r e f e r e n c e s t a t i o n n u m b e r .
( v i ) T h e h o u r l y h e i g h t d a t a c a r d s c o n t a i n t h e f o l l o w i n g i n f o r m a t i o n i n t h e f o r m a t ( I 1 , 1 X , I 5 , 7 X ,
3 I 2 , 1 2 A 4 )
K O L I = 1 o r 2 i n d i c a t e s w h e t h e r t h i s s p e c i c c a r d i s t h e r s t o r s e c o n d
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
13/66
5
o n e f o r t h a t d a y ,
= o t h e r w i s e i n d i c a t e s a n o n - d a t a c a r d w h i c h i s i g n o r e d
J S T N = t i d a l s t a t i o n n u m b e r
I D , I M , I Y = d a y , m o n t h a n d y e a r o f t h e h e i g h t s o n t h i s c a r d .
( K A R D ( J ) , J = 1 , 1 2 ) = h o u r l y h e i g h t s i n i n t e g e r f o r m . T h e n a l c o n s t i t u e n t a m p l i t u d e s
u n l e s s a a r e i n u n i t s 1 / 1 0 0 o f t h o s e f o r t h e h o u r l y h e i g h t n o n z e r o
f o r O B S F A C i s r e a d ( s e e ( i ) ) . M i s s i n g v a l u e s s h o u l d b e s p e c i e d a s
a b l a n k e l d o r 9 9 9 9 .
W h e n K O L I = 1 , t h e r s t h o u r l y h e i g h t o n t h e d a t a c a r d i s a s s u m e d t o b e a t 0 1 0 0 h a n d
w h e n K O L I = 2 , i t i s a s s u m e d t o b e a t 1 3 0 0 h . T h e t i m e z o n e o f t h e s e o b s e r v a t i o n s d e t e r m i n e s
t h e n a t u r e o f t h e G r e e n w i c h p h a s e l a g ( s e e S e c t i o n 2 . 3 . 1 ) .
A f t e r t h e i n i t i a l a n a l y s i s o f a c o m p u t e r r u n i s c o m p l e t e d , c o n t r o l r e t u r n s t o i n p u t ( i v ) .
S u c c e s s i v e c a r d s a r e r e a d t h e n u n t i l e i t h e r a 0 o r 8 v a l u e i s f o u n d f o r I N D Y
T h e h o u r l y h e i g h t d a t a c a r d s n e e d n o t b e g i n a n d e n d s o a s t o i n c l u d e e x a c t l y t h e a n a l y s i s
p e r i o d . T h e p r o g r a m i g n o r e s d a t a o u t s i d e t h i s r a n g e . H o w e v e r i f m o r e t h a n o n e a n a l y s i s i s
d e s i r e d f r o m a s i n g l e j o b s u b m i s s i o n a n d h o u r l y h e i g h t d a t a c a r d s d o e x t e n d b e y o n d t h e r s t
a n a l y s i s p e r i o d , c a r e s h o u l d b e t a k e n t o e n s u r e t h a t o n e o f t h e s e c a r d s d o e s n o t h a v e K O L I = 0
o r b l a n k , o t h e r w i s e t h e j o b w i l l b e t e r m i n a t e d . T h i s i s b e c a u s e a l l s u c c e s s i v e c a r d s a f t e r t h e
o n e c o n t a i n i n g t h e l a s t h o u r o f t h e d e s i r e d a n a l y s i s p e r i o d a r e r e a d i n i n p u t ( i v ) f o r m a t .
1 . 4 O u t p u t
T h r e e l o g i c a l u n i t s a r e u s e d f o r t h e o u t p u t o f r e s u l t s f r o m t h e t i d a l h e i g h t s a n a l y s i s p r o g r a m .
D e v i c e n u m b e r 6 i s t h e l i n e p r i n t e r , 2 i s u s e d f o r a n a l y s i s r e s u l t s a n d 1 0 c o n t a i n s h o u r l y
s y n t h e s i z e d v a l u e s b a s e d o n t h e a n a l y s i s r e s u l t s 6 i s r e q u i r e d f o r a l l p r o g r a m r u n s w h e r e a s t h e
u s e o f 2 a n d 1 0 i s c o n t r o l l e d b y t h e i n p u t v a r i a b l e s I O U T 1 a n d I N D P R w h i c h a r e r e a d f r o m
d e v i c e 4 .
R e c o m m e n d a t i o n s f o r t h e u s e o f m o v i n g a v e r a g e l t e r s o n t h e e l e v a t i o n d a t a p r i o r t o s u b m i s -
s i o n f o r a n a l y s i s , a n d t h e s c a l i n g c o m p e n s a t i o n m e t h o d u s e d i n t h e i m p r o v e d a n a l y s i s p r o g r a m
a r e f o u n d i n F o r e m a n ( 1 9 7 8 ) o r G o d i n ( 1 9 7 2 ) .
W h e n I O U T 1 i s 6 , I N D P R i s o t h e r t h a n 1 , a n d t h e r e a r e n o i n f e r r e d c o n s t i t u e n t s , t h e o n l y
o u t p u t i s t w o p a g e s o n t h e l i n e p r i n t e r . T h e r s t o f t h e s e l i s t s t h e c o n s t i t u e n t s i n c l u d e d i n
t h e l e a s t s q u a r e s t , t h e i r f r e q u e n c i e s i n c y c l e s / h ( a l t h o u g h e i g h t d e c i m a l p l a c e s a r e g i v e n ,
d e p e n d i n g o n c o m p u t e r a c c u r a c y , l e s s t h a n t h i s n u m b e r m a y b e s i g n i c a n t ) , t h e C a n d S
c o e c i e n t v a l u e s ( s e e S e c t i o n 2 . 2 . 1 ) m e a s u r e d i n u n i t s O B S F A C t i m e s t h o s e f o r t h e h o u r l y h e i g h t s ,
a n d t h e i r r e s p e c t i v e s t a n d a r d d e v i a t i o n e s t i m a t e s . I t a l s o s p e c i e s t h e n u m b e r o f h o u r l y h e i g h t
o b s e r v a t i o n s ( e x c l u d i n g g a p s ) w i t h i n t h e a n a l y s i s p e r i o d , t h e a v e r a g e a n d s t a n d a r d d e v i a t i o n o f
t h e o r i g i n a l o b s e r v a t i o n s , t h e r o o t m e a n s q u a r e r e s i d u a l e r r o r , a n d t h e m a t r i x c o n d i t i o n n u m b e r .
I n t h e c o l u m n s t i t l e d A L , G L , A , a n d G , t h e s e c o n d p a g e r e s p e c t i v e l y l i s t s t h e a m p l i t u d e s a n d
p h a s e s ( d e g r e e s ) o b t a i n e d f o r e a c h c o n s t i t u e n t f r o m t h e C a n d S c o e c i e n t v a l u e s , a n d t h e
s a m e a m p l i t u d e s a n d p h a s e s a f t e r n o d a l m o d u l a t i o n a n d a s t r o n o m i c a l a r g u m e n t a d j u s t m e n t s .
T h e i n i t i a l a n d n a l h o u r o f t h e a n a l y s i s a r e a l s o s p e c i e d a l o n g w i t h t h e R a y l e i g h c r i t e r i o n
c o n s t a n t ( ` s e p a r a t i o n ' ) , t h e m i d p o i n t o f t h e a n a l y s i s p e r i o d , t h e t o t a l n u m b e r o f p o s s i b l e h o u r l y
o b s e r v a t i o n s i n t h e a n a l y s i s p e r i o d , a n d t h e t o t a l n u m b e r o f p o s s i b l e o b s e r v a t i o n s u s e d i n t h e
a n a l y s i s . T h i s l a s t v a l u e i n c l u d e s g a p s i n t h e r e c o r d a n d i s t h e l a r g e s t o d d n u m b e r l e s s t h a n
o r e q u a l t o t h e t o t a l n u m b e r o f p o s s i b l e h o u r l y o b s e r v a t i o n s ( i f t h e t o t a l n u m b e r o f p o s s i b l e
h o u r l y o b s e r v a t i o n s i s a n e v e n n u m b e r , t h e l a s t h o u r i s i g n o r e d ) . I f t h e r e i s a t l e a s t o n e
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
14/66
6
i n f e r r e d c o n s t i t u e n t , p a g e 2 r e s u l t s a r e r e p e a t e d w i t h t h e i n c l u s i o n o f i n f e r r e d c o n s t i t u e n t s a n d
a p p r o p r i a t e a d j u s t m e n t s t o t h e c o n s t i t u e n t s f r o m w h i c h t h e i n f e r e n c e s w e r e m a d e . A p p e n d i x 7 . 3
l i s t s t h e n a l p a g e o f r e s u l t s o b t a i n e d f r o m t h e i n p u t v a l u e o f A p p e n d i x 7 . 2 .
T h e o n l y e e c t o f c h a n g i n g t h e v a l u e o f I O U T 1 t o 2 ( r e g a r d l e s s o f I N D P R ' s v a l u e ) i s t o
s t o r e o n l e 2 , t h e s a m e i n f o r m a t i o n a s t h e s e c o n d ( a n d t h i r d ) p a g e ( s ) o f t h e l i n e p r i n t e r .
T h e l i s t o f c o n s t i t u e n t n a m e s , a m p l i t u d e s a n d G r e e n w i c h p h a s e l a g s b e g i n s o n l i n e 5 o f t h i s
l e a n d i s i n t h e c o r r e c t f o r m a t f o r i n p u t t o t h e t i d a l h e i g h t s p r e d i c t i o n p r o g r a m , n a m e l y
( 5 X , A 5 , 2 8 X , F 8 . 4 , F 7 . 2 )
W h e n I N D P R e q u a l s 4 , d e v i c e 1 0 w i l l c o n t a i n h o u r l y p r e d i c t i o n s c a l c u l a t e d f r o m t h e a n a l y s i s
r e s u l t s . V a l u e s a r e s p e c i e d o n l y f o r t h e a n a l y s i s p e r i o d , i n c l u d i n g t h o s e i n t e r v a l s w h e r e t h e r e
w e r e g a p s i n t h e o r i g i n a l r e c o r d , a n d a r e i n t h e s a m e m e a s u r e m e n t u n i t s a n d s c a l i n g a s t h e
o r i g i n a l d a t a . T h e f o r m a t u s e d i s t h e s a m e a s f o r i n p u t t y p e ( v i ) o f l o g i c a l u n i t 4 .
1 . 5 P r o g r a m C o n v e r s i o n , M o d i c a t i o n s , S t o r a g e a n d D i m e n s i o n G u i d e l i n e s
T h e s o u r c e p r o g r a m a n d c o n s t i t u e n t d a t a p a c k a g e d e s c r i b e d i n t h i s m a n u a l h a v e b e e n
t e s t e d o n v a r i o u s m a i n f r a m e , P C a n d w o r k s t a t i o n c o m p u t e r s a t t h e I n s t i t u t e o f O c e a n S c i e n c e s ,
P a t r i c i a B a y . A l t h o u g h a s m u c h o f t h e p r o g r a m a s p o s s i b l e w a s w r i t t e n i n b a s i c F O R T R A N ,
s o m e c h a n g e s m a y b e r e q u i r e d b e f o r e t h e p r o g r a m a n d d a t a p a c k a g e c a n b e u s e d o n o t h e r
i n s t a l l a t i o n s . P l e a s e w r i t e o r c a l l t h e a u t h o r i f a n y p r o b l e m s a r e e n c o u n t e r e d .
T h e p r o g r a m i n i t s p r e s e n t f o r m r e q u i r e s a p p r o x i m a t e l y 6 8 , 0 0 0 b y t e s f o r t h e s t o r a g e o f i t s
i n s t r u c t i o n s a n d a r r a y s .
C h a n g i n g t h e n u m b e r o r t y p e o f c o n s t i t u e n t s i n t h e s t a n d a r d d a t a p a c k a g e m a y r e q u i r e s o m e
a l t e r a t i o n s t o t h e a n a l y s i s p r o g r a m . I f c o n s t i t u e n t s a r e a d d e d t o t h e s t a n d a r d d a t a p a c k a g e , t h e
d i m e n s i o n s o f s e v e r a l a r r a y s m a y h a v e t o b e a l t e r e d . R e s t r i c t i o n s o n t h e m i n i m u m d i m e n s i o n
o f s u c h a r r a y s a r e n o w g i v e n .
L e t
M T O T b e t h e t o t a l n u m b e r o f p o s s i b l e c o n s t i t u e n t s c o n t a i n e d i n t h e d a t a p a c k a g e
( p r e s e n t l y 1 4 6 ) ,
M b e t h e n u m b e r o f c o n s t i t u e n t s c o n s i d e r e d f o r i n c l u s i o n i n t h e a n a l y s i s ( p r e s e n t l y
6 9 p l u s t h e n u m b e r o f s h a l l o w w a t e r c o n s t i t u e n t s s p e c i c a l l y d e s i g n a t e d f o r
i n c l u s i o n ) ,
M C O N b e t h e n u m b e r o f m a i n c o n s t i t u e n t s i n t h e s t a n d a r d d a t a p a c k a g e ( p r e s e n t l y 4 5 ) ,
M S A T b e t h e s u m o f t h e t o t a l n u m b e r o f s a t e l l i t e s f o r t h e s e m a i n c o n s t i t u e n t s a n d t h e
n u m b e r o f m a i n c o n s t i t u e n t s w i t h n o s a t e l l i t e s ( p r e s e n t l y 1 6 2 p l u s 8 f o r t h e
v e r s i o n o f t h e c o n s t i t u e n t d a t a p a c k a g e , l i s t e d i n A p p e n d i x 7 . 1 , t h a t c o n t a i n s n o
t h i r d - o r d e r s a t e l l i t e s f o r b o t h N
2
a n d L
2
) ,
M S H A L b e t h e s u m f o r a l l s h a l l o w w a t e r c o n s t i t u e n t s , o f t h e n u m b e r o f m a i n c o n s t i t u e n t s
f r o m w h i c h e a c h i s d e r i v e d ( p r e s e n t l y 2 5 1 ) .
T h e n i n t h e m a i n p r o g r a m , a r r a y s K O N T A B , F R E Q a n d K M P R s h o u l d h a v e m i n i m u m d i m e n s i o n
M T O T + 1 a r r a y s K O N , C , S , S I G , E R C , E R S , A , E P S , K O , A A a n d G D s h o u l d h a v e m i n i m u m d i m e n s i o n
M a r r a y N K O N s h o u l d h a v e d i m e n s i o n a t l e a s t a s l a r g e a s t h e n u m b e r o f e x t r a s h a l l o w w a t e r
c o n s t i t u e n t s s p e c i c a l l y d e s i g n a t e d f o r a n a l y s i s i n c l u s i o n ( i t s p r e s e n t m a x i m u m i s 1 5 ) a n d a r r a y s
Z a n d X P s h o u l d b e l a r g e e n o u g h t o c o n t a i n t h e h o u r l y h e i g h t s ( a n d g a p s ) i n t h e a n a l y s i s p e r i o d
( i t s p r e s e n t m a x i m u m i s 3 7 5 d a y s ) .
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
15/66
7
I n s u b r o u t i n e I N P U T a r r a y Z s h o u l d b e d i m e n s i o n e d t h e s a m e a s i n t h e m a i n p r o g r a m ,
w h i l e K A R D a n d I H T s h o u l d b e d i m e n s i o n e d 1 2 .
I n t h e o t h e r s u b r o u t i n e O U T P U T , Z i s i n a c o m m o n b l o c k a n d s h o u l d b e d i m e n s i o n e d a s
i n t h e m a i n p r o g r a m , X P i s i n t h e a r g u m e n t l i s t a n d n e e d o n l y h a v e d i m e n s i o n 2 , a n d a r r a y s
M O N T H a n d I H T s h o u l d h a v e d i m e n s i o n 1 2 a n d 2 4 r e s p e c t i v e l y .
I n s u b r o u t i n e V U F , a r r a y s V U a n d F s h o u l d h a v e m i n i m u m d i m e n s i o n M T O T a r r a y s K O N
a n d N J s h o u l d h a v e m i n i m u m d i m e n s i o n M T O T + 1 a r r a y s I I , J J , K K , L L , M M , N N a n d S E M I s h o u l d
h a v e m i n i m u m d i m e n s i o n M C O N + 1 a r r a y s E E , L D E L , M D E L , N D E L , I R a n d P H s h o u l d h a v e m i n i m u m
d i m e n s i o n M S A T a n d K O N C O , C O E F s h o u l d h a v e m i n i m u m d i m e n s i o n M S H A L + 4
I n s u b r o u t i n e I N F E R , a r r a y s K O N A N , K O N I N , S I G A N , S I G I N , R a n d Z E T A c a n p r e s e n t l y a c c o m -
m o d a t e a m a x i m u m o f n i n e i n f e r r e d c o n s t i t u e n t s .
I n s u b r o u t i n e S C F I T 2 , a r r a y s X , X P , C , S , E R C , E R S a n d F s h o u l d h a v e t h e s a m e d i m e n s i o n a s
Z , X P , C , S , E R C , E R S a n d S I G i n t h e m a i n p r o g r a m a n d a r r a y s R H S a n d A s h o u l d h a v e m i n i m u m
d i m e n s i o n 2 M - 1 a n d M ( 2 M - 1 ) r e s p e c t i v e l y . A C a n d A S s h o u l d h a v e t h e s i z e o f A a n d c a r e s h o u l d
b e t a k e n t h a t t h r o u g h t h e i r e q u i v a l e n c e r e l a t i o n s h i p s , n e i t h e r A C a n d A S , n o r R H S C a n d R H S S
o v e r l a p .
F i n a l l y , i n s u b r o u t i n e C H L S K Y , a r r a y s A a n d F s h o u l d h a v e m i n i m u m d i m e n s i o n s M ( 2 M - 1 )
a n d 2 M - 1 r e s p e c t i v e l y .
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
16/66
8
2 T I D A L H E I G H T S A N A L Y S I S P R O G R A M D E T A I L S
2 . 1 C o n s t i t u e n t D a t a P a c k a g e
2 . 1 . 1 A s t r o n o m i c a l v a r i a b l e s
T h e a s t r o n o m i c a l v a r i a b l e s r e q u i r e d b y t h e t i d a l a n a l y s i s p r o g r a m w e r e u s e d b y D o o d s o n
( 1 9 2 1 ) i n h i s d e v e l o p m e n t o f t h e t i d a l p o t e n t i a l . F r o m t h e m o n e c a n c a l c u l a t e t h e p o s i t i o n o f
t h e s u n o r m o o n , a n d h e n c e t h e t i d e g e n e r a t i n g f o r c e s , a t a n y t i m e . T h e s e v a r a i b l e s a r e :
S ( t ) = m e a n l o n g i t u d e o f t h e m o o n
H ( t ) = m e a n l o n g i t u d e o f t h e s u n
P ( t ) = m e a n l o n g i t u d e o f t h e l u n a r p e r i g e e
N
0
( t ) = n e g a t i v e o f t h e l o n g i t u d e o f t h e m e a n a s c e n d i n g n o d e
P
0
( t ) = m e a n l o n g i t u d e o f t h e s o l a r p e r i g e e ( p e r i h e l i o n ) .
F o r H , N
0
a n d P
0
t h e s e l o n g i t u d e s a r e m e a s u r e d a l o n g t h e e c l i p t i c e a s t w a r d f r o m t h e
m e a n v e r n a l e q u i n o x p o s i t i o n a t t i m e t w h i l e f o r S a n d P t h e y a r e m e a s u r e d i n t h e e c l i p t i c
e a s t w a r d f r o m t h e m e a n v e r n a l e q u i n o x p o s i t i o n a t t i m e t t o t h e m e a n a s c e n d i n g m o d e o f t h e
l u n a r o r b i t , a n d t h e n a l o n g t h i s o r b i t . T o g e t h e r w i t h t h e r a t e s o f c h a n g e o f t h e s e v a r i a b l e s ,
t h e l o c a l m e a n l u n a r t i m e , a n d t h e D o o d s o n n u m b e r s f o r e a c h t i d a l c o n s t i t u e n t , o n e c a n
c a l c u l a t e t h e c o n s t i t u e n t f r e q u e n c i e s , t h e i r a s t r o n o m i c a l a r g u m e n t p h a s e a n g l e s , V , a n d t h e i r
n o d a l m o d u l a t i o n p h a s e , u , a n d a m p l i t u d e , f , c o r r e c t i o n s .
T h e v a l u e s o f t h e a s t r o n o m i c a l v a r i a b l e s a n d c o n s t i t u e n t f r e q u e n c i e s i n t h e p r o g r a m a r e
c a l c u l a t e d u s i n g t h e p o w e r s e r i e s e x p a n s i o n f o r m u l a e g i v e n o n p a g e s 9 8 a n d 1 0 7 o f t h e E x -
p l a n a t o r y S u p p l e m e n t t o t h e A s t r o n o m i c a l E p h e m e r i s a n d t h e A m e r i c a n E p h e m e r i s a n d N a u t i c a l
A l m a n a c ( 1 9 6 1 ) . T h e s e f o r m u l a e w e r e d e r i v e d f r o m N e w c o m b ' s T a b l e s o f t h e S u n a n d a r e v i s i o n
o f B r o w n ' s l u n a r t h e o r y ( u s e d i n t h e d e v e l o p m e n t o f h i s T a b l e s o f M o t i o n o f t h e M o o n ) s o t h a t
i t i s i n a c c o r d w i t h N e w c o m b ' s .
( F o r t h o s e i n t e r e s t e d , e v e n h i g h e r o r d e r e d a p p r o x i m a t i o n s c a n b e f o u n d i n A s t r o n o m i c a l
F o r m u l a e f o r C a l c u l a t o r s b y J e a n M e e u s . ) I n p a r t i c u l a r , t h e a s t r o n o m i c a l v a r i a b l e s a n d f r e -
q u e n c i e s a r e c a l c u l a t e d a t t h e c e n t r a l h o u r o f t h e a n a l y s i s p e r i o d a n d i n o r d e r t o g a i n p r e c i s i o n
t
0
, t h e r e f e r e n c e t i m e o r i g i n , i s t a k e n t o b e 0 0 0 0 E T .
1
T h i s l a t t e r d a t e , i t w a s f e l t , w o u l d b e
c l o s e r t o t h e a n a l y s i s p e r i o d o f m o s t r e c o r d s t h a n t h e p r e v i o u s c h o i c e o f 0 0 0 0 E T J a n u a r y 1 ,
1 9 0 1 , a n d h e n c e w o u l d y i e l d m o r e a c c u r a t e r e s u l t s v i a t h e l i n e a r a p p r o x i m a t i o n .
I n k e e p i n g w i t h t h e c h o i c e o f r e f e r e n c e t i m e o r i g i n a n d a s t r o n o m i c a l v a r i a b l e s p e c i c a t i o n s , t
s h o u l d b e m e a s u r e d i n E p h e m e r i s t i m e . H o w e v e r , t h e c o r r e c t i o n f r o m U n i v e r s a l t i m e i s i r r e g u l a r
a n d i n m o s t c a s e s s m a l l , s o i t h a s b e e n a s s u m e d f o r c o m p u t a t i o n a l p u r p o s e s t h a t a l l o b s e r v a t i o n s
a r e r e c o r d e d i n E T .
2 . 1 . 2 C h o i c e o f c o n s t i t u e n t s a n d R a y l e i g h c o m p a r i s o n p a i r s
T h e r e i s a m a x i m u m o f 1 4 6 p o s s i b l e t i d a l c o n s t i t u e n t s t h a t c a n b e i n c l u d e d i n t h e t i d a l
a n a l y s i s , 4 5 o f t h e s e a r e a s t r o n o m i c a l i n o r i g i n ( m a i n c o n s t i t u e n t s ) w h i l e t h e r e m a i n i n g 1 0 1 a r e
1
E p h e m e r i s T i m e ( E T ) i s t h e u n i f o r m m e a s u r e o f t i m e d e n e d b y t h e l a w s o f d y n a m i c s a n d d e t e r m i n e d
i n p r i n c i p l e f r o m t h e o r b i t a l m o t i o n o f t h e E a r t h a s r e p r e s e n t e d b y N e w c o m b ' s T a b l e s o f t h e S u n . U n i v e r s a l
o r G r e e n w i c h M e a n T i m e i s d e n e d b y t h e r o t a t i o n a l m o t i o n o f t h e E a r t h a n d i s n o t r i g o r o u s l y u n i f o r m .
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
17/66
9
s h a l l o w w a t e r c o n s t i t u e n t s .
2
B e c a u s e c o m p u t a t i o n t i m e ( a n d c o s t ) o f t h e c o m p u t e r p r o g r a m
i n c r e a s e s a p p r o x i m a t e l y a s t h e s q u a r e o f t h e n u m b e r o f c o n s t i t u e n t s i n c l u d e d i n t h e a n a l y s i s ,
a n d b e c a u s e f o r m a n y t i d a l s t a t i o n s , m o s t o f t h e s h a l l o w w a t e r c o n s t i t u e n t s a r e i n s i g n i c a n t , a
s m a l l e r s t a n d a r d p a c k a g e w a s s e e n a s a d e q u a t e f o r g e n e r a l u s e . B a s e d o n t h e s u g g e s t i o n s o f G .
G o d i n , i t w a s d e c i d e d t h a t t h i s p a c k a g e c o n t a i n a l l t h e m a i n c o n s t i t u e n t s a n d 2 4 o f t h e s h a l l o w
w a t e r . H o w e v e r , p r o v i s i o n w a s m a d e s o t h a t o t h e r s h a l l o w w a t e r c o n s t i t u e n t s a m o n g t h e 7 7
r e m a i n i n g c o u l d b e i n c l u d e d i f d e s i r e d .
T h e R a y l e i g h c o m p a r i s o n c o n s t i t u e n t i s u s e d f o r t h e p u r p o s e o f d e c i d i n g w h e t h e r o r n o t
a s p e c i c c o n s t i t u e n t s h o u l d b e i n c l u d e d i n t h e a n a l y s i s . I f F
0
i s t h e f r e q u e n c y o f s u c h a
c o n s t i t u e n t , F
1
i s t h e f r e q u e n c y o f i t s R a y l e i g h c o m p a r i s o n c o n s t i t u e n t a n d T i s t h e t i m e s p a n
o f t h e p r o p o s e d r e c o r d t o b e a n a l y s e d , t h e n t h e c o n s t i t u e n t w i l l b e i n c l u d e d i n t h e a n a l y s i s
o n l y i f F
0
; F
1
T R A Y R A Y i s c o m m o n l y g i v e n t h e v a l u e 1 a l t h o u g h i t c a n b e s p e c i e d
d i e r e n t l y i n t h e p r o g r a m .
I n o r d e r t o d e t e r m i n e t h e s e t o f R a y l e i g h c o m p a r i s o n p a i r s , i t i s i m p o r t a n t t o c o n s i d e r ,
w i t h i n a g i v e n c o n s t i t u e n t g r o u p ( e . g . d i u r n a l o r s e m i d i u r n a l ) , t h e o r d e r o f c o n s t i t u e n t i n c l u s i o n
i n t h e a n a l y s i s a s T ( t h e t i m e s p a n o f t h e r e c o r d t o b e a n a l y s e d ) i n c r e a s e s . A s s u m i n g t h i s p o i n t
o f v i e w , t h e s p e c i c o b j e c t i v e s u s e d w h e n c o n s t r u c t i n g t h e s e t l i s t e d i n A p p e n d i x 7 . 1 w e r e :
( i ) w i t h i n e a c h c o n s t i t u e n t g r o u p , w h e n p o s s i b l e , h a v e t h e o r d e r o f c o n s t i t u e n t s e l e c t i o n c o r r e -
s p o n d w i t h d e c r e a s i n g m a g n i t u d e o f t i d a l p o t e n t i a l a m p l i t u d e ( a s c a l c u l a t e d b y C a r t w r i g h t
a n d E d d e n ( 1 9 7 3 ) ) ,
( i i ) w h e n p o s s i b l e , c o m p a r e a c a n d i d a t e c o n s t i t u e n t w i t h w h i c h e v e r o f t h e n e i g h b o u r i n g , a l r e a d y
s e l e c t e d c o n s t i t u e n t s , t h a t i s n e a r e s t i n f r e q u e n c y ,
( i i i ) w h e n t h e r e a r e t w o n e i g h b o u r i n g c o n s t i t u e n t s o f r e l a t i v e l y e q u a l t i d a l p o t e n t i a l a m p l i t u d e ,
r a t h e r t h a n w a i t i n g u n t i l t h e r e c o r d l e n g t h i s s u c i e n t t o p e r m i t t h e s e l e c t i o n o f b o t h a t
t h e s a m e t i m e ( i . e . b y c o m p a r i n g t h e m t o e a c h o t h e r ) , c h o o s e a r e p r e s e n t a t i v e o f t h e p a i r
w h o s e i n c l u s i o n w i l l b e a s e a r l y a s p o s s i b l e . T h i s w i l l g i v e i n f o r m a t i o n s o o n e r a b o u t t h a t
f r e q u e n c y r a n g e , a n d v i a i n f e r e n c e , s t i l l e n a b l e s o m e i n f o r m a t i o n t o b e o b t a i n e d o n b o t h
c o n s t i t u e n t s .
T h e R a y l e i g h c o m p a r i s o n p a i r s c h o s e n f o r t h e l o w f r e q u e n c y , d i u r n a l , s e m i d i u r n a l a n d t e r -
d i u r n a l c o n s t i t u e n t g r o u p s a r e g i v e n i n T a b l e s 1 , 2 , 3 a n d 4 r e s p e c t i v e l y . F i g u r e s g i v e n f o r t h e
l e n g t h o f r e c o r d r e q u i r e d f o r c o n s t i t u e n t i n c l u s i o n a s s u m e a R a y l e i g h c r i t e r i o n c o n s t a n t v a l u e
( i n p u t v a r i a b l e R A Y O P T ) o f 1 . 0 .
2 Q
1
a n d S I G
1
p r o v i d e a n e x a m p l e o f o b j e c t i v e ( i i i ) . B e c a u s e 2 Q
1
h a s a g r e a t e r f r e q u e n c y
s e p a r a t i o n f o r Q
1
a n d h e n c e w o u l d a p p e a r i n a n a n a l y s i s o f s h o r t e r r e c o r d l e n g t h t h a n S I G
1
,
i t w a s c h o s e n a s t h e r e p r e s e n t a t i v e .
H o w e v e r , i t c a n b e s e e n i n s e v e r a l c a s e s , t h a t i t w a s n o t p o s s i b l e o r f e a s i b l e t o a d h e r e
t o a l l t h e o b j e c t i v e s j u s t o u t l i n e d . C h o o s i n g a R a y l e i g h c o m p a r i s o n c o n s t i t u e n t f r o m t h e l i s t
o f t h o s e c o n s t i t u e n t s a l r e a d y i n c l u d e d i n t h e a n a l y s i s p r o v e d t o b e d i c u l t n e a r t h e f r e q u e n c y
e d g e s o f c o n s t i t u e n t g r o u p s . U p w a r d a r r o w s i n d i c a t e f a i l u r e t o u p h o l d t h i s o b j e c t i v e . O O
1
i s
s u c h a c a s e . F o r i t , t h e p o t e n t i a l c o m p a r i s o n p a i r s w e r e S O
1
, K
1
a n d J
1
. T h e r s t o f t h e s e
w o u l d r e s u l t i n b o t h S O
1
a n d O O
1
a p p e a r i n g a t t h e s a m e l a t e r t i m e t h a n h a d J
1
o r K
1
b e e n
2
T h e c r i t e r i o n f o r s e l e c t i n g t h e s e m a i n c o n s t i t u e n t s w a s t o i n c l u d e a l l t h e d i u r n a l a n d s e m i d i u r n a l
c o n s t i t u e n t s w i t h C a r t w r i g h t a n d E d d e n ( 1 9 7 3 ) t i d a l p o t e n t i a l a m p l i t u d e s g r e a t e r t h a n 0 . 0 0 2 5 0 , a l o n g w i t h
M
3
a n d t h e m o s t i m p o r t a n t l o w f r e q u e n c y c o n s t i t u e n t s . S e c t i o n 2 . 1 . 3 g i v e s t h e a n a l o g o u s s h a l l o w w a t e r
c o n s t i t u e n t c r i t e r i o n .
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
18/66
1 0
T
a
b
l
e
1
O
r
d
e
r
o
f
S
l
o
w
e
r
-
t
h
a
n
-
D
i
u
r
n
a
l
C
o
n
s
t
i
t
u
e
n
t
S
e
l
e
c
t
i
o
n
i
n
A
c
c
o
r
d
a
n
c
e
w
i
t
h
t
h
e
R
a
y
l
e
i
g
h
C
r
i
t
e
r
i
o
n
T
i
d
a
l
P
o
t
e
n
t
i
a
l
A
m
p
l
i
t
u
d
e
f
o
r
M
a
i
n
C
o
n
s
t
i
t
u
e
n
t
s
S
h
o
w
n
w
i
t
h
i
n
B
r
a
c
k
e
t
s
L
i
n
e
s
w
i
t
h
A
r
r
o
w
s
D
e
n
o
t
e
L
i
n
k
s
w
i
t
h
R
a
y
l
e
i
g
h
C
o
m
p
a
r
i
s
o
n
P
a
i
r
s
L
e
n
g
t
h
o
f
F
r
e
q
u
e
n
c
y
D
i
e
r
e
n
c
e
s
(
c
y
c
l
e
s
/
h
)
1
0
3
b
e
t
w
e
e
n
N
e
i
g
h
b
o
u
r
i
n
g
C
o
n
s
t
i
t
u
e
n
t
s
R
e
c
o
r
d
(
h
)
R
e
q
u
i
r
e
d
f
o
r
C
o
n
s
t
i
t
u
e
n
t
I
n
c
l
u
s
i
o
n
Z
O
S
A
S
S
A
M
S
M
M
M
M
S
F
M
F
1
3
Z
O
3
5
5
(
1
3
6
9
)
M
S
F
7
6
4
(
8
2
5
4
)
M
M
4
3
8
3
(
7
2
8
1
)
(
1
5
6
4
7
)
S
S
A
M
F
4
9
4
2
(
1
5
7
9
)
M
S
M
8
7
6
6
(
1
1
5
6
)
S
A
0.114
07
0.114
08
1.081
62
0.202
37
1.309
78
0.228
16
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
19/66
1 1
T
a
b
l
e
2
O
r
d
e
r
o
f
C
o
n
s
t
i
t
u
e
n
t
S
e
l
e
c
t
i
o
n
i
n
A
c
c
o
r
d
a
n
c
e
w
i
t
h
t
h
e
R
a
y
l
e
i
g
h
C
r
i
t
e
r
i
o
n
T
i
d
a
l
P
o
t
e
n
t
i
a
l
A
m
p
l
i
t
u
d
e
f
o
r
M
a
i
n
C
o
n
s
t
i
t
u
e
n
t
s
i
s
S
h
o
w
n
w
i
t
h
i
n
B
r
a
c
k
e
t
s
L
i
n
e
s
w
i
t
h
A
r
r
o
w
s
D
e
n
o
t
e
L
i
n
k
s
w
i
t
h
R
a
y
l
e
i
g
h
C
o
m
p
a
r
i
s
o
n
P
a
i
r
s
L
e
n
g
t
h
o
f
F
r
e
q
u
e
n
c
y
D
e
r
e
n
c
e
s
(
c
y
c
l
e
s
/
h
)
1
0
3
b
e
t
w
e
e
n
N
e
g
h
b
o
u
r
n
g
C
o
n
s
t
t
u
e
n
t
s
R
e
c
o
r
d
(
h
)
R
e
q
u
r
e
d
f
o
r
C
o
n
s
t
t
u
e
n
t
I
n
c
l
u
s
o
n
A
L
P
1
2
Q
1
S
I
G
1
Q
1
R
H
O
1
O
1
T
A
U
1
B
E
T
1
N
O
1
C
H
I
1
P
I
1
P
1
S
1
K
1
P
S
I
1
P
H
I
1
T
H
E
1
J
1
S
O
1
O
O
1
U
P
S
1
2
4
(
5
3
0
1
1
)
K
1
3
2
8
(
3
7
6
9
4
)
O
1
6
5
1
(
1
6
2
4
)
O
O
1
6
6
2
(
9
5
5
)
(
7
2
1
7
)
(
2
9
6
4
)
(
2
9
6
4
)
(
3
1
1
)
2
Q
1
Q
1
N
O
1
J
1
U
P
S
1
7
6
4
(
2
7
8
)
A
L
P
1
4
3
8
3
(
4
9
3
)
(
2
7
8
)
(
1
7
5
4
3
)
(
7
5
5
)
T
A
U
1
B
E
T
1
P
1
P
H
I
1
S
O
1
4
9
4
2
(
1
1
5
2
)
(
1
3
7
1
)
(
5
6
7
)
(
5
6
7
)
S
I
G
1
R
H
O
1
C
H
I
1
T
H
E
1
8
7
6
7
(
1
0
2
8
)
(
4
1
6
)
(
4
2
2
)
P
I
1
S
1
P
S
I
1
1.309
78
0.20237
1.309
78
0.202
37
1.309
78
0.228
16
1.081
62
0.228
16
0.202
37
0.967
58
0.114
07
0.114
08
0.114
07
0.114
07
0.114
08
1.081
62
0.202
37
1.30
978
0.228
16
1512
15
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
20/66
1 2
T
a
b
l
e
3
O
r
d
e
r
o
f
S
e
m
i
d
i
u
r
n
a
l
C
o
n
s
t
i
t
u
e
n
t
S
e
l
e
c
t
i
o
n
i
n
A
c
c
o
r
d
a
n
c
e
w
i
t
h
t
h
e
R
a
y
l
e
i
g
h
C
r
i
t
e
r
i
o
n
T
i
d
a
l
P
o
t
e
n
t
i
a
l
A
m
p
l
i
t
u
d
e
f
o
r
M
a
i
n
C
o
n
s
t
i
t
u
e
n
t
s
i
s
S
h
o
w
n
w
i
t
h
i
n
B
r
a
c
k
e
t
s
L
i
n
e
s
w
i
t
h
A
r
r
o
w
s
D
e
n
o
t
e
L
i
n
k
s
w
i
t
h
R
a
y
l
e
i
g
h
C
o
m
p
a
r
i
s
o
n
P
a
i
r
s
L
e
n
g
t
h
o
f
F
r
e
q
u
e
n
c
y
D
e
r
e
n
c
e
s
(
c
y
c
l
e
s
/
h
)
1
0
3
b
e
t
w
e
e
n
N
e
g
h
b
o
u
r
n
g
C
o
n
s
t
t
u
e
n
t
s
R
e
c
o
r
d
(
h
)
R
e
q
u
r
e
d
f
o
r
C
o
n
s
t
t
u
e
n
t
I
n
c
l
u
s
o
n
O
Q
2
E
P
S
2
2
N
2
M
U
2
N
2
N
U
2
G
A
M
2
H
1
M
2
H
2
M
K
S
2
L
D
A
2
L
2
T
2
S
2
R
2
K
2
M
S
N
2
E
T
A
2
1
3
(
9
0
8
0
9
)
M
2
3
5
5
(
4
2
2
4
8
)
S
2
6
6
2
(
1
7
3
8
6
)
(
6
4
3
)
N
2
E
T
A
2
7
6
4
(
6
7
1
)
(
2
7
7
6
)
(
2
5
9
7
)
E
P
S
2
M
U
2
L
2
4
3
8
3
(
1
1
4
9
8
)
M
K
S
2
K
2
M
S
N
2
4
9
4
2
(
2
5
9
)
(
2
3
0
1
)
(
3
3
0
2
)
(
6
7
0
)
O
Q
2
2
N
2
N
U
2
L
D
A
2
8
7
6
7
(
3
1
3
)
(
2
7
7
)
(
2
4
7
6
)
(
3
5
5
)
H
1
H
2
T
2
R
2
1
1
3
2
6
(
2
7
3
)
G
A
M
2
0.228
16
1.283
99
0.114
08
0.114
07
0.114
07
1.195
70
0.202
37
1.081
62
0.114
08
0.114
08
0.114
08
0
.0883
0
1.107
41
0.202
37
1.309
78
0.202
37
1.30978
0.202
37
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
21/66
1 3
T
a
b
l
e
4
O
r
d
e
r
o
f
T
e
r
d
i
u
r
n
a
l
C
o
n
s
t
i
t
u
e
n
t
S
e
l
e
c
t
i
o
n
i
n
A
c
c
o
r
d
a
n
c
e
w
i
t
h
t
h
e
R
a
y
l
e
i
g
h
C
r
i
t
e
r
i
o
n
T
i
d
a
l
P
o
t
e
n
t
i
a
l
A
m
p
l
i
t
u
d
e
f
o
r
M
a
i
n
C
o
n
s
t
i
t
u
e
n
t
s
i
s
S
h
o
w
n
w
i
t
h
i
n
B
r
a
c
k
e
t
s
L
i
n
e
s
w
i
t
h
A
r
r
o
w
s
D
e
n
o
t
e
L
i
n
k
s
w
i
t
h
R
a
y
l
e
i
g
h
C
o
m
p
a
r
i
s
o
n
P
a
i
r
s
L
e
n
g
t
h
o
f
F
r
e
q
u
e
n
c
y
D
i
e
r
e
n
c
e
s
(
c
y
c
l
e
s
/
h
)
1
0
3
b
e
t
w
e
e
n
N
e
i
g
h
b
o
u
r
i
n
g
C
o
n
s
t
i
t
u
e
n
t
s
R
e
c
o
r
d
(
h
)
R
e
q
u
i
r
e
d
f
o
r
C
o
n
s
t
i
t
u
e
n
t
I
n
c
l
u
s
i
o
n
M
O
3
M
3
S
O
3
M
K
3
S
K
3
2
5
(
1
1
8
8
)
M
3
3
5
5
S
K
3
6
5
6
M
O
3
M
K
3
4
3
8
3
S
O
3
2.821
93
0.228
16
1.296
89
1.525
05
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
22/66
1 4
c h o s e n . H e n c e , i n f o r m a t i o n a b o u t O O
1
w o u l d b e u n n e c e s s a r i l y d e l a y e d . A l t h o u g h , d u e t o t h e
t i d a l p o t e n t i a l a m p l i t u d e o f J
1
, o b j e c t i v e ( i ) i s v i o l a t e d w i t h b o t h t h e s e c o n d a n d t h i r d c h o i c e s ,
i t w a s f e l t t h a t t h e t h i r d w a s a b e t t e r c o m p r o m i s e . W i t h i t , O O
1
o n l y a p p e a r s 1 1 h s o o n e r
t h a n J
1
K
2
i s a n e x a m p l e o f a n u n a v o i d a b l e v i o l a t i o n o f o b j e c t i v e ( i ) . B e c a u s e i t i s s o c l o s e i n
f r e q u e n c y t o S
2
, i t s i m p o r t a n c e a s a m a j o r s e m i d i u r n a l c o n s t i t u e n t d o e s n o t i n s u r e i t a n e a r l y
i n c l u s i o n i n t h e a n a l y s i s p a c k a g e .
B e c a u s e s h a l l o w w a t e r c o n s t i t u e n t s d o n o t h a v e a t i d a l p o t e n t i a l a m p l i t u d e , o b j e c t i v e ( i )
d o e s n o t a p p l y t o t h e m . H o w e v e r , b a s e d o n h i s e x p e r i e n c e , G o d i n w a s a b l e t o s u g g e s t a
h i e r a r c h y o f t h e i r r e l a t i v e i m p o r t a n c e . A f u r t h e r c r i t e r i a u s e d w h e n s e l e c t i n g c o m p a r i s o n p a i r s
f o r t h e m w a s t h a t n o s h a l l o w w a t e r c o n s t i t u e n t s h o u l d a p p e a r i n a n a n a l y s i s b e f o r e a l l t h e m a i n
c o n s t i t u e n t s , f r o m w h i c h i t i s d e r i v e d , h a v e a l s o b e e n s e l e c t e d . T a b l e 5 s h o w s t h a t t h i s h a s
T a b l e 5 S h a l l o w W a t e r C o n s t i t u e n t s i n t h e S t a n d a r d D a t a P a c k a g e .
S h a l l o w R e c o r d L e n g t h ( h ) C o m p o n e n t M a i n C o n s t i t u e n t s a n d R e c o r d
W a t e r R e q u i r e d f o r L e n g t h s ( h ) R e q u i r e d f o r T h e i r I n c l u s i o n
C o n s t i t u e n t C o n s t i t u e n t I n c l u s i o n i n t h e A n a l y s i s
S O
1
4 3 8 3 S
2
3 5 5 O
1
3 2 8
M K S
2
4 3 8 3 M
2
1 3 K
2
4 3 8 3 S
2
3 5 6
M S N
2
4 3 8 3 M
2
1 3 S
2
3 5 5 N
2
6 6 2
M O
3
6 5 6 M
2
1 3 O
1
3 2 8
S O
3
4 3 8 3 S
2
3 5 5 O
1
3 2 8
M K
3
6 5 6 M
2
1 3 K
1
2 4
S K
3
3 5 5 S
2
3 5 5 K
1
2 4
M N
4
6 6 2 M
2
1 3 N
2
6 6 2
M
4
2 5 M
2
1 3
S N
4
7 6 4 S
2
3 5 5 N
2
6 6 2
M S
4
3 5 5 M
2
1 3 S
2
3 5 5
M K
4
4 3 8 3 M
2
1 3 K
2
4 3 8 3
S
4
3 5 5 S
2
3 5 5
S K
4
4 3 8 3 S
2
3 5 5 K
2
4 3 8 3
2 M K
5
2 4 M
2
1 3 K
1
2 4
2 S K
5
1 7 8 S
2
3 5 5 K
1
2 4
2 M N
6
6 6 2 M
2
1 3 N
2
6 6 2
M
6
2 6 M
2
1 3
2 M S
6
3 5 5 M
2
1 3 S
2
3 5 5
2 M K
6
4 3 8 3 M
2
1 3 K
2
4 3 8 3
2 S M
6
3 5 5 S
2
3 5 5 M
2
1 3
M S K
6
4 3 8 3 M
2
1 3 S
2
3 5 5 K
2
4 3 8 3
3 M K
7
2 4 M
2
1 3 K
1
2 4
M
8
2 6 M
2
1 3
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
23/66
1 5
b e e n u p h e l d f o r a l l s h a l l o w w a t e r c o n s t i t u e n t s i n t h e s t a n d a r d 6 9 c o n s t i t u e n t d a t a p a c k a g e .
W e r e c o m m e n d t h a t t h e o b j e c t i v e s o u t l i n e d h e r e b e e m p l o y e d w h e n c h o o s i n g t h e R a y l e i g h
c o m p a r i s o n c o n s t i t u e n t f o r a n y a d d i t i o n s t o t h e l i s t o f p o s s i b l e c o n s t i t u e n t s t o b e i n c l u d e d i n
t h e a n a l y s i s .
2 . 1 . 3 S a t e l l i t e c o n s t i t u e n t s a n d n o d a l m o d u l a t i o n
D o o d s o n ' s ( 1 9 2 1 ) d e v e l o p m e n t o f t h e t i d a l p o t e n t i a l c o n t a i n s a v e r y l a r g e n u m b e r o f c o n -
s t i t u e n t s . D u e t o t h e g r e a t l e n g t h o f r e c o r d r e q u i r e d f o r t h e i r s e p a r a t i o n , s e v e r a l o f t h e s e c a n b e
c o n s i d e r e d , f o r a l l i n t e n t s a n d p u r p o s e s , u n a n a l y s a b l e . T h e s t a n d a r d a p p r o a c h t o t h i s p r o b l e m
i s t o f o r m c l u s t e r s c o n s i s t i n g o f a l l c o n s t i t u e n t s w i t h t h e s a m e r s t t h r e e D o o d s o n n u m b e r s .
T h e m a j o r c o n t r i b u t o r i n t e r m s o f t i d a l p o t e n t i a l a m p l i t u d e l e n d s i t s n a m e t o t h e c l u s t e r a n d
t h e l e s s e r c o n s t i t u e n t s a r e c a l l e d s a t e l l i t e s .
T h e m e t h o d o f a n a l y s i s u s e s t h i s m a i n a n d s a t e l l i t e c o n s t i t u e n t a p p r o a c h i n t h e f o l l o w i n g
m a n n e r . T h e R a y l e i g h c r i t e r i a i s a p p l i e d t o t h e m a i n c o n s t i t u e n t f r e q u e n c i e s t o d e t e r m i n e
w h e t h e r o r n o t t h e y a r e t o b e i n c l u d e d i n t h e a n a l y s i s . F o r e a c h o f t h o s e s o c h o s e n , w e a n a l y s e
a t i t s f r e q u e n c y a n d o b t a i n a n a p p a r e n t a m p l i t u d e a n d p h a s e . H o w e v e r , b e c a u s e t h e s e r e s u l t s
a r e a c t u a l l y d u e t o t h e c u m u l a t i v e e e c t o f a l l t h e c o n s t i t u e n t s i n t h a t c l u s t e r , a n a d j u s t m e n t
i s m a d e s o t h a t o n l y t h e c o n t r i b u t i o n d u e t o t h e m a i n c o n s t i t u e n t i s f o u n d . T h i s a d j u s t m e n t
i s c a l l e d t h e n o d a l m o d u l a t i o n .
I n o r d e r t o m a k e t h e n o d a l m o d u l a t i o n c o r r e c t i o n t o t h e a m p l i t u d e a n d p h a s e o f a m a i n
c o n s t i t u e n t , i t i s n e c e s s a r y t o k n o w t h e r e l a t i v e a m p l i t u d e s a n d p h a s e s o f t h e s a t e l l i t e s . A s
i s c o m m o n l y d o n e , i t i s a s s u m e d i n t h i s p r o g r a m t h a t t h e s a m e r e l a t i o n s h i p a s i s f o u n d w i t h
t h e e q u i l i b r i u m t i d e ( t i d a l p o t e n t i a l ) , h o l d s w i t h t h e a c t u a l t i d e . T h a t i s , t h e t i d a l p o t e n t i a l
a m p l i t u d e r a t i o o f a s a t e l l i t e t o i t s m a i n c o n s t i t u e n t i s a s s u m e d t o b e e q u a l t o t h e c o r r e s p o n d i n g
t i d a l h e i g h t s a m p l i t u d e r a t i o , a n d t h e d i e r e n c e i n t i d a l p o t e n t i a l p h a s e e q u a l s t h e d i e r e n c e i n
t i d a l h e i g h t p h a s e .
T h e s o u r c e o f t h e t i d a l p o t e n t i a l a m p l i t u d e r a t i o , a s f o u n d i n t h e c o n s t i t u e n t d a t a p a c k a g e
o f A p p e n d i x 7 . 1 , i s C a r t w r i g h t a n d T a y l e r ( 1 9 7 1 ) a n d C a r t w r i g h t a n d E d d e n ( 1 9 7 3 ) . U s i n g
n e w c o m p u t a t i o n m e t h o d s a n d t h e l a t e s t v a l u e s f o r t h e a s t r o n o m i c a l c o n s t a n t s , t h e y o b t a i n e d
m o r e a c c u r a t e r e s u l t s t h a n t h o s e f r o m t h e p r e v i o u s l y u s e d D o o d s o n c o m p u t a t i o n s . I t s h o u l d b e
n o t e d t h a t i n s e v e r a l c a s e s ( w h e n e v e r t h e s a t e l l i t e a r i s e s v i a t h e t h i r d - o r d e r t e r m ) , t h i s v e r s i o n
o f t h e c o n s t i t u e n t d a t a p a c k a g e r e q u i r e s t h a t t h e a m p l i t u d e r a t i o b e m u l t i p l i e d b y a l a t i t u d e
c o r r e c t i o n f a c t o r .
P h a s e d i e r e n c e s b e t w e e n s a t e l l i t e s a n d m a i n c o n s t i t u e n t s a r i s e w h e n t h e t i d a l p o t e n t i a l
d e v e l o p m e n t y i e l d s d i e r e n t t r i g o n o m e t r i c t e r m s f o r t h e s e c o n s t i t u e n t s . T h e c o m m o n c o n v e n t i o n
i s t o e x p r e s s a l l t e r m s i n c o s i n e f o r m a n d s o a n e x t r a ;
1
4
c y c l e p h a s e s h i f t i s i n t r o d u c e d i f t h e
t e r m w a s o r i g i n a l l y a s i n e . S a t e l l i t e s r e q u i r i n g s u c h a s h i f t a r e c a l l e d t h i r d o r d e r . A f u r t h e r
1
2
c y c l e c h a n g e i s a l s o i n t r o d u c e d w h e n a l l n e g a t i v e a m p l i t u d e s a r e m a d e p o s i t i v e .
B e c a u s e s e v e r a l t e s t a n a l y s e s i n d i c a t e l e s s c o n s i s t e n t r e s u l t s w h e n t h i r d - o r d e r s a t e l l i t e s a r e
i n c l u d e d i n t h e N
2
a n d L
2
n o d a l m o d u l a t i o n , G o d i n h a s d e c i d e d t o d e l e t e t h e s e f r o m t h e p r e s e n t
s t a n d a r d c o n s t i t u e n t d a t a p a c k a g e . I n s t e a d h e s u g g e s t s t h a t t h e r e s u l t s o f a n a l y s e s w i t h t h i s
p a c k a g e s h o u l d b e c o m p a r e d w i t h t h o s e o f p r e v i o u s a n a l y s e s i n o r d e r t o n d t h e m o s t s u i t a b l e
a d j u s t m e n t f o r t h e s e c o n s t i t u e n t s .
T h e o n l y o t h e r m a i n c o n s t i t u e n t s t h a t d o n o t h a v e a l l t h e i r s a t e l l i t e s i n c l u d e d f o r n o d a l
m o d u l a t i o n a r e t h e s l o w f r e q u e n c y c o n s t i t u e n t s . F o r t h e m , n o s a t e l l i t e s a r e s p e c i e d . B e c a u s e
l o w f r e q u e n c y n o i s e m a y b e a s m u c h a s a n o r d e r o f m a g n i t u d e g r e a t e r t h a n t h e s a t e l l i t e c o n -
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
24/66
1 6
t r i b u t i o n s , a n d M
m
, M
s f
a n d M
f
w h e n t h e y a r e d e t e c t a b l e a r e o f t e n o f s h a l l o w w a t e r o r i g i n ,
t h e e e c t o f m a k i n g c o r r e c t i o n s f o r t h e e x p e c t e d s a t e l l i t e s w o u l d b e t o o b s c u r e f u r t h e r , r a t h e r
t h a n c l a r i f y t h e a c t u a l l o w f r e q u e n c y p e r i o d i c s i g n a l .
S e c t i o n 2 . 3 . 2 g i v e s f u r t h e r d e t a i l s o n t h e n o d a l m o d u l a t i o n c o r r e c t i o n .
2 . 1 . 4 S h a l l o w w a t e r c o n s t i t u e n t s
S h a l l o w w a t e r t i d a l c o n s t i t u e n t s a r i s e f r o m t h e d i s t o r t i o n o f m a i n c o n s t i t u e n t t i d a l o s c i l l a -
t i o n s i n s h a l l o w w a t e r . B e c a u s e t h e s p e e d o f p r o p a g a t i o n o f a p r o g r e s s i v e w a v e i s a p p r o x i m a t e l y
p r o p o r t i o n a l t o t h e s q u a r e r o o t o f t h e d e p t h o f w a t e r i n w h i c h i t i s t r a v e l l i n g , s h a l l o w w a t e r
h a s t h e e e c t o f r e t a r d i n g t h e t r o u g h o f a w a v e m o r e t h a n t h e c r e s t . T h i s d i s t o r t s t h e o r i g i n a l
s i n u s o i d a l w a v e s h a p e a n d i n t r o d u c e s h a r m o n i c s i g n a l s t h a t a r e n o t p r e d i c t e d i n t i d a l p o t e n t i a l
d e v e l o p m e n t . T h e f r e q u e n c i e s o f t h e s e d e r i v e d h a r m o n i c s c a n b e f o u n d b y c a l c u l a t i n g t h e e e c t
o f n o n - l i n e a r t e r m s i n t h e h y d r o d y n a m i c e q u a t i o n s o f m o t i o n o n a s i g n a l d u e t o o n e o r m o r e
m a i n c o n s t i t u e n t s ( s e e G o d i n ( 1 9 7 2 ) , p p . 1 5 4 { 1 6 4 f o r f u r t h e r d e t a i l s ) .
T h e s h a l l o w w a t e r c o n s t i t u e n t s c h o s e n f o r i n c l u s i o n i n t h e s t a n d a r d 6 9 c o n s t i t u e n t d a t a
p a c k a g e w e r e s u g g e s t e d b y G . G o d i n . T h e y a r e l i s t e d i n T a b l e 5 a n d a r e d e r i v e d o n l y f r o m
t h e l a r g e s t m a i n c o n s t i t u e n t s , n a m e l y M
2
, S
2
, N
2
, K
2
, K
1
a n d O
1
, u s i n g t h e l o w e s t t y p e s o f
p o s s i b l e i n t e r a c t i o n . T h e 7 7 a d d i t i o n a l s h a l l o w w a t e r c o n s t i t u e n t s t h a t c a n b e i n c l u d e d i n t h e
a n a l y s i s i f s o d e s i r e d a r e d e r i v e d f r o m l e s s e r m a i n c o n s t i t u e n t s a n d h i g h e r t y p e s o f i n t e r a c t i o n .
I n t h e c o n s t i t u e n t d a t a p a c k a g e l i s t i n g o f A p p e n d i x 7 . 1 , t h e y c a n b e s p o t t e d b y t h e i r l a c k o f a
R a y l e i g h c o m p a r i s o n c o n s t i t u e n t .
W h e n s h a l l o w w a t e r e e c t s a r e n o t i c e a b l e , m a i n c o n s t i t u e n t s , i f t h e y a r e c l o s e i n f r e q u e n c y ,
m a y c o e x i s t o r b e m a s k e d b y c o n s t i t u e n t s o f n o n - l i n e a r o r i g i n . T h e r e s u l t a n t n o d a l m o d u l a t i o n
w i l l b e d u e t o t h e p a i r a n d t h u s w i l l n o t c o i n c i d e t o t h e c a l c u l a t e d m o d u l a t i o n o f t h e m a i n
c o n s t i t u e n t . I n s u s p e c t e d c a s e s , t h e e e c t i v e n e s s o f n o d a l c o r r e c t i o n s i n a s e r i e s o f s u c c e s s i v e
a n a l y s e s w i l l i n d i c a t e t h e p r e s e n c e o f p a i r s o r e m p h a s i z e t h e p r e d o m i n a n c e o f o n e c o n s t i t u e n t
o v e r t h e o t h e r . T a b l e 6 ( t a k e n f r o m u n p u b l i s h e d n o t e s o f G o d i n ) l i s t s c o m p o u n d c o n s t i t u e n t s
w h i c h m a y c o e x i s t w i t h o r m a s k c o n s t i t u e n t s o f d i r e c t a s t r o n o m i c a l o r i g i n . I n a l l c a s e s e x c e p t
S O
1
a n d M O
3
, t h e m a i n r a t h e r t h a n t h e c o m p o u n d c o n s t i t u e n t i s i n c l u d e d i n t h e s t a n d a r d
c o n s t i t u e n t d a t a p a c k a g e .
2 . 2 T h e L e a s t S q u a r e s M e t h o d o f A n a l y s i s
2 . 2 . 1 F o r m u l a t i o n o f t h e p r o b l e m
T h e r s t s t a g e i n t h e a c t u a l a n a l y s i s o f t i d a l r e c o r d s i s t h e l e a s t s q u a r e s t f o r c o n s t i t u e n t
a m p l i t u d e a n d p h a s e . I f t h e t i d a l r e c o r d i s o f m i n i m u m l e n g t h 1 3 h , t h e p r e s e n t p r o g r a m
a n d d a t a p a c k a g e i n s u r e t h a t t h e c o n s t a n t c o n s t i t u e n t s Z
0
a n d M
2
a r e a l w a y s i n c l u d e d i n t h e
a n a l y s i s . I f
j
f o r j = 1 , M a r e t h e f r e q u e n c i e s ( c y c l e s / h ) o f t h e o t h e r t i d a l c o n s t i t u e n t s
c h o s e n f o r i n c l u s i o n i n t h e a n a l y s i s b y t h e R a y l e i g h c r i t e r i o n , t h e n t h e p r o b l e m i s t o n d t h e
a m p l i t u d e s , A
j
, a n d p h a s e s ,
j
, o f t h e f u n c t i o n C
0
+
P
M
j = 1
A
j
c o s 2 (
j
t
i
;
j
) ] t h a t b e s t t t h e
s e r i e s o f o b s e r v a t i o n s y ( t
i
) , i = 1 N
3
A s s u m i n g N > 2 M + 1 w e s e e t h a t i t i s i m p o s s i b l e t o
3
I n o r d e r t o m i n i m i z e t h e l o s s o f a c c u r a c y d u e t o r o u n d o , t h e a v e r a g e o f t h e h o u r l y h e i g h t s o b s e r v a t i o n s
i s s u b t r a c t e d f r o m a l l o r i g i n a l v a l u e s . T h e y ( t
i
) v a l u e s m e n t i o n e d i n a l l c o m p u t a t i o n s h e n c e f o r t h a r e a c t u a l l y
t h e r e s u l t a n t d e v i a t i o n s . A t t h e e n d o f a l l c a l c u l a t i o n s , C
0
i s a d j u s t e d b y t h i s m e a n v a l u e .
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
25/66
1 7
T a b l e 6 S h a l l o w W a t e r C o n s t i t u e n t s t h a t M a y M a s k M a i n C o n s t i t u e n t s .
M a i n C o n s t i t u e n t C o m p o n e n t C o n s t i t u e n t
w h i c h M a y C o e x i s t a t o r
N e a r i t s F r e q u e n c y
Q
1
N K
1
O
1
N K
1
* *
T A U
1
M P
1
* *
N O
1
* N O
1
* *
P
1
S K
1
* *
K
1
M O
1
J
1
M Q
1
S O
1
S O
1
O Q
2
O Q
2
* *
E P S
2
M N S
2
2 N
2
O
2
* *
M U
2
2 M S
2
N
2
K Q
2
* *
G A M
2
O P
2
* *
M
2
K O
2
* *
L
2
2 M N
2
* *
S
2
K P
2
K
2
K
2
M O
3
M O
3
* *
M
3
N K
3
* *
W i t h M
1
a s a s a t e l l i t e .
T h e m o d u l a t i o n o r f r e q u e n c y o f t h e c o m p o u n d c o n s t i t u e n t
i s s u c i e n t l y d i e r e n t t h a t t h e p a i r c o u l d b e s e p a r a t e d i f a
l o n g e n o u g h r e c o r d o f h i g h p r e c i s i o n w e r e a v a i l a b l e .
s o l v e t h e s y s t e m y ( t
i
) = C
0
+
P
M
j = 1
A
j
c o s 2 (
j
t
i
;
j
) ] e x a c t l y b e c a u s e i t i s o v e r d e t e r m i n e d .
H e n c e , i t i s n e c e s s a r y t o a d o p t a c r i t e r i o n w h i c h w i l l e n a b l e u n i q u e o p t i m u m v a l u e s f o r t h e
p a r a m e t e r s A
j
a n d
j
t o b e f o u n d . T h e m o s t c o m m o n o p t i m i z a t i o n c r i t e r i o n u s e d , a n d t h e o n e
c h o s e n h e r e , i s t h e l e a s t s q u a r e s t e c h n i q u e .
R e - e x p r e s s i n g
P
M
j = 1
A
j
c o s 2 (
j
t
i
;
j
) ] a s
M
X
j = 1
C
j
c o s ( 2
j
t
i
) + S
j
s i n ( 2
j
t
i
)
w h e r e A
j
= ( C
2
j
+ S
2
j
)
1 = 2
a n d 2 j = a r c t a n S
j
= C
j
, s o t h a t t h e t t i n g f u n c t i o n i s l i n e a r i n t h e
p a r a m e t e r s S
j
a n d C
j
a n d h e n c e m o r e e a s i l y s o l v e d , a n d r e w r i t i n g y ( t
i
) a s y
i
, t h e o b j e c t i v e o f
t h e l e a s t s q u a r e s t e c h n i q u e i s t o m i n i m i z e
T =
M
X
i = 1
y
i
; C
0
;
M
X
j = 1
( C
j
c o s 2
j
t
i
+ S
j
s i n 2
j
t
i
)
2
-
8/12/2019 1977 - Manual for Tidal Heights Analysis and Prediction - Foreman
26/66
1 8
F i g u r e 1 T h e m a t r i x e q u a t i o n B x = y r e s u l t i n g f r o m t h e l e a s t s q u a r e s t f o r c o n s t i t u e n t
a m p l i t u d e s a n d p h a s e s .
f o r C
0
a n d a l l C
j
S
j
j = 1 M . T h i s i s d o n e b y s o l v i n g t h e f o l l o w i n g 2 M + 1 s i m u l t a n e o u s
e q u a t i o n s f o r j = 1 M :
0 =
@ T
@ C
0
= 2
N
X
i = 1
y
1
; C
0
;
M
X
j = 1
C
j
c o s 2
j
t
i
;
M
X
j = 1
S
j
s i n 2
j
t
i
( ; 1 )
0 =
@ T
@ C
0
= 2
N
X
i = 1
y
1
; C
0
;
M
X
j = 1
C
j
c o s 2
j
t
i
;
M
X
j = 1
S
j
s i n 2
j
t
i
( ; c o s 2
j
t
i
)
0 =
@ T
@ C
0
= 2
N
X
i = 1
y
1
; C
0
;
M
X
j = 1
C
j
c o s 2
j
t
i
;
M
X
j = 1
S
j
s i n 2
j
t
i
( ; s i n 2
j
t
i
)
T h i s r e s u l t s i n t h e m a t r i x e q u a t i o n B x = y o f F i g u r e 1 .
C
k
=
N
X
i = 1
c o s 2