斗室望远分享 http://blog.sciencenet.cn/u/slrseer 乱涂乱画,不正不经

博文

生成JPL/DE421星历的方法

已有 6006 次阅读 2020-3-11 14:34 |个人分类:职业病|系统分类:科研笔记| JPL, 行星星历, DE421

首先获取下列文件,附md5sum

5d16eed1d78d00d40b1bb06e7f3c5fa8  asc2eph.f
1967a44569ce50710c7403afe150fc2d  ascp1900.421
041e82d6fdd5ba971d74968b6af25078  ascp2050.421
c05f534ca15326e745bc310649529bc9  header.421
da6e477366358b8007a26dec5b2c14a8  README.txt
00bf84539ba3c65a5ac5ab3eaab85b24  testeph1.f
84cf0bb3bbe830500e65f029a6b4db9d  testeph.f
0cb939eedc498d815f748491a17e4250  testpo.421
13b5be47f5bfdfe8404b5547ea1b28fc  userguide.txt


获取方法是JPL的FTP。

首先阅读README.txt和userguide.txt,做以下动作。

动作分为几个步骤:
1,转换。把文本方式的系数表ascp1900.421、ascp2050.421,带头文件header.421转换为二进制数据文件JPLEPH。
2,测试。使用testeph.f和testpo.421测试生成的文件,验证数据正确性。


在类unix平台上(我用的Ubuntu 13.04),
1,转换。
修改asc2eph.f 中的74行,取消注释,使得PARAMETER ( NRECL = 4 )生效。
将文件保存为asc2eph-nrecl_4.f
编译:gfortran asc2eph-nrecl_4.f -o asc2eph-nrecl_4.exe
运行转换:cat header.421 ascp1900.421 ascp2050.421 | ./asc2eph-nrecl_4.exe
此时,输出头文件信息,可以看到 KSIZE =  2036 字样。
运行结束后,同目录下就生成JPLEPH二进制数据文件。
2,测试。
按照testeph.f中的注释,一处处修改代码,取消注释,使某些语句生效。主要有:输入文件名JPLEPH,NRECL=4,KSIZE=2036,CALL FSIZER3,等等。
编译:gfortran testeph-ksize_2036-fsizer3.f -o testeph-ksize_2036-fsizer3.exe
运行测试:cat testpo.405 | ./testeph-ksize_2036-fsizer3.exe > testeph.out.txt

在testeph.out.txt中查看difference非常小,于是测试通过。

【完】


附录:所有文件的列表,及MD5校验码。

5d16eed1d78d00d40b1bb06e7f3c5fa8  asc2eph.f
4225b85c900097a4d53c697b3c451441  asc2eph-nrecl_4.exe
24df47b6cc8f8b40b6c5d73f718c1d04 
asc2eph-nrecl_4.f

1967a44569ce50710c7403afe150fc2d  ascp1900.421
041e82d6fdd5ba971d74968b6af25078  ascp2050.421
c05f534ca15326e745bc310649529bc9  header.421
0b034c16a1295d968542075df837b9bb  JPLEPH
da6e477366358b8007a26dec5b2c14a8 
README.txt 

00bf84539ba3c65a5ac5ab3eaab85b24  testeph1.f
84cf0bb3bbe830500e65f029a6b4db9d  testeph.f
4c40c38d4727fc6e7efadd3fe4e58955  testeph-ksize_2036-fsizer3.exe
3fd9e3e25225c17cfcc396b8b5e4fba8 
testeph-ksize_2036-fsizer3.f

284914c22227649b05f89281525e016b  testeph.out.txt

0cb939eedc498d815f748491a17e4250  testpo.421
13b5be47f5bfdfe8404b5547ea1b28fc 
userguide.txt




https://blog.sciencenet.cn/blog-47397-1222961.html

上一篇:ITRF用什么椭球?
下一篇:MLRS_Lunar_Code v1.0 测试笔记
收藏 IP: 210.72.87.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-12-22 09:31

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部