[daip] HGEOM fails on a big cube
Andrea Petric
apetric at aoc.nrao.edu
Mon Oct 16 13:40:20 EDT 2000
Hi DAIP,
HGEOM run on a 1024x1024x63 cube dies on channel 63, line 950, as
follows:
> inp hgeom
AIPS 1: HGEOM Task to make an image consistent with another image
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'N1058 C+CS+D' Input image (name)
AIPS 1: INCLASS 'ICLN 3' (class)
AIPS 1: INSEQ 1 (seq. #)
AIPS 1: INDISK 2 (disk unit #)
AIPS 1: IN2NAME 'N1058 C+CS+D' Second input image (name)
AIPS 1: IN2CLASS 'MOM1 ' (class)
AIPS 1: IN2SEQ 2 (seq. #)
AIPS 1: IN2DISK 2 (disk unit #)
AIPS 1: OUTNAME 'N1058 C+CS+D' Output image (name)
AIPS 1: OUTCLASS 'TMP ' (class)
AIPS 1: OUTSEQ 0 (seq #)
AIPS 1: OUTDISK 5 (disk unit #)
AIPS 1: BLC *all 0 Bottom left corner of image
AIPS 1: 0=>entire image
AIPS 1: TRC *all 0 Top right corner of image
AIPS 1: 0=>entire image
AIPS 1: IMSIZE 0 0 Output image size in pixels
AIPS 1: APARM *all 0 (1)=terp.ord., (2)=zero flag
AIPS 1: (3,4) = new reference pixel
AIPS 1: = (0,0) => center in output
>
> prta 'hgeom' ; prtmsg
HGEOM1 10:46:56 Task HGEOM (release of 31DEC99) begins
HGEOM1 10:46:56 You are using a non-standard program
HGEOM1 10:48:56 Create N1058 C+CS+D.TMP . 1 (MA) on disk 5 cno 41
HGEOM1 10:48:56 Using 699924 words of buffer
HGEOM1 10:54:00 ERROR EXIT AT 980 IN GEOSUB: 5085 17 9
HGEOM1 10:54:00 SUBROUTINE MAPIO RETURNED ERROR 9 AT 85
HGEOM1 10:54:00 LOGICAL END OF FILE
HGEOM1 10:54:00 DISK 2 SLOT 2 USER 2154 (MA)FILE: N1058 C+CS+D.ICLN 3.1
HGEOM1 10:54:00 Destroyed MA image file: catno= 41 disk= 5
HGEOM1 10:54:00 Purports to die of UNNATURAL causes
HGEOM1 10:54:00 spectra 31DEC99 TST: Cpu= 251.8 Real= 424
We looked at a TVMOVIE of the cube & ran IMEAN on it; both worked without
any problems (apart from the IMEAN bug just reported), and certainly channel
63 itself both looks fine & gives a reasonable IMEAN. There are no blanked
pixels in the input cube, and the header of the output file (which we checked
before the thing was destroyed, while HGEOM was running) looked fine:
10:47:18 Got(1) disk= 5 user=2154 type=MA N1058 C+CS+D.TMP.1
10:47:19 imh
10:47:19 Image=N1058 (MA) Filename=N1058 C+CS+D.TMP . 1
10:47:19 Telescope=VLA Receiver=VLA
10:47:19 Observer=AR296 User #= 2154
10:47:19 Observ. date=14-JUN-1993 Map date=16-OCT-2000
10:47:19 Minimum=-3.28992074E-03 Maximum= 1.29607338E-02 JY/BEAM
10:47:19 ----------------------------------------------------------------
10:47:19 Type Pixels Coord value at Pixel Coord incr Rotat
10:47:19 RA---SIN 1016 02 40 23.200 508.00 -4.000 0.00
10:47:19 DEC--SIN 1016 37 07 48.000 509.00 4.000 0.00
10:47:19 FELO-HEL 63 5.1800000E+05 33.00 -2.5852092E+03 0.00
10:47:19 STOKES 1 1.0000000E+00 1.00 1.0000000E+00 0.00
10:47:19 ----------------------------------------------------------------
10:47:19 Coordinate equinox 1950.00
10:47:19 Map type=NORMAL Number of iterations= 13922
10:47:19 Conv size= 15.43 X 13.52 Position angle= 87.00
10:47:19 Rest freq 1420.406 Vel type: OPTICAL wrt SUN
10:47:19 Alt ref. value 1.41726E+09 wrt pixel -30.00
For reference, IMHEAD on the original image gives
>tget hgeom
>imh
AIPS 1: Image=N1058 (MA) Filename=N1058 C+CS+D.ICLN 3. 1
AIPS 1: Telescope=VLA Receiver=VLA
AIPS 1: Observer=AR296 User #= 2154
AIPS 1: Observ. date=14-JUN-1993 Map date=31-MAY-1996
AIPS 1: Minimum=-3.28992074E-03 Maximum= 1.29607338E-02 JY/BEAM
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: RA---SIN 1024 02 40 23.200 512.00 -4.000 0.00
AIPS 1: DEC--SIN 1024 37 07 48.000 513.00 4.000 0.00
AIPS 1: FELO-HEL 63 5.1800000E+05 33.00 -2.5852092E+03 0.00
AIPS 1: STOKES 1 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 1950.00
AIPS 1: Map type=NORMAL Number of iterations= 13922
AIPS 1: Conv size= 15.43 X 13.52 Position angle= 87.00
AIPS 1: Rest freq 1420.406 Vel type: OPTICAL wrt SUN
AIPS 1: Alt ref. value 1.41726E+09 wrt pixel -30.00
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type CC is 126
AIPS 1: Maximum version number of extension files of type PL is 15
AIPS 1: Keyword = 'WTNOISE ' value = 1.184647E+00
AIPS 1: Keyword = 'PARANGLE' value = -2.717787E+01
AIPS 1: Keyword = 'ZENANGLE' value = 1.026279E+02
and IMHEAD on the "reference" image (#2) gives
>inna in2na
>incla in2cla
>inse in2se
>indi in2di
>imh
AIPS 1: Image=N1058 (MA) Filename=N1058 C+CS+D.MOM1 . 2
AIPS 1: Telescope=VLA Receiver=VLA
AIPS 1: Observer=AR296 User #= 2154
AIPS 1: Observ. date=14-JUN-1993 Map date=23-SEP-2000
AIPS 1: Pixel type: FLOATING Magic value blanking
AIPS 1: Minimum= 4.43028938E+05 Maximum= 5.98141500E+05 M/S
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: RA---SIN 1016 02 40 23.200 508.00 -4.000 0.00
AIPS 1: DEC--SIN 1016 37 07 48.000 509.00 4.000 0.00
AIPS 1: STOKES 1 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 1950.00
AIPS 1: Map type=NORMAL Number of iterations= 13922
AIPS 1: Conv size= 15.43 X 13.52 Position angle= 87.00
AIPS 1: Rest freq 1420.406 Vel type: OPTICAL wrt SUN
AIPS 1: Alt ref. value 1.41726E+09 wrt pixel -30.00
AIPS 1: Maximum version number of extension files of type HI is 1
>
Thanks,
Andreea & Michael R.
More information about the Daip
mailing list