add code to reproduce the temperature anomaly result