added function to compute model forecast error based on WRF forecasts