% Leo Baldreich 2023-02-22 % Time series of output values of gmf.m for station GRAZ % Lat: 47.0671 deg, Lon: 15.4935 deg, Ell. height: 538.30 m % Mapping functions calculated for 5 deg elevation % % Example code for recreating this table without the header: % % % defining station coordinates, epochs for calculation and elevation angle % GRAZ = [47.0671*pi/180,15.4935*pi/180,538.30]; % mjds = 58849:1/4:58859; % zd = 85*pi/180; % % % allocating output array % gmfout = zeros(length(mjds),2); % % % looping through epochs and filling the output array % for i = 1:length(mjds) % [gmfout(i,1), gmfout(i,2)] = gmf(mjds(i),GRAZ(1),GRAZ(2),GRAZ(3),zd); % end % % % writing output to file % file = fopen('gmf_series.txt','w'); % for i = 1:length(mjds) % fprintf(file,'%8.2f %7.4f %7.4f\n',[mjds(i),gmfout(i,:)]); % end % columns: % 1: MJD for which the values were computed % 2: gmfh, hydrostatic mapping function % 3: gmfw, wet mapping function 58849.00 10.1636 10.8001 58849.25 10.1637 10.8002 58849.50 10.1637 10.8002 58849.75 10.1637 10.8002 58850.00 10.1638 10.8003 58850.25 10.1638 10.8003 58850.50 10.1639 10.8004 58850.75 10.1639 10.8004 58851.00 10.1639 10.8004 58851.25 10.1640 10.8005 58851.50 10.1640 10.8005 58851.75 10.1641 10.8005 58852.00 10.1641 10.8006 58852.25 10.1641 10.8006 58852.50 10.1642 10.8006 58852.75 10.1642 10.8007 58853.00 10.1642 10.8007 58853.25 10.1643 10.8007 58853.50 10.1643 10.8008 58853.75 10.1643 10.8008 58854.00 10.1644 10.8008 58854.25 10.1644 10.8009 58854.50 10.1644 10.8009 58854.75 10.1645 10.8009 58855.00 10.1645 10.8010 58855.25 10.1645 10.8010 58855.50 10.1646 10.8010 58855.75 10.1646 10.8011 58856.00 10.1646 10.8011 58856.25 10.1647 10.8011 58856.50 10.1647 10.8012 58856.75 10.1647 10.8012 58857.00 10.1648 10.8012 58857.25 10.1648 10.8012 58857.50 10.1648 10.8013 58857.75 10.1648 10.8013 58858.00 10.1649 10.8013 58858.25 10.1649 10.8013 58858.50 10.1649 10.8014 58858.75 10.1650 10.8014 58859.00 10.1650 10.8014