|||
首先获取下列文件,附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
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-12-22 09:31
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社