ReadMe
The current IRAF version is poorly optimized. The script itself was made as a prototype to verify whether the skewness method can work in real applications or not.
At the moment, I recommend you to use the IDL version posted on this site.
<< Issues of the IRAF version >>
1. use imcalc.cl
- contisub.cl has to run under "stsdas"
- contisub.cl writes each residual image to HDD
- "_tempskew.fits" is written and deleted for each process. This is a huge drawback of this iraf-based routine.
Though this can work for most cases, this writing to HDD should be fixed by some persons familiar
with IRAF programming.
2. iterstatnew.cl
- iterstat.cl has a clipping procedure to remove anomalies, such as cosmic rays and saturated stars. But sometimes it does not work great. Best way is to choose a section to avoid anomalies. Again, please use the IDL version until a better GUI treatment in iraf is available.
<<Script details>>
The script name is `contisub.cl'. This was made based on `msky.cl' and `iterstat.cl'. The contisub.cl was written by modifying msky.cl. And I added some lines to iterstat.cl to call skewness value too and named it iterstatnew.cl; the original iterstat.cl does not have skewness calls.
msky.cl needs iterstat.cl. My contisub.cl needs msky.cl and iterstatnew.cl. And implicitly I used STSDAS/imcalc.
**If you have any questions, please contact me : shongScience(at)gmail(dot)com.