【GPS&Matlab笔记】根据广播星历（.12N文件）计算卫星位置代码

最新推荐文章于 2023-12-14 16:11:18 发布

一滴色彩雫

最新推荐文章于 2023-12-14 16:11:18 发布

阅读量4.2k

点赞数 9

分类专栏： gnss 文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/qq_46922074/article/details/125562878

版权

gnss 专栏收录该内容

3 篇文章 1 订阅

订阅专栏

2022年6月、GNSS与Matlab课程结合大作业测试代码。
matlab源代码：

%测试文件
clear,clc;
%%
%操作部分：
filename='XS013100.12N';
pathout = 'GPS_sat_Position7126.txt';
%若没有就创建，若已存在就清空txt文件
fid=fopen(pathout,'w');
fclose(fid);
%从PRN1号卫星开始
for startLine = 4:8:68 %对于这个文件,第5行开始隔8行到69行一个卫星。
    XYZ=Culpos(startLine,filename);
    outPut(startLine,XYZ)
end
%%
%读取卫星星历文件并计算函数部分：
function [XYZ]=Culpos(startLine,filename)
    fid=fopen(filename);
    %按行读取星历文件的每一个元素
    data=textscan(fid,'%f %f %f %f %f %f %f %f %f%f\n','HeaderLines',startLine);%第startLine+1行开扫
    %将读取的数组转换成数值型矩阵
    B1=cell2mat(data);
    data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
    B2=cell2mat(data);
    data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
    B3=cell2mat(data);
    data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
    B4=cell2mat(data);
    data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
    B5=cell2mat(data);
    data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
    B6=cell2mat(data);
    data=textscan(fid,'%f %f %f %f\n','HeaderLines',1);
    B7=cell2mat(data);
    data=textscan(fid,'%f\n','HeaderLines',1);
    B8=cell2mat(data);
    %%
    %定义常量GM,Ve
    GM=3.986004e14;
    Ve=7.2921151467e-5;
    %求出轨道长半轴a
    sqrta=B3(4);
    a=sqrta^2;
    e=B3(2);
    %求出平均角速度n0
    n0=sqrt(GM/(a^3));
    XYZ=[];
    %每1分钟的卫星坐标
    for i=0:1:20
        %计算当前时刻的儒略日
        JulianDay=fix(365.25*(2004-1))+fix(30.6001*(1+13))+30+(16+i/60+0/3600)/24+1720981.5;
        %将儒略日转换成GPS时
    t2=(JulianDay-2444244.5)/7-fix((JulianDay-2444244.5)/7);
    t=t2*604800;
    %计算卫星钟差改正后的当前时刻的GPS时
    t=t-(B1(8)+B1(9)*(t-B4(1))+B1(10)*(t-B4(1))^2);
    %计算参考时刻的儒略日
    JulianDayoe=fix(365.25*(2004-1))+fix(30.6001*(1+13))+30+(8+0/60+0/3600)/24+1720981.5;
    %将儒略日转换成GPS时
    t2oe=(JulianDayoe-2444244.5)/7-fix((JulianDayoe-2444244.5)/7);
    toe=t2oe*604800;
    %计算时间差tk
    tk=t-toe;
    %计算改正平角速度
    n=n0+B2(3);
    %计算平近点角Ek
    Mk=B2(4)+n*tk;
    %迭代计算，求出偏近点角Ek
    Ek=Mk;
    Ek0=inf;
    while abs(Ek0-Ek)>1e-100
        Ek0=Ek
        Ek=Mk+e*sin(Ek0)
    end
    %计算真近点角
    Sinfk=(sqrt(1-e^2)*sin(Ek))/(1-e*cos(Ek));
    Cosfk=(cos(Ek)-e)/(1-e*cos(Ek));
    fk=atan2(Sinfk,Cosfk);
    %计算升交角距
    phik=fk+B5(3);
    %从卫星星历文件中读取各个方向的振幅Crs,Cuc,Cus,Cic,Cis,Crc
    Crs=B2(2);Cuc=B3(1);Cus=B3(3);Cic=B4(2);Cis=B4(4);Crc=B5(2);
    %计算卫星轨道摄动项改正数
    deltauk=Cus*sin(2*phik)+Cuc*cos(2*phik);
    deltark=Crs*sin(2*phik)+Crc*cos(2*phik);
    deltaik=Cis*sin(2*phik)+Cic*cos(2*phik);
    %计算改正后的升交角距uk
    uk=phik+deltauk;
    %计算改正后的向径rk
    rk=a*(1-e*cos(Ek))+deltark;
    %计算改正后的倾角ik
    ik=B5(1)+deltaik+B6(1)*tk;
    %计算观测瞬间升交点的经度Lk
    Lk=B4(3)+(B5(4)-Ve)*tk-Ve*B4(1);
    %计算卫星在轨道平面内的坐标x,y,z
    xyz=[rk*cos(uk);rk*sin(uk);0];
    %计算旋转矩阵的值
    RzRx=[cos(Lk),-sin(Lk)*cos(ik),sin(Lk)*sin(ik);sin(Lk),cos(Lk)*cos(ik),-cos(Lk)*sin(ik);0,sin(ik),cos(ik)];
    %计算卫星在协议地球坐标系中的位置
    XYZ=vpa([XYZ;1e-3*(RzRx*xyz)'],10);
    end
end
%%
%输出到txt文件中部分：
function []=outPut(startLine,XYZ)
    pathout = 'GPS_sat_Location7126.txt';
    fid=fopen(pathout,'a');
    fprintf(fid,'提取开始的行号为：'); 
    fprintf(fid,'%d',startLine+1);  
    fprintf(fid,'----------------------Its7126testCode_sep_Line_  ^_^  ----------------------------一滴色彩雫\n'); 
    fprintf(fid,[repmat('%f\t', 1, size(XYZ,2)), '\n'], XYZ'); %以矩阵原样输出
    fclose(fid);
end
%%

XS013100.12N
文件：

     2.10           N: GPS NAV DATA                         RINEX VERSION / TYPE
HuaceNav                                                    PGM / RUN BY / DATE 
     GPS RECEIVER NAVIGATION DATA                           COMMENT             
                                                            END OF HEADER       
01 12 11  5  2  0  0.0  2.80956272036D-04  1.81898940355D-12  0.00000000000D+00
     0.00000000000D+00  1.14062500000D+01  4.75805533496D-09  5.13439931108D-01
     5.42029738426D-07  1.38616340701D-03  6.88061118126D-06  5.15370963478D+03
     9.36000000000D+04  1.49011611938D-08  2.16953569760D+00  4.65661287308D-08
     9.59938584427D-01  2.47906250000D+02  5.05762664302D-01 -8.10033741157D-09
     1.99651173416D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00  8.38190317154D-09  3.30000000000D+01
     9.36000000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
07 12 11  5  1 59 44.0  1.41716562212D-04  4.54747350886D-12  0.00000000000D+00
     0.00000000000D+00  2.80937500000D+01  4.37125350883D-09  1.42192075099D-01
     1.32992863655D-06  6.06139330193D-03  3.64892184734D-06  5.15371718979D+03
     9.35840000000D+04  9.31322574615D-08 -9.52836101060D-01  7.45058059692D-09
     9.75746915227D-01  3.14562500000D+02 -2.94195569443D+00 -8.25927260328D-09
    -4.64305054455D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.11758708954D-08  3.70000000000D+01
     9.35840000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
08 12 11  5  1 59 28.0  2.69571319222D-06  2.27373675443D-13  0.00000000000D+00
     0.00000000000D+00  1.90625000000D+01  3.86658963041D-09 -4.89043260306D-01
     9.16421413422D-07  1.27701068996D-02  3.81842255592D-06  5.15363484001D+03
     9.35680000000D+04  2.25380063057D-07 -8.66009063291D-01  1.86264514923D-08
     9.98214108264D-01  3.23968750000D+02 -2.90504734069D+00 -8.06283584948D-09
    -3.67872466222D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -4.19095158577D-09  1.00000000000D+00
     9.35680000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
11 12 11  5  2  0  0.0 -3.19930724800D-04 -3.86535248253D-12  0.00000000000D+00
     0.00000000000D+00  1.63125000000D+01  6.41991027217D-09  3.18685790199D-01
     8.30739736557D-07  1.37224610662D-02  7.13206827641D-06  5.15367549706D+03
     9.36000000000D+04 -2.47731804848D-07  1.88761796145D+00  1.52736902237D-07
     8.88075336314D-01  2.05562500000D+02  1.11484555642D+00 -9.11002232614D-09
     2.32509684962D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.16415321827D-08  6.60000000000D+01
     9.36000000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
17 12 11  5  2  0  0.0  8.15535895526D-05 -3.41060513165D-12  0.00000000000D+00
     0.00000000000D+00 -4.35625000000D+01  4.16017328792D-09  2.95444425157D+00
    -2.19047069550D-06  7.67974182963D-03  1.01346522570D-05  5.15375967216D+03
     9.36000000000D+04  1.67638063431D-08  1.16111145250D+00 -6.33299350739D-08
     9.63675520885D-01  1.83156250000D+02 -2.29242800650D+00 -7.68103423162D-09
     3.28942273195D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.07102096081D-08  5.50000000000D+01
     9.36000000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
19 12 11  5  2  0  0.0 -3.19271814078D-04 -3.52429196937D-12  0.00000000000D+00
     0.00000000000D+00 -2.56562500000D+01  4.39304013061D-09  2.15798610565D+00
    -1.34110450745D-06  8.66045872681D-03  9.43988561630D-06  5.15363704109D+03
     9.36000000000D+04 -1.32247805595D-07  1.21367844688D+00  5.02914190292D-08
     9.60956829538D-01  1.98375000000D+02  2.69247945406D-01 -7.89925760721D-09
     2.13580325049D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.44354999065D-08  2.40000000000D+01
     9.36000000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
28 12 11  5  2  0  0.0  2.01629940420D-04  3.63797880709D-12  0.00000000000D+00
     0.00000000000D+00 -9.02500000000D+01  3.85194616330D-09 -2.64756628469D+00
    -4.60818409920D-06  1.82792989071D-02  6.35161995888D-06  5.15363742638D+03
     9.36000000000D+04  1.89989805222D-07  1.28945528905D-01 -1.86264514923D-09
     9.81203909801D-01  2.64437500000D+02 -1.81516324532D+00 -7.57031533402D-09
    -7.28601777760D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.07102096081D-08  2.80000000000D+01
     9.36000000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
20 12 11  5  2  0  0.0  7.24382698536D-05  1.13686837722D-12  0.00000000000D+00
     0.00000000000D+00  4.71875000000D+00  5.44701260422D-09 -1.36930725337D+00
     3.83704900742D-07  5.21686091088D-03  7.55302608013D-06  5.15360239983D+03
     9.36000000000D+04  4.28408384323D-08 -3.11362311552D+00  1.35973095894D-07
     9.29551660541D-01  2.18937500000D+02  1.27302040506D+00 -8.40499295884D-09
    -2.37152735506D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -7.91624188423D-09  5.30000000000D+01
     9.36000000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
32 12 11  5  2  0  0.0 -5.04904892296D-04 -3.52429196937D-12  0.00000000000D+00
     0.00000000000D+00  8.65625000000D+00  4.96092092798D-09  7.59699554022D-01
     4.06056642532D-07  1.17990599247D-02  7.80634582043D-06  5.15353764153D+03
     9.36000000000D+04  2.34693288803D-07 -2.97808027373D+00 -8.00937414169D-08
     9.50909680794D-01  2.25593750000D+02 -4.77852427077D-01 -8.23141430002D-09
    -2.43581574722D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -3.25962901115D-09  7.90000000000D+01
     9.36000000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
04 12 11  5  2  0  0.0  3.67774628103D-05  1.12549969344D-11  0.00000000000D+00
     0.00000000000D+00  9.84375000000D+00  5.20771692231D-09 -1.48256272694D+00
     6.53788447380D-07  1.03886879515D-02  6.39818608761D-06  5.15367927933D+03
     9.36000000000D+04  2.45869159698D-07  2.16901742229D+00  1.30385160446D-08
     9.37512244840D-01  2.35343750000D+02  9.26380911415D-01 -8.27427322812D-09
     2.88583449230D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -6.51925802231D-09  8.50000000000D+01
     9.36000000000D+04  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
07 12 11  5  4  0  0.0  1.41749158502D-04  4.54747350886D-12  0.00000000000D+00
     0.00000000000D+00  3.06562500000D+01  4.41196949052D-09  1.19459476392D+00
     1.77510082722D-06  6.06159924064D-03  3.93949449062D-06  5.15371810722D+03
     1.00800000000D+05  5.96046447754D-08 -9.52895296577D-01  6.70552253723D-08
     9.75746584608D-01  3.10031250000D+02 -2.94187137915D+00 -8.21677083292D-09
    -7.53602819154D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.11758708954D-08  3.80000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
08 12 11  5  3 59 28.0  2.69711017609D-06  2.27373675443D-13  0.00000000000D+00
     0.00000000000D+00  1.85312500000D+01  3.90730561211D-09  5.61153092039D-01
     1.01514160633D-06  1.27697419375D-02  3.81469726563D-06  5.15363414574D+03
     1.00768000000D+05  1.32247805595D-07 -8.66066952422D-01  1.62050127983D-07
     9.98213825921D-01  3.21250000000D+02 -2.90504334546D+00 -8.03247744207D-09
    -8.46463830045D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -4.19095158577D-09  2.00000000000D+00
     1.00768000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
17 12 11  5  4  0  0.0  8.15289095044D-05 -3.41060513165D-12  0.00000000000D+00
     0.00000000000D+00 -3.98125000000D+01  4.18731727571D-09 -2.27871218105D+00
    -2.00606882572D-06  7.68005801365D-03  9.78820025921D-06  5.15375864792D+03
     1.00800000000D+05 -2.04890966415D-08  1.16105644824D+00 -1.11758708954D-07
     9.63677811815D-01  1.92156250000D+02 -2.29233167334D+00 -7.66674792225D-09
     2.82154610015D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     3.00000000000D+00  0.00000000000D+00 -1.07102096081D-08  5.60000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
01 12 11  5  4  0  0.0  2.80969310552D-04  1.81898940355D-12  0.00000000000D+00
     0.00000000000D+00  1.24062500000D+01  4.67233747875D-09  1.56274466707D+00
     5.71832060814D-07  1.38690031599D-03  7.22706317902D-06  5.15371227837D+03
     1.00800000000D+05 -1.49011611938D-08  2.16947752905D+00  2.04890966415D-08
     9.59939915682D-01  2.38093750000D+02  5.06617343469D-01 -8.09176562595D-09
     1.77150236161D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00  8.38190317154D-09  7.50000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
11 12 11  5  4  0  0.0 -3.19958664477D-04 -3.86535248253D-12  0.00000000000D+00
     0.00000000000D+00  4.43750000000D+00  6.33347810050D-09  1.36875533388D+00
     1.41561031342D-07  1.37237855233D-02  7.62008130550D-06  5.15367959976D+03
     1.00800000000D+05 -3.07336449623D-07  1.88755241687D+00 -8.00937414169D-08
     8.88076756807D-01  1.96375000000D+02  1.11496815188D+00 -9.04716256492D-09
     2.02508435289D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.16415321827D-08  6.80000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
04 12 11  5  4  0  0.0  3.68584878743D-05  1.12549969344D-11  0.00000000000D+00
     0.00000000000D+00  1.17812500000D+01  5.24486132667D-09 -4.32393757496D-01
     8.45640897751D-07  1.03887419682D-02  6.87688589096D-06  5.15367910576D+03
     1.00800000000D+05  1.86264514923D-08  2.16895701255D+00  1.58324837685D-07
     9.37513049444D-01  2.35000000000D+02  9.26394652604D-01 -8.52964100807D-09
     2.53581991279D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -6.51925802231D-09  8.70000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
28 12 11  5  4  0  0.0  2.01656017452D-04  3.63797880709D-12  0.00000000000D+00
     0.00000000000D+00 -6.93750000000D+01  4.04481133977D-09 -1.59750011011D+00
    -3.47383320332D-06  1.82807632955D-02  6.35720789433D-06  5.15363212395D+03
     1.00800000000D+05  3.05473804474D-07  1.28891501877D-01  3.16649675369D-07
     9.81202228908D-01  2.64718750000D+02 -1.81502964286D+00 -7.79496754883D-09
    -2.10365905441D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.07102096081D-08  3.10000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
32 12 11  5  4  0  0.0 -5.04930503666D-04 -3.52429196937D-12  0.00000000000D+00
     0.00000000000D+00  1.33125000000D+01  4.82377235805D-09  1.80990799300D+00
     5.97909092903D-07  1.17991642328D-02  7.81938433647D-06  5.15353895378D+03
     1.00800000000D+05  1.00582838058D-07 -2.97813788344D+00  1.28522515297D-07
     9.50908014530D-01  2.23843750000D+02 -4.77794931471D-01 -8.05390690612D-09
    -3.69301097159D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -3.25962901115D-09  8.00000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
20 12 11  5  4  0  0.0  7.24461860955D-05  1.13686837722D-12  0.00000000000D+00
     0.00000000000D+00  7.15625000000D+00  5.43665502993D-09 -3.18975123211D-01
     3.03611159325D-07  5.21720887627D-03  7.51949846745D-06  5.15360359383D+03
     1.00800000000D+05 -4.28408384323D-08 -3.11368443665D+00  3.16649675369D-08
     9.29549561253D-01  2.22843750000D+02  1.27291915504D+00 -8.52928385034D-09
    -2.05365697163D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -7.91624188423D-09  5.40000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
10 12 11  5  3 59 44.0 -5.74588775635D-05 -1.25055521494D-12  0.00000000000D+00
     0.00000000000D+00  9.40625000000D+00  4.98449333844D-09 -1.30116519633D+00
     5.28991222382D-07  1.18976483354D-02  7.52881169319D-06  5.15372942352D+03
     1.00784000000D+05  2.90572643280D-07 -3.04640315476D+00  5.21540641785D-08
     9.45168902056D-01  2.20656250000D+02  7.69841476450D-01 -8.16355433052D-09
    -2.08580116771D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -2.79396772385D-09  1.50000000000D+01
     1.00784000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
02 12 11  5  4  0  0.0  4.08003572375D-04  1.13686837722D-12  0.00000000000D+00
     0.00000000000D+00  9.90625000000D+00  5.10449833713D-09  2.56434347184D+00
     5.79282641411D-07  1.15479802480D-02  7.23451375961D-06  5.15377860451D+03
     1.00800000000D+05 -1.88127160072D-07  2.15226815868D+00 -1.17346644402D-07
     9.38452109726D-01  2.25843750000D+02 -2.72560346695D+00 -8.17891211309D-09
     2.70368404786D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.76951289177D-08  5.80000000000D+01
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
09 12 11  5  4  0  0.0  1.97479967028D-04  2.72848410532D-12  0.00000000000D+00
     0.00000000000D+00  3.69375000000D+01  4.32303721471D-09  6.27033024225D-01
     1.83098018169D-06  1.72197232023D-02  4.60632145405D-06  5.15365665054D+03
     1.00800000000D+05 -1.97440385818D-07 -9.77472585097D-01 -9.68575477600D-08
     9.84409924030D-01  3.06343750000D+02  1.64824231461D+00 -8.41999358367D-09
     3.10727228751D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     3.00000000000D+00  0.00000000000D+00 -5.58793544769D-09  1.13000000000D+02
     1.00800000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
10 12 11  5  6  0  0.0 -5.74677251279D-05 -1.25055521494D-12  0.00000000000D+00
     0.00000000000D+00  8.25000000000D+00  5.02913805522D-09 -2.48606338080D-01
     5.42029738426D-07  1.18984634755D-02  7.43009150028D-06  5.15373173904D+03
     1.08000000000D+05  2.04890966415D-08 -3.04646305540D+00  1.24797224998D-07
     9.45166049365D-01  2.29218750000D+02  7.69765774829D-01 -8.45070914881D-09
    -2.23580741607D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -2.79396772385D-09  1.60000000000D+01
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
02 12 11  5  6  0  0.0  4.08011488616D-04  1.13686837722D-12  0.00000000000D+00
     0.00000000000D+00  1.30625000000D+01  5.03128100162D-09 -2.66871932232D+00
     6.66826963425D-07  1.15474594058D-02  7.19539821148D-06  5.15377955627D+03
     1.08000000000D+05 -6.51925802231D-08  2.15220959075D+00 -1.39698386192D-07
     9.38454315806D-01  2.28250000000D+02 -2.72560555600D+00 -8.12712424163D-09
     2.84297556420D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.76951289177D-08  5.90000000000D+01
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
17 12 11  5  6  0  0.0  8.15037637949D-05 -3.41060513165D-12  0.00000000000D+00
     0.00000000000D+00 -3.29687500000D+01  4.38982571100D-09 -1.22873984096D+00
    -1.77882611752D-06  7.68092204817D-03  9.30391252041D-06  5.15375612259D+03
     1.08000000000D+05  2.51457095146D-07  1.16100193260D+00  8.00937414169D-08
     9.63678437944D-01  2.03718750000D+02 -2.29217592084D+00 -7.88639992878D-09
     7.78603860547D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     3.00000000000D+00  0.00000000000D+00 -1.07102096081D-08  7.30000000000D+01
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
08 12 11  5  6  0  0.0  2.69804149866D-06  2.27373675443D-13  0.00000000000D+00
     0.00000000000D+00  2.27500000000D+01  3.92909223389D-09  1.61589580317D+00
     1.39325857162D-06  1.27704449696D-02  4.77395951748D-06  5.15363746071D+03
     1.08000000000D+05 -5.58793544769D-09 -8.66124796203D-01  2.58907675743D-07
     9.98213088610D-01  3.09687500000D+02 -2.90491868436D+00 -7.98604693662D-09
    -7.82175437889D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -4.19095158577D-09  9.00000000000D+00
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
04 12 11  5  6  0  0.0  3.69395129383D-05  1.12549969344D-11  0.00000000000D+00
     0.00000000000D+00  1.58750000000D+01  5.17200114889D-09  6.17773516427D-01
     1.04121863842D-06  1.03894958738D-02  7.14696943760D-06  5.15368058777D+03
     1.08000000000D+05 -2.10478901863D-07  2.16889682372D+00  9.12696123123D-08
     9.37514665969D-01  2.31531250000D+02  9.26409012608D-01 -8.43249410437D-09
     8.75036448780D-11  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -6.51925802231D-09  8.80000000000D+01
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
28 12 11  5  6  0  0.0  2.01682094485D-04  3.63797880709D-12  0.00000000000D+00
     0.00000000000D+00 -6.00000000000D+01  4.22981904608D-09 -5.47227892295D-01
    -3.08081507683D-06  1.82817617897D-02  6.18025660515D-06  5.15363513565D+03
     1.08000000000D+05 -1.93715095520D-07  1.28832647220D-01  1.04308128357D-07
     9.81198554058D-01  2.72875000000D+02 -1.81509841610D+00 -8.33070415012D-09
    -1.37505727665D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.07102096081D-08  4.00000000000D+01
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
20 12 11  5  6  0  0.0  7.24541023374D-05  1.13686837722D-12  0.00000000000D+00
     0.00000000000D+00  7.93750000000D+00  5.36022327481D-09  7.31268501868D-01
     2.57045030594D-07  5.21715916693D-03  7.18608498573D-06  5.15360347176D+03
     1.08000000000D+05 -4.65661287308D-08 -3.11374573585D+00  1.86264514923D-09
     9.29548016412D-01  2.23437500000D+02  1.27290549577D+00 -8.52821237713D-09
    -2.26438003480D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -7.91624188423D-09  9.20000000000D+01
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
09 12 11  5  6  0  0.0  1.97499990463D-04  2.72848410532D-12  0.00000000000D+00
     0.00000000000D+00  3.70937500000D+01  4.14660129402D-09  1.67719764149D+00
     1.58883631229D-06  1.72197560314D-02  4.11644577980D-06  5.15365742111D+03
     1.08000000000D+05  2.16066837311D-07 -9.77532493055D-01 -3.94880771637D-07
     9.84412612874D-01  3.09937500000D+02  1.64826748265D+00 -7.94104506212D-09
     1.26790995640D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     3.00000000000D+00  0.00000000000D+00 -5.58793544769D-09  1.41000000000D+02
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
05 12 11  5  6  0  0.0 -3.64434439689D-04 -2.61479726760D-12  0.00000000000D+00
     0.00000000000D+00  5.68750000000D+00  4.81305762603D-09 -5.41120460936D-01
     3.78116965294D-07  3.03652207367D-03  7.62753188610D-06  5.15377803230D+03
     1.08000000000D+05  4.28408384323D-08 -3.06707891244D+00  1.11758708954D-08
     9.50308418537D-01  2.22968750000D+02  2.54467043917D-01 -8.16534011919D-09
    -2.21437795202D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.02445483208D-08  4.80000000000D+01
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00
12 12 11  5  6  0  0.0  9.99365001917D-05  2.04636307899D-12  0.00000000000D+00
     0.00000000000D+00 -7.67187500000D+01  4.10909973193D-09  1.92499591043D+00
    -3.97488474846D-06  4.01161285117D-03  6.05732202530D-06  5.15368676567D+03
     1.08000000000D+05 -5.58793544769D-08  1.04069837281D-01  1.37835741043D-07
     9.82037335694D-01  2.72625000000D+02  1.43975287389D-01 -7.97533220460D-09
    -1.93579491934D-10  0.00000000000D+00  1.71300000000D+03  0.00000000000D+00
     2.00000000000D+00  0.00000000000D+00 -1.21071934700D-08  6.00000000000D+01
     1.08000000000D+05  0.00000000000D+00  0.00000000000D+00  0.00000000000D+00