% Leo Baldreich 2023-02-23 % 10 deg grid of output values of gmf.m at MJD 58849 for an elevation of 5 deg % The grid spans from 90 deg Lat / 0 deg Lon in the top left to -90 deg Lat / 350 deg Lon in the bottom right % All values calculated at an ellipsoidal height taken from orography_ell_10x10.txt % % Example code for recreating the grids without the header: % % % reading ellipsoidal heights from file % hell_10x10 = readmatrix('orography_ell_10x10.txt'); % % % defining lat/lon vectors and epoch for calculation % lat = [90:-10:-90]*pi/180; % lon = [0:10:350]*pi/180; % time = 58849; % % % allocating output array % gmfout_grid = zeros(19,36,2); % % % looping through lat/lon filling the output array % for i = 1:19 % for j = 1:36 % [gmfout_grid(i,j,1),gmfout_grid(i,j,2)] = gmf(time,lat(i),lon(j),hell_10x10(i,j),zd); % end % end % % % writing output to file % file = fopen('gmf_grid.txt','w'); % gmf_ = ["%% gmfh, hydrostatic mapping function:\n",... % "\n%% gmfw, wet mapping function:\n"]; % for k = 1:2 % fprintf(file,gmf_(k)); % for i = 1:19 % for j = 1:36 % fprintf(file,'%7.4f ',gmfout_grid(i,j,k)); % end % fprintf(file,'\n'); % end % end % gmfh, hydrostatic mapping function: 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2089 10.2010 10.2030 10.2089 10.2022 10.2025 10.2039 10.2066 10.2049 10.2063 10.2068 10.2101 10.2072 10.2067 10.2067 10.2059 10.2053 10.2047 10.2036 10.2034 10.2030 10.2027 10.2031 10.2036 10.2044 10.2052 10.2061 10.2065 10.2343 10.2364 10.2171 10.2314 10.2508 10.2526 10.2363 10.2158 10.2008 10.1803 10.1817 10.1898 10.1894 10.1878 10.1901 10.1929 10.1962 10.1979 10.2115 10.2130 10.2145 10.2083 10.2189 10.2116 10.2050 10.1971 10.1960 10.1905 10.1889 10.1910 10.1884 10.1902 10.1943 10.1929 10.2018 10.1986 10.2005 10.2032 10.2225 10.1995 10.2230 10.2552 10.2373 10.1839 10.1811 10.1615 10.1775 10.1671 10.1720 10.1766 10.1791 10.1862 10.1820 10.1846 10.1907 10.2002 10.2073 10.2070 10.2087 10.2190 10.2101 10.1910 10.1875 10.1806 10.1765 10.1875 10.1876 10.2046 10.2016 10.1865 10.1847 10.1900 10.1869 10.1897 10.1921 10.1857 10.1792 10.1732 10.1668 10.1619 10.1609 10.1510 10.1630 10.1716 10.1641 10.1638 10.1598 10.1671 10.1710 10.1734 10.2294 10.2235 10.2179 10.2043 10.1951 10.1998 10.1827 10.1776 10.1716 10.1665 10.1617 10.1582 10.1541 10.1518 10.1504 10.1854 10.1775 10.1740 10.1769 10.1736 10.1818 10.1655 10.1592 10.1540 10.1494 10.1477 10.1481 10.1551 10.1462 10.1617 10.1707 10.1997 10.1421 10.1443 10.1940 10.1743 10.1664 10.1850 10.1851 10.1624 10.1601 10.1622 10.1554 10.1505 10.1461 10.1430 10.1406 10.1383 10.1365 10.1351 10.1342 10.1794 10.1886 10.1564 10.1462 10.1528 10.1418 10.1393 10.1373 10.1353 10.1347 10.1349 10.1367 10.1322 10.1323 10.1258 10.1278 10.1415 10.1224 10.1550 10.1471 10.2035 10.2469 10.2305 10.1578 10.1359 10.1269 10.1241 10.1216 10.1193 10.1180 10.1180 10.1180 10.1182 10.1188 10.1197 10.1197 10.1201 10.1451 10.1346 10.1184 10.1175 10.1166 10.1166 10.1178 10.1188 10.1193 10.1196 10.1247 10.1138 10.1156 10.1161 10.1126 10.1062 10.1111 10.1023 10.1009 10.1083 10.1050 10.1244 10.1073 10.1029 10.1021 10.1016 10.1008 10.1003 10.1008 10.1024 10.1038 10.1057 10.1065 10.1077 10.1083 10.1066 10.1051 10.1646 10.1044 10.1031 10.1020 10.1033 10.1053 10.1070 10.1072 10.1057 10.1144 10.1019 10.1106 10.1055 10.1061 10.1448 10.1119 10.0953 10.0938 10.0942 10.0963 10.0973 10.0983 10.0982 10.0979 10.0984 10.0970 10.0965 10.0969 10.0980 10.0992 10.1002 10.1012 10.1020 10.1015 10.0999 10.0986 10.0986 10.0988 10.0990 10.1278 10.0992 10.1008 10.1018 10.1010 10.0989 10.1084 10.0979 10.1013 10.1051 10.1363 10.1022 10.0973 10.0968 10.0954 10.0948 10.0958 10.1149 10.1025 10.1040 10.1025 10.1006 10.0996 10.0985 10.0980 10.0980 10.0984 10.0994 10.1004 10.1006 10.1003 10.0996 10.0993 10.0994 10.0992 10.1042 10.1038 10.1026 10.0995 10.0999 10.0999 10.0990 10.0983 10.1000 10.0993 10.1250 10.1307 10.0989 10.0989 10.0990 10.0984 10.0975 10.0973 10.0971 10.0974 10.1028 10.0985 10.0984 10.1100 10.1041 10.0988 10.0983 10.0983 10.0990 10.1003 10.1010 10.1011 10.1016 10.1024 10.1020 10.1004 10.1009 10.1033 10.1050 10.1032 10.1116 10.0999 10.1004 10.1005 10.1004 10.0993 10.1251 10.1257 10.0982 10.0998 10.1001 10.1003 10.1001 10.0994 10.0988 10.0973 10.0967 10.1035 10.1018 10.0978 10.1000 10.1020 10.1007 10.1002 10.1007 10.1018 10.1023 10.1026 10.1038 10.1043 10.1036 10.1015 10.0993 10.1326 10.1030 10.1084 10.1015 10.1011 10.1018 10.1016 10.1052 10.1036 10.1241 10.1339 10.1036 10.1042 10.1050 10.1056 10.1063 10.1068 10.1065 10.1045 10.1127 10.1051 10.1038 10.1149 10.1063 10.1077 10.1072 10.1062 10.1063 10.1067 10.1067 10.1065 10.1069 10.1081 10.1081 10.1074 10.1052 10.2077 10.1050 10.1069 10.1048 10.1067 10.1074 10.1068 10.1201 10.1190 10.1184 10.1185 10.1186 10.1183 10.1181 10.1186 10.1195 10.1207 10.1210 10.1200 10.1181 10.1162 10.1160 10.1170 10.1180 10.1175 10.1172 10.1167 10.1165 10.1167 10.1164 10.1161 10.1169 10.1185 10.1199 10.1205 10.1198 10.1415 10.1175 10.1163 10.1172 10.1189 10.1209 10.1214 10.1413 10.1411 10.1407 10.1404 10.1399 10.1391 10.1380 10.1392 10.1379 10.1385 10.1387 10.1383 10.1372 10.1357 10.1346 10.1337 10.1326 10.1316 10.1310 10.1304 10.1307 10.1312 10.1315 10.1318 10.1332 10.1349 10.1368 10.1380 10.1381 10.1478 10.1353 10.1343 10.1350 10.1369 10.1392 10.1408 10.1597 10.1607 10.1609 10.1606 10.1604 10.1591 10.1586 10.1578 10.1573 10.1571 10.1567 10.1562 10.1554 10.1542 10.1528 10.1513 10.1499 10.1485 10.1479 10.1476 10.1481 10.1481 10.1487 10.1492 10.1504 10.1514 10.1525 10.1532 10.1532 10.1524 10.1522 10.1519 10.1527 10.1544 10.1566 10.1582 10.1719 10.1758 10.1759 10.1841 10.2024 10.2269 10.2309 10.1725 10.2205 10.2315 10.2283 10.2269 10.2253 10.2228 10.2225 10.2057 10.1866 10.1655 10.1653 10.1651 10.1642 10.1642 10.1643 10.1643 10.1647 10.1647 10.1648 10.1644 10.1640 10.1853 10.1654 10.1644 10.1652 10.1669 10.1687 10.1700 10.2388 10.2513 10.2564 10.2591 10.2606 10.2575 10.2570 10.2629 10.2660 10.2672 10.2600 10.2524 10.2482 10.2409 10.2322 10.2252 10.2074 10.1794 10.1785 10.1761 10.1781 10.1825 10.1904 10.2031 10.2085 10.2160 10.2219 10.2178 10.1886 10.1741 10.1770 10.1945 10.1760 10.1867 10.2239 10.2389 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 10.2443 % gmfw, wet mapping function: 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.7557 10.8233 10.8234 10.8217 10.8181 10.8126 10.8054 10.7969 10.7873 10.7773 10.7674 10.7579 10.7492 10.7416 10.7353 10.7303 10.7267 10.7243 10.7230 10.7225 10.7227 10.7235 10.7247 10.7264 10.7290 10.7327 10.7378 10.7446 10.7531 10.7631 10.7739 10.7849 10.7953 10.8044 10.8119 10.8176 10.8213 10.8396 10.8402 10.8397 10.8365 10.8293 10.8187 10.8062 10.7937 10.7826 10.7732 10.7653 10.7591 10.7546 10.7523 10.7521 10.7537 10.7567 10.7606 10.7648 10.7684 10.7703 10.7698 10.7664 10.7607 10.7544 10.7498 10.7494 10.7550 10.7669 10.7836 10.8018 10.8181 10.8297 10.8361 10.8385 10.8391 10.8049 10.8086 10.8122 10.8101 10.7996 10.7842 10.7713 10.7671 10.7712 10.7781 10.7825 10.7837 10.7846 10.7881 10.7941 10.8007 10.8067 10.8123 10.8182 10.8241 10.8285 10.8294 10.8245 10.8130 10.7957 10.7757 10.7580 10.7485 10.7511 10.7656 10.7862 10.8048 10.8149 10.8154 10.8102 10.8053 10.7920 10.7960 10.7980 10.7905 10.7700 10.7453 10.7351 10.7515 10.7860 10.8155 10.8232 10.8134 10.8033 10.8038 10.8119 10.8184 10.8188 10.8159 10.8135 10.8139 10.8171 10.8218 10.8257 10.8254 10.8173 10.8000 10.7753 10.7508 10.7370 10.7407 10.7582 10.7782 10.7905 10.7935 10.7918 10.7904 10.8147 10.8093 10.8052 10.7978 10.7778 10.7521 10.7488 10.7879 10.8503 10.8882 10.8730 10.8257 10.7901 10.7861 10.7989 10.8056 10.7998 10.7888 10.7788 10.7707 10.7666 10.7698 10.7812 10.7973 10.8124 10.8190 10.8096 10.7857 10.7615 10.7533 10.7634 10.7794 10.7912 10.7996 10.8088 10.8155 10.8035 10.7882 10.7912 10.8110 10.8212 10.8112 10.8103 10.8499 10.9107 10.9318 10.8818 10.8014 10.7547 10.7599 10.7840 10.7934 10.7880 10.7827 10.7791 10.7703 10.7586 10.7539 10.7599 10.7744 10.7963 10.8212 10.8350 10.8270 10.8074 10.7975 10.8026 10.8088 10.8069 10.8057 10.8121 10.8153 10.7275 10.7176 10.7473 10.8095 10.8595 10.8666 10.8532 10.8603 10.8847 10.8785 10.8180 10.7454 10.7184 10.7425 10.7743 10.7837 10.7835 10.7936 10.8063 10.8043 10.7905 10.7797 10.7748 10.7715 10.7754 10.7936 10.8144 10.8190 10.8094 10.8051 10.8114 10.8117 10.7960 10.7760 10.7632 10.7493 10.6588 10.6666 10.7170 10.7907 10.8437 10.8473 10.8168 10.7884 10.7749 10.7569 10.7206 10.6872 10.6870 10.7163 10.7424 10.7490 10.7542 10.7741 10.7954 10.7989 10.7895 10.7826 10.7781 10.7672 10.7559 10.7569 10.7658 10.7654 10.7538 10.7462 10.7501 10.7546 10.7482 10.7313 10.7080 10.6803 10.6720 10.6909 10.7321 10.7665 10.7759 10.7614 10.7331 10.7035 10.6832 10.6745 10.6699 10.6653 10.6659 10.6753 10.6866 10.6941 10.7029 10.7181 10.7326 10.7376 10.7376 10.7430 10.7526 10.7574 10.7561 10.7539 10.7501 10.7374 10.7140 10.6896 10.6777 10.6857 10.7074 10.7227 10.7134 10.6865 10.7430 10.7601 10.7748 10.7582 10.7226 10.6992 10.6940 10.6923 10.6900 10.6957 10.7047 10.6992 10.6768 10.6577 10.6579 10.6704 10.6791 10.6801 10.6806 10.6859 10.6959 10.7110 10.7323 10.7570 10.7771 10.7859 10.7828 10.7699 10.7453 10.7083 10.6750 10.6733 10.7108 10.7559 10.7686 10.7515 10.7884 10.8021 10.8051 10.7741 10.7318 10.7158 10.7296 10.7487 10.7605 10.7715 10.7778 10.7618 10.7238 10.6929 10.6916 10.7073 10.7122 10.7005 10.6898 10.6929 10.7042 10.7162 10.7324 10.7576 10.7844 10.8012 10.8082 10.8115 10.8049 10.7754 10.7331 10.7140 10.7391 10.7824 10.8015 10.7925 10.7765 10.7912 10.8013 10.7894 10.7681 10.7612 10.7718 10.7850 10.7938 10.8022 10.8064 10.7936 10.7657 10.7456 10.7492 10.7642 10.7666 10.7524 10.7382 10.7350 10.7362 10.7341 10.7335 10.7423 10.7578 10.7721 10.7863 10.8049 10.8185 10.8085 10.7759 10.7491 10.7513 10.7714 10.7813 10.7762 10.7521 10.7639 10.7766 10.7793 10.7728 10.7657 10.7617 10.7596 10.7604 10.7668 10.7746 10.7749 10.7673 10.7631 10.7703 10.7824 10.7868 10.7807 10.7718 10.7649 10.7570 10.7451 10.7341 10.7305 10.7343 10.7421 10.7546 10.7727 10.7889 10.7899 10.7734 10.7541 10.7469 10.7504 10.7525 10.7500 10.7606 10.7640 10.7677 10.7676 10.7619 10.7529 10.7438 10.7380 10.7383 10.7449 10.7536 10.7589 10.7593 10.7588 10.7618 10.7678 10.7729 10.7751 10.7750 10.7730 10.7682 10.7612 10.7550 10.7527 10.7543 10.7579 10.7630 10.7692 10.7739 10.7727 10.7656 10.7578 10.7552 10.7572 10.7594 10.7598 10.7991 10.7995 10.7987 10.7966 10.7933 10.7898 10.7876 10.7880 10.7912 10.7959 10.7992 10.7983 10.7931 10.7855 10.7786 10.7745 10.7735 10.7749 10.7777 10.7809 10.7839 10.7867 10.7898 10.7934 10.7965 10.7980 10.7971 10.7941 10.7897 10.7847 10.7808 10.7798 10.7825 10.7879 10.7934 10.7973 10.8402 10.8458 10.8506 10.8550 10.8594 10.8638 10.8684 10.8725 10.8755 10.8761 10.8732 10.8663 10.8559 10.8434 10.8307 10.8195 10.8111 10.8059 10.8037 10.8040 10.8061 10.8093 10.8129 10.8160 10.8178 10.8179 10.8160 10.8127 10.8091 10.8063 10.8055 10.8075 10.8121 10.8188 10.8263 10.8336 10.8652 10.8722 10.8791 10.8855 10.8914 10.8964 10.9003 10.9028 10.9035 10.9022 10.8989 10.8936 10.8868 10.8788 10.8703 10.8619 10.8541 10.8474 10.8420 10.8378 10.8348 10.8327 10.8313 10.8303 10.8296 10.8290 10.8285 10.8283 10.8286 10.8298 10.8319 10.8352 10.8396 10.8451 10.8514 10.8582 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666 10.8666