Skip to content

Padang, Sumatera Barat

Mixed, mainly semidiurnal
Period
01 Jan 202601 Agu 2026212 days
Datum
Nol sensor stasiun (bukan MSL)
gaps
1333 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

The observed record, the fitted prediction overlaid on it, and the residual in its own band on the same time axis. The residual holds weather, surge, and everything the harmonic model does not explain.
  • 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
67%

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.

Tide chart for Padang, Sumatera Barat over 212 days: a black line for the recorded water level, a blue line over it computed from 8 constituents, and an ochre band beneath for the difference, on the same time axis. Residual RMS 0.0913 m.This part was not seen during fittingZ₀ ioc-sensor-zero7.007.508.008.5002/0130/0127/0227/0324/0422/0519/0617/070.00Residual
ObservedPredictedResidual
Condition number κ
1.44good
Residual RMS
0.0913 m
Mean level Z₀
7.5777 m
Record length
141.3 days

Unexplained, over the fitted part

0.0913m

Mean distance between the black line and the blue one, over the part used to fit.

Unexplained, over the held-out part

0.1005m

Over the part deliberately not seen while fitting. This is the real test.

Constituents

Harmonic constituents

Harmonic constants per constituent
ConstituentAmplitude H (m)Phase g (°)Lag (h)Speed (°/h)Period (h)Nodal factor fNodal correction u (°)
M2255.5550.3565±0.0023337.5±0.411.6428.98410412.4210.96590.79
S2273.5550.1651±0.002221.8±0.80.7330.00000012.0001.00000.00
K1165.5550.1188±0.00201.6±1.00.1115.04106923.9341.10682.86
O1145.5550.0763±0.0019347.9±1.424.9513.94303625.8191.1728-3.23
N2245.6550.0680±0.0023322.5±1.911.3428.43973012.6580.96590.79
M4455.5550.0162±0.002433.8±8.40.5857.9682086.2100.93291.57
Q1135.6550.0161±0.0019331.0±6.824.7013.39866126.8681.1728-3.23
MS4473.5550.0071±0.0023156.2±18.62.6558.9841046.1030.96590.79

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.

  1. Stretch 1 01 Jan 202623 Feb 2026
  2. Stretch 2 23 Feb 202617 Apr 2026
  3. Stretch 3 17 Apr 202609 Jun 2026
  4. Stretch 4 09 Jun 202601 Agu 2026
The same numbers, four different stretches of the record
ConstituentMean H (m)Lowest — highestSwingPhase swing
M20.35900.3580 0.36001%±0.7°
S20.14950.1093 0.194857%±20.3°
K10.12590.0971 0.138333%±12.9°
O10.07720.0748 0.08007%±4.5°
N20.07520.0630 0.086131%±7.0°
M40.01690.0140 0.018526%±4.2°
Q10.01620.0122 0.017734%±16.1°
MS40.00690.0060 0.007522%±34.1°

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.

  1. 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.

  2. 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 06:00 UTC332.600.8941-0.4479+0.34347.8820
    12 Mar 2026 07:00 UTC1.590.99910.0414+0.31237.8580
    12 Mar 2026 08:00 UTC30.570.85390.5204+0.20307.7570
    12 Mar 2026 09:00 UTC59.550.49480.8690+0.04287.6490
  3. 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.3181 and b = -0.1319. The last column above is a·cos + b·sin: what M2 alone contributes at that hour.
  4. Step 4 — Two lines of trigonometry. The pair (a, b) is amplitude and phase in another form: H = √(a² + b²) ÷ f = 0.3565 m, and g = atan2(b, a) = 337.5°. 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.

Method