Sabang, Aceh
Semidiurnal- Period
- 01 Jan 2026 — 01 Agu 2026212 days
- Datum
- Nol sensor stasiun (bukan MSL)
- gaps
- 6479 jam
- Source
- IOC Sea Level Station Monitoring Facility (UNESCO/IOC & VLIZ)
Source, licence and datum notes
IOC menyajikan data mentah terhadap nol sensor masing-masing stasiun. Tinggi di sini tidak merujuk MSL, LAT, maupun chart datum.
Licence: Akses terbuka dengan kewajiban sitasi (VLIZ/IOC, DOI 10.14284/482)
Flanders Marine Institute (VLIZ); Intergovernmental Oceanographic Commission (IOC) (2026): Sea level station monitoring facility. https://www.ioc-sealevelmonitoring.org — DOI 10.14284/482
Not for navigationThe official Indonesian tide tables are published by Pushidrosal.
The chart
The record
- Observed the water level the gauge actually recorded
- Predicted what the fitted constituents say it should be
- Residual the difference — weather, wind and instrument trouble
The rest is held out and predicted. Shorten the window and constituents start dropping out — the same Rayleigh criterion, moved by the window rather than by the length of the record.
Result
Request refused
S2 and K2 cannot be separated on a record of 141.3 days.
- S2 and K2 differ by only 0.0821°/h: separating them needs a record of 182.6 days, and 141.3 are available.
- K1 and P1 differ by only 0.0821°/h: separating them needs a record of 182.6 days, and 141.3 are available.
Record length: 141.3 days · Length required: 182.6 days
What follows is the largest set this window can still support: M2, S2, N2, K1, O1, Q1, M4, MS4.
- Condition number κ
- 1.45good
- Residual RMS
- 0.1050 m
- Mean level Z₀
- 8.0249 m
- Record length
- 141.3 days
Unexplained, over the fitted part
0.1050m
Mean distance between the black line and the blue one, over the part used to fit.
Unexplained, over the held-out part
0.1326m
Over the part deliberately not seen while fitting. This is the real test.
Constituents
Harmonic constituents
| Constituent | Amplitude H (m) | Phase g (°) | Lag (h) | Speed (°/h) | Period (h) | Nodal factor f | Nodal correction u (°) |
|---|---|---|---|---|---|---|---|
| M2255.555 | 0.5016±0.0027 | 86.0±0.3 | 2.97 | 28.984104 | 12.421 | 0.9659 | 0.79 |
| S2273.555 | 0.2713±0.0026 | 123.3±0.5 | 4.11 | 30.000000 | 12.000 | 1.0000 | 0.00 |
| K1165.555 | 0.1010±0.0023 | 37.5±1.3 | 2.49 | 15.041069 | 23.934 | 1.1068 | 2.86 |
| N2245.655 | 0.0912±0.0027 | 82.1±1.7 | 2.89 | 28.439730 | 12.658 | 0.9659 | 0.79 |
| O1145.555 | 0.0488±0.0022 | 4.2±2.6 | 0.30 | 13.943036 | 25.819 | 1.1728 | -3.23 |
| MS4473.555 | 0.0115±0.0027 | 182.6±13.2 | 3.10 | 58.984104 | 6.103 | 0.9659 | 0.79 |
| M4455.555 | 0.0095±0.0028 | 148.4±16.6 | 2.56 | 57.968208 | 6.210 | 0.9329 | 1.57 |
| Q1135.655 | 0.0064±0.0022 | 305.9±19.8 | 22.83 | 13.398661 | 26.868 | 1.1728 | -3.23 |
H is half the height of that constituent’s wave; g is how late it arrives behind the Moon or Sun, and the column beside it says the same thing in hours — g divided by the constituent’s speed. f and u are the 18.6-year corrections, applied and still shown.
The ± figure is one standard error, taken from the diagonal of the solve’s covariance matrix. It assumes the record’s noise is independent from hour to hour, and tidal residuals are not: weather lasts for days, so neighbouring hours are wrong in the same direction. The true uncertainty is larger than what is printed here — by how much, this site does not compute.
A property of the place, or of these months?
The same numbers, four different stretches of the record
This whole site rests on one claim: that amplitude and phase are properties of the place — what a coastline does to the same forcing. But every number above came from one particular stretch of 2026. Here the record is cut into 4 equal, non-overlapping pieces and each is fitted on its own. Where an amplitude barely moves between them it is behaving like a property of the harbour. Where it swings, the record is doing the work rather than the place.
- Stretch 1 01 Jan 2026 — 23 Feb 2026
- Stretch 2 23 Feb 2026 — 17 Apr 2026
- Stretch 3 17 Apr 2026 — 09 Jun 2026
- Stretch 4 09 Jun 2026 — 01 Agu 2026
| Constituent | Mean H (m) | Lowest — highest | Swing | Phase swing |
|---|---|---|---|---|
| M2 | 0.5111 | 0.4874 — 0.5337 | 9% | ±1.4° |
| S2 | 0.2405 | 0.1622 — 0.3239 | 67% | ±19.2° |
| K1 | 0.1084 | 0.0842 — 0.1203 | 33% | ±11.1° |
| N2 | 0.0963 | 0.0827 — 0.1071 | 25% | ±11.2° |
| O1 | 0.0475 | 0.0454 — 0.0511 | 12% | ±2.0° |
| M4 | 0.0126 | 0.0054 — 0.0190 | 108% | ±41.0° |
| MS4 | 0.0111 | 0.0054 — 0.0161 | 97% | ±49.0° |
| Q1 | 0.0063 | 0.0049 — 0.0084 | 55% | ±38.4° |
One constituent, from the start
How M2 came out of this record
The table above is the answer. This is the route to it, in this record’s own numbers — nothing simplified or rounded for the example.
- Step 1 — Where the Moon and Sun are. Every constituent has six Doodson numbers saying how many times each astronomical element enters its angle. For M2 they are 2 0 0 0 0 0, so its equilibrium argument is V(t) = 2τ. The element polynomials are from Meeus, Astronomical Algorithms (2nd ed.), chapters 22 and 47.
τ = mean lunar time · s = the Moon’s mean longitude · h = the Sun’s mean longitude · p = lunar perigee · N = the Moon’s ascending node · p′ = solar perigee.
Step 2 — Two columns, hour by hour. cos(V + u) and sin(V + u) are evaluated at every sample in the record. These are M2’s two columns of the design matrix; four consecutive hours of this record:
How M2 came out of this record Time (UTC) V(t)° cos(V+u) sin(V+u) M2’s share (m) Recorded (m) 12 Mar 2026 03:00 UTC 245.65 -0.3998 -0.9166 -0.4566 7.8330 12 Mar 2026 04:00 UTC 274.63 0.0945 -0.9955 -0.4779 7.9540 12 Mar 2026 05:00 UTC 303.62 0.5650 -0.8251 -0.3794 8.0680 12 Mar 2026 06:00 UTC 332.60 0.8941 -0.4479 -0.1859 8.1730 - Step 3 — Least squares picks one pair of numbers. The solve looks for the a and b that bring Σ (a·cos + b·sin) as close to the record as it can get, for every constituent at once. For M2 it returned a = 0.0342 and b = 0.4832. The last column above is a·cos + b·sin: what M2 alone contributes at that hour.
- Step 4 — Two lines of trigonometry. The pair (a, b) is amplitude and phase in another form: H = √(a² + b²) ÷ f = 0.5016 m, and g = atan2(b, a) = 86.0°. Exactly the numbers on M2’s row above. The division by f takes out the 18.6-year nodal correction, so the constant does not depend on which year you happened to observe.