What's Old in MIRIAD

This document describes changes and developments in the MIRIAD software older than 6 months.
March 27, 1996
new subroutines
  • align.for - check for alignment of two images. Extracted from CLEAN.
  • planet.for - planetary ephemeris, with only moderate accuracy.
new programs
  • imblr - replace blanked values with a new value
  • imcfn - computes confusion noise for an interferometer image
  • imsad - image "search and destroy"
  • gethd - print out header varialbles
  • planets - Print ephemerides of planets and the Sun
  • solar - report some parameters on solar system objects
  • mostess - Maximum Entropy deconvolution for a mosaiced observation
  • tvclip - Interactive editing of a UV data set on a TV device (see tvflag)
  • closure - plots closure phase and triple product amplitude
  • uvfstats - print flagging statistics
  • uvtcor - redo the Tsys correction, using a median smoothed Tsys curve
  • uvtplot - plot Tsys for each field of a mosaiced observation

October 5, 1995
RCS Archive With the departure of Mark Stupar (who kept the miriad distribution alive) we wanted to come up with a less work-intensive method of source code distribution. Doug Roberts (droberts@lai.ncsa.uiuc.edu) has set up an remotely accessible RCS system which will allow programmers to check new or modified source code into the system. Distribution of source code, in the form of tar file, for outside users remains the same.

January 31, 1995
PGPLOT 5.0 Jim Morgan has updated just about every Miriad site within BIMA for this new major update of our graphics subroutine library PGPLOT. He has tried to account for every device that is usually used. The following changes apply:
  • The driver /xdisp has been removed (along with pgdisp).
  • The /xw driver has been completely re-written and is about 3-4 times faster than before. It now uses a client-server model which will automatically start up the first time you call any PGPLOT routine (this will appear as an icon on your screen).
  • A new driver /xs is just like /xw except it stays up after you are finished with your plot.
  • Multiple X-window plot windows may be used together by specifying them as device=#/xs (or #/xw) where # is any integer number.
  • There is a new image plotting routine which both Neil and Jim are incorporating into their plotting routines (cg* and WIP) which will permit color palette loading both on the screen and on hardcopy devices. This will remove the need for my pscolor post processor routine (indeed, it does not work with PGPLOT 5.0's /cps drivers).
  • GIF and PPM hardcopy drivers have been added and are available.
  • C-bindings for each Fortran callable PGPLOT routine are available. The include file and library are in $MIRLIB. Use of these routines requires an ANSI-level compiler.
  • Programmers should not use the new rubber band routine PGBAND. Jim believes it is only currently functional with the new /xw and /xs drivers.
  • Online documentation for subroutines is now availble in full HTML format.
Further information (and X resources for the /xw and /xs drivers), can be found in the file $MIR/borrow/pgplot/aaaread.me. Any problems should be sent to morgan@astro.umd.edu .

LINUX An experimental version of miriad running under the LINUX operating system is now available and will be distributed officially in the next upcoming release. There are some patch files available and be sure to read the install.LINUX file.

mosaic 3.0 The mosaic program has had a major upgrade. Starting with this new version (3.0) the old huge single source code has been broken into smaller pieces (see $MIR/src/spec/mosaic) and also supports the CM5.

December 6, 1994
Linear Mosaicing: Large ATNF update ( 06dec94 ) to allow linear mosaicing in INVERT ( options=mosaic. Two new programs, MOSMEM (maximum entropy deconvolve) and YYY (mosaic Steer CLEAN) can then be used to deconvolve. Finally RESTOR does what is has always done. Not much documentation is present yet.
Primary Beam: Programs that deal with primary beams, should have their keyword pbfwhm replaced by the pbtype . For example (from ELLINT) "pbfwhm=120" has become "pbtype=gaus(120)" (or "pbtype=hatcreek").
  • new programs:
    • mosmem - Maximum Entropy deconvolution for a mosaiced image (rjs)
    • mospsf - Determine approximate PSF of a mosaic experiment
    • mossdi - Mosaic Steer CLEAN algorithm (rjs)
    • mossen - Determine the rms and gain in a mosaiced image (rjs)
    • pbplot - Plot primary beam shapes (rjs)
  • new subroutines:
    • hdtab.for - A set of routines, used by INVERT, to build up the header of a datase
    • mapper.for - A set of subroutines to form an image from visibility data
    • mc.for - A set of routines to perform mosaic "convolution" operations
    • mostab.for - A set of routines to handle multi-pointing vis data sets, mosaic tables, and to mosaic images
    • rest.for - A set of routines to aid restoring mosaiced images
    • packALPHA.c - data conversion routines for little endian machines (PC, DEC alpha)
  • A new dataset item, mostable , is created and used by the new multi-pointing mosaicing software. Programs that enumerate items for copy (using hdcopy), need to have their item-list (sometimes erroneously called keywords) extended. Program like IMLIST and PRTHD know about this new dataset-item.
  • New pb* routines to handle primary beams. This due to the new mosaicing software.

November 15, 1994
  • programs
    • gaufit : Fix erroneous 'Center outside range' for cdelt3<0 (bpw)
    • imspec, imstat : new option 'eformat' (bpw)
    • bee : Added flags to BELIST, and bypass questions in EDIT (mchw)
  • new subroutines:
    • tabflux : alias installed, for example with this patch calflux source=pk0420 device=/xw now plots whole history.

October 24, 1994
Major ATNF extensions ( 23sep94 and 30sep94 ) to better handle coordinates into Miriad. Primarely to handle RA and DEC accurately, but also to a lesser extent to better support images with an optical velocity convention (FELO axis).
  • new programs
    • uvsub : Subtract point source models from visibilities (nebk)
    • uvclip : determine data to flag from channel-0 uv dataset (njt)
    • uvpnt : create closely spaced uv-data from multiple pointings (mchw)
  • new subroutines:
    • co.for, cosubs.for : convert between coordinate systems (rjs)
    • keyline.for : linetype parser for user interface
    • lsqfit.for : least squares
    • nswc.for : borrowed numerical routines
    • lsqfit.for : least squares
  • passfit and calapply cannot handle the new corellator modes, the problem has been found and will be distributed shortly to all sites.

June 23, 1994
  • A duplicate version of uvpflag (version 8-Mar-93) has been removed from the system. The correct version is dated 12-Jan-94 or later.
  • fixed subroutines:
        boxes    - potential I/O closing error
        cgpgsubs - box routines improved
        fitsio   - tolerant RA/DEC labeling in non-AIPS uv fits files
        ephem    - added subroutine sunradec 
        obspar   - added WSRT/evector as 90 degrees
        uvdat    - less verbosities
        cgsubs   - moved routines
        ofm      - Add fixed zero colour contours & hardcopy device capability
        sort     - fixed sorting 1 element 
        pghline  - add pgbbuf/pgebuf calls
  • Some unlisted programs output format was polished up, and for many programs the description of the region= has been clarified if they allow pixel blanking and/or non-rectangular regions.
  • gapply (Teuben; 21-jun-94). Added options=jyperk for baseline based gains (jyperk=Jy/K) to improve the noise estimates from INVERT.
  • contsen (Sault; 10-mar-94). Added xrange= keyword.
  • uvflux (Sault; 05-Nov-93) Use double precision, better doc, print number visibs
  • uvlin (Sault; 14-mar-94) Fixed shift option, and added sun and twofit options.
  • uvplt (Killeen; 17-mar-94) Added 2pass option and changed nocolor option.
  • uvspec (Killeen; 22-Mar-94) Added unflagged option.
  • cgcurs (Killeen; 11-Mar-94) Add spatial binning.
  • cgdisp (Killeen; 11-Mar-94) Add spatial binning and OPTIONS=SINGLE.
  • cgslice (Killeen; 11-Mar-94) Add spatial binning.
  • cgspec (Killeen; 18-Mar-94) Add "dspectrum" image type.
  • mfcal (Killeen; 30-Mar-94) Add options=interpolate.
  • gpcal (Sault; 23-Mar-94) Closure checking.
  • gpplt (Sault/Killeen; 29-Mar-94) Fix bugs in listing of bandpass table.
  • copyhd and delhd (Sault; 25-apr-94) Corrected call to closing the dataset which could potentially have faulted.
  • fits (Sault: 6-Apr-94) Check antenna table ref freq for validity

May 11, 1994
New Programs:
  • imfit fits a model, M(x,y,z) to a Miriad image, I(x,y,z). The model is specified by up to 6 parameters, and an initial guess. The convergence and required accuracy of the initial guess depend on the model option. Other models may be added by changing a subroutine. This program is currently under development. The author, wright@ral.berkeley.edu solicits comments.
  • implot (Mel Wright; 29-apr-94) fixed a bug in beam position, courtesy D.Wilner.
  • corset3 (Jim/Mel; 26-apr-94) keyword if was changed to iffreq . Various format changes, tables in the doc page were corrected
  • uvindex (Neil Killeen; 28-mar-94) A new keyword interval was added.
  • uvlist (Peter/Mel; 15-apr-94) A efficient new keyword scale was added to optionally rescale the output amplitudes.
  • gfiddle (Peter; 3-may-94) now doesn't allow you to fit or flag if the output file already exists, in accordance with normal miriad practice.
  • MAXWIN related cosmetic changes to uvindex, uvcat, uvcal, uvspect, uvgains, zeesim, and calapply. A new global variable, MAXWIN, is now defined in the include file maxdim.h instead of private ones. This was necessary to align software for the new correlator software, which is able to handle up to 16 different spectral windows.

March 10, 1994
  • corset3 (Mel Wright; 10-Mar-94) Corrected the plot and channel numbers for new modes 6,7 and 8.
  • mosaic (Bart Wakker; 6-Jan-94) Many many changes, mostly to make the code run on the Illinois CM5.
March 8, 1994
  • uvlin (Bob Sault; 23-Dec-93) The shift option was not shifting data back to the original position.
  • uvfft (Mel Wright; 07-Mar-94) Documentation improved and various extra checks were added to the code.
March 6, 1994
  • tvflag (Bob Sault; 4-Mar-94) should now be able to handle data that is not in time order.
  • uvgen (Mel Wright; 03-Mar-94) has an added elevation limit and uses standard keywords dec= and elev=. Note that, in order to make it's behavior more similar to uvtrack, the default behavior of the program has now slightly changed. Also handles the new generalized hybrid correlator model. (see also corset3).

February 20, 1994
New Programs:
  • imheq (Neil Killeen). IMHEQ writes a new image to which histogram equalization has been applied.
  • tvflag (Bob Sault). If you had mysterious results, running it again with this new version may resolve them.
  • uvindex (Bob Sault) now generates messages when there is a long gap in the observing time.
  • gapply (Peter Teuben) now has an option to allow extrapolation of amplitudes. This allows you to amplitude calibrate with a source that was not part of the regular observing time interval. Be sure a low order polynomial (0) was used in gfiddle, since the nearest edge value is used for extrapolation.

February 6, 1994
New Programs:
  • cgslice (Neil Killeen). CGSLICE displays an image via a contour plot or a grey scale on a PGPLOT device. The cursor (or a text file with slice positions) is then used to define the end points of 1-D slices which are marked on the image, and then plotted.
  • visibility datasets can now be read by by anybody, you don't need write permission to just use UVLIST or any program that does not attempt to write.

Home New General Proposal Status U Info U Tools Memos Learn

Please direct comments, criticisms, corrections and contributions to: teuben@astro.umd.edu.
Page last modified: Tuesday, 20-Oct-1998 14:08:22 EDT.