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

博文

测光-作光变曲线图总结

已有 6430 次阅读 2009-11-7 22:48 |个人分类:IRAF测光|系统分类:科研笔记| 数据, 光变曲线, 作图

根据已知的数据画光变曲线的总结

第一阶段:

1、把相应要用到的五个数据源(expt、flist、res.dat、utime、coord.dat)整理出    
    来;
 
2、修改coord.dat。根据XX网络已知星的参数,选取与目标星参数比较接近的星作比较星,删除coord.dat文件前面
    的信息,将比较星提前(可选9颗提前)。

3、修改expt。根据程序中读取的数值,需更改格式如(EXPTIME =       20.000000这之间的空格应与程序的读出位
    置一致)

4、修改res.dat。根据程序中读取的数值,需更改格式。如
  (QXAnd_B_0000.fit              284.282  186.516  14.493  0.009)

第二阶段:

1、建立相应的数据存放文件,将an文件下的(sofe pag)中的fortran程序拷入其中;

2、得出的五个数据源拷入刚建立的目录中;

3、修改程序,编译读出数据(主要用到fmtlst.f 和 getlcc.f,要会读其中的变量,学会改数据)。

4、编译plot.f。应与pgplot建立起相应的环境变量操作如下:

sudo gfortran -o plot plot.f -L/usr/local/pgplot -L/usr/local/X11R6/lib -lpgplot -lX11

;编译plot.f文件,使其与pgplot和X11R6建立关联  

6、运行下列程序画出光变曲线:

   export PGPLOT_DIR=/usr/local/pgplot;导出的目录
   export PGPLOT_DEV=/Xserve;显示的设备
   ./plot;在运行另一文件夹中编译好的plot文件

 
7、根据画出的光变曲线图选取较好的比较星作光变曲线图。此时会产生一个tmp临时文件,用于后面kw程序执行。

8、执行 ./kw 程序,取出minimum times。会让输入文件 (tmp),然后有问是否擦除bad point,(一般选NO),用鼠标选取最小值的范围,会等到某一范围内的极小值。依次对第一天的数据这样做,将得到的极小值放入一个mini的文件中。

9、编译epoch,运行 ./epo 程序,会让输入数据文件(将mini输入),还会输入周期p,T0。最后产生一个tmp临时文件。

10、编译oc1.f,将tmp和oc拷入pgplot目录下,运行oc会得到EPOCH-L(o-c)图。


                                                                                  2009-11-24于河北兴隆[2013-11-14修改于BNU]



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

上一篇:pgplot与f77 or gfortran的关联
下一篇:详论fortran格式化输出
收藏 IP: 114.255.218.*| 热度|

0

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

数据加载中...

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

GMT+8, 2024-11-21 23:15

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部