TidalTimeLag

[manuflo97]

%Add a function to calculate the tidal time lag
Q = 35600; % Jupiter quality factor from Lainey (Nature)
Qio = 50; % Io quality factor
omega = 1.76e-4; % Jupiter spin rotation in rad/s
n = 2*pi/152853.5; % Mean motion

TimeLagJup = (atan(1/Q))/(2abs(n-omega)) % on the planet
TimeLagIo = (152853.5/(2pi))*atan(1/Qio) % on the satellite