[daip] SAD usage

[Eric Greisen]

On 06/22/2016 02:43 PM, Stockdale, Christopher wrote:
> Eric,
>
> I’ve got things up and working well, so now I have a technical question
> or two about SAD/TVSAD.  I’ve got a lot of diffuse emission in the
> spiral arms of M83 and SAD is doing a fair job of getting brighter
> sources, although missing a few and doing some odd very large grouping.
> I could send a screen capture of the TV with the star file outputs to
> show, but I’d like to figure out how to reduce some of the extraneous
> boxes and try to get some autofits.
>
> On a related note, is it possible to run jmfit or imfit on a text file
> with source positions in RA and DEC on all of our sources that we detect
> in the full band image?  I’ve imaged the C band and L band fully and
> partially and I’d like to use my source list for the full image to batch
> the imfit/jmfit for the subband images.
>

I wrote a small proc that seems to work - RUN CSPROC to compile the 
procedure.  Then set up the adverbs for JMFIT or IMFIT (some 
instructions in CSPROC.001).  Edit the file CSSCRIPT.001 for your 
coordinates and the number you want to fit (4 in the example).  Then
RIN CSSCRIPT (note VERSION must point to the directory containing CSPROC 
and CSSCRIPT).

Eric


-------------- next part --------------
$   proc to read script file for coordinates, do JMFIT
$   all adverbs for JMFIT/IMFIT including TASK must be set
$   in advance:
$      NGAUSS=1 OR 2 IF ctype(2)>1
$      GMAX=GPOS=GWIDTH=0
$      select OUTVERS to save results in MF file
PROCEDURE  COFIT(II)
imsize=11
dowait=1
for i = 1:ii; read coordina; print ii, coordina;
   cowindow; go; print fmax; print fpos; print fwidth
   end
dowait=-1
finish
-------------- next part --------------
$
COFIT(4)
17 48 32.0  70 05 51
17 51 14.5  69 48 50
17 51 29.7  69 54 27
17 51 35.7  69 58 29
$