stn_id = {'bbry';'ldes';'jplm';'widc'};

istn = 4;

[mat_date, year, doy, Un, Ue, Uu, Sn, Se, Su] = load_gps_data(stn_id{istn});

subplot(3,1,1)
plot(mat_date,Un,'.');
datetick('x',10)
xlabel('Date')
ylabel('Un, m')
title(stn_id{istn})

subplot(3,1,2)
plot(mat_date,Ue,'.');
datetick('x',10)
xlabel('Date')
ylabel('Ue, m')
title(stn_id{istn})

subplot(3,1,3)
plot(mat_date,Uu,'.');
datetick('x',10)
xlabel('Date')
ylabel('Uu, m')
title(stn_id{istn})