我需要知道FINDCLUMPS输出的云块的线宽是怎么定义的。于是查看findclumps.c中的语句
cupidStoreClumps( "OUTCAT", xloc, ndfs, nsig, deconv, backoff, ishape,
                        velax, beamcorr, "Output from CUPID:FINDCLUMPS",
                        usewcs, gotwcs ? iwcs : NULL, dataunits,
                        confgrp, logfile, &nclumps, status );
于是查找cupidstoreclumps.c，其中提到
cpars = cupidClumpDesc( indf, deconv, wcsmap, wcsfrm, dataunits,
                                 beamcorr, backoff, stccol, velax, cpars,
                                 &names, &units, &ncpar, &ok, &stcptr,
                                 &region, status );
于是再查找cupidclumpdesc.c，其中关于宽度有
v = v0 - beamcorr[ 1 ]*beamcorr[ 1 ]/5.5451774;
其中beamcorr[ 1 ]乃是波束的半高全宽，5.5451774就是8ln2。所以从这里看，最后输出的是线宽的方均根值。于是我如果想换算成半高全宽，还得乘以。

转载本文请联系原作者获取授权，同时请注明本文来自钱磊科学网博客。
链接地址：https://m.sciencenet.cn/blog-117333-378395.html

上一篇：编程点滴（二）数据怎么存？
下一篇：用IDL从Fits文件中读二进制表