tides

代码说明：

   T_XTIDE Tidal prediction YOUT=T_XTIDE(STATION) makes a tidal prediction for the current day using the harmonics file XTIDE. Predictions at other times can be generated using [...]=T_XTIDE(STATION,TIM) where TIM is a vector of matlab-format decimal days ( DATENUM). If TIM is a scalar then it is taken as the number of days present over which to make predictions. If TIM>1000 then it is taken to be the starting time in matlab-format. at the matlab time JD (UTC, but see code for details) where tau = lunar time s = mean longitude of the moon h = mean longitude of the sun p = mean longitude of the lunar perigee np = negative of the longitude of the mean ascending node pp = mean longitude of the perihelion (solar perigee)