国际参考地磁场(IGRF)是国际商定地球磁场的数学模型。这个程序是一个FORTRAN子例程的计算转化为MATLAB。它不使用FORTRAN墨西哥人编译文件,这可能使它慢但更容易使用的优势(如不需要编译)。事实上,我写程序的目的就是提供一个IGRF实现MATLAB以最小的“小题大做”。Another motivation was a vectorized IGRF function, which this function is (with a separate routine adapted directly from the FORTRAN code that is faster for scalars implemented as well).

提供以下文件:
-igrf。m:地球磁场计算点(年代)。
-igrfline。m:给出了坐标沿着磁场线从一个给定的点。
-getigrfcoefs。m:从.dat文件提取系数提供了IGRF网站.mat文件并保存它们。
-igrfcoefs。12 IGRF代垫:IGRF系数(最近的2015)。
-loadigrfcoefs。m:加载适当的IGRF系数在给定的时间(做必要的插值)。
- *平*。dat: 12代IGRF系数数据文件。
-plotbline:情节磁场线。
-plotbearth:情节磁场线。
-geod2ecef:大地坐标转换为ECEF坐标。
-ecef2geod: ECEF坐标转换为大地。

运行的唯一先决条件IGRF函数或函数IGRFLINE igrfcoefs是将文件。垫在MATLAB的搜索路径中。程序被设计为可伸缩与时间:新IGRF代被释放,只是将旧.dat文件替换为新版本在一个子文件夹称为“datfiles”getigrfcoefs函数在相同的目录中。m是位于getigrfcoefs运行,然后替换它生成的文件(igrfcoefs.mat)老.mat文件。每五年更新发生,持续更新发生在2015年。新的.dat文件将希望继续上传到以下网站(虽然似乎是2007年以来的文件没有更新):
https://ccmc.gsfc.nasa.gov/pub/modelweb/geomagnetic/igrf/fortran_code/

最后,我有两个示例脚本显示IGRFLINE函数是如何工作的:plotbline。m和plotbearth.m。这些脚本都使用的映射工具箱情节地球仪磁场线绘制,但如果用户没有包,一个全球原油只有显示经度和纬度线。

我做了一些粗略的比较与在线IGRF计算器http://ccmc.gsfc.nasa.gov/modelweb/models/igrf_vitmo.php,发现这个函数精确到1元。我不清楚为什么两者之间有差异,但我猜想舍入误差。