学 点 新 知--做 点 事 情分享 http://blog.sciencenet.cn/u/dabing 平时多记点,东西就多了...

博文

利用IRAF计算日心儒略日的方法

已有 5276 次阅读 2013-12-4 19:05 |个人分类:IRAF测光|系统分类:科研笔记

具体处理过程:

1.进入IRAF setjd设置:ecl>noao>astutil>epar setjd

注: (observa=     )输入观测站的名称,名称包含在obsdb.dat中,在终端查找

astro@astro-desktop:~$ locate obsdb.dat

/iraf/iraf/noao/lib/obsdb.dat

astro@astro-desktop:~$ sudo gedit /iraf/iraf/noao/lib/obsdb.dat    (字符编码选择: 西欧)

打开可以看到如下的信息(列出部分),查看你的观测站名称是否包含其中,如果包含可以直接在(observa=    )中输入简称,没有可以手动添加,下面的第一条‘konkoly’的信息就是我手动加入的,然后保存。

---------------------------------------------------------------------

# Observatory Parameters.
# Some were taken from the Almanac and others have been submitted.
# Observatories wishing to be added or make changes in the default
# distributed database should send information to iraf@noao.edu.

observatory = "konkoly"
   name = "Hungarian Academy of Sciences, Konkoly Observatory"
   longitude = -19:00.7            #观测站赤经
   latitude = 46:10.8              #观测站赤纬
   altitude = 160.                   #观测站海拔高度
   timezone = -1                     #观测站所处时区(注:这里是匈牙利时区 +1,输入为 timezone= -1)

#timezone只是当地时间和世界时间的转换,如果前面的utday 和uttime为UT时间(yes),timezone不起转换的作用

# Submitted by Biwei JIANG 11/28/95

observatory = "BAO"

   name = "Beijing XingLong Observatory"

   longitude = 242:25.5

   latitude = 40:23.6

   altitude = 950.

   timezone = -8                     #观测站所处时区(注:这里是中国时区 +8,输入为 timezone= -8)


# Submitted by Pat van Heerden 11/20/06
observatory = "boyden"
   name = Boyden Observatory
   longitude = 332:35:40.0
   latitude = -29:02:20
   altitude = 1387
   timezone = -2
------------------------------------------------------------------------------

2、运行计算

astutil> setjd kasc*.fits > time.dat    #将数据输入time.dat文件中

打开time.dat得到如下结果:

#              Image             jd            hjd       ljd
# SETJD: Observatory parameters for Hungarian Academy of Sciences, Konkoly Observatory
#    timezone = -1
kasc-f2_V_0001.fits  2454522.27453  2454522.27606   2454522
kasc-f2_V_0002.fits  2454522.27614  2454522.27767   2454522
kasc-f2_V_0003.fits  2454522.27775  2454522.27929   2454522

 .....

参考IDL计算程序: http://idlastro.gsfc.nasa.gov/ftp/pro/astro/observatory.pro

检验计算的HJD: http://www.physics.sfasu.edu/astro/javascript/hjd.html

                                                        2013-12-05于BNU



https://blog.sciencenet.cn/blog-287062-747053.html

上一篇:基础物理-高能物理[Hyperphysics]
下一篇:Eclipsing Binary Stars
收藏 IP: 210.31.76.*| 热度|

0

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

数据加载中...

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

GMT+8, 2024-10-19 23:00

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部