! @(#)necshow.prg 17.1.1.1 (ES0-DMD) 01/25/02 17:51:12 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++ ! !.COPYRIGHT (C) 1991 European Southern Observatory !.IDENT echshow.prg !.AUTHOR Pascal Ballester, ESO - Garching !.KEYWORDS Spectroscopy, Echelle, !.PURPOSE Display echelle session parameters !.VERSION 1.0 Creation 21-AUG-1991 PB ! !------------------------------------------------------- ! SET/FORMAT I1 F12.5 ECHO/OFF SESSLINE = 0 ! Initialize the line counter IF P8(1:1) .EQ. "H" GOTO HELP ! IF P1(1:1) .EQ. "G" GOTO GENERAL IF P1(1:1) .EQ. "O" GOTO ORDER IF P1(1:1) .EQ. "B" GOTO BKG IF P1(1:1) .EQ. "C" GOTO FILT IF P1(1:1) .EQ. "E" GOTO EXTR IF P1(1:1) .EQ. "W" GOTO WLC IF P1(1:1) .EQ. "R" GOTO REBIN IF P1(1:1) .EQ. "F" THEN IF P1(2:2) .EQ. "F" THEN GOTO FLAT ELSE GOTO FLUX ENDIF ENDIF IF P1(1:1) .EQ. "M" GOTO MERGE ! GENERAL: ! SHOW/SESS "**************** Echelle Parameters *******************" SHOW/SESS " " SHOW/SESS "--- (G) General Information ----" SHOW/SESS " Instrument : INSTR = {INSTRUM}" SHOW/SESS " Dispersers : GRATING = {GRATING} CROSS = {CROSS}" SHOW/SESS " Detector : CCD = {CCD}" SHOW/SESS " Binning factor : CCDBIN = {CCDBIN(1)},{CCDBIN(2)}" SHOW/SESS " Image size : IMSIZE = {IMSIZE(1)},{IMSIZE(2)}" SHOW/SESS " Scan limits : SCAN = {SCAN(1)},{SCAN(2)}" SHOW/SESS " Orders : {ECHORD(1)} orders, from {ECHORD(2)} to {ECHORD(3)}" SHOW/SESS " Wavel. range : {WLRANGE(1)} - {WLRANGE(2)} wavel. units" SHOW/SESS " " ! IF P1(1:1) .NE. "?" GOTO RET ! ORDER: ! SHOW/SESS "--- (O) Order definition ----" SHOW/SESS " Method (STD/COM/HOUGH) : DEFMTD = {DEFMTD}" SHOW/SESS " Polynomial Degree(X,Y) : DEFPOL = {DEFPOL(1)},{DEFPOL(2)}" SHOW/SESS " Order Reference Frame : ORDREF = {ORDREF}" IF DEFMTD(1:1) .EQ. "H" THEN SHOW/SESS " Parameters for the method Hough:" SHOW/SESS " Number of orders : NBORDI = {NBORDI}" SHOW/SESS " Width of the orders (pixels) : WIDTHI = {WIDTHI}" SHOW/SESS " Threshold : THRESI = {THRESI}" SHOW/SESS " Parameter SCAN is also taken into account during Hough order definition" SHOW/SESS " Values of NBORDI, WIDTHI and THRESI will be estimated automatically" SHOW/SESS " if these parameters are set to the null value." ELSE SHOW/SESS " Parameters for the methods STD and COM:" SHOW/SESS " Width of the orders (pixels) : WIDTH1 = {WIDTH1}" SHOW/SESS " Threshold : THRES1 = {THRES1}" SHOW/SESS " Approximate slope (unitless) : SLOPE = {SLOPE}" ENDIF ! SHOW/SESS " " IF P1(1:1) .NE. "?" GOTO RET ! BKG: SHOW/SESS "--- (B) Background ---" SHOW/SESS " Method (POLY,SPLINE,SMOOTH) : BKGMTD = {BKGMTD}" SHOW/SESS " Visualisation flag (YES,NO) : BKGVISU = {BKGVISU}" IF BKGMTD(1:1) .EQ. "P" THEN SHOW/SESS " Parameters for method POLY" SHOW/SESS " Step (in pixels) : BKGSTEP = {BKGSTEP}" SHOW/SESS " Polynomial degress (X,Y) : BKGPOL = {BKGPOL(1)},{BKGPOL(2)}" ELSE SHOW/SESS " Radius (X,Y in pixels) : BKGRAD = {BKGRAD(1)},{BKGRAD(2)} " IF BKGMTD(1:2) .EQ. "SP" THEN SHOW/SESS " Step (in pixels) : BKGSTEP = {BKGSTEP}" SHOW/SESS " Spline degreee : BKGDEG = {BKGDEG}" SHOW/SESS " Smoothing factor : BKGSMO = {BKGSMO}" ENDIF IF BKGMTD(1:2) .EQ. "SM" THEN SHOW/SESS " Number of iterations : BKGNIT = {BKGNIT}" ENDIF ENDIF ! SHOW/SESS " " IF P1(1:1) .NE. "?" GOTO RET ! FILT: SHOW/SESS "--- (C) Cosmic Rays Filtering ---" SHOW/SESS " Median filter (X,Y,THRES): MEDFILT= {MEDFILT(1)},{MEDFILT(2)},{MEDFILT(3)}" SHOW/SESS " Cosmics filter (X,Y,Nit) : CRFILT = {CRFILT(1)},{CRFILT(2)},{CRFILT(3)}" SHOW/SESS " CCD Params (ron,g,thres) : CCDFILT= {CCDFILT(1)},{CCDFILT(2)},{CCDFILT(3)}" ! SHOW/SESS " " IF P1(1:1) .NE. "?" GOTO RET ! EXTR: SHOW/SESS "--- (E) Spectrum and Sky Extraction ---" SHOW/SESS " " SHOW/SESS " -- Spectrum Extraction" SHOW/SESS " Extraction method (LINEAR, AVERAGE, OPTIMAL) : EXTMTD = {EXTMTD} SHOW/SESS " Slit Length (in pixels) : SLIT = {SLIT}" IF EXTMTD(1:1) .NE. "O" THEN SHOW/SESS " Offset (in pixels) : OFFSET = {OFFSET}" ELSE SHOW/SESS " Read-Out Noise : RON = {RON}" SHOW/SESS " Gain : GAIN = {GAIN}" SHOW/SESS " Threshold (in std. dev.) : EXTSIGMA = {EXTSIGMA}" ENDIF SHOW/SESS " " SHOW/SESS " -- Sky Definition and Extraction:" SHOW/SESS " Number of sky windows : NSKY = {NSKY}" SHOW/SESS " Sky feature position (,auto,cursor): POSSKY = {POSSKY}" SHOW/SESS " Offset limits of 1st window : SKYWIND(1-2) = {SKYWIND(1)},{SKYWIND(2)}" IF NSKY .EQ. 2 THEN SHOW/SESS " Offset limits of 2nd window : SKYWIND(3-4) = {SKYWIND(3)},{SKYWIND(4)}" ENDIF ! SHOW/SESS " " IF P1(1:1) .NE. "?" GOTO RET ! WLC: SHOW/SESS "--- (W) Wavelength Calibration ---" SHOW/SESS " Search Method (GRAVITY, GAUSSIAN) : SEAMTD = {SEAMTD}" SHOW/SESS " Analysis window width : WIDTH2 = {WIDTH2}" SHOW/SESS " Threshold above the background : THRES2 = {THRES2}" SHOW/SESS " Wavelength calibr. meth. : WLCMTD = {WLCMTD}" SHOW/SESS " (Possible methods: PAIR, ANGLE, GUESS, LINEAR, ORDER, RESTART) SHOW/SESS " Wavelength calibr. option (1D/2D) : WLCOPT = {WLCOPT}" SHOW/SESS " Regress. method (Standard,Robust) : WLCREG = {WLCREG}" SHOW/SESS " Visualisation flag (Yes,No) : WLCVISU = {WLCVISU}" IF WLCMTD(1:1) .EQ. "G" THEN SHOW/SESS " Method Guess. Guess Session : GUESS = {GUESS}" ENDIF SHOW/SESS " Arc frame : WLC = {WLC}" SHOW/SESS " Line Catalog : LINCAT = {LINCAT}" SHOW/SESS " Polynomial Degree : DC = {DC}" SHOW/SESS " Tolerance : TOL = {TOL}" SHOW/SESS "NB: Parameter CCDBIN is also involved" SHOW/SESS " Mini, Maxi number of iterations : WLCNITER = {WLCNITER(1)},{WLCNITER(2)}" SHOW/SESS " Iteration Loop parameters:" SHOW/SESS " WLCITER(1-3) = {WLCITER(1)}, {WLCITER(2)}, {WLCITER(3)} SHOW/SESS " WLCITER(4-5) = {WLCITER(4)}, {WLCITER(5)}" ! SHOW/SESS " " IF P1(1:1) .NE. "?" GOTO RET ! REBIN: SHOW/SESS "--- (R) Spectrum Resampling---" SHOW/SESS " Resampling Option (Yes/No) : REBOPT = {REBOPT}" SHOW/SESS " Resampling Method : REBMTD = {REBMTD}" SHOW/SESS " Resampling step or reference frame: SAMPLE = {SAMPLE}" ! SHOW/SESS " " IF P1(1:1) .NE. "?" GOTO RET ! FLAT: SHOW/SESS "--- (F) Flat Field correction ---" SHOW/SESS " Flat-Field option (Yes,No) : FFOPT = {FFOPT}" IF FFOPT(1:1) .EQ. "Y" THEN SHOW/SESS " Flat-Field frame : FLAT = {FLAT}" SHOW/SESS " Background corrected flat-field : CORRECT = {CORRECT}" SHOW/SESS " Blaze function : BLAZE = {BLAZE}" ENDIF SHOW/SESS " " ! IF P1(1:1) .NE. "?" GOTO RET ! FLUX: SHOW/SESS "--- (F) Flux Calibration ---" SHOW/SESS " Response Correction Option (Yes, No) : RESPOPT = {RESPOPT}" IF RESPOPT(1:1) .EQ. "Y" THEN SHOW/SESS " Response Corection Method : RESPMTD = {RESPMTD}" SHOW/SESS " (Possible methods: STD, IUE, FLAT)" IF RESPMTD(1:1) .EQ. "S" THEN SHOW/SESS " Standard Star : STD = {STD}" SHOW/SESS " Response Table : RESPONSE = {RESPONSE}" SHOW/SESS " Flux Table : FLUXTAB = {FLUXTAB}" SHOW/SESS " Median Filter param. (X,Y,thres): FILTMED = {FILTMED(1)}," SHOW/SESS " Smooth filter param. (X,Y,THRES): FILTSMO = {FILTSMO(1)}" SHOW/SESS " Valid pixel range : PIXNUL = {PIXNUL(1)}, {PIXNUL(2)}" ENDIF IF RESPMTD(1:1) .EQ. "I" THEN SHOW/SESS " RIPMTD = {RIPMTD} RIPFRZ = {RIPFRZ}" IF RIPMTD(1:1) .NE. "S" THEN SHOW/SESS " LAMBDA1 = {LAMBDA1} LAMBDA2 = {LAMBDA2}" ENDIF SHOW/SESS " RIPK = {RIPK} ALPHA = {ALPHA}" ENDIF IF RESPMTD(1:1) .EQ. "M" THEN SHOW/SESS " OBSTAB = {OBSTAB} RESPONSE = {RESPONSE}" SHOW/SESS " RIPORD = {RIPORD} RIPRAD = {RIPRAD(1)},{RIPRAD(2)}" SHOW/SESS " RIPPIX = {RIPPIX} RIPREF = {RIPREF}" SHOW/SESS " RIPBIN = {RIPBIN}" ENDIF ENDIF ! SHOW/SESS " " IF P1(1:1) .NE. "?" GOTO RET ! MERGE: SHOW/SESS "--- (M) Merging ---" SHOW/SESS " Orders merging option (Yes, No) : MGOPT = {MGOPT}" show/sess " Order merging method (Aver, Noapp) : MRGMTD = {MRGMTD}" IF MGOPT(1:1) .EQ. "Y" THEN IF MRGMTD(1:1) .EQ. "N" THEN SHOW/SESS " First, last order number : MRGORD = {MRGORD(1)},{MRGORD(2)}" ELSE SHOW/SESS " Zeroed interval (wavel. units) : DELTA = {DELTA}" ENDIF ENDIF ! SHOW/SESS " " IF P1(1:1) .NE. "?" GOTO RET ! RET: SET/FORMAT RETURN HELP: SHOW/SESS " **** Echelle Main Help Message ****" SHOW/SESS " " SHOW/SESS " This is the list of high level commands useful for an echelle" SHOW/SESS " reduction in the sequential order you will need them. Optional" SHOW/SESS " steps are indicated by (opt.). Entries in brackets [] indicate" SHOW/SESS " the entries of SHOW/ECHELLE which must be checked before running" SHOW/SESS " the command. " SHOW/SESS " More help is available for each command with HELP " SHOW/SESS " " SHOW/SESS " INTAPE/FITS" SHOW/SESS " ROTATE/ECHELLE" SHOW/SESS " PREPARE/ECHELLE (opt.)" SHOW/SESS " SHOW/ECHELLE" SHOW/SESS " HELP/ECHELLE" SHOW/SESS " SET/ECHELLE" SHOW/SESS " CALIBRATE/ECHELLE [Ord. Def.] [Extr.] [Search.] [Wav. cal.]" SHOW/SESS " SAVE/ECHELLE" SHOW/SESS " INIT/ECHELLE" SHOW/SESS " SELECT/BACKGROUND (opt.)" SHOW/SESS " FILTER/ECHELLE (opt.) [Background] [Filtering]" SHOW/SESS " FLAT/ECHELLE (opt.) [Flat-field] [Backgr.] [Extr.] [Rebin]" SHOW/SESS " RESPONSE/ECHELLE (opt.) [Backgr.] [Extr.] [Rebin] [Response]" SHOW/SESS " REPEAT/ECHELLE (opt.)" SHOW/SESS " EDIT/ECHELLE (opt.)" SHOW/SESS " REDUCE/ECHELLE [Backgr.] [Extr.] [Rebin] [Response] [Merging]" SHOW/SESS " " SHOW/SESS "TUTORIAL/ECH demonstrates the use of the different commands" SHOW/SESS "SHOW/ECHELLE display a list of the parameters to be set (SET/ECH" SHOW/SESS " or direct assignment. The list depends on the chosen method" SHOW/SESS "HELP/ECH method provide a list of all methods available" SHOW/SESS " for the different steps of the calibration and reduction" SHOW/SESS "HELP/ECH option provide a list of all available options" RETURN