Improving Madden–Julian oscillation simulation in atmospheric general circulation models by coupling with a one-dimensional snow–ice–thermocline ocean model