[daip] DOFARS problems
Eric Greisen
egreisen at nrao.edu
Sun Oct 16 16:25:50 EDT 2011
Frazer Owen wrote:
> DOFARS seems to have problems with keeping track of where it puts the
> weight file.
>
> Script was run on pilabo under user # 2256.
>
> AIPS 4: DOFARS: Procedure to create Faraday rotation measure images
> AIPS 4: Adverbs Values Comments
> AIPS 4: ----------------------------------------------------------------
> AIPS 4: INNAME 'A2256ALL' First image name
> AIPS 4: INCLASS 'QFLATC' First image class
> AIPS 4: INSEQ 2000 First image seq. #
> AIPS 4: INDISK 7 First image disk drive #
> AIPS 4: IN2NAME 'A2256ALL' Second image name
> AIPS 4: IN2CLASS 'QFLATC' Second image class
> AIPS 4: IN2SEQ 2000 Second image seq. #
> AIPS 4: IN2DISK 7 Second image disk drive #
> AIPS 4: OUTNAME 'A2256C55W' First output image name
> AIPS 4: OUTCLASS 'RMTFA' First output image class
> AIPS 4: OUTSEQ 14 First output image seq. #
> AIPS 4: OUTDISK 7 First output image disk #
> AIPS 4: OUT2NAME 'A2256C55W' Second output image
> AIPS 4: OUT2CLAS 'RMTFP' Second output image class
> AIPS 4: OUT2SEQ 14 Second output image seq #
> AIPS 4: OUT2DISK 7 Second output image disk #
> AIPS 4: INFILE *all ' ' Input file of weights
> AIPS 4: DOALIGN 1 Should images be coincident?
> AIPS 4: (See HELP.)
> AIPS 4: BLC 900 900 Bottom left corner
> AIPS 4: 14 *rest 0
> AIPS 4: TRC 1100 1100 Top right corner
> AIPS 4: *rest 0 See help for the 1st axis
> AIPS 4: (must be frequencies!!!)
> AIPS 4: APARM 200 2 Parameters for algorithm:
> AIPS 4: 1 1 1 number of pixels at
> AIPS 4: 2 0 half of the Fourier
> AIPS 4: 0 0 transform output
> AIPS 4: 100 0 The whole number is
> AIPS 4: 2*APARM(1)+1
> AIPS 4: 2 cell size in 1/m^2
> AIPS 4: 0 =>
> AIPS 4: PI/(4*(Lmax^2-Lmin^2))
> AIPS 4: Lmax,Lmin-max,min lambda
> AIPS 4: at the data
> AIPS 4: 3 0 => regular output
> AIPS 4: 1 => output is RMTF
> AIPS 4: 4 0=> CLEANed Fourier
> AIPS 4: transform
> AIPS 4: 1=> unCLEANed Fourier
> AIPS 4: transform
> AIPS 4: 5 0=>original(shifted back)
> AIPS 4: RE/IM are sent out
> AIPS 4: 1=>the shifted RE/IM are
> AIPS 4: sent out
> AIPS 4: the data are sent out
> AIPS 4: 6 is not used
> AIPS 4: 7 0=> convolve the clean
> AIPS 4: components
> AIPS 4: 1=> no convolve
> AIPS 4: 8 0=> use the Gaussian as
> AIPS 4: the convolve function
> AIPS 4: 1=> use the Re of RMTF as
> AIPS 4: the convolve function
> AIPS 4: 9 full width of Gaussian
> AIPS 4: convolve function, at 0.5
> AIPS 4: level, in 1/m^2
> AIPS 4: 0 => 1
> AIPS 4: 10 what send to output?
> AIPS 4: 0 => sum of CLEAN and
> AIPS 4: residual
> AIPS 4: 1 => CLEAN result
> AIPS 4: 2 => residual
> AIPS 4: GAIN 0.1 Gain in the CLEAN
> AIPS 4: NITER 0 Maximum number of clean
> AIPS 4: components
> AIPS 4: FLUX 8.000E-05 Minimum flux of clean
> AIPS 4: component (Jy)
> AIPS 4: DOHIST -2 -2 => copy 1st HI only
> AIPS 4:
> go dofars
> AIPS 4: 'FITS:DOFARS.WEIGHTS' DOES NOT EXIST
> RSPEC4: Task RSPEC (release of 31DEC11) begins
> RSPEC4: Created plot file version 3
> RSPEC4: GFINIS: number records used 6
> RSPEC4: Appears to have ended successfully
> RSPEC4: pilabo 31DEC11 TST: Cpu= 6.3 Real= 7 IO= 1203
> AIPS 4: Resumes
> TRANS4: Task TRANS (release of 31DEC11) begins
> TRANS4: Create A2256ALL .TRANS1. 1 (MA) on disk 7 cno 2289
> TRANS4: Copied CC file from vol/cno/vers 7 2189 1 to 7 2289 1
> TRANS4: Copied CC file from vol/cno/vers 7 2189 2 to 7 2289 2
> ...
> TRANS4: Copied CC file from vol/cno/vers 7 2189 75 to 7 2289 75
> TRANS4: Copied CC file from vol/cno/vers 7 2189 76 to 7 2289 76
> TRANS4: Copied CC file from vol/cno/vers 7 2189 77 to 7 2289 77
> TRANS4: Begin row swap with subimaging
> TRANS4: Begin in-core x-y transpose
> TRANS4: Appears to have ended successfully
> TRANS4: pilabo 31DEC11 TST: Cpu= 0.6 Real= 13 IO= 185
> AIPS 4: Resumes
> TRANS4: Task TRANS (release of 31DEC11) begins
> TRANS4: Create A2256ALL .TRANS2. 1 (MA) on disk 7 cno 2290
> TRANS4: Copied CC file from vol/cno/vers 7 2189 1 to 7 2290 1
> TRANS4: Copied CC file from vol/cno/vers 7 2189 2 to 7 2290 2
> TRANS4: Copied CC file from vol/cno/vers 7 2189 3 to 7 2290 3
> TRANS4: Copied CC file from vol/cno/vers 7 2189 4 to 7 2290 4
> ...
> TRANS4: Copied CC file from vol/cno/vers 7 2189 75 to 7 2290 75
> TRANS4: Copied CC file from vol/cno/vers 7 2189 76 to 7 2290 76
> TRANS4: Copied CC file from vol/cno/vers 7 2189 77 to 7 2290 77
> TRANS4: Begin row swap with subimaging
> TRANS4: Begin in-core x-y transpose
> TRANS4: Appears to have ended successfully
> TRANS4: pilabo 31DEC11 TST: Cpu= 0.8 Real= 1 IO= 185
> AIPS 4: Resumes
> FARS 4: Task FARS (release of 31DEC11) begins
> FARS 4: ZTXOP2: using translated file name =
> FARS 4: ZTXOP2: /DATA/PILABO_1/FITS/DOFARS.WEIGHTS
> FARS 4: ZTXOP2: FILE TO BE READ DOES NOT EXIST
> FARS 4: ERROR
> FARS 4: MINIT: FILE FOR LUN 21 NOT OPEN
> FARS 4: COULD NOT INITIALIZE OUTPUT FILE 3 1 IER= 1
> FARS 4: Purports to die of UNNATURAL causes
> FARS 4: pilabo 31DEC11 TST: Cpu= 0.0 Real= 0 IO= 1
> AIPS 4: Resumes
> AIPS 4: RETURN CODE 16 RECEIVED: STOPPING
> inp rspec
> AIPS 4: RSPEC : Task to plot spectrum of the rms of an image
> AIPS 4: Adverbs Values Comments
> AIPS 4: ----------------------------------------------------------------
> AIPS 4: USERID 0 User ID. ignored
> AIPS 4: INNAME 'A2256ALL' Image name (name)
> AIPS 4: INCLASS 'TRANS1' Image name (class)
> AIPS 4: INSEQ 0 Image name (seq. #)
> AIPS 4: INDISK 7 Disk drive #
> AIPS 4: BLC 14 *rest 0 Bottom left corner 0->use all
> AIPS 4: (1,2) -> exclude area
> AIPS 4: TRC *all 0 Top right corner 0 -> use all
> AIPS 4: (1,2) -> exclude area
> AIPS 4: OPTYPE ' ' 'MEDI' median window
> AIPS 4: 'LOW' use tight limits in rms
> AIPS 4: else use less tight limits
> AIPS 4: OPCODE 'WEIT' ' ' => RMS
> AIPS 4: 'WEIT' => WEIGHT=1/(RMS**2)
> AIPS 4: DOOUTPUT -1 > 0 => write S/N data cube
> AIPS 4: OUTNAME 'A2256C55W' Output name(name).
> AIPS 4: OUTCLASS 'RMTFA' Output name(class).
> AIPS 4: OUTSEQ 14 Output name(seq. #).
> AIPS 4: 0 => highest unique
> AIPS 4: OUTDISK 7 Output image disk drive #
> AIPS 4: 0 => highest with room
> AIPS 4: PIXRANGE 0 0 Range of intensities to plot
> AIPS 4: ZINC 1 Increment on freq axis.
> AIPS 4: < 0 => invert Z axis in plot
> AIPS 4: SMOOTH *all 0 Frequency smoothing function
> AIPS 4: LTYPE 3 Type of labeling: 1 border,
> AIPS 4: 2 no ticks, 3 standard, 4 rel
> AIPS 4: to center, 5 rel to subim cen
> AIPS 4: 6 pixels, 7-10 as 3-6 with
> AIPS 4: only tick labels
> AIPS 4: <0 -> no date/time
> AIPS 4: DOCENTER 1 > 0 => plot with box-like
> AIPS 4: lines, else plot line to
> AIPS 4: points
> AIPS 4: DOTV -1 > 0 Do plot on the TV, else
> AIPS 4: make a plot file
> AIPS 4: GRCHAN 0 Graphics channel 0 => 1.
> AIPS 4: DOPRINT 0 > 0 => use the terminal, 0 =>
> AIPS 4: no output, < 0 => printer or
> AIPS 4: file, > 72 => terminal width
> AIPS 4: OUTPRINT 'FITS:DOFARS.WEIGHTS'
> AIPS 4: Printer disk file to save
> AIPS 4: DOSLICE -1 > 0 => save as SLice file
> AIPS 4: XYRATIO 0 Stretch the X axis wrt Y axis
> AIPS 4: 0 -> 1.33
> >exit
> AIPS 3: Begin check for any 'standard' scratch files
> AIPS 3: Scratch files -- destroyed: 0 still active: 0
> AIPS 3: Destroyed empty user catalog on disk 3
> AIPS 3: Destroyed empty user catalog on disk 5
> AIPS 3: Destroyed empty user catalog on disk 6
> AIPS 3: User 2256: 899 messages, oldest written 13-OCT-2011 08:45:06
> AIPS 3: Deleted 0 messages
> AIPS 3: Saved POPS environment in area named 'LASTEXIT'
> AIPS 3: Session command-line history saved.
> AIPS 3: pilabo 31DEC11 TST: Cpu= 0.0 Real= 18 IO= 2
RSPEC changed adverbs from docrt to doprint and the procedure did not
catch up. Fixed now - RUN DOFARS and it should work.
ERic Greisen
More information about the Daip
mailing list