(Bernese GNSS 软件5.2版本)制作并安装JPL星历(行星和月球)文件
-----------------------------------------
来自JPL的行星和月球星历用于计算来自太阳、月亮和其他行星的万有引力,进而用于确定卫星轨道。
-----------------------------------------
星历文件来源
JPL星历由JPL的M.E. Standish等人完成,具体信息可在
http://ssd.jpl.nasa.gov?planet_eph_export上查看。(注:打不开的话,如果能上谷歌,可以在http://1t.click/aKXp查看内容)
星历文件可以在网址http://www.willbell.com/software/jpl.htm上的“Download the ephemerides”下载(注:没有找到)。这个网站上也提供含二进制和文件格式的JPL星历文件CD光盘售卖,价格24.95美元(注:人民币约176元)。
-----------------------------------------
下载星历
星历文件可以用匿名从JPL的ftp服务器ssd.jpl.nasa.gov上下载,具体在
pub/eph/planets/ascii/de405目录下(注:正常打不开,可以从https://github.com/SGL-UT/GPSTk/tree/master/ext/apps/geomatics/JPLeph/JPL上下载)。
下载文件如下:
header.405
testpo.405
ascp1980.405 : 用于1980 and 2000年间的
ascp2000.405 : 用于2000 and 2020年间的
ascp2020.405 : 用于2020 and 2040年间的
-----------------------------------------
制作步骤
安装JPL星历需要Fortran程序,这个程序已经在Bernese GNSS软件中了,这个程序原本是JPL的程序,这里进行了调整便于使用。
需要用“ASC2EPH”程序把下载的文本格式星历转换为二进制格式。
JPL文本格式的星历每20年一个。
注意:请确认Bernese软件的环境变量(例如 ${XG})可以被程序找到。
——LINUX 平台:
确认环境变量已经加载。
依赖与你使用的shell
“. BERN52/GPS/LOADGPS.setvar”或者
“source BERN52/GPS/LOADGPS.setvar”
(如果有必要把路径加入到LOADGPS.setvar中)
可以加入到登录脚本中去。
——Windows平台
在安装完成Bernese软件后需要重启(win9x系统)或者重新登录即可,否则可能提示找不到命令。
1)合并文本文件
合并头文件和文本星历文件成为一个临时文件:
——Linux系统:
cat header.405 > temp.405
cat ascp1980.405 >> temp.405
cat ascp2000.405 >> temp.405
cat ascp2020.405 >> temp.405
——Windows系统:
copy header.405+ascp1980.405+ascp2000.405+ascp2020.405 temp.405
2)把文本文件转换为二进制
准备一个输入文件ASC2EPH.INP,包括:
——第一行:文本文件的文件名
——第二行:结果二进制文件的文件名(例如:JPLEPH)
例如:
echo temp.405 >ASC2EPH.INP
echo JPLEPH >>ASC2EPH.INP
在“temp.405”文件夹中运行程序:
——Linux系统:使用命令 ${XG}/ASC2EPH
——Windows系统:使用命令 %XG%\ASC2EPH
如果提示“ok FORTRAN STOP”则程序运行正常结束,生成了JPLEPH文件将包含二进制的DE405星历。
另外,可以运行以上2步骤,分别的把ascp????.405进行转换,然后再用程序“BINMERGE”把二进制文件进行合并。
BINMERGE.INP文件示例:
temp1.405
temp2.405
JPLEPH
DE405
3)检查生成的文件可读性和与“testpo.405”比较数值的正确性
——准备输入文件TESTEPH.INP,内容包括:
——第一行:准备比较的二进制文件的文件名(例如:JPLEPH)
——第二行:检核文件名(testpo.405)
——第三行:星历的DE号(DE405)
TESTEPH.INP文件样例如下:
JPLEPH
testpo.405
DE405
——运行程序
Linux系统:使用命令 ${XG}/TESTEPH
Windows系统:使用命令 %XG%\TESTEPH
正常情况下,应该没有警告和错误提示。
“TESTEPH”程序会给出新生成的JPLEPH星历文件的时间覆盖范围,时间格式儒略日(2444208.50 to 2466160.50)。
最后面部分最后一列给出与检核文件的差值。
输出示例:
line - jed - ... --- jpl value --- --- user value --- -- diff. --
100 2447221.5 ... -0.4321082450831 -0.4321082450831 0.36138E-13
200 2450265.5 ... 2269.9483911098296 2269.9483911098296 -0.00000E+00
300 2453310.5 ... -0.0104722325301 -0.0104722325301 0.15380E-13
400 2456352.5 ... 0.1312929459907 0.1312929459907 0.30143E-13
-----------------------------------------
JPL星历文件归位
Bernese GNSS软件中的JPL星历文件
——位于通用文件目录:${X}/GEN (Linux系统)或者 %X%\GEN (Windows 系统)
——具备.EPH后缀
——默认文件名DE405.EPH
把新生成的星历文件放到Bernese GNSS软件中,很容易的把“JPLEPH”文件复制为DE405.EPH文件同时放到指定目录(Linux系统:${X}/GEN (UNIX) ,Windows系统:%X%\GEN )即可:
- Linux系统: cp JPLEPH $X/GEN/DE405.EPH
- Windows系统: copy JPLEPH %X%\GEN\DE405.EPH
制作好的文件DE405.EPH下载地址:
https://download.csdn.net/download/zzh_my/12177638
https://download.csdn.net/download/zzh_my/12172256(不可用)