Fatal error with wsrtfits op=uvout

20 views
Skip to first unread message

Sarrvesh Sridhar

unread,
Aug 5, 2019, 4:23:21 PM8/5/19
to miriad
Hi,

I am trying to convert a Westerbork (single source single band) miriad data to uvfits using the command "wsrtfits op=uvout in=input out=output". However, the task fails with the following error message:

wsrtfits: version 1.3.1 21-Jul-04
### Warning [wsrtfits]:  Error writing to scratch file; check your $TMPDIR
### Fatal Error [wsrtfits]:  Invalid argument

The variable TMPDIR is set to /tmp but changing it does not seem to make a difference. Any idea what is going on? 

Best regards,
Sarrvesh

Peter Teuben

unread,
Aug 5, 2019, 4:26:52 PM8/5/19
to mir...@googlegroups.com

Hi Sarrvesh

    is it possible the disk where it wants to write to is too full for the output ?

wsrtfits is modeled after the fits command, which has the same argument structure. If problems remain, check to see if the source code in fits.for changed how to deal with $TMPDIR

peter

--
You received this message because you are subscribed to the Google Groups "miriad" group.
To unsubscribe from this group and stop receiving emails from it, send an email to miriad+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/miriad/58c085e9-7720-46ff-8277-797bf42301e3%40googlegroups.com.

melvyn wright

unread,
Aug 5, 2019, 4:45:08 PM8/5/19
to mir...@googlegroups.com

Hi Peter,

 grep did not find TMPDIR  --  can you offer some more clues ?

Melvyn.

/usr/local/miriad/src/subs/fitsio.for  /usr/local/miriad/src/subs/fitsio.h
morph>grep TMPDIR  $MIRSUBS/fit*
morph>grep TMPDIR $MIRPROG/convert/fit*

morph>ls $MIRPROG/convert/fit*
/usr/local/miriad/src/prog/convert/fits.for
/usr/local/miriad/src/prog/convert/fits.h
morph>grep TMPDIR $MIRSUBS/*.for
morph>grep TMPDIR $MIRSUBS/*.h
morph>grep TMPDIR $MIRPROG/convert/fit*

Peter Teuben

unread,
Aug 5, 2019, 4:48:50 PM8/5/19
to mir...@googlegroups.com

Mel,

  ah, but you forgot we have C files in $MIRSUBS (also f2c files). But in dio.c and scrio.c we use TMPDIR


peter

Sarrvesh Sridhar

unread,
Aug 6, 2019, 12:56:03 AM8/6/19
to miriad
Hi Peter,

I can confirm that there is enough space on the disk. I tried running wsrtfits with miriad_2012.7 and it does not produce the same error. I'll compare the source codes later today to see if something was changed. 

Best regards,
Sarrvesh

melvyn wright

unread,
Aug 6, 2019, 11:47:58 AM8/6/19
to mir...@googlegroups.com
Re Peter's note, i suppose that it may be  the usage of scrio in the fits routines that is likely the problem.
Are scrio and dio also updated with the updated code ?

Melvyn

Mel,

  ah, but you forgot we have C files in $MIRSUBS (also f2c files). But in dio.c and scrio.c we use TMPDIR

peter
--
You received this message because you are subscribed to the Google Groups "miriad" group.
To unsubscribe from this group and stop receiving emails from it, send an email to miriad+un...@googlegroups.com.

Peter Teuben

unread,
Aug 6, 2019, 11:55:49 AM8/6/19
to mir...@googlegroups.com

Nothing really changed in those C routines since 2011. Only a MAXPATH patch in dio.c in 2014.


I tried my current version of miriad, using exportuvfits() in casa to write a UVFITS file, then converted it via fits and wsrtfits and they both "worked".

casa -c "exportuvfits('aver_7.ms','aver_7.fits')"
fits in=aver_7.fits  out=aver_7.mir op=uvin
wsrtfits in=aver_7.fits  out=aver_7.mir2 op=uvin

Peter Teuben

unread,
Aug 6, 2019, 12:00:28 PM8/6/19
to mir...@googlegroups.com

another trick to figure out what the system is thinking that might point you to the error is the linux strace command, just prepend it to your command and watch the waterfall.

On 8/5/19 4:23 PM, Sarrvesh Sridhar wrote:
--
You received this message because you are subscribed to the Google Groups "miriad" group.
To unsubscribe from this group and stop receiving emails from it, send an email to miriad+un...@googlegroups.com.
Reply all
Reply to author
Forward
0 new messages