Konrad R. W. Tristram

Staff Astronomer
European Southern Observatory
Main
Research
Publications
Downloads
About me
Other

Reports

There are currently four reports for download. They are all of rather technical nature and address the positioning of the source in MIDI, the MID-infrared Interferometric instrument at the Very Large Telescope Interferometer. The reports are ordered chronologically.

 

IDL programmes

Below are a few IDL (Interactive Data Language by ITT Visual Information Solutions) programmes I wrote for my work. Most of the programmes are for reducing, visualising and modelling interferometric data. Feel free to download and use any of the routines. Many programmes require the AstroLib package (http://idlastro.gsfc.nasa.gov/) and/or a recent version of the MIA+EWS package (http://www.strw.leidenuniv.nl/~nevec/MIDI/) to be installed. The first is actually contained in the latter. A more detailed description of every routine is contained along with the code in the individual files. If you have any suggestions for improving the code, please contact me.

ctableplot.pro   Plotting tools
CTABLEPLOT–Make a plot of a color table (2009-05-04).
ctableplot, num, /ps
midi_2bbfit.pro   MIDI two component black body fit
MIDI_2BBPRE–Example routine: Create the variables needed for a fit (2009-07-22).
transf = midi_2bbpre()
MIDI_2BBMDL–Create a 2 component elliptical Gaussian black body brightness distribution (2009-07-22).
gmodel = midi_2bbmdl(arrsiz, wavele, params, parepl, pscale=2.0, tausil=tausil, /no_uvp, fluxes=fluxes, /clumps, /reclmp, sizclp=sizclp, ranvec=ranvec)
MIDI_2BBDAT–Calculate and return the "visibility data" for the model (2009-07-22).
mdlvis = midi_2bbdat(transf, params, /clumps, /reclmp, sizclp=sizclp, ranvec=ranvec)
MIDI_2BBFIT–Run the fit for the 2 component Gaussian black body model (2009-07-22).
midi_2bbfit, transf, result
MIDI_2BBCMP–Compare a 2 component Gaussian black body model to the data (2009-08-18).
midi_2bbcmp, transf, params, wavele, wavran=wavran, /corflx, /visamp, /visphi, /normal, pmulti=[4,6], psfile=psfile, /compact
midi_comp.pro   MIDI model comparison
MIDI_PREPHT–Prefix the photometry to a MIDI result structure (2012-11-12).
outdat = midi_prepht(orgdat, flxsel, /noweight, /totflx)
MIDI_EXTVIS–Extract visibilities from a model at specified UV coordinates (2010-09-23).
result = midi_extvis(uvplan, uvcoor, rotang, /uvplot, /direct)
MIDI_DATCOM–Interpolate MIDI data at certain wavelengths (2012-10-17).
result = modl_datcom(caldat, wavele, smtval, wavwid=0.0)
MIDI_DPHASE–Determine the differential phase from input phase data (2009-08-18).
result = modl_dphase(mdlvis, method, wavran=[8,13], smtval=3, /output, /silent)
MIDI_PLTCMP–Plot a comparison of measured data and modelled data (2012-11-19).
modl_pltcmp, datvis, mdlvis, mdlvi1, mdlvi2, wavran=[8,13], /corflx, /visamp, /visphi, /normal, pmulti=[4,6], ewidth=0.0, psfile='psname', /compact
midi_redu.pro   MIDI data reduction
MIDI_DATSET__DEFINE–Automatically define a structure of type MIDI_DATSET (2010-02-05).
datset = {midi_datset}
MIDI_GETKEY–Get FITS header keywords for a list of MIDI files (2008-08-26).
result = midi_getkey(filist, keywrd, direct='/data/mydata/')
MIDI_FILLIN–Automatically fill in information in a MIDI data set (2012-02-07).
midi_fillin, datset, /rmpath, /identi, /redpar, /caldat
MIDI_GETDAT–Read raw and calibrated MIDI data from EWS into a struct (2012-11-17).
result = midi_getdat(datset, /rawdat, /fixbli, /silent)
MIDI_CALINF–Get calibrator information for a list of calibrated MIDI results (2008-08-26).
calinf = midi_calinf(caldat, datset)
MIDI_RVBCAL–Recalibrate fluxes with templates by Roy van Boekel - OBSOLETE (2010-07-13).
result = midi_rvbcal(caldat, calnam, calflx, /output)
result = midi_rvbcal(caldat, calinf.calnam, calinf.calflx)
MIDI_REDPLT–Plot the fluxes and visibilities of a MIDI result structure (2012-11-01).
midi_redplt, result, ctable, colour, xrange=xrange, /rawdat, /plterr, psfile='psname'
MIDI_MODPHT–Modify the sky bands and recompile EWS photometry programmes. (2010-07-12).
midi_modpht, skypos, /verbose
MIDI_REDSET–Fully reduce the data in a MIDI data set (2012-08-30).
midi_redset, datset, 'prefix', /noredu , /shimsk, /fringe, maxopd=200, /slopef
MIDI_PLTPRP–Plot a FITS header keyword for several MIDI data files (2009-10-09).
midi_pltprp, fileli, keywrd, ytitle, psfile='psname', yrange=[0,0], titles='Title', labels=labels, dprcat=dprcat, xoffse=xoffse, yoffse=yoffse, angles=angles
MIDI_PLTPDS–Plot common FITS header keywords for fringe tracks in a MIDI data set (2008-10-14).
midi_pltpds, datset, titles, pspath='psname'
MIDI_A2JAVE–Calculate the average and scatter of conversion factors (2012-08-29).
a2jave = midi_a2jave(adu2jy, mjdobs, nummax, sigclp=sigclp, flgdat=flgdat, /silent)
MIDI_ADU2JY–Plot the ADU to JY conversion factor for a set of calibrators (2012-10-11).
adu2jy = midi_adu2jy(datset, 'Title', pspath='psname', /noplot, /plterr, caldatabase=caldatabase, sigclp=1.5, /silent)
MIDI_CALERR–Estimate and add the errors from the variation of the transfer function (2012-10-10).
newdat = midi_calerr(caldat, datset, calibs, mincal=mincal, maxcal=maxcal)
MIDI_STRC2T–Print out an IDL structure as an ASCII table (2009-06-26).
midi_strc2t, datset, tagnam, maxwid=40
MIDI_DATA2T–Print out MIDI data as an ASCII file (2009-06-23).
midi_data2t, caldat, filnam, wavele
MIDI_DATMRG–Merge several data sets using Gaussian error propagation (2011-01-26).
result = midi_datmrg(caldat)
MIDI_PLOT3D–Make a 3D plot of the UV plane using iPlot (2009-07-01).
midi_plot3d, caldat, wavele, ranges
midi_utils.pro   MIDI data tools
MIDI_CHPDAT–Reduce chopped data from MIDI (2010-11-03).
chpres = midi_chpdat(rawdat, err_stat)
MIDI_WAVAXI–Helper function: Plot a wavelength axis on top of the previous plot (2008-03-15).
midi_wavaxi, disper, xtickl=0.02, /pltlab
MIDI_OUTDIS–Helper function: Plot one detector window for dispersed MIDI data (2009-06-30).
midi_outdis, detwin, positi, offset, ctable=3, disper='prism', /llabel, /rlabel, /tlabel, /blabel
MIDI_OUTFIE–Helper function: Plot one detector window for undispersed MIDI data (2010-12-06).
midi_outfie, detwin, positi, offset, ctable=3, /llabel, /rlabel, /tlabel
MIDI_PLTWIN–Make a plot of the two detector windows of MIDI data (2012-11-01).
midi_pltwin, rawdat, ctable, psfile='psname', metafile=metafile, mindat=mindat, maxdat=maxdat, /linear, /sqroot, /logari, /invers, factor=3, titles='Title', /encaps
MIDI_PHOTOM–Chop a MIDI photometry and extract the spectrum (2012-11-01).
photom = midi_photom(phtdat, mskfil=mskfil, sigma=3.0, /reset, /chop, /output, /nosky, porder=porder, phtbnd=[12,21], skydis=3, skywid=3, skysmo=10, psfile='psname', /invers, titstr='Title'])
MIDI_PHTFIL–Reduce a MIDI photometry data file (2012-04-03).
photom = midi_phtfil(phtfil, mskfil=mskfil, /nosky, phtbnd=[12,21], skydis=4, psfile='psname', /invers)
MIDI_PLTFRI–Show a video of the delay line position and the fringe (2012-04-03).
midi_pltfri, ifdata, frames, delay, msknam=msknam, tsmoot=20, /davera, factor=3
MIDI_GAUDIS–Helper function: Fit a Gaussian to every channel in a dispersed MIDI window (2011-07-13).
fitpar = midi_gaudis(detwin, rrange=[9,26], smtval=[0,0], /output, /quiet)
MIDI_MSKFIT–Create a MIDI mask from A and B photometry files (2012-08-30).
mskdat = midi_mskfit(infile, msknam, factor=1.25, additi=2.0, smtval=[0,0], /output, psfile='psname', skymsk=skymsk, skydis=1.0, skywid=4.0, yshift=0.0, linpar=1.0)
MIDI_FLAGOK–Determine which frames where flagged as good in a fringe track (2012-04-03).
flagok = midi_flagok(indata, frames)
MIDI_UNWRAP–Unwrap phases assuming a smooth phase change (2012-04-03).
newphi = midi_unwrap(oldphi, smooth, maxval=40)
MIDI_GETPHI–Get the water vapour phases from a FITS file (2012-07-19).
phases = midi_getphi('/tmp/tmp.ungroupdelay.fits')
MIDI_PLTPHI–Plot the water vapour phases for a MIDI fringe track (2012-08-03).
midi_pltphi, filnam, xrange=xrange, yrange=yrange, /noflag, titles='Title', /symbol, refphi=refphi, smtval=40, psfile='psname'
MIDI_PLTOPD–Plot the Optical Path Difference for a MIDI fringe track (2012-09-17).
midi_pltopd, filnam, xrange=xrange, yrange=yrange, /noflag, titles='Title', /symbol, /pltimg, gsmooth=0, asmooth=0, minmax=minmax, /invers, /phases, /noinst, /nogdel, refdel=refdel, getdel=getdel, maxopd=150, pltpos=pltpos, /noxlab, /flip, psfile='psname'
MIDI_PLTOVL–Plot the beam overlap for the photometry, fringes and masks (2012-11-01).
midi_pltovl, dattag, wavran, mskfil='msknam.fits', /normal, pltpos=pltpos, psfile='psname', ierror=ierror
MIDI_PLTRES–Plot a page of diagnostic plots for reduced MIDI data (2012-12-12).
midi_pltres, dattag, psfile, adu2jy, mskfil='msknam.fits', yrange=yrange
MIDI_PLTRAD–Plot the correlated flux or visibility as a function of baseline length (2010-10-05).
midi_pltrad, caldat, wavele, wavdel, /corflx, /visamp, xrange=xrange, yrange=yrange, /ylog, titles='Title', colour=[50], psfile='psname'
MIDI_PLTUVP–Plot points in a UV plane with values from a MIDI result struct (2012-03-29).
midi_pltuvp, caldat, wavele, tagnam='corflx', zrange=zrange, declin=declin, blname=['U1-U2'], minalt=21, pblran=[-100.,+100.], titles='Title', /sqroot, barlab='Color bar label', /noline, /labels, /oldwin, pltsiz=0, psfile='psname'
modl_geom.pro   MIDI models - geometrical models
MODL_EGAUSS–Generate a generalised elliptical gaussian or uniform disk (2009-07-13).
result = modl_egauss(arrsiz, fwhmax, iangle, pangle, boxyns=0.0, xoffse=0.0, yoffse=0.0, expone=2.0, /ufdisk, /silent)
MODL_GAUSSM–Generate an set of generalised elliptical Gaussians (2009-07-13).
result = modl_gaussm(arrsiz, wavele, fwhmax, iangle, pangle, boxyns=0.0, pscale=1.0, xoffse=0.0, yoffse=0.0, /ufdisk, expone=2.0, /silent, /output)
MODL_POLARM–Generate a polar model for a single wavelength (2009-07-13).
result = modl_polarm(arrsiz, fwhmax, iangle, pangle, expone=2.0, porder=1.0, xoffse=0.0, yoffse=0.0, /gaussd, pscale=1.0, /silent, /output)
MODL_TDISKM–Generate a temperature disk model with T(r) = r^(-EXPONE) (2009-07-13).
result = modl_tdiskm(arrsiz, wavele, innrad, iangle, pangle, xoffse=0.0, yoffse=0.0, expone=0.9, t_maxi=500.0, silabs=silabs, sildep=0.5, pscale=1.0, /fillit, /behave, /output, /silent)
modl_utils.pro   MIDI models - model utilities
MODL_MKHEAD–Generate a default FITS header for the models (2009-07-13).
header = modl_mkhead(arrsiz, functi, object, pscale)
MODL_FITS2M–Read FITS files into a model struct (2011-02-21).
result = modl_fits2m(patter, rotang, /silent)
MODL_M2FITS–Write a model into a FITS file or into several FITS files (2009-07-13).
modl_m2fits, models, filnam, /cube
MODL_M2VISI–Calculate the UV-plane, i.e. the complex visibilities for a model (2009-07-13).
result = modl_m2visi(models, /output, /silent)
MODL_M2PLOT–Plot a flux distribution and/or UV-plane to screen or PS file (2012-11-01).
modl_m2plot, uvplan, psfile, metafile=metafile, minval=1.0e-20, /flux, /fftd, delay=0.0, factor=1.0, maxpbl=100, /overpl
MODL_M2GAUS–Determine the gaussian width of a brightness distribution (2009-07-13).
modl_m2gaus, models, psfile='psname'
plotcbar.pro   Plotting tools
PLOTCBAR–Helper function: Add a vertical colour bar to a plot (2010-08-11).
plotcbar, values, colour, barpos, ctable, barlab=barlab, tickle=0.29, majval=majval, minval=minval, format='', /endlab
plotfits.pro   Plotting tools
PLOTFITS–Make a plot of FITS imaging data (2012-11-01).
plotfits, fitsim, header, minmax=minmax, /linear, /sqroot, /logari, /asinhy, ctable=3b, /invers, factor=1.0, /pixels, /decima, /sexage, xrange=xrange, yrange=yrange, /offset, psfile='psname', /encaps, pltset=pltset, titles='Title', positi=positi, barpos=barpos, bargap=0.0, colbar='Intensity', colmap=colmap, colinv=colinv, /dfonts
setplot.pro   Plotting tools
SETPLOT–Set the plot environment for plots on screen or to PS or METAFILE (2012-11-01).
setplot, pltset, pltpos, /small, /normal, /large, /longer, psfile='psname', /encaps, metafile=metafile, /close, /oldwin
sirs_data.pro   Spitzer IRS data
SIRS_READAT–Read Spitzer IRS data from a pipeline reduced file (2009-07-02).
result = sirs_readat(filnam, /nocrop)
SIRS_COMPRE–Compress the overlapping orders of an IRS spectrum (2009-07-02).
result = sirs_compre(data)
strctxt.pro   General
GETTAGS–Helper function: Recursively get the tag names of an IDL struct (2009-06-26).
gettags, struct, vartyp=vartyp, /expand
REFPTRA–Helper function: Dereference an array of pointers (2009-06-26).
refptra, arrptr, arrmin
STRCTXT–Print out an IDL structure as an ASCII table (2012-08-27).
strctxt, struct, tagnam, maxwid=40, filnam='filnam.txt'
vlti_utils.pro   VLTI utilities
VLTI_HRANGE–Calculate the hour angle range where an object is visible (2009-07-06).
hrange = vlti_hrange(declin, latitu, minalt, /silent)
VLTI_BLPROJ–Calculate projected baselines positions or tracks for the VLTI (2012-11-17).
bltrck = vlti_blproj(declin, basein, hangin, hrange=[-3.0,+3.0], hdelta=0.1, numele=512)
VLTI_UVTRCK–Calculate a set of UV tracks for a certain declination (2012-11-17).
uvtrck = vlti_uvtrck(declin, hrs, /print)
VLTI_UVCHCK–Check the UV coordinates saved in the fits header (2012-11-17).
result = vlti_uvchck(filnam, /silent, bltrak=bltrak)
VLTI_PLTUVP–Plot a UV plane and plot UV tracks for a certain declination (2012-11-01).
vlti_pltuvp, pltset, declin=declin, blname=['U1-U2'], minalt=21, pblran=[-100.,+100.], titles='Title', pltpos=pltpos, barpos=barpos, /noline, /noarrow, /nolabel, /oldwin, pltsiz=pltsiz, psfile='psname', /close
© 2005-2024 Konrad R. W. Tristram Last modified: 15/04/2014, 01:40. Valid HTML 4.01! Valid CSS!
This webpage is best viewed with at least a resolution of 1024×786, 24 bit colour, Mozilla 5.0, IE 8.0 or Opera 8.0, as well as Java(Script) enabled.