|||
根据已知的数据画光变曲线的总结
第一阶段:
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]
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-21 23:15
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社