GMP和GPT模型的矢量化版本的Matlab代码。原始Fortran和Matlab代码可在IERS 2010公约中找到(http://iers-conventions.obspm.fr/conventions_versions.php#official_target),
维也纳地图功能开放存取数据(http://vmf.geo.tuwien.ac.at/codes/）

测试示例(Matlab代码):

N = 86400;
% % GMF
DMJD = 55055*ones(n,1);
dat = 0.6708665767*ones(n,1);
Dlon = -1.393397187*ones(n,1);
DHGT = 844.715*ones(n,1);
Zd = 1.278564131*ones(n,1);
抽搐
[gmfh,gmfw] = gmf(dmjd, dat,dlon,dhgt,zd);
toc
测试用例:
%给定输入:DMJD = 55055D0
% DLAT = 0.6708665767D0弧度(NRAO, Green Bank, WV)
% DLON = -1.393397187D0弧度
% DHGT = 844.715D0米
% ZD = 1.278564131D0弧度
％
%预期产量:GMFH = 3.425245519339138678 11.064702723565059
% GMFW = 3.449589116182419257 11.313101285827178
流(GMF % .15f \ t %。15f\n'， gmfh(1)， gmfw(1));
流(' \ n ');

抽搐
[gmfh,gmfw] = gmf_deriv(dmjd, dat,dlon,dhgt,zd);
toc
流(“GMF_deriv % .15f \ t % .15f \ t % .15f \ %。15f\n'， gmfh(1，:)， gmfw(1，:));
流(' \ n ');
% % GPT
DMJD = 55055*ones(n,1);
dat = 0.6708665767*ones(n,1);
Dlon = -1.393397187*ones(n,1);
DHGT = 812.546*ones(n,1);
抽搐
[pres,temp,undu] = gpt(dmjd, dat, dlon, dhgt);
toc
测试用例:
%给定输入:DMJD = 55055D0
% DLAT = 0.6708665767D0弧度(NRAO, Green Bank, WV)
% DLON = -1.393397187D0弧度
% DHGT = 812.546米
%预期输出:PRES = 918.0710638757363995D0 hPa
% TEMP = 19.31914181012882992D0摄氏度
% UNDU = -42.19185643717770517D0米
流(“GPT % .15f \ t % .15f \ %。15f\n'， pres(1)， temp(1)， undu(1));
流(' \ n ');

享受吧!