Investigating the natural variability of sea-surface heat fluxes with an intermediate-complexity coupled model