function [rmse, corr,bias] = calc_rmse_r (soil_stn, soil_wrf, stn_number)


stn  = soil_stn(1,:, stn_number);
for irun =1:10
    runopt = soil_wrf(1,:,stn_number,irun);

    rmse(irun) = sqrt(mean( runopt - stn ).^2);
    
    corr_temp = corrcoef(stn, runopt);
    corr(irun) = corr_temp(1,2);
    bias(irun) = mean(abs(runopt - stn));
end


end 
