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

soil_stn_mean = squeeze(mean(soil_stn,3));
soil_wrf_mean = squeeze(mean(soil_wrf,3));

stn  = soil_stn_mean;
    for irun =1:13
        runopt = soil_wrf_mean(1,:,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 
