write/run scripts to compute 22year anomalies