矢量化版本的Matlab代码GMP和GPT模型。原始Fortran和Matlab代码可以在日前发表2010公约(http://iers-conventions.obspm.fr/conventions_versions.php official_target),
维也纳映射功能开放获取数据(http://vmf.geo.tuwien.ac.at/codes/)

测试示例(Matlab代码):

n = 86400;
% % GMF
dmjd = 55055 *的(n - 1);
dlat = 0.6708665767 *的(n - 1);
dlon = -1.393397187 *的(n - 1);
dhgt = 844.715 *的(n - 1);
zd = 1.278564131 *的(n - 1);
抽搐
[gmfh, gmfw] = gmf (dmjd、dlat dlon, dhgt, zd);
toc
%测试用例:
%给定的输入:DMJD = 55055 d0
% DLAT = 0.6708665767 d0弧度(西弗吉尼亚州NRAO,绿色银行的)
% DLON = -1.393397187 d0弧度
% DHGT = 844.715 d0米
% ZD = 1.278564131 d0弧度
%
%的预期输出:GMFH = 3.425245519339138678 - 11.064702723565059
% GMFW = 3.449589116182419257 - 11.313101285827178
流(GMF % .15f \ t %。15 f \ n ', gmfh (1) gmfw (1));
流(' \ n ');

抽搐
[gmfh, gmfw] = gmf_deriv (dmjd、dlat dlon, dhgt, zd);
toc
流(“GMF_deriv % .15f \ t % .15f \ t % .15f \ %。15 f \ n ', gmfh (1:), gmfw (: 1));
流(' \ n ');
% % GPT
dmjd = 55055 *的(n - 1);
dlat = 0.6708665767 *的(n - 1);
dlon = -1.393397187 *的(n - 1);
dhgt = 812.546 *的(n - 1);
抽搐
(临时总统,过分]= gpt (dmjd、dlat dlon, dhgt);
toc
%测试用例:
%给定的输入:DMJD = 55055 d0
% DLAT = 0.6708665767 d0弧度(西弗吉尼亚州NRAO,绿色银行的)
% DLON = -1.393397187 d0弧度
% DHGT = 812.546米
%的预期输出:总统= 918.0710638757363995 d0 hPa
% TEMP = 19.31914181012882992 d0摄氏度
%过分= -42.19185643717770517 d0米
流(“GPT % .15f \ t % .15f \ %。15 f \ n”,总统(1),临时(1),过分(1));
流(' \ n ');

享受吧!