[daip] GAL
Jessica Rosenberg
jlrosenberg at cfa.harvard.edu
Mon Jul 10 15:41:24 EDT 2006
Hello,
I am trying to run GAL and am experiencing some trouble... I am running
the routine on a blanked cube to get the rotation curve. The problem is
that I am getting an error message that I don't understand and the
routine does not appear to be doing any fitting (doesn't seem to even
be going through an iterations) so I am hoping you can offer a
suggestion about what is going wrong (I have tried this over a range of
radii and using 'ex' as the function as well as 'cc' and the results
are the same). Here is the information that is returned by the routine:
>go
AIPS 2: Waiting for returned adverbs
GAL 2: Task GAL (release of 31DEC06) begins
GAL 2: # pixels: 70747, # non-blanked : 0
GAL 2: cc curve, Rmin= 20 Rmax= 50 arcsec
GAL 2: COSmin = -1.00 COSmax = 1.00
GAL 2: SINmin = -1.00 SINmax = 1.00
GAL 2: *** no special weighting, all weights unity ***
GAL 2: GAL : information parameter # 0
GAL 2: after 0 iterations the solution vector is
GAL 2: 254.00 254.00 162.00 54.00 1225.00 80.00 15.00 1.00
GAL 2: 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
GAL 2: R.A. : 11 26 47.00 Dec : 59 08 50.00
GAL 2: 0.00 0.00
GAL 2: Create IC691_3_B5 .GAL . 14 (MA) on disk 1 cno 30
GAL 2: ZTXOP2: using translated file name =
GAL 2: ZTXOP2: /home/jlrosenb/IC691/newincl_rotcur.txt
GAL 2: X, Y RANGES 0.0000E+00 0.0000E+00 IMPROPER
GAL 2: PLAXES: IERR = -1
GAL 2: ZCLOSE: LUN = 16 ALREADY CLOSED IN FTAB ON SEARCH
GAL 2: returns adverbs to AIPS
GAL 2: Appears to have ended successfully
GAL 2: rglinux5 31DEC06 TST: Cpu= 0.1 Real= 0
AIPS 2: Resumes
The parameters are set as follows:
AIPS 2: GAL Task to determine parameters from a velocity field.
AIPS 2: Adverbs Values Comments
AIPS 2: ----------------------------------------------------------------
AIPS 2: USERID 0 User ID. 0=>current user
AIPS 2: 32000=>all users
AIPS 2: INNAME 'IC691_3_B5' Image name
AIPS 2: INCLASS 'BLANK' Image class
AIPS 2: INSEQ 1 Image seq. #
AIPS 2: INDISK 1 Disk drive #
AIPS 2: IN2NAME ' ' Weight image name
AIPS 2: IN2CLASS ' ' Weight image class
AIPS 2: IN2SEQ 0 Weight image seq. #
AIPS 2: IN2DISK 0 Weight image disk drive #
AIPS 2: OUTNAME ' ' Output image name (name)
AIPS 2: OUTCLASS ' ' Output image name (class)
AIPS 2: OUTSEQ 0 Output image name (seq. #)
AIPS 2: OUTDISK 1 Output disk drive #
AIPS 2: BLC *all 0 Bottom left corner of image
AIPS 2: 0=>entire image
AIPS 2: TRC 263 269 Top right corner of image
AIPS 2: 38 *rest 0 0=>entire image
AIPS 2: SOLCON 0.001 Convergence criterion.
AIPS 2: 0=> 0.001
AIPS 2: FUNCTYPE 'cc' Type of rotation curve: 'BR':
AIPS 2: Brandt, 'EX' : exponential,
AIPS 2: 'CC' : constant, 'SB' : solid
AIPS 2: body.
AIPS 2: APARM 254 254 First guesses to parameters:
AIPS 2: 162 54 1: x - coordinate centre, 2:
AIPS 2: 1225 80 y - coordinate centre, 3: po-
AIPS 2: 15 1 sition angle major axis, 4:
AIPS 2: 0 0 inclination, 5: systemic ve-
AIPS 2: locity, 6: maximum velocity,
AIPS 2: 7: radius maximum velocity
AIPS 2: ("), 8: exponent (0=>1).
AIPS 2: CPARM 2 1 (1) : residual map; 0 : none,
AIPS 2: *rest 0 1 : uses whole field, 2 :
AIPS 2: uses selected field, 3: uses
AIPS 2: values from INFILE. (2) :
AIPS 2: 0/1 : plot (no/yes). (3) :
AIPS 2: fixed parameters. (4) : 0/1:
AIPS 2: use IN2C map for weighting
AIPS 2: (no/yes) 5: 0/1: copy old
AIPS 2: history (no/yes) 6: 0: write
AIPS 2: one line to output;-1/1 write
AIPS 2: position to output file.
AIPS 2: DPARM 20 50 1 : lower limit radius; 2 :
AIPS 2: *rest 0 upper limit radius; 3 : lower
AIPS 2: limit cos (pos. angle in ga-
AIPS 2: laxy); 4 : upper limit; 5 and
AIPS 2: 6 : idem for sine; 7 : max.
AIPS 2: plotted radius; 8: max. plot-
AIPS 2: ted rot. velocity; 9 :0 (def)
AIPS 2: => fit parameters, 1=> just
AIPS 2: plot using inputs. 10: width
AIPS 2: annuli (in pixel separation).
AIPS 2: INFILE ' '
AIPS 2: Text file containing rotation
AIPS 2: curve information.
AIPS 2: OUTFILE '/home/jlrosenb/IC691/newincl_rotcur.txt'
AIPS 2: Text file to which one line
AIPS 2: of output results is written.
AIPS 2: PIXSTD 0 Estimated rms uncertainty in
AIPS 2: the observed radial velocity
AIPS 2: (km/s) at one pixel. 0=>10
AIPS 2: FACTOR 0 Scale plot symbol by FACTOR
AIPS 2: DOTV 1 > 0 Do plot on the TV, else
AIPS 2: make a plot file
AIPS 2: GRCHAN 1 Graphics channel 0 => 1.
The output file gives:
20.0 50.0 254.0 254.0 162.0 54.0 1225.0 80.0 0.0 0.0
Thanks!
Jessica
---------------------------------------------------------------------
Jessica Rosenberg
Harvard-Smithsonian Center for Astrophysics
60 Garden Street, MS 65
Cambridge, MA 02138-1516
Phone: 617-496-7740
E-mail: jlrosenberg at cfa.harvard.edu
---------------------------------------------------------------------
More information about the Daip
mailing list