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

fname = [stn_id 'CleanUnf.neu'];

a = load(fname); %file from SOPAC

decimal_year = a(:,1);
year         = a(:,2);
doy          = a(:,3); %day of year

Un = a(:,4);
Ue = a(:,5);
Uu = a(:,6);
Sn = a(:,7);
Se = a(:,8);
Su = a(:,9);

ndays = length(Uu);


for i = 1:ndays
  mat_date(i) = yearday(year(i),doy(i)); %internal matlab date format
end