电离层-计算(修正)地磁纬度

1.下载XCFRAME

见MathWorks

2.导入数据

T=repmat([2020,05,11,0,0,0],[73,1]);
lon=(-180:5:180)';
lat=(-90:2.5:90)';
height=ones(numel(lon),1)*350;

3.计算地磁纬度

temp=xcframe(T,[lat,lon,height], [], [], 'GEO', 'MAG', 'geod84', 'sph' );
mlat=temp(:,2);

4.计算修正地磁纬度

[~,~,D,I]=igrfmagm(height,lat,lon,decyear(T),13);
dip=rad2deg(atan(deg2rad(I)./sqrt(cosd(lat))));

① 地磁纬度

② 修正地磁纬度