BIPM_ocean_tides_2007.pdf

IERS Workshop on Conventions: 20-21 September 2007 Ocean and atmospheric tides standards (used for EIGEN gravity field modeling) Richard Biancale (CNES/GRGS) The height of ocean tides is expressed by a sum over n waves : ( ) ( )( )λ?ψθλ?λ?ξ ,)t(cos),(Zt,, nn n n ?=∑ Zn is the amplitude of the wave n, ψ n is the phase nθ is the Doodson argument which is expressed in linear combination of 6 variables : ( ) ( ) ( ) ( ) ( ) sn pnNnpnhnsnn)t( 55555 654321 ?+′?+?+?+?+= τθ These 6 variables with decreasing frequencies represent the fundamental arguments according to the Sun and Moon motions : τ : angle of the mean lunar day (1.03505 d) s : angle of the mean tropic month (27.32158 d) h : angle of the mean tropic year (365.2422 d) p : angle of the mean lunar perigee (8.8473 y) N‘ : angle of the mean lunar node (18.6129 y) ps : angle of the perihelion (20940.28 y) n1 (= 0, 1, 2, 3...) defines the specie (long period, diurnal, semi-diurnal, ter-diurnal...), n2 the group (in general :1 ≤ n2 ≤ 9) and n3 the constituent (1 ≤ n3 ≤ 9). Ocean tides modeling The amplitude (Zn) and phase (ψn) of the different waves of tides represented by cotidal maps can be expanded in spherical harmonic functions of Zncosψn and Znsinψn : Then we have : ( ) ( )∑ += n nnnnnn ZZt θψθψλ?ξ sinsincoscos,, ( ) ( ) ( ) ( ) ( )? θλθλ θλθλ θλθλ θλθλ sinP cosmsinsinmcos bc cosmsinsinmcos bc sinmsincosmcos da sinmsincosmcos da lm n l m nn lm,nlm,n nn lm,nlm,n nn lm,nlm,n nn lm,nlm,n ∑∑∑ ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? + + + + + + + ? ? = 2 2 2 2 Writing : , 2 bc S, 2 bc Sand 2 da C, 2 da C lm,nlm,nlm,n lm,nlm,n lm,n lm,nlm,n lm,n lm,nlm,n lm,n ? = + = + = ? = ?+?+ ( ) ( ) ( ) ( )? ? ? ? ? += += ∑∑ ∑∑ ?λλψ ?λλψ sinsincossin sinsincoscos ,, ,, lm l m lmnlmnnn lm l m lmnlmnnn PmdmcZ PmbmaZ the height of tide becomes : ( ) ( ) ( ) ( )[ ]∑∑∑∑ ? + ±± ±++±+= λχθλχθ?λ?ξ msinSmcosCsinPt,, nlm,nnnlm,n n l m lm Procedure # reading the point grid by lec_grille(_quart)_degre # converting by gpgm2; p