#!/bin/bash -xv # the script processes a set of DARK frames. # The images are overscan corrected and stacked without # any rescaling. This script is the parallel version # and uses the 'preprocess' program to perform the # reduction. In contrary to the script version based # on FLIPS, the individual preprocessed images are NOT # saved to disk but only the final coadded calibration # frames. # It assumes that all images from ONE chip reside # in the same physical directory # $1: master directory # $2: subdirectory with the DARK images # $3: chips to work on # The result images are in the $1/$2 directory and have # the same basic name as the subdirectory, f.I. # The images are in .../DARK. Then the stacked # images have the names .../DARK/DARK_i.fits, where i=1..NCHIPS. # # 27.04.2005: # I rewrote the script to use the MAXIMAGE parameter in preprocess # and the imcat-based imcombflat for the actual coaddition. # preliminary work: INSTDIR=instruments_professional + INSTDIR=instruments_professional if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=instruments_commercial fi + '[' '!' -f instruments_professional/Moravian10z.ini ']' + INSTDIR=instruments_commercial if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=~/.theli/instruments_user/ fi + '[' '!' -f instruments_commercial/Moravian10z.ini ']' + INSTDIR=/home/albireo/.theli/instruments_user/ . ${INSTDIR}/${INSTRUMENT:?}.ini + . /home/albireo/.theli/instruments_user//Moravian10z.ini # ************************************************************** # Instrument config file created by THELI # ************************************************************** # Name of the instrument (must be unique!) INSTRUMENT=Moravian10z ++ INSTRUMENT=Moravian10z # Number of chips NCHIPS=1 ++ NCHIPS=1 # We need to export these export INSTRUMENT ++ export INSTRUMENT export NCHIPS ++ export NCHIPS # Load directories for programs . progs.ini ++ . progs.ini # If we include this script in bash verbose and trace mode we do not want # all the (verbose) output from this script: USE_X=0 +++ USE_X=0 if [[ ${-/x} != $- ]] ; then USE_X=1 set +xv fi +++ [[ hvB != hvxB ]] +++ USE_X=1 +++ set +xv # # this is needed for the GUI # P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh +++ P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh P_GETPMMOBJECTS=${BIN}/get_pmm_objects +++ P_GETPMMOBJECTS=/home/albireo/THELI/theli//bin/Linux_64//get_pmm_objects P_LDACDESC=${BIN}/ldacdesc +++ P_LDACDESC=/home/albireo/THELI/theli//bin/Linux_64//ldacdesc P_LDACDELTAB=${BIN}/ldacdeltab +++ P_LDACDELTAB=/home/albireo/THELI/theli//bin/Linux_64//ldacdeltab P_LDACDELKEY=${BIN}/ldacdelkey +++ P_LDACDELKEY=/home/albireo/THELI/theli//bin/Linux_64//ldacdelkey P_LDACPUTXY=${BIN}/ldacputxy +++ P_LDACPUTXY=/home/albireo/THELI/theli//bin/Linux_64//ldacputxy P_FITSFLIP=${BIN}/fitsflip +++ P_FITSFLIP=/home/albireo/THELI/theli//bin/Linux_64//fitsflip P_DEBLOOM=${BIN}/debloom +++ P_DEBLOOM=/home/albireo/THELI/theli//bin/Linux_64//debloom P_EQUALISEBAYERFLAT=${BIN}/equalisebayerflat +++ P_EQUALISEBAYERFLAT=/home/albireo/THELI/theli//bin/Linux_64//equalisebayerflat P_FITSCOLLAPSE=${BIN}/fitscollapse +++ P_FITSCOLLAPSE=/home/albireo/THELI/theli//bin/Linux_64//fitscollapse P_FITSWRAP=${BIN}/fitswrap +++ P_FITSWRAP=/home/albireo/THELI/theli//bin/Linux_64//fitswrap P_FITSCUBESLICE=${BIN}/fitscubeslice +++ P_FITSCUBESLICE=/home/albireo/THELI/theli//bin/Linux_64//fitscubeslice P_FITSADDKEY=${BIN}/fitsaddkey +++ P_FITSADDKEY=/home/albireo/THELI/theli//bin/Linux_64//fitsaddkey P_FITSCUT=${BIN}/fitscut +++ P_FITSCUT=/home/albireo/THELI/theli//bin/Linux_64//fitscut P_FITSBIN=${BIN}/fitsbin +++ P_FITSBIN=/home/albireo/THELI/theli//bin/Linux_64//fitsbin P_FITSBLOOMDETECT=${BIN}/fitsbloomdetect +++ P_FITSBLOOMDETECT=/home/albireo/THELI/theli//bin/Linux_64//fitsbloomdetect P_FITSFILL=${BIN}/fitsfill +++ P_FITSFILL=/home/albireo/THELI/theli//bin/Linux_64//fitsfill P_FITSGRAD=${BIN}/fitsgrad +++ P_FITSGRAD=/home/albireo/THELI/theli//bin/Linux_64//fitsgrad P_FITSDEMOSAICBAYER=${BIN}/fitsdemosaicbayer +++ P_FITSDEMOSAICBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitsdemosaicbayer P_FITS3MAX=${BIN}/fits3max +++ P_FITS3MAX=/home/albireo/THELI/theli//bin/Linux_64//fits3max P_FITSGAUSS=${BIN}/fitsgauss +++ P_FITSGAUSS=/home/albireo/THELI/theli//bin/Linux_64//fitsgauss P_FITS2TIFF=${BIN}/fits2tiff +++ P_FITS2TIFF=/home/albireo/THELI/theli//bin/Linux_64//fits2tiff P_TIFF2FITS=${BIN}/tiff2fits +++ P_TIFF2FITS=/home/albireo/THELI/theli//bin/Linux_64//tiff2fits P_CORRLIRIS=${BIN}/corrliris +++ P_CORRLIRIS=/home/albireo/THELI/theli//bin/Linux_64//corrliris P_CROSSTALK=${BIN}/crosstalk +++ P_CROSSTALK=/home/albireo/THELI/theli//bin/Linux_64//crosstalk P_FINDGSC=${BIN}/findgsc +++ P_FINDGSC=/home/albireo/THELI/theli//bin/Linux_64//findgsc P_FINDUSNOB1=${BIN}/findusnob1 +++ P_FINDUSNOB1=/home/albireo/THELI/theli//bin/Linux_64//findusnob1 P_FIND2MASS=${BIN}/find2mass +++ P_FIND2MASS=/home/albireo/THELI/theli//bin/Linux_64//find2mass P_FINDSDSS7=${BIN}/findsdss7 +++ P_FINDSDSS7=/home/albireo/THELI/theli//bin/Linux_64//findsdss7 P_FINDUCAC3=${BIN}/finducac3 +++ P_FINDUCAC3=/home/albireo/THELI/theli//bin/Linux_64//finducac3 P_FINDNOMAD1=${BIN}/findnomad1 +++ P_FINDNOMAD1=/home/albireo/THELI/theli//bin/Linux_64//findnomad1 P_VIZQUERY=${BIN}/vizquery +++ P_VIZQUERY=/home/albireo/THELI/theli//bin/Linux_64//vizquery P_SESAME=${BIN}/sesame +++ P_SESAME=/home/albireo/THELI/theli//bin/Linux_64//sesame P_ACLIENT=${BIN}/aclient +++ P_ACLIENT=/home/albireo/THELI/theli//bin/Linux_64//aclient P_SCAMP=${BIN}/scamp +++ P_SCAMP=/home/albireo/THELI/theli//bin/Linux_64//scamp P_RENAME=${BIN}/rename +++ P_RENAME=/home/albireo/THELI/theli//bin/Linux_64//rename P_MISSFITS=${BIN}/missfits +++ P_MISSFITS=/home/albireo/THELI/theli//bin/Linux_64//missfits P_SWARPFILTER=${BIN}/swarpfilter +++ P_SWARPFILTER=/home/albireo/THELI/theli//bin/Linux_64//swarpfilter P_FITSBLOCKEDIT=${BIN}/fitsblockedit +++ P_FITSBLOCKEDIT=/home/albireo/THELI/theli//bin/Linux_64//fitsblockedit P_FITSMASK=${BIN}/fitsmask +++ P_FITSMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsmask P_FITSMEDIAN=${BIN}/fitsmedian +++ P_FITSMEDIAN=/home/albireo/THELI/theli//bin/Linux_64//fitsmedian P_FITSNANMASK=${BIN}/fitsnanmask +++ P_FITSNANMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsnanmask P_FITSSTAT=${BIN}/fitsstat +++ P_FITSSTAT=/home/albireo/THELI/theli//bin/Linux_64//fitsstat P_FITSSPLITBAYER=${BIN}/fitssplitbayer +++ P_FITSSPLITBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitssplitbayer P_FITSGUESSGEOM=${BIN}/fitsguessgeom +++ P_FITSGUESSGEOM=/home/albireo/THELI/theli//bin/Linux_64//fitsguessgeom P_SPLITSUPRIMECAM=${BIN}/splitsuprimecam +++ P_SPLITSUPRIMECAM=/home/albireo/THELI/theli//bin/Linux_64//splitsuprimecam P_PASTEDUALCHANNEL=${BIN}/paste_dualchannel +++ P_PASTEDUALCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_dualchannel P_PASTEQUADCHANNEL=${BIN}/paste_quadchannel +++ P_PASTEQUADCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_quadchannel P_SPLITGMOSMULTIPORT=${BIN}/split_gmos_multiport +++ P_SPLITGMOSMULTIPORT=/home/albireo/THELI/theli//bin/Linux_64//split_gmos_multiport P_SPLITGROND=${BIN}/split_grond +++ P_SPLITGROND=/home/albireo/THELI/theli//bin/Linux_64//split_grond P_GETSHIFT=${BIN}/get_shift +++ P_GETSHIFT=/home/albireo/THELI/theli//bin/Linux_64//get_shift P_FITSNONLINEARITY=${BIN}/fitsnonlinearity +++ P_FITSNONLINEARITY=/home/albireo/THELI/theli//bin/Linux_64//fitsnonlinearity P_FITSPIC=${BIN}/fitspic +++ P_FITSPIC=/home/albireo/THELI/theli//bin/Linux_64//fitspic P_FITSPOLYGON=${BIN}/fitspolygon +++ P_FITSPOLYGON=/home/albireo/THELI/theli//bin/Linux_64//fitspolygon P_EXPANDSEXMASK=${BIN}/expand_sexmask +++ P_EXPANDSEXMASK=/home/albireo/THELI/theli//bin/Linux_64//expand_sexmask P_DCRAW=${BIN}/dcraw +++ P_DCRAW=/home/albireo/THELI/theli//bin/Linux_64//dcraw P_DECSEXCONV=${BIN}/decsexconv +++ P_DECSEXCONV=/home/albireo/THELI/theli//bin/Linux_64//decsexconv P_STIFF=${BIN}/stiff +++ P_STIFF=/home/albireo/THELI/theli//bin/Linux_64//stiff P_SUBSKY=${BIN}/subsky +++ P_SUBSKY=/home/albireo/THELI/theli//bin/Linux_64//subsky P_GIFSICLE=${BIN}/gifsicle +++ P_GIFSICLE=/home/albireo/THELI/theli//bin/Linux_64//gifsicle P_XY2SKY=${BIN}/xy2sky +++ P_XY2SKY=/home/albireo/THELI/theli//bin/Linux_64//xy2sky P_SKY2XY=${BIN}/sky2xy +++ P_SKY2XY=/home/albireo/THELI/theli//bin/Linux_64//sky2xy P_ISFITS=${BIN}/isfits +++ P_ISFITS=/home/albireo/THELI/theli//bin/Linux_64//isfits P_GETKEY=${BIN}/getkey +++ P_GETKEY=/home/albireo/THELI/theli//bin/Linux_64//getkey P_SKY2XY_HELPER=${BIN}/sky2xy_helper +++ P_SKY2XY_HELPER=/home/albireo/THELI/theli//bin/Linux_64//sky2xy_helper P_GET_ROTIMSIZE=${BIN}/get_rotimsize +++ P_GET_ROTIMSIZE=/home/albireo/THELI/theli//bin/Linux_64//get_rotimsize P_CLEANSTRING=${BIN}/cleanstring +++ P_CLEANSTRING=/home/albireo/THELI/theli//bin/Linux_64//cleanstring P_NUMCPU=${BIN}/numcpu +++ P_NUMCPU=/home/albireo/THELI/theli//bin/Linux_64//numcpu P_GETSTATS=${BIN}/getstats +++ P_GETSTATS=/home/albireo/THELI/theli//bin/Linux_64//getstats P_CDELT2CD=${BIN}/cdelt2cd +++ P_CDELT2CD=/home/albireo/THELI/theli//bin/Linux_64//cdelt2cd P_GETCDMATRIX=${BIN}/get_cdmatrix +++ P_GETCDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//get_cdmatrix P_GETPIXSCALE=${BIN}/get_pixscale +++ P_GETPIXSCALE=/home/albireo/THELI/theli//bin/Linux_64//get_pixscale P_GETPOSANGLE=${BIN}/get_posangle +++ P_GETPOSANGLE=/home/albireo/THELI/theli//bin/Linux_64//get_posangle P_GETREFCATRADIUS=${BIN}/get_refcat_radius +++ P_GETREFCATRADIUS=/home/albireo/THELI/theli//bin/Linux_64//get_refcat_radius P_MERGEDAT=${BIN}/mergedat +++ P_MERGEDAT=/home/albireo/THELI/theli//bin/Linux_64//mergedat P_ROTATECDMATRIX=${BIN}/rotate_cdmatrix +++ P_ROTATECDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//rotate_cdmatrix P_BACKGROUNDDYNAMIC=${BIN}/background_dynamic +++ P_BACKGROUNDDYNAMIC=/home/albireo/THELI/theli//bin/Linux_64//background_dynamic P_BACKGROUNDSTATIC=${BIN}/background_static +++ P_BACKGROUNDSTATIC=/home/albireo/THELI/theli//bin/Linux_64//background_static P_MJD=${BIN}/mjd +++ P_MJD=/home/albireo/THELI/theli//bin/Linux_64//mjd P_LST=${BIN}/lst +++ P_LST=/home/albireo/THELI/theli//bin/Linux_64//lst P_IMALYZER=${BIN}/imalyzer +++ P_IMALYZER=/home/albireo/THELI/theli//bin/Linux_64//imalyzer P_E2ELL=${BIN}/e2ell +++ P_E2ELL=/home/albireo/THELI/theli//bin/Linux_64//e2ell P_FITSSMOOTHEDGE=${BIN}/fitssmoothedge +++ P_FITSSMOOTHEDGE=/home/albireo/THELI/theli//bin/Linux_64//fitssmoothedge P_FITSEXTRACTEXT=${BIN}/fitsextractext +++ P_FITSEXTRACTEXT=/home/albireo/THELI/theli//bin/Linux_64//fitsextractext P_ERRTEST=${BIN}/errtest +++ P_ERRTEST=/home/albireo/THELI/theli//bin/Linux_64//errtest P_BUILDINDEX=${BIN_ANET}/build-astrometry-index +++ P_BUILDINDEX=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//build-astrometry-index P_SOLVEFIELD=${BIN_ANET}/solve-field +++ P_SOLVEFIELD=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//solve-field P_XCORR=${BIN}/xcorr +++ P_XCORR=/home/albireo/THELI/theli//bin/Linux_64//xcorr # # load the parameter files # . ~/.theli/param_set1.ini +++ . /home/albireo/.theli/param_set1.ini PROJECTNAME=Pferdekopf ++++ PROJECTNAME=Pferdekopf NPARA=1 ++++ NPARA=1 GUIVERSION=2.9.5 ++++ GUIVERSION=2.9.5 KERNEL=Linux ++++ KERNEL=Linux NFRAMES=47 ++++ NFRAMES=47 V_PRE_SPLITMIRCUBE=N ++++ V_PRE_SPLITMIRCUBE=N V_DO_BIAS=N ++++ V_DO_BIAS=N V_DO_FLAT=Y ++++ V_DO_FLAT=Y V_NONLINCORR=N ++++ V_NONLINCORR=N V_AP_SERVER=vizier.u-strasbg.fr ++++ V_AP_SERVER=vizier.u-strasbg.fr V_AP_MAGLIM=20 ++++ V_AP_MAGLIM=20 V_AP_RADIUS=30 ++++ V_AP_RADIUS=30 V_AP_REFRA=header ++++ V_AP_REFRA=header V_AP_REFDEC=header ++++ V_AP_REFDEC=header . ~/.theli/param_set2.ini +++ . /home/albireo/.theli/param_set2.ini V_PRE_RENAME_CHECKED=0 ++++ V_PRE_RENAME_CHECKED=0 V_RENAME_FITSKEY=ARCFILE ++++ V_RENAME_FITSKEY=ARCFILE V_SORT_FITSKEY=OBJECT ++++ V_SORT_FITSKEY=OBJECT V_SORT_BIAS= ++++ V_SORT_BIAS= V_SORT_DARK= ++++ V_SORT_DARK= V_SORT_DOMEFLAT= ++++ V_SORT_DOMEFLAT= V_SORT_SKYFLAT= ++++ V_SORT_SKYFLAT= V_SORT_STD= ++++ V_SORT_STD= V_PRE_XTALK_NOR_CHECKED=0 ++++ V_PRE_XTALK_NOR_CHECKED=0 V_PRE_XTALK_ROW_CHECKED=0 ++++ V_PRE_XTALK_ROW_CHECKED=0 V_PRE_XTALK_MULTI_CHECKED=0 ++++ V_PRE_XTALK_MULTI_CHECKED=0 V_PRE_XTALK_NOR_AMPLITUDE= ++++ V_PRE_XTALK_NOR_AMPLITUDE= V_PRE_XTALK_ROW_AMPLITUDE= ++++ V_PRE_XTALK_ROW_AMPLITUDE= V_PRE_XTALK_NOR_BUTTONID=-1 ++++ V_PRE_XTALK_NOR_BUTTONID=-1 V_PRE_XTALK_ROW_BUTTONID=-1 ++++ V_PRE_XTALK_ROW_BUTTONID=-1 V_PRE_XTALK_MULTI_BUTTONID=-1 ++++ V_PRE_XTALK_MULTI_BUTTONID=-1 V_PRE_XTALK_MULTI_NSECTION= ++++ V_PRE_XTALK_MULTI_NSECTION= V_CAL_OVSCANNLOW=0 ++++ V_CAL_OVSCANNLOW=0 V_CAL_OVSCANNHIGH=1 ++++ V_CAL_OVSCANNHIGH=1 V_CAL_BIASNLOW=0 ++++ V_CAL_BIASNLOW=0 V_CAL_BIASNHIGH=1 ++++ V_CAL_BIASNHIGH=1 V_CAL_DARKNLOW=0 ++++ V_CAL_DARKNLOW=0 V_CAL_DARKNHIGH=1 ++++ V_CAL_DARKNHIGH=1 V_CAL_FLATNLOW=0 ++++ V_CAL_FLATNLOW=0 V_CAL_FLATNHIGH=1 ++++ V_CAL_FLATNHIGH=1 V_CAL_USEDARK=FALSE ++++ V_CAL_USEDARK=FALSE V_BACK_NLOW=0 ++++ V_BACK_NLOW=0 V_BACK_NHIGH=1 ++++ V_BACK_NHIGH=1 V_BACK_DETECTTHRESH=1.3 ++++ V_BACK_DETECTTHRESH=1.3 V_BACK_DETECTMINAREA=5 ++++ V_BACK_DETECTMINAREA=5 V_BACK_MASKEXPAND= ++++ V_BACK_MASKEXPAND= V_BACK_ILLUMSMOOTH= ++++ V_BACK_ILLUMSMOOTH= V_BACK_FRINGESMOOTH= ++++ V_BACK_FRINGESMOOTH= V_BACK_WINDOWSIZE=0 ++++ V_BACK_WINDOWSIZE=0 V_BACK_GAPSIZE=1.0 ++++ V_BACK_GAPSIZE=1.0 V_BACK_APPLYMODE=0 ++++ V_BACK_APPLYMODE=0 V_BACK_COMBINEMETHOD=0 ++++ V_BACK_COMBINEMETHOD=0 V_BACK_SEXFILTER=Y ++++ V_BACK_SEXFILTER=Y V_BACK_ADJUSTGAINS=N ++++ V_BACK_ADJUSTGAINS=N V_BACK_FRINGESCALE=Y ++++ V_BACK_FRINGESCALE=Y V_BACK_SKY2PASS=Y ++++ V_BACK_SKY2PASS=Y V_COLLDETECTTHRESH= ++++ V_COLLDETECTTHRESH= V_COLLDETECTMINAREA= ++++ V_COLLDETECTMINAREA= V_COLLMASKEXPAND= ++++ V_COLLMASKEXPAND= V_COLLREJECTTHRESH=1.5 ++++ V_COLLREJECTTHRESH=1.5 V_COLLDIRECTION=x ++++ V_COLLDIRECTION=x V_COLLXMIN= ++++ V_COLLXMIN= V_COLLXMAX= ++++ V_COLLXMAX= V_COLLYMIN= ++++ V_COLLYMIN= V_COLLYMAX= ++++ V_COLLYMAX= V_COLLMASKACTION=1 ++++ V_COLLMASKACTION=1 V_COLLAUTOTHRESHOLD=0 ++++ V_COLLAUTOTHRESHOLD=0 V_WEIGHTBINMIN=-100 ++++ V_WEIGHTBINMIN=-100 V_WEIGHTBINMAX=500 ++++ V_WEIGHTBINMAX=500 V_WEIGHTBINSIZE=4 ++++ V_WEIGHTBINSIZE=4 V_WEIGHTLOWTHRESHOLD= ++++ V_WEIGHTLOWTHRESHOLD= V_WEIGHTHIGHTHRESHOLD= ++++ V_WEIGHTHIGHTHRESHOLD= V_GLOBWFLATMIN=0.5 ++++ V_GLOBWFLATMIN=0.5 V_GLOBWFLATMAX=1.6 ++++ V_GLOBWFLATMAX=1.6 V_GLOBWDARKMIN= ++++ V_GLOBWDARKMIN= V_GLOBWDARKMAX= ++++ V_GLOBWDARKMAX= V_GLOBWDARKDIR= ++++ V_GLOBWDARKDIR= V_COSMICSTHRESHOLD=0.1 ++++ V_COSMICSTHRESHOLD=0.1 V_COSMICDT=6 ++++ V_COSMICDT=6 V_COSMICDMIN=1 ++++ V_COSMICDMIN=1 V_DEFECT_KERNELSIZE= ++++ V_DEFECT_KERNELSIZE= V_DEFECT_ROWTOL= ++++ V_DEFECT_ROWTOL= V_DEFECT_COLTOL= ++++ V_DEFECT_COLTOL= V_DEFECT_CLUSTOL= ++++ V_DEFECT_CLUSTOL= V_DEFECT_KERNELSIZE_SF= ++++ V_DEFECT_KERNELSIZE_SF= V_DEFECT_ROWTOL_SF= ++++ V_DEFECT_ROWTOL_SF= V_DEFECT_COLTOL_SF= ++++ V_DEFECT_COLTOL_SF= V_DEFECT_CLUSTOL_SF= ++++ V_DEFECT_CLUSTOL_SF= V_BLOOMLOWLIMIT=20000 ++++ V_BLOOMLOWLIMIT=20000 V_BLOOMMINAREA=100 ++++ V_BLOOMMINAREA=100 V_BLOOMWIDTH=0 ++++ V_BLOOMWIDTH=0 V_WEIGHT_BINOUTLIER=FALSE ++++ V_WEIGHT_BINOUTLIER=FALSE V_GLOBW_UNIFORMWEIGHT=FALSE ++++ V_GLOBW_UNIFORMWEIGHT=FALSE V_MASKBLOOMSPIKE=0 ++++ V_MASKBLOOMSPIKE=0 V_AP_DETECTTHRESH=5.0 ++++ V_AP_DETECTTHRESH=5.0 V_AP_DETECTMINAREA=5 ++++ V_AP_DETECTMINAREA=5 V_DEBLENDMINCONT=0.0005 ++++ V_DEBLENDMINCONT=0.0005 V_AP_LOWNUM=1 ++++ V_AP_LOWNUM=1 V_SEXCATMINFWHM=1.5 ++++ V_SEXCATMINFWHM=1.5 V_SEXCATFLAG=0 ++++ V_SEXCATFLAG=0 V_SATURATION= ++++ V_SATURATION= V_SEXBACKLEVEL= ++++ V_SEXBACKLEVEL= V_AP_FILTER=N ++++ V_AP_FILTER=N V_SCAMP_MAXPOSANGLE=2.0 ++++ V_SCAMP_MAXPOSANGLE=2.0 V_SCAMP_MAXOFFSET=2.0 ++++ V_SCAMP_MAXOFFSET=2.0 V_SCAMP_MAXSCALE=1.05 ++++ V_SCAMP_MAXSCALE=1.05 V_SCAMP_SNLOW=5 ++++ V_SCAMP_SNLOW=5 V_SCAMP_SNHIGH=20 ++++ V_SCAMP_SNHIGH=20 V_SCAMP_FGROUPRADIUS=1.0 ++++ V_SCAMP_FGROUPRADIUS=1.0 V_SCAMP_CROSSIDRADIUS=0 ++++ V_SCAMP_CROSSIDRADIUS=0 V_SCAMP_ASTREFWEIGHT=1.0 ++++ V_SCAMP_ASTREFWEIGHT=1.0 V_SCAMP_ASTRINSTRUKEY=FILTER ++++ V_SCAMP_ASTRINSTRUKEY=FILTER V_SCAMP_PHOTINSTRUKEY=FILTER ++++ V_SCAMP_PHOTINSTRUKEY=FILTER V_SCAMP_POLYORDER=3 ++++ V_SCAMP_POLYORDER=3 V_SCAMP_POLYORDER2= ++++ V_SCAMP_POLYORDER2= V_SCAMP_DISTORTGROUPS= ++++ V_SCAMP_DISTORTGROUPS= V_SCAMP_DISTORTKEYS= ++++ V_SCAMP_DISTORTKEYS= V_SCAMP_STABILITY=INSTRUMENT ++++ V_SCAMP_STABILITY=INSTRUMENT V_SCAMP_MOSAICTYPE=UNCHANGED ++++ V_SCAMP_MOSAICTYPE=UNCHANGED V_SCAMP_FOCALPLANE=DEFAULT ++++ V_SCAMP_FOCALPLANE=DEFAULT V_SCAMP_MATCHFLIPPED=N ++++ V_SCAMP_MATCHFLIPPED=N V_SCAMP_MATCHWCS=Y ++++ V_SCAMP_MATCHWCS=Y V_ANET_MAXSCALE=1.05 ++++ V_ANET_MAXSCALE=1.05 V_ANET_POLYORDER=3 ++++ V_ANET_POLYORDER=3 V_ANET_SCAMPDISTORT=N ++++ V_ANET_SCAMPDISTORT=N V_ABSPHOT_FILTER=B ++++ V_ABSPHOT_FILTER=B V_ABSPHOT_COLORFIXED=0.1 ++++ V_ABSPHOT_COLORFIXED=0.1 V_ABSPHOT_EXTINCTION=0 ++++ V_ABSPHOT_EXTINCTION=0 V_ABSPHOT_ZPMIN=24 ++++ V_ABSPHOT_ZPMIN=24 V_ABSPHOT_ZPMAX=25 ++++ V_ABSPHOT_ZPMAX=25 V_ABSPHOT_STDMINMAG=0.0 ++++ V_ABSPHOT_STDMINMAG=0.0 V_ABSPHOT_MAXPHOTERR=0.05 ++++ V_ABSPHOT_MAXPHOTERR=0.05 V_ABSPHOT_WCSERR=10 ++++ V_ABSPHOT_WCSERR=10 V_ABSPHOT_WCSERRCHECKBOX=Y ++++ V_ABSPHOT_WCSERRCHECKBOX=Y V_ABSPHOT_STDCAT_DIRECT=SDSS ++++ V_ABSPHOT_STDCAT_DIRECT=SDSS V_ABSPHOT_STDFILTERM2=u ++++ V_ABSPHOT_STDFILTERM2=u V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI ++++ V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI V_ABSPHOT_STDFILTER=U ++++ V_ABSPHOT_STDFILTER=U V_ABSPHOT_STDCOLOR=UmB ++++ V_ABSPHOT_STDCOLOR=UmB V_ABSPHOT_CALMODE=RUNCALIB ++++ V_ABSPHOT_CALMODE=RUNCALIB V_ABSPHOT_FITTINGMETHODM2=CHIP ++++ V_ABSPHOT_FITTINGMETHODM2=CHIP V_COADD_REFRA= ++++ V_COADD_REFRA= V_COADD_REFDEC= ++++ V_COADD_REFDEC= V_COADD_IDENT=(null) ++++ V_COADD_IDENT=(null) V_COADD_SEEING= ++++ V_COADD_SEEING= V_COADD_RZP= ++++ V_COADD_RZP= V_COADD_PIXSCALE=0.89 ++++ V_COADD_PIXSCALE=0.89 V_COADD_SIZEX= ++++ V_COADD_SIZEX= V_COADD_SIZEY= ++++ V_COADD_SIZEY= V_COADD_SKYPOSANGLE= ++++ V_COADD_SKYPOSANGLE= V_COADD_CHIPS="" ++++ V_COADD_CHIPS= V_COADD_FILTERTHRESHOLD= ++++ V_COADD_FILTERTHRESHOLD= V_COADD_FILTERCLUSTERSIZE= ++++ V_COADD_FILTERCLUSTERSIZE= V_COADD_FILTERBORDERWIDTH= ++++ V_COADD_FILTERBORDERWIDTH= V_COADD_SMOOTHEDGE= ++++ V_COADD_SMOOTHEDGE= V_COADD_PROPMOTIONRA= ++++ V_COADD_PROPMOTIONRA= V_COADD_PROPMOTIONDEC= ++++ V_COADD_PROPMOTIONDEC= V_COADD_KERNEL=LANCZOS3 ++++ V_COADD_KERNEL=LANCZOS3 V_COADD_FILTER= ++++ V_COADD_FILTER= V_COADD_PROJECTION=TAN ++++ V_COADD_PROJECTION=TAN V_COADD_CELESTIALTYPE=EQUATORIAL ++++ V_COADD_CELESTIALTYPE=EQUATORIAL V_COADD_COMBINETYPE=WEIGHTED ++++ V_COADD_COMBINETYPE=WEIGHTED . ~/.theli/param_set3.ini +++ . /home/albireo/.theli/param_set3.ini # Geographical latitude of the observatory OBSLAT= ++ OBSLAT= # Pixel scale in x and y [degrees] (CDELT1/2 keywords) PIXSCX=-0 ++ PIXSCX=-0 PIXSCY=0 ++ PIXSCY=0 # Pixel scale in arcsec PIXSCALE= ++ PIXSCALE= # Detector gain GAIN=1.0 ++ GAIN=1.0 # Overscan regions. The values are the x-range. # The min is given by OVSCANX1, the max by OVSCANX2. # We need separate regions for each chip. OVSCANX1=([1]=0) ++ OVSCANX1=([1]=0) OVSCANX2=([1]=0) ++ OVSCANX2=([1]=0) # Sections for the cutting. Select those pixels that receive light. # Give min values for x and y, and the extent in x and y. CUTX=([1]=1) ++ CUTX=([1]=1) CUTY=([1]=1) ++ CUTY=([1]=1) SIZEX=([1]=3358) ++ SIZEX=([1]=3358) SIZEY=([1]=2536) ++ SIZEY=([1]=2536) # Reference pixels for each chip. REFPIXX=([1]=100) ++ REFPIXX=([1]=100) REFPIXY=([1]=100) ++ REFPIXY=([1]=100) # Camera type TYPE=OPT ++ TYPE=OPT # # the resultdir is where the output coadded images # will go. If ONE image of the corresponding chip # is a link the image will go to THIS directory for CHIP in $3 do RESULTDIR[${CHIP}]="$1/$2" done + for CHIP in '$3' + RESULTDIR[${CHIP}]=/home/albireo/14-17.02.2015/rgb//dark . bash.include + . bash.include function get_wcs { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do if [ "${NCHIPS}" = "1" ]; then ext=0 else ext=${k} fi # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good" ] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_wcsext { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do ext=1 # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good" ] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_statminsection { local xmin local ymin local xmax local ymax if [ "${STATXMIN[${CHIP}]}"_A != "_A" ] && \ [ "${STATXMAX[${CHIP}]}"_A != "_A" ] && \ [ "${STATYMIN[${CHIP}]}"_A != "_A" ] && \ [ "${STATYMAX[${CHIP}]}"_A != "_A" ]; then xmin=${STATXMIN[${CHIP}]} xmax=${STATXMAX[${CHIP}]} ymin=${STATYMIN[${CHIP}]} ymax=${STATYMAX[${CHIP}]} else k=1 statminsection="" xmax=${SIZEX[1]} ymax=${SIZEY[1]} while [ "${k}" -le "${NCHIPS}" ] do if [ ${SIZEX[k]} -lt ${xmax} ]; then xmax=${SIZEX[k]} fi if [ ${SIZEY[k]} -lt ${ymax} ]; then ymax=${SIZEY[k]} fi k=$((${k}+1)) done xmin=2 ymin=2 xmax=$((${xmax}-1)) ymax=$((${ymax}-1)) fi statminsection="${xmin},${xmax},${ymin},${ymax}" } function p_ls { ls -l | grep ^- | ${P_GAWK} '{print $NF}' } # the following extracts keywords; # leading and trailing whitespace is removed (2nd call to gawk) # and blanks inside the keyword value are replaced with underscores (e.g. for OBJECT keyword) function get_key { ${P_DFITS} $1 | ${P_FITSORT} -d "$2" | ${P_GAWK} 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' | ${P_GAWK} '{$1=$1}1' | sed 's/ /_/g' } function get_keyext { ${P_DFITS} -x 1 $1 | ${P_FITSORT} -d "$2" | ${P_GAWK} '{if ($2!="KEY_N/A") print $2}' } function get_region { # merge global and local region files if they exist regionfile="" region="" test -f ${TEMPDIR}/regionfile_${CHIP}.reg && rm ${TEMPDIR}/regionfile_${CHIP}.reg if [ -f ${INSTDIR}/${INSTRUMENT}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}.reg >> \ ${TEMPDIR}/regionfile_${CHIP}.reg regionfile="${TEMPDIR}/regionfile_${CHIP}.reg" region="-d ${regionfile}" fi if [ -f ${INSTDIR}/${INSTRUMENT}_${CHIP}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}_${CHIP}.reg >> \ ${TEMPDIR}/regionfile_${CHIP}.reg regionfile="${TEMPDIR}/regionfile_${CHIP}.reg" region="-d ${regionfile}" fi } function get_region_chip { # merge global and local region files if they exist # for a different chip than the currently processed one # (constant sky subtraction) regionfile="" region="" test -f ${TEMPDIR}/regionfile_$1.reg_$$ && rm ${TEMPDIR}/regionfile_$1.reg_$$ if [ -f ${INSTDIR}/${INSTRUMENT}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}.reg >> \ ${TEMPDIR}/regionfile_$1.reg_$$ regionfile="${TEMPDIR}/regionfile_$1.reg_$$" region="-d ${regionfile}" fi if [ -f ${INSTDIR}/${INSTRUMENT}_$1.reg ]; then cat ${INSTDIR}/${INSTRUMENT}_$1.reg >> \ ${TEMPDIR}/regionfile_$1.reg_$$ regionfile="${TEMPDIR}/regionfile_$1.reg_$$" region="-d ${regionfile}" fi } function prepare_sex_for_region { backstring="" weightstring="" if [ -f ${regionfile} ] && [ "${regionfile}_A" != "_A" ]; then # estimate the background, ignoring blanked pixels background=`${P_FITSSTAT} -i ${FILE} -m median -q` backstring="-BACK_TYPE MANUAL -BACK_VALUE ${background}" weightstring="-WEIGHT_TYPE MAP_WEIGHT -WEIGHT_IMAGE /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits" # check if the dummy weight exists if [ ! -e /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits ]; then mkdir -p /$1/WEIGHTS/ ${P_IC} -p -32 -c ${SIZEX[${CHIP}]} ${SIZEY[${CHIP}]} '1' > ${TEMPDIR}/dummy_${CHIP}.fits ${P_FITSPOLYGON} -i ${TEMPDIR}/dummy_${CHIP}.fits -p ${regionfile} \ -o /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits rm ${TEMPDIR}/dummy_${CHIP}.fits fi fi } function prepare_imred_for_backgroundmodel { # The relevant settings depending on the apply mode of the background model flatstring="" fringestring="" #################################################### # subtract background model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 0 ]; then fringestring="-FLAT N -FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE ${FRINGEIMAGE}" fi #################################################### # divide background model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 1 ]; then flatstring="-FRINGE N -FLAT Y -FLAT_THRESHHOLD 0.1 -FLAT_SCALE ${V_BACK_ADJUSTGAINS}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FLATIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FLATIMAGE=/$1/$2/$2_${CHIP}.fits fi flatstring="${flatstring} -FLAT_IMAGE ${FLATIMAGE} -FLAT_SCALEIMAGE ${FLATIMAGE}" fi #################################################### # subtract fringing model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 2 ]; then fringestring="-FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_FRINGESMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe_smooth.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE /$1/$2/$2_${CHIP}_illum.fits" fi #################################################### # divide background model, subtract fringing model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 3 ]; then flatstring="-FLAT Y -FLAT_THRESHHOLD 0.1 -FLAT_SCALE ${V_BACK_ADJUSTGAINS}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FLATIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FLATIMAGE=/$1/$2/$2_${CHIP}.fits fi flatstring="${flatstring} -FLAT_IMAGE ${FLATIMAGE} -FLAT_SCALEIMAGE ${FLATIMAGE}" fringestring="-FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_FRINGESMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe_smooth.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE /$1/$2/$2_${CHIP}_illum.fits" fi } for CHIP in $3 do # check whether the master calibration file # exists already. if yes, then delete it. test -f $1/$2/$2_${CHIP}.fits && rm $1/$2/$2_${CHIP}.fits FILES=`ls $1/$2/*_${CHIP}.fits` image=`ls $1/$2/*_${CHIP}.fits | ${P_GAWK} '(NR==1)'` EXPTIME=`get_key ${image} EXPTIME` for FILE in ${FILES} do if [ -L ${FILE} ]; then LINK=`${P_READLINK} ${FILE}` BASE=`basename ${LINK} .fits` DIR=`dirname ${LINK}` ln -s ${DIR}/${BASE}OC.fits $1/$2/${BASE}OC.fits RESULTDIR[${CHIP}]=`dirname ${LINK}` fi done MAXX=$(( ${CUTX[${CHIP}]} + ${SIZEX[${CHIP}]} - 1 )) MAXY=$(( ${CUTY[${CHIP}]} + ${SIZEY[${CHIP}]} - 1 )) do_overscan=Y if [ "${OVSCANX1[${CHIP}]}" = "0" ] && \ [ "${OVSCANX2[${CHIP}]}" = "0" ]; then do_overscan=N fi # check the number of exposures numexp=`ls $1/$2/*_${CHIP}.fits | wc | ${P_GAWK} '{print $1}'` # if less than 3 exposures are available: if [ "${numexp}" -lt 3 ]; then echo PREPROCESS_FAILURE: Need at least three images of chip ${CHIP} to create master dark continue fi # if numexp <= NFRAMES combine the images right away, # otherwise write them down and then combine them. xmin=2 xmax=$(( ${SIZEX[${CHIP}]}-1 )) ymin=2 ymax=$(( ${SIZEY[${CHIP}]}-1 )) if [ "${numexp}" -gt "${NFRAMES}" ]; then # overscan correct and trim DARK frames ${P_IMRED_ECL:?} `ls $1/$2/*_${CHIP}.fits` \ -c ${CONF}/imred.conf \ -MAXIMAGES ${NFRAMES} \ -OVERSCAN ${do_overscan} \ -OVERSCAN_REGION ${OVSCANX1[${CHIP}]},${OVSCANX2[${CHIP}]} \ -OVERSCAN_REJECT ${V_CAL_OVSCANNLOW},${V_CAL_OVSCANNHIGH} \ -OUTPUT Y \ -OUTPUT_DIR ${RESULTDIR[${CHIP}]}/ \ -TRIM Y \ -OUTPUT_SUFFIX OC.fits \ -TRIM_REGION ${CUTX[${CHIP}]},${MAXX},${CUTY[${CHIP}]},${MAXY} \ -STATSSEC ${xmin},${xmax},${ymin},${ymax} # and combine them ls -1 $1/$2/*_${CHIP}OC.fits > ${TEMPDIR}/dark_images_$$ ${P_IMCOMBFLAT_IMCAT} -i ${TEMPDIR}/dark_images_$$\ -o ${RESULTDIR[${CHIP}]}/$2_${CHIP}.fits \ -e ${V_CAL_DARKNLOW} ${V_CAL_DARKNHIGH} fi if [ "${numexp}" -le "${NFRAMES}" ]; then # overscan correct and trim DARK frames ${P_IMRED_ECL:?} `ls $1/$2/*_${CHIP}.fits` \ -c ${CONF}/imred.conf \ -MAXIMAGES ${NFRAMES} \ -OVERSCAN ${do_overscan} \ -OVERSCAN_REGION ${OVSCANX1[${CHIP}]},${OVSCANX2[${CHIP}]} \ -OVERSCAN_REJECT ${V_CAL_OVSCANNLOW},${V_CAL_OVSCANNHIGH} \ -OUTPUT N \ -OUTPUT_DIR ${RESULTDIR[${CHIP}]}/ \ -TRIM Y \ -OUTPUT_SUFFIX OC.fits \ -TRIM_REGION ${CUTX[${CHIP}]},${MAXX},${CUTY[${CHIP}]},${MAXY} \ -STATSSEC ${xmin},${xmax},${ymin},${ymax} \ -COMBINE Y \ -COMBINE_MAXVAL 1e8 \ -COMBINE_RESCALE N \ -COMBINE_OUTPUT ${RESULTDIR[${CHIP}]}/$2_${CHIP}.fits \ -COMBINE_REJECT ${V_CAL_DARKNLOW},${V_CAL_DARKNHIGH} fi MD5KEY=`get_key ${RESULTDIR[${CHIP}]}/$2_${CHIP}.fits DATAMD5` if [ "${MD5KEY}"_A != "_A" ]; then ${P_REPLACEKEY} ${RESULTDIR[${CHIP}]}/$2_${CHIP}.fits "EXPTIME = ${EXPTIME}" "DATAMD5" else ${P_FITSADDKEY} -i ${RESULTDIR[${CHIP}]}/$2_${CHIP}.fits -k EXPTIME -v ${EXPTIME} fi if [ "${RESULTDIR[${CHIP}]}" != "$1/$2" ]; then ln -s ${RESULTDIR[${CHIP}]}/$2_${CHIP}.fits $1/$2/$2_${CHIP}.fits fi \rm -f /$1/$2/*_${CHIP}OC.fits if [ "${RESULTDIR[${CHIP}]}" != "$1/$2" ]; then \rm -f ${RESULTDIR[${CHIP}]}/*_${CHIP}OC.fits fi done + for CHIP in '$3' + test -f /home/albireo/14-17.02.2015/rgb//dark/dark_1.fits ls $1/$2/*_${CHIP}.fits ++ ls /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0001_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0002_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0003_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0004_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0005_1.fits + FILES='/home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0001_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0002_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0003_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0004_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0005_1.fits' ls $1/$2/*_${CHIP}.fits | ${P_GAWK} '(NR==1)' ++ /usr/bin/gawk '(NR==1)' ++ ls /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0001_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0002_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0003_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0004_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0005_1.fits + image=/home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits get_key ${image} EXPTIME ++ get_key /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits EXPTIME ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d EXPTIME ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + EXPTIME=600.00 + for FILE in '${FILES}' + '[' -L /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits ']' + for FILE in '${FILES}' + '[' -L /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0001_1.fits ']' + for FILE in '${FILES}' + '[' -L /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0002_1.fits ']' + for FILE in '${FILES}' + '[' -L /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0003_1.fits ']' + for FILE in '${FILES}' + '[' -L /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0004_1.fits ']' + for FILE in '${FILES}' + '[' -L /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0005_1.fits ']' + MAXX=3358 + MAXY=2536 + do_overscan=Y + '[' 0 = 0 ']' + '[' 0 = 0 ']' + do_overscan=N ls $1/$2/*_${CHIP}.fits | wc | ${P_GAWK} '{print $1}' ++ wc ++ /usr/bin/gawk '{print $1}' ++ ls /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0001_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0002_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0003_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0004_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0005_1.fits + numexp=6 + '[' 6 -lt 3 ']' + xmin=2 + xmax=3357 + ymin=2 + ymax=2535 + '[' 6 -gt 47 ']' + '[' 6 -le 47 ']' ls $1/$2/*_${CHIP}.fits ++ ls /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0001_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0002_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0003_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0004_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0005_1.fits + /home/albireo/THELI/theli//bin/Linux_64//preprocess /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0000_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0001_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0002_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0003_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0004_1.fits /home/albireo/14-17.02.2015/rgb//dark/dark600sek_0005_1.fits -c /home/albireo/THELI/theli//gui/reduction/imred.conf -MAXIMAGES 47 -OVERSCAN N -OVERSCAN_REGION 0,0 -OVERSCAN_REJECT 0,1 -OUTPUT N -OUTPUT_DIR /home/albireo/14-17.02.2015/rgb//dark/ -TRIM Y -OUTPUT_SUFFIX OC.fits -TRIM_REGION 1,3358,1,2536 -STATSSEC 2,3357,2,2535 -COMBINE Y -COMBINE_MAXVAL 1e8 -COMBINE_RESCALE N -COMBINE_OUTPUT /home/albireo/14-17.02.2015/rgb//dark/dark_1.fits -COMBINE_REJECT 0,1 -> Part 1: Preliminary work: -> Part 2: performing run processing -> 2.1: loading input images for run 1 loading framelist...: 1 out of 6 loading framelist...: 2 out of 6 loading framelist...: 3 out of 6 loading framelist...: 4 out of 6 loading framelist...: 5 out of 6 loading framelist...: 6 out of 6 -> 2.2: trimming images -> Part 3: stacking cube without rescaling get_key ${RESULTDIR[${CHIP}]}/$2_${CHIP}.fits DATAMD5 ++ get_key /home/albireo/14-17.02.2015/rgb//dark/dark_1.fits DATAMD5 ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//dark/dark_1.fits ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d DATAMD5 ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' + MD5KEY=6a8dfba3ee01a449348dc51dc1966c20 + '[' 6a8dfba3ee01a449348dc51dc1966c20_A '!=' _A ']' + /home/albireo/THELI/theli//bin/Linux_64//replacekey /home/albireo/14-17.02.2015/rgb//dark/dark_1.fits 'EXPTIME = 600.00' DATAMD5 + '[' /home/albireo/14-17.02.2015/rgb//dark '!=' /home/albireo/14-17.02.2015/rgb//dark ']' + rm -f '//home/albireo/14-17.02.2015/rgb//dark/*_1OC.fits' + '[' /home/albireo/14-17.02.2015/rgb//dark '!=' /home/albireo/14-17.02.2015/rgb//dark ']' #!/bin/bash -xv # the script processes a set of FLAT frames # the images are overscan corrected, debiased and stacked with # mode rescaling. #$1: main directory (filter) #$2: Bias directory #$3: Flat directory #$4: chips to process # preliminary work: INSTDIR=instruments_professional + INSTDIR=instruments_professional if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=instruments_commercial fi + '[' '!' -f instruments_professional/Moravian10z.ini ']' + INSTDIR=instruments_commercial if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=~/.theli/instruments_user/ fi + '[' '!' -f instruments_commercial/Moravian10z.ini ']' + INSTDIR=/home/albireo/.theli/instruments_user/ . ${INSTDIR}/${INSTRUMENT:?}.ini + . /home/albireo/.theli/instruments_user//Moravian10z.ini # ************************************************************** # Instrument config file created by THELI # ************************************************************** # Name of the instrument (must be unique!) INSTRUMENT=Moravian10z ++ INSTRUMENT=Moravian10z # Number of chips NCHIPS=1 ++ NCHIPS=1 # We need to export these export INSTRUMENT ++ export INSTRUMENT export NCHIPS ++ export NCHIPS # Load directories for programs . progs.ini ++ . progs.ini # If we include this script in bash verbose and trace mode we do not want # all the (verbose) output from this script: USE_X=0 +++ USE_X=0 if [[ ${-/x} != $- ]] ; then USE_X=1 set +xv fi +++ [[ hvB != hvxB ]] +++ USE_X=1 +++ set +xv # # this is needed for the GUI # P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh +++ P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh P_GETPMMOBJECTS=${BIN}/get_pmm_objects +++ P_GETPMMOBJECTS=/home/albireo/THELI/theli//bin/Linux_64//get_pmm_objects P_LDACDESC=${BIN}/ldacdesc +++ P_LDACDESC=/home/albireo/THELI/theli//bin/Linux_64//ldacdesc P_LDACDELTAB=${BIN}/ldacdeltab +++ P_LDACDELTAB=/home/albireo/THELI/theli//bin/Linux_64//ldacdeltab P_LDACDELKEY=${BIN}/ldacdelkey +++ P_LDACDELKEY=/home/albireo/THELI/theli//bin/Linux_64//ldacdelkey P_LDACPUTXY=${BIN}/ldacputxy +++ P_LDACPUTXY=/home/albireo/THELI/theli//bin/Linux_64//ldacputxy P_FITSFLIP=${BIN}/fitsflip +++ P_FITSFLIP=/home/albireo/THELI/theli//bin/Linux_64//fitsflip P_DEBLOOM=${BIN}/debloom +++ P_DEBLOOM=/home/albireo/THELI/theli//bin/Linux_64//debloom P_EQUALISEBAYERFLAT=${BIN}/equalisebayerflat +++ P_EQUALISEBAYERFLAT=/home/albireo/THELI/theli//bin/Linux_64//equalisebayerflat P_FITSCOLLAPSE=${BIN}/fitscollapse +++ P_FITSCOLLAPSE=/home/albireo/THELI/theli//bin/Linux_64//fitscollapse P_FITSWRAP=${BIN}/fitswrap +++ P_FITSWRAP=/home/albireo/THELI/theli//bin/Linux_64//fitswrap P_FITSCUBESLICE=${BIN}/fitscubeslice +++ P_FITSCUBESLICE=/home/albireo/THELI/theli//bin/Linux_64//fitscubeslice P_FITSADDKEY=${BIN}/fitsaddkey +++ P_FITSADDKEY=/home/albireo/THELI/theli//bin/Linux_64//fitsaddkey P_FITSCUT=${BIN}/fitscut +++ P_FITSCUT=/home/albireo/THELI/theli//bin/Linux_64//fitscut P_FITSBIN=${BIN}/fitsbin +++ P_FITSBIN=/home/albireo/THELI/theli//bin/Linux_64//fitsbin P_FITSBLOOMDETECT=${BIN}/fitsbloomdetect +++ P_FITSBLOOMDETECT=/home/albireo/THELI/theli//bin/Linux_64//fitsbloomdetect P_FITSFILL=${BIN}/fitsfill +++ P_FITSFILL=/home/albireo/THELI/theli//bin/Linux_64//fitsfill P_FITSGRAD=${BIN}/fitsgrad +++ P_FITSGRAD=/home/albireo/THELI/theli//bin/Linux_64//fitsgrad P_FITSDEMOSAICBAYER=${BIN}/fitsdemosaicbayer +++ P_FITSDEMOSAICBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitsdemosaicbayer P_FITS3MAX=${BIN}/fits3max +++ P_FITS3MAX=/home/albireo/THELI/theli//bin/Linux_64//fits3max P_FITSGAUSS=${BIN}/fitsgauss +++ P_FITSGAUSS=/home/albireo/THELI/theli//bin/Linux_64//fitsgauss P_FITS2TIFF=${BIN}/fits2tiff +++ P_FITS2TIFF=/home/albireo/THELI/theli//bin/Linux_64//fits2tiff P_TIFF2FITS=${BIN}/tiff2fits +++ P_TIFF2FITS=/home/albireo/THELI/theli//bin/Linux_64//tiff2fits P_CORRLIRIS=${BIN}/corrliris +++ P_CORRLIRIS=/home/albireo/THELI/theli//bin/Linux_64//corrliris P_CROSSTALK=${BIN}/crosstalk +++ P_CROSSTALK=/home/albireo/THELI/theli//bin/Linux_64//crosstalk P_FINDGSC=${BIN}/findgsc +++ P_FINDGSC=/home/albireo/THELI/theli//bin/Linux_64//findgsc P_FINDUSNOB1=${BIN}/findusnob1 +++ P_FINDUSNOB1=/home/albireo/THELI/theli//bin/Linux_64//findusnob1 P_FIND2MASS=${BIN}/find2mass +++ P_FIND2MASS=/home/albireo/THELI/theli//bin/Linux_64//find2mass P_FINDSDSS7=${BIN}/findsdss7 +++ P_FINDSDSS7=/home/albireo/THELI/theli//bin/Linux_64//findsdss7 P_FINDUCAC3=${BIN}/finducac3 +++ P_FINDUCAC3=/home/albireo/THELI/theli//bin/Linux_64//finducac3 P_FINDNOMAD1=${BIN}/findnomad1 +++ P_FINDNOMAD1=/home/albireo/THELI/theli//bin/Linux_64//findnomad1 P_VIZQUERY=${BIN}/vizquery +++ P_VIZQUERY=/home/albireo/THELI/theli//bin/Linux_64//vizquery P_SESAME=${BIN}/sesame +++ P_SESAME=/home/albireo/THELI/theli//bin/Linux_64//sesame P_ACLIENT=${BIN}/aclient +++ P_ACLIENT=/home/albireo/THELI/theli//bin/Linux_64//aclient P_SCAMP=${BIN}/scamp +++ P_SCAMP=/home/albireo/THELI/theli//bin/Linux_64//scamp P_RENAME=${BIN}/rename +++ P_RENAME=/home/albireo/THELI/theli//bin/Linux_64//rename P_MISSFITS=${BIN}/missfits +++ P_MISSFITS=/home/albireo/THELI/theli//bin/Linux_64//missfits P_SWARPFILTER=${BIN}/swarpfilter +++ P_SWARPFILTER=/home/albireo/THELI/theli//bin/Linux_64//swarpfilter P_FITSBLOCKEDIT=${BIN}/fitsblockedit +++ P_FITSBLOCKEDIT=/home/albireo/THELI/theli//bin/Linux_64//fitsblockedit P_FITSMASK=${BIN}/fitsmask +++ P_FITSMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsmask P_FITSMEDIAN=${BIN}/fitsmedian +++ P_FITSMEDIAN=/home/albireo/THELI/theli//bin/Linux_64//fitsmedian P_FITSNANMASK=${BIN}/fitsnanmask +++ P_FITSNANMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsnanmask P_FITSSTAT=${BIN}/fitsstat +++ P_FITSSTAT=/home/albireo/THELI/theli//bin/Linux_64//fitsstat P_FITSSPLITBAYER=${BIN}/fitssplitbayer +++ P_FITSSPLITBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitssplitbayer P_FITSGUESSGEOM=${BIN}/fitsguessgeom +++ P_FITSGUESSGEOM=/home/albireo/THELI/theli//bin/Linux_64//fitsguessgeom P_SPLITSUPRIMECAM=${BIN}/splitsuprimecam +++ P_SPLITSUPRIMECAM=/home/albireo/THELI/theli//bin/Linux_64//splitsuprimecam P_PASTEDUALCHANNEL=${BIN}/paste_dualchannel +++ P_PASTEDUALCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_dualchannel P_PASTEQUADCHANNEL=${BIN}/paste_quadchannel +++ P_PASTEQUADCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_quadchannel P_SPLITGMOSMULTIPORT=${BIN}/split_gmos_multiport +++ P_SPLITGMOSMULTIPORT=/home/albireo/THELI/theli//bin/Linux_64//split_gmos_multiport P_SPLITGROND=${BIN}/split_grond +++ P_SPLITGROND=/home/albireo/THELI/theli//bin/Linux_64//split_grond P_GETSHIFT=${BIN}/get_shift +++ P_GETSHIFT=/home/albireo/THELI/theli//bin/Linux_64//get_shift P_FITSNONLINEARITY=${BIN}/fitsnonlinearity +++ P_FITSNONLINEARITY=/home/albireo/THELI/theli//bin/Linux_64//fitsnonlinearity P_FITSPIC=${BIN}/fitspic +++ P_FITSPIC=/home/albireo/THELI/theli//bin/Linux_64//fitspic P_FITSPOLYGON=${BIN}/fitspolygon +++ P_FITSPOLYGON=/home/albireo/THELI/theli//bin/Linux_64//fitspolygon P_EXPANDSEXMASK=${BIN}/expand_sexmask +++ P_EXPANDSEXMASK=/home/albireo/THELI/theli//bin/Linux_64//expand_sexmask P_DCRAW=${BIN}/dcraw +++ P_DCRAW=/home/albireo/THELI/theli//bin/Linux_64//dcraw P_DECSEXCONV=${BIN}/decsexconv +++ P_DECSEXCONV=/home/albireo/THELI/theli//bin/Linux_64//decsexconv P_STIFF=${BIN}/stiff +++ P_STIFF=/home/albireo/THELI/theli//bin/Linux_64//stiff P_SUBSKY=${BIN}/subsky +++ P_SUBSKY=/home/albireo/THELI/theli//bin/Linux_64//subsky P_GIFSICLE=${BIN}/gifsicle +++ P_GIFSICLE=/home/albireo/THELI/theli//bin/Linux_64//gifsicle P_XY2SKY=${BIN}/xy2sky +++ P_XY2SKY=/home/albireo/THELI/theli//bin/Linux_64//xy2sky P_SKY2XY=${BIN}/sky2xy +++ P_SKY2XY=/home/albireo/THELI/theli//bin/Linux_64//sky2xy P_ISFITS=${BIN}/isfits +++ P_ISFITS=/home/albireo/THELI/theli//bin/Linux_64//isfits P_GETKEY=${BIN}/getkey +++ P_GETKEY=/home/albireo/THELI/theli//bin/Linux_64//getkey P_SKY2XY_HELPER=${BIN}/sky2xy_helper +++ P_SKY2XY_HELPER=/home/albireo/THELI/theli//bin/Linux_64//sky2xy_helper P_GET_ROTIMSIZE=${BIN}/get_rotimsize +++ P_GET_ROTIMSIZE=/home/albireo/THELI/theli//bin/Linux_64//get_rotimsize P_CLEANSTRING=${BIN}/cleanstring +++ P_CLEANSTRING=/home/albireo/THELI/theli//bin/Linux_64//cleanstring P_NUMCPU=${BIN}/numcpu +++ P_NUMCPU=/home/albireo/THELI/theli//bin/Linux_64//numcpu P_GETSTATS=${BIN}/getstats +++ P_GETSTATS=/home/albireo/THELI/theli//bin/Linux_64//getstats P_CDELT2CD=${BIN}/cdelt2cd +++ P_CDELT2CD=/home/albireo/THELI/theli//bin/Linux_64//cdelt2cd P_GETCDMATRIX=${BIN}/get_cdmatrix +++ P_GETCDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//get_cdmatrix P_GETPIXSCALE=${BIN}/get_pixscale +++ P_GETPIXSCALE=/home/albireo/THELI/theli//bin/Linux_64//get_pixscale P_GETPOSANGLE=${BIN}/get_posangle +++ P_GETPOSANGLE=/home/albireo/THELI/theli//bin/Linux_64//get_posangle P_GETREFCATRADIUS=${BIN}/get_refcat_radius +++ P_GETREFCATRADIUS=/home/albireo/THELI/theli//bin/Linux_64//get_refcat_radius P_MERGEDAT=${BIN}/mergedat +++ P_MERGEDAT=/home/albireo/THELI/theli//bin/Linux_64//mergedat P_ROTATECDMATRIX=${BIN}/rotate_cdmatrix +++ P_ROTATECDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//rotate_cdmatrix P_BACKGROUNDDYNAMIC=${BIN}/background_dynamic +++ P_BACKGROUNDDYNAMIC=/home/albireo/THELI/theli//bin/Linux_64//background_dynamic P_BACKGROUNDSTATIC=${BIN}/background_static +++ P_BACKGROUNDSTATIC=/home/albireo/THELI/theli//bin/Linux_64//background_static P_MJD=${BIN}/mjd +++ P_MJD=/home/albireo/THELI/theli//bin/Linux_64//mjd P_LST=${BIN}/lst +++ P_LST=/home/albireo/THELI/theli//bin/Linux_64//lst P_IMALYZER=${BIN}/imalyzer +++ P_IMALYZER=/home/albireo/THELI/theli//bin/Linux_64//imalyzer P_E2ELL=${BIN}/e2ell +++ P_E2ELL=/home/albireo/THELI/theli//bin/Linux_64//e2ell P_FITSSMOOTHEDGE=${BIN}/fitssmoothedge +++ P_FITSSMOOTHEDGE=/home/albireo/THELI/theli//bin/Linux_64//fitssmoothedge P_FITSEXTRACTEXT=${BIN}/fitsextractext +++ P_FITSEXTRACTEXT=/home/albireo/THELI/theli//bin/Linux_64//fitsextractext P_ERRTEST=${BIN}/errtest +++ P_ERRTEST=/home/albireo/THELI/theli//bin/Linux_64//errtest P_BUILDINDEX=${BIN_ANET}/build-astrometry-index +++ P_BUILDINDEX=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//build-astrometry-index P_SOLVEFIELD=${BIN_ANET}/solve-field +++ P_SOLVEFIELD=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//solve-field P_XCORR=${BIN}/xcorr +++ P_XCORR=/home/albireo/THELI/theli//bin/Linux_64//xcorr # # load the parameter files # . ~/.theli/param_set1.ini +++ . /home/albireo/.theli/param_set1.ini PROJECTNAME=Pferdekopf ++++ PROJECTNAME=Pferdekopf NPARA=1 ++++ NPARA=1 GUIVERSION=2.9.5 ++++ GUIVERSION=2.9.5 KERNEL=Linux ++++ KERNEL=Linux NFRAMES=47 ++++ NFRAMES=47 V_PRE_SPLITMIRCUBE=N ++++ V_PRE_SPLITMIRCUBE=N V_DO_BIAS=N ++++ V_DO_BIAS=N V_DO_FLAT=Y ++++ V_DO_FLAT=Y V_NONLINCORR=N ++++ V_NONLINCORR=N V_AP_SERVER=vizier.u-strasbg.fr ++++ V_AP_SERVER=vizier.u-strasbg.fr V_AP_MAGLIM=20 ++++ V_AP_MAGLIM=20 V_AP_RADIUS=30 ++++ V_AP_RADIUS=30 V_AP_REFRA=header ++++ V_AP_REFRA=header V_AP_REFDEC=header ++++ V_AP_REFDEC=header . ~/.theli/param_set2.ini +++ . /home/albireo/.theli/param_set2.ini V_PRE_RENAME_CHECKED=0 ++++ V_PRE_RENAME_CHECKED=0 V_RENAME_FITSKEY=ARCFILE ++++ V_RENAME_FITSKEY=ARCFILE V_SORT_FITSKEY=OBJECT ++++ V_SORT_FITSKEY=OBJECT V_SORT_BIAS= ++++ V_SORT_BIAS= V_SORT_DARK= ++++ V_SORT_DARK= V_SORT_DOMEFLAT= ++++ V_SORT_DOMEFLAT= V_SORT_SKYFLAT= ++++ V_SORT_SKYFLAT= V_SORT_STD= ++++ V_SORT_STD= V_PRE_XTALK_NOR_CHECKED=0 ++++ V_PRE_XTALK_NOR_CHECKED=0 V_PRE_XTALK_ROW_CHECKED=0 ++++ V_PRE_XTALK_ROW_CHECKED=0 V_PRE_XTALK_MULTI_CHECKED=0 ++++ V_PRE_XTALK_MULTI_CHECKED=0 V_PRE_XTALK_NOR_AMPLITUDE= ++++ V_PRE_XTALK_NOR_AMPLITUDE= V_PRE_XTALK_ROW_AMPLITUDE= ++++ V_PRE_XTALK_ROW_AMPLITUDE= V_PRE_XTALK_NOR_BUTTONID=-1 ++++ V_PRE_XTALK_NOR_BUTTONID=-1 V_PRE_XTALK_ROW_BUTTONID=-1 ++++ V_PRE_XTALK_ROW_BUTTONID=-1 V_PRE_XTALK_MULTI_BUTTONID=-1 ++++ V_PRE_XTALK_MULTI_BUTTONID=-1 V_PRE_XTALK_MULTI_NSECTION= ++++ V_PRE_XTALK_MULTI_NSECTION= V_CAL_OVSCANNLOW=0 ++++ V_CAL_OVSCANNLOW=0 V_CAL_OVSCANNHIGH=1 ++++ V_CAL_OVSCANNHIGH=1 V_CAL_BIASNLOW=0 ++++ V_CAL_BIASNLOW=0 V_CAL_BIASNHIGH=1 ++++ V_CAL_BIASNHIGH=1 V_CAL_DARKNLOW=0 ++++ V_CAL_DARKNLOW=0 V_CAL_DARKNHIGH=1 ++++ V_CAL_DARKNHIGH=1 V_CAL_FLATNLOW=0 ++++ V_CAL_FLATNLOW=0 V_CAL_FLATNHIGH=1 ++++ V_CAL_FLATNHIGH=1 V_CAL_USEDARK=FALSE ++++ V_CAL_USEDARK=FALSE V_BACK_NLOW=0 ++++ V_BACK_NLOW=0 V_BACK_NHIGH=1 ++++ V_BACK_NHIGH=1 V_BACK_DETECTTHRESH=1.3 ++++ V_BACK_DETECTTHRESH=1.3 V_BACK_DETECTMINAREA=5 ++++ V_BACK_DETECTMINAREA=5 V_BACK_MASKEXPAND= ++++ V_BACK_MASKEXPAND= V_BACK_ILLUMSMOOTH= ++++ V_BACK_ILLUMSMOOTH= V_BACK_FRINGESMOOTH= ++++ V_BACK_FRINGESMOOTH= V_BACK_WINDOWSIZE=0 ++++ V_BACK_WINDOWSIZE=0 V_BACK_GAPSIZE=1.0 ++++ V_BACK_GAPSIZE=1.0 V_BACK_APPLYMODE=0 ++++ V_BACK_APPLYMODE=0 V_BACK_COMBINEMETHOD=0 ++++ V_BACK_COMBINEMETHOD=0 V_BACK_SEXFILTER=Y ++++ V_BACK_SEXFILTER=Y V_BACK_ADJUSTGAINS=N ++++ V_BACK_ADJUSTGAINS=N V_BACK_FRINGESCALE=Y ++++ V_BACK_FRINGESCALE=Y V_BACK_SKY2PASS=Y ++++ V_BACK_SKY2PASS=Y V_COLLDETECTTHRESH= ++++ V_COLLDETECTTHRESH= V_COLLDETECTMINAREA= ++++ V_COLLDETECTMINAREA= V_COLLMASKEXPAND= ++++ V_COLLMASKEXPAND= V_COLLREJECTTHRESH=1.5 ++++ V_COLLREJECTTHRESH=1.5 V_COLLDIRECTION=x ++++ V_COLLDIRECTION=x V_COLLXMIN= ++++ V_COLLXMIN= V_COLLXMAX= ++++ V_COLLXMAX= V_COLLYMIN= ++++ V_COLLYMIN= V_COLLYMAX= ++++ V_COLLYMAX= V_COLLMASKACTION=1 ++++ V_COLLMASKACTION=1 V_COLLAUTOTHRESHOLD=0 ++++ V_COLLAUTOTHRESHOLD=0 V_WEIGHTBINMIN=-100 ++++ V_WEIGHTBINMIN=-100 V_WEIGHTBINMAX=500 ++++ V_WEIGHTBINMAX=500 V_WEIGHTBINSIZE=4 ++++ V_WEIGHTBINSIZE=4 V_WEIGHTLOWTHRESHOLD= ++++ V_WEIGHTLOWTHRESHOLD= V_WEIGHTHIGHTHRESHOLD= ++++ V_WEIGHTHIGHTHRESHOLD= V_GLOBWFLATMIN=0.5 ++++ V_GLOBWFLATMIN=0.5 V_GLOBWFLATMAX=1.6 ++++ V_GLOBWFLATMAX=1.6 V_GLOBWDARKMIN= ++++ V_GLOBWDARKMIN= V_GLOBWDARKMAX= ++++ V_GLOBWDARKMAX= V_GLOBWDARKDIR= ++++ V_GLOBWDARKDIR= V_COSMICSTHRESHOLD=0.1 ++++ V_COSMICSTHRESHOLD=0.1 V_COSMICDT=6 ++++ V_COSMICDT=6 V_COSMICDMIN=1 ++++ V_COSMICDMIN=1 V_DEFECT_KERNELSIZE= ++++ V_DEFECT_KERNELSIZE= V_DEFECT_ROWTOL= ++++ V_DEFECT_ROWTOL= V_DEFECT_COLTOL= ++++ V_DEFECT_COLTOL= V_DEFECT_CLUSTOL= ++++ V_DEFECT_CLUSTOL= V_DEFECT_KERNELSIZE_SF= ++++ V_DEFECT_KERNELSIZE_SF= V_DEFECT_ROWTOL_SF= ++++ V_DEFECT_ROWTOL_SF= V_DEFECT_COLTOL_SF= ++++ V_DEFECT_COLTOL_SF= V_DEFECT_CLUSTOL_SF= ++++ V_DEFECT_CLUSTOL_SF= V_BLOOMLOWLIMIT=20000 ++++ V_BLOOMLOWLIMIT=20000 V_BLOOMMINAREA=100 ++++ V_BLOOMMINAREA=100 V_BLOOMWIDTH=0 ++++ V_BLOOMWIDTH=0 V_WEIGHT_BINOUTLIER=FALSE ++++ V_WEIGHT_BINOUTLIER=FALSE V_GLOBW_UNIFORMWEIGHT=FALSE ++++ V_GLOBW_UNIFORMWEIGHT=FALSE V_MASKBLOOMSPIKE=0 ++++ V_MASKBLOOMSPIKE=0 V_AP_DETECTTHRESH=5.0 ++++ V_AP_DETECTTHRESH=5.0 V_AP_DETECTMINAREA=5 ++++ V_AP_DETECTMINAREA=5 V_DEBLENDMINCONT=0.0005 ++++ V_DEBLENDMINCONT=0.0005 V_AP_LOWNUM=1 ++++ V_AP_LOWNUM=1 V_SEXCATMINFWHM=1.5 ++++ V_SEXCATMINFWHM=1.5 V_SEXCATFLAG=0 ++++ V_SEXCATFLAG=0 V_SATURATION= ++++ V_SATURATION= V_SEXBACKLEVEL= ++++ V_SEXBACKLEVEL= V_AP_FILTER=N ++++ V_AP_FILTER=N V_SCAMP_MAXPOSANGLE=2.0 ++++ V_SCAMP_MAXPOSANGLE=2.0 V_SCAMP_MAXOFFSET=2.0 ++++ V_SCAMP_MAXOFFSET=2.0 V_SCAMP_MAXSCALE=1.05 ++++ V_SCAMP_MAXSCALE=1.05 V_SCAMP_SNLOW=5 ++++ V_SCAMP_SNLOW=5 V_SCAMP_SNHIGH=20 ++++ V_SCAMP_SNHIGH=20 V_SCAMP_FGROUPRADIUS=1.0 ++++ V_SCAMP_FGROUPRADIUS=1.0 V_SCAMP_CROSSIDRADIUS=0 ++++ V_SCAMP_CROSSIDRADIUS=0 V_SCAMP_ASTREFWEIGHT=1.0 ++++ V_SCAMP_ASTREFWEIGHT=1.0 V_SCAMP_ASTRINSTRUKEY=FILTER ++++ V_SCAMP_ASTRINSTRUKEY=FILTER V_SCAMP_PHOTINSTRUKEY=FILTER ++++ V_SCAMP_PHOTINSTRUKEY=FILTER V_SCAMP_POLYORDER=3 ++++ V_SCAMP_POLYORDER=3 V_SCAMP_POLYORDER2= ++++ V_SCAMP_POLYORDER2= V_SCAMP_DISTORTGROUPS= ++++ V_SCAMP_DISTORTGROUPS= V_SCAMP_DISTORTKEYS= ++++ V_SCAMP_DISTORTKEYS= V_SCAMP_STABILITY=INSTRUMENT ++++ V_SCAMP_STABILITY=INSTRUMENT V_SCAMP_MOSAICTYPE=UNCHANGED ++++ V_SCAMP_MOSAICTYPE=UNCHANGED V_SCAMP_FOCALPLANE=DEFAULT ++++ V_SCAMP_FOCALPLANE=DEFAULT V_SCAMP_MATCHFLIPPED=N ++++ V_SCAMP_MATCHFLIPPED=N V_SCAMP_MATCHWCS=Y ++++ V_SCAMP_MATCHWCS=Y V_ANET_MAXSCALE=1.05 ++++ V_ANET_MAXSCALE=1.05 V_ANET_POLYORDER=3 ++++ V_ANET_POLYORDER=3 V_ANET_SCAMPDISTORT=N ++++ V_ANET_SCAMPDISTORT=N V_ABSPHOT_FILTER=B ++++ V_ABSPHOT_FILTER=B V_ABSPHOT_COLORFIXED=0.1 ++++ V_ABSPHOT_COLORFIXED=0.1 V_ABSPHOT_EXTINCTION=0 ++++ V_ABSPHOT_EXTINCTION=0 V_ABSPHOT_ZPMIN=24 ++++ V_ABSPHOT_ZPMIN=24 V_ABSPHOT_ZPMAX=25 ++++ V_ABSPHOT_ZPMAX=25 V_ABSPHOT_STDMINMAG=0.0 ++++ V_ABSPHOT_STDMINMAG=0.0 V_ABSPHOT_MAXPHOTERR=0.05 ++++ V_ABSPHOT_MAXPHOTERR=0.05 V_ABSPHOT_WCSERR=10 ++++ V_ABSPHOT_WCSERR=10 V_ABSPHOT_WCSERRCHECKBOX=Y ++++ V_ABSPHOT_WCSERRCHECKBOX=Y V_ABSPHOT_STDCAT_DIRECT=SDSS ++++ V_ABSPHOT_STDCAT_DIRECT=SDSS V_ABSPHOT_STDFILTERM2=u ++++ V_ABSPHOT_STDFILTERM2=u V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI ++++ V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI V_ABSPHOT_STDFILTER=U ++++ V_ABSPHOT_STDFILTER=U V_ABSPHOT_STDCOLOR=UmB ++++ V_ABSPHOT_STDCOLOR=UmB V_ABSPHOT_CALMODE=RUNCALIB ++++ V_ABSPHOT_CALMODE=RUNCALIB V_ABSPHOT_FITTINGMETHODM2=CHIP ++++ V_ABSPHOT_FITTINGMETHODM2=CHIP V_COADD_REFRA= ++++ V_COADD_REFRA= V_COADD_REFDEC= ++++ V_COADD_REFDEC= V_COADD_IDENT=(null) ++++ V_COADD_IDENT=(null) V_COADD_SEEING= ++++ V_COADD_SEEING= V_COADD_RZP= ++++ V_COADD_RZP= V_COADD_PIXSCALE=0.89 ++++ V_COADD_PIXSCALE=0.89 V_COADD_SIZEX= ++++ V_COADD_SIZEX= V_COADD_SIZEY= ++++ V_COADD_SIZEY= V_COADD_SKYPOSANGLE= ++++ V_COADD_SKYPOSANGLE= V_COADD_CHIPS="" ++++ V_COADD_CHIPS= V_COADD_FILTERTHRESHOLD= ++++ V_COADD_FILTERTHRESHOLD= V_COADD_FILTERCLUSTERSIZE= ++++ V_COADD_FILTERCLUSTERSIZE= V_COADD_FILTERBORDERWIDTH= ++++ V_COADD_FILTERBORDERWIDTH= V_COADD_SMOOTHEDGE= ++++ V_COADD_SMOOTHEDGE= V_COADD_PROPMOTIONRA= ++++ V_COADD_PROPMOTIONRA= V_COADD_PROPMOTIONDEC= ++++ V_COADD_PROPMOTIONDEC= V_COADD_KERNEL=LANCZOS3 ++++ V_COADD_KERNEL=LANCZOS3 V_COADD_FILTER= ++++ V_COADD_FILTER= V_COADD_PROJECTION=TAN ++++ V_COADD_PROJECTION=TAN V_COADD_CELESTIALTYPE=EQUATORIAL ++++ V_COADD_CELESTIALTYPE=EQUATORIAL V_COADD_COMBINETYPE=WEIGHTED ++++ V_COADD_COMBINETYPE=WEIGHTED . ~/.theli/param_set3.ini +++ . /home/albireo/.theli/param_set3.ini # Geographical latitude of the observatory OBSLAT= ++ OBSLAT= # Pixel scale in x and y [degrees] (CDELT1/2 keywords) PIXSCX=-0 ++ PIXSCX=-0 PIXSCY=0 ++ PIXSCY=0 # Pixel scale in arcsec PIXSCALE= ++ PIXSCALE= # Detector gain GAIN=1.0 ++ GAIN=1.0 # Overscan regions. The values are the x-range. # The min is given by OVSCANX1, the max by OVSCANX2. # We need separate regions for each chip. OVSCANX1=([1]=0) ++ OVSCANX1=([1]=0) OVSCANX2=([1]=0) ++ OVSCANX2=([1]=0) # Sections for the cutting. Select those pixels that receive light. # Give min values for x and y, and the extent in x and y. CUTX=([1]=1) ++ CUTX=([1]=1) CUTY=([1]=1) ++ CUTY=([1]=1) SIZEX=([1]=3358) ++ SIZEX=([1]=3358) SIZEY=([1]=2536) ++ SIZEY=([1]=2536) # Reference pixels for each chip. REFPIXX=([1]=100) ++ REFPIXX=([1]=100) REFPIXY=([1]=100) ++ REFPIXY=([1]=100) # Camera type TYPE=OPT ++ TYPE=OPT . bash.include + . bash.include function get_wcs { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do if [ "${NCHIPS}" = "1" ]; then ext=0 else ext=${k} fi # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good" ] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_wcsext { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do ext=1 # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good" ] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_statminsection { local xmin local ymin local xmax local ymax if [ "${STATXMIN[${CHIP}]}"_A != "_A" ] && \ [ "${STATXMAX[${CHIP}]}"_A != "_A" ] && \ [ "${STATYMIN[${CHIP}]}"_A != "_A" ] && \ [ "${STATYMAX[${CHIP}]}"_A != "_A" ]; then xmin=${STATXMIN[${CHIP}]} xmax=${STATXMAX[${CHIP}]} ymin=${STATYMIN[${CHIP}]} ymax=${STATYMAX[${CHIP}]} else k=1 statminsection="" xmax=${SIZEX[1]} ymax=${SIZEY[1]} while [ "${k}" -le "${NCHIPS}" ] do if [ ${SIZEX[k]} -lt ${xmax} ]; then xmax=${SIZEX[k]} fi if [ ${SIZEY[k]} -lt ${ymax} ]; then ymax=${SIZEY[k]} fi k=$((${k}+1)) done xmin=2 ymin=2 xmax=$((${xmax}-1)) ymax=$((${ymax}-1)) fi statminsection="${xmin},${xmax},${ymin},${ymax}" } function p_ls { ls -l | grep ^- | ${P_GAWK} '{print $NF}' } # the following extracts keywords; # leading and trailing whitespace is removed (2nd call to gawk) # and blanks inside the keyword value are replaced with underscores (e.g. for OBJECT keyword) function get_key { ${P_DFITS} $1 | ${P_FITSORT} -d "$2" | ${P_GAWK} 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' | ${P_GAWK} '{$1=$1}1' | sed 's/ /_/g' } function get_keyext { ${P_DFITS} -x 1 $1 | ${P_FITSORT} -d "$2" | ${P_GAWK} '{if ($2!="KEY_N/A") print $2}' } function get_region { # merge global and local region files if they exist regionfile="" region="" test -f ${TEMPDIR}/regionfile_${CHIP}.reg && rm ${TEMPDIR}/regionfile_${CHIP}.reg if [ -f ${INSTDIR}/${INSTRUMENT}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}.reg >> \ ${TEMPDIR}/regionfile_${CHIP}.reg regionfile="${TEMPDIR}/regionfile_${CHIP}.reg" region="-d ${regionfile}" fi if [ -f ${INSTDIR}/${INSTRUMENT}_${CHIP}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}_${CHIP}.reg >> \ ${TEMPDIR}/regionfile_${CHIP}.reg regionfile="${TEMPDIR}/regionfile_${CHIP}.reg" region="-d ${regionfile}" fi } function get_region_chip { # merge global and local region files if they exist # for a different chip than the currently processed one # (constant sky subtraction) regionfile="" region="" test -f ${TEMPDIR}/regionfile_$1.reg_$$ && rm ${TEMPDIR}/regionfile_$1.reg_$$ if [ -f ${INSTDIR}/${INSTRUMENT}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}.reg >> \ ${TEMPDIR}/regionfile_$1.reg_$$ regionfile="${TEMPDIR}/regionfile_$1.reg_$$" region="-d ${regionfile}" fi if [ -f ${INSTDIR}/${INSTRUMENT}_$1.reg ]; then cat ${INSTDIR}/${INSTRUMENT}_$1.reg >> \ ${TEMPDIR}/regionfile_$1.reg_$$ regionfile="${TEMPDIR}/regionfile_$1.reg_$$" region="-d ${regionfile}" fi } function prepare_sex_for_region { backstring="" weightstring="" if [ -f ${regionfile} ] && [ "${regionfile}_A" != "_A" ]; then # estimate the background, ignoring blanked pixels background=`${P_FITSSTAT} -i ${FILE} -m median -q` backstring="-BACK_TYPE MANUAL -BACK_VALUE ${background}" weightstring="-WEIGHT_TYPE MAP_WEIGHT -WEIGHT_IMAGE /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits" # check if the dummy weight exists if [ ! -e /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits ]; then mkdir -p /$1/WEIGHTS/ ${P_IC} -p -32 -c ${SIZEX[${CHIP}]} ${SIZEY[${CHIP}]} '1' > ${TEMPDIR}/dummy_${CHIP}.fits ${P_FITSPOLYGON} -i ${TEMPDIR}/dummy_${CHIP}.fits -p ${regionfile} \ -o /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits rm ${TEMPDIR}/dummy_${CHIP}.fits fi fi } function prepare_imred_for_backgroundmodel { # The relevant settings depending on the apply mode of the background model flatstring="" fringestring="" #################################################### # subtract background model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 0 ]; then fringestring="-FLAT N -FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE ${FRINGEIMAGE}" fi #################################################### # divide background model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 1 ]; then flatstring="-FRINGE N -FLAT Y -FLAT_THRESHHOLD 0.1 -FLAT_SCALE ${V_BACK_ADJUSTGAINS}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FLATIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FLATIMAGE=/$1/$2/$2_${CHIP}.fits fi flatstring="${flatstring} -FLAT_IMAGE ${FLATIMAGE} -FLAT_SCALEIMAGE ${FLATIMAGE}" fi #################################################### # subtract fringing model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 2 ]; then fringestring="-FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_FRINGESMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe_smooth.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE /$1/$2/$2_${CHIP}_illum.fits" fi #################################################### # divide background model, subtract fringing model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 3 ]; then flatstring="-FLAT Y -FLAT_THRESHHOLD 0.1 -FLAT_SCALE ${V_BACK_ADJUSTGAINS}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FLATIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FLATIMAGE=/$1/$2/$2_${CHIP}.fits fi flatstring="${flatstring} -FLAT_IMAGE ${FLATIMAGE} -FLAT_SCALEIMAGE ${FLATIMAGE}" fringestring="-FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_FRINGESMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe_smooth.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE /$1/$2/$2_${CHIP}_illum.fits" fi } # # the resultdir is where the output coadded images # will go. If ONE image of the corresponding chip # is a link the image will go to THIS directory for CHIP in $4 do RESULTDIR[${CHIP}]="$1/$3" done + for CHIP in '$4' + RESULTDIR[${CHIP}]=/home/albireo/14-17.02.2015/rgb//flat # Correction overscan; create config file on the fly # we do not assume that the overscan is the same for all chips for CHIP in $4 do # check whether the master calibration file # exists already. if yes, then delete it. test -f $1/$3/$3_${CHIP}.fits && rm $1/$3/$3_${CHIP}.fits FILES=`ls $1/$3/*_${CHIP}.fits` for FILE in ${FILES} do BAYERPATTERN=`get_key ${FILE} BAYERPAT` if [ -L ${FILE} ]; then LINK=`${P_READLINK} ${FILE}` BASE=`basename ${LINK} .fits` DIR=`dirname ${LINK}` ln -s ${DIR}/${BASE}OC.fits $1/$2/${BASE}OC.fits RESULTDIR[${CHIP}]=`dirname ${LINK}` fi done MAXX=$(( ${CUTX[${CHIP}]} + ${SIZEX[${CHIP}]} - 1 )) MAXY=$(( ${CUTY[${CHIP}]} + ${SIZEY[${CHIP}]} - 1 )) do_overscan=Y if [ "${OVSCANX1[${CHIP}]}" = "0" ] && \ [ "${OVSCANX2[${CHIP}]}" = "0" ]; then do_overscan=N fi # check the number of exposures numexp=`ls $1/$3/*_${CHIP}.fits | wc | ${P_GAWK} '{print $1}'` # if less than 3 exposures are available: if [ "${numexp}" -lt 3 ]; then echo PREPOROCESS_FAILURE: Need at least three images of chip ${CHIP} to create master flat continue fi get_statminsection # if numexp <= NFRAMES combine the images right away, # otherwise write them down and then combine them. if [ "${numexp}" -gt "${NFRAMES}" ] || [ "${V_NONLINCORR}" = "Y" ]; then # overscan correct, trim and bias correct FLAT frames ${P_IMRED_ECL:?} `ls /$1/$3/*_${CHIP}.fits` \ -c ${CONF}/imred.conf \ -MAXIMAGES ${NFRAMES} \ -OVERSCAN ${do_overscan} \ -OVERSCAN_REGION ${OVSCANX1[${CHIP}]},${OVSCANX2[${CHIP}]} \ -OVERSCAN_REJECT ${V_CAL_OVSCANNLOW},${V_CAL_OVSCANNHIGH} \ -BIAS ${V_DO_BIAS} \ -BIAS_IMAGE /$1/$2/$2_${CHIP}.fits \ -OUTPUT Y \ -OUTPUT_DIR /$1/$3/ \ -TRIM Y \ -OUTPUT_SUFFIX OC.fits \ -TRIM_REGION ${CUTX[${CHIP}]},${MAXX},${CUTY[${CHIP}]},${MAXY} \ -STATSSEC ${statminsection} # correct for nonlinearity if [ "${V_NONLINCORR}" = "Y" ] && [ -f ${INSTDIR}/${INSTRUMENT}".nldat" ]; then ls /$1/$3/*_${CHIP}OC.fits > ${TEMPDIR}/nonlinlist_$$ cat ${TEMPDIR}/nonlinlist_$$ |\ { while read file do BASE=`basename ${file} _${CHIP}OC.fits` ${P_FITSNONLINEARITY} -i ${file} -o ${RESULTDIR[${CHIP}]}/${BASE}_${CHIP}OCnl.fits \ -c ${CHIP} -f ${INSTDIR}/${INSTRUMENT}".nldat" mv ${RESULTDIR[${CHIP}]}/${BASE}_${CHIP}OCnl.fits \ ${RESULTDIR[${CHIP}]}/${BASE}_${CHIP}OC.fits done } fi # and combine them ${P_IMSTATS} -d 6 `ls /$1/$3/*_${CHIP}OC.fits` \ -o ${TEMPDIR}/flat_images_$$ -t -1e8 1e8 ${P_IMCOMBFLAT_IMCAT} -i ${TEMPDIR}/flat_images_$$\ -o ${RESULTDIR[${CHIP}]}/$3_${CHIP}.fits \ -s 1 -e ${V_CAL_FLATNLOW} ${V_CAL_FLATNHIGH} fi if [ "${numexp}" -le "${NFRAMES}" ] && [ "${V_NONLINCORR}" = "N" ]; then # overscan correct, trim and bias correct FLAT frames ${P_IMRED_ECL:?} `ls /$1/$3/*_${CHIP}.fits` \ -c ${CONF}/imred.conf \ -MAXIMAGES ${NFRAMES} \ -OVERSCAN ${do_overscan} \ -OVERSCAN_REGION ${OVSCANX1[${CHIP}]},${OVSCANX2[${CHIP}]} \ -OVERSCAN_REJECT ${V_CAL_OVSCANNLOW},${V_CAL_OVSCANNHIGH} \ -BIAS ${V_DO_BIAS} \ -BIAS_IMAGE /$1/$2/$2_${CHIP}.fits \ -OUTPUT N \ -OUTPUT_DIR ${RESULTDIR[${CHIP}]}/ \ -TRIM Y \ -OUTPUT_SUFFIX OC.fits \ -TRIM_REGION ${CUTX[${CHIP}]},${MAXX},${CUTY[${CHIP}]},${MAXY} \ -COMBINE Y \ -COMBINE_RESCALE Y \ -COMBINE_OUTPUT ${RESULTDIR[${CHIP}]}/$3_${CHIP}.fits \ -COMBINE_REJECT ${V_CAL_FLATNLOW},${V_CAL_FLATNHIGH} \ -STATSSEC ${statminsection} fi # merge global and individual region files, if they exist get_region # mask the flat if [ "${region}"_A != "_A" ]; then # we mask it with a very small negative value, so that the # corresponding pixels in the flat-fielded data will be VERY # negative, and thus outside the range where statistics is probed ${P_FITSPOLYGON} \ -i ${RESULTDIR[${CHIP}]}/$3_${CHIP}.fits \ -o ${RESULTDIR[${CHIP}]}/$3_${CHIP}.masked.fits \ -p ${regionfile} \ -v -1e-9 mv ${RESULTDIR[${CHIP}]}/$3_${CHIP}.masked.fits ${RESULTDIR[${CHIP}]}/$3_${CHIP}.fits fi if [ ! "${BAYERPATTERN}_A" = "_A" ]; then FILE=$1/$3/$3_${CHIP}.fits BASE=`basename ${FILE} _${CHIP}.fits` ${P_EQUALISEBAYERFLAT} -i ${FILE} -o ${BASE}_tmp.fits mv ${BASE}_tmp.fits ${FILE} fi if [ "${RESULTDIR[${CHIP}]}" != "$1/$3" ]; then ln -s ${RESULTDIR[${CHIP}]}/$3_${CHIP}.fits $1/$3/$3_${CHIP}.fits fi \rm -f /$1/$3/*_${CHIP}OC.fits if [ "${RESULTDIR[${CHIP}]}" != "$1/$3" ]; then \rm -f ${RESULTDIR[${CHIP}]}/*_${CHIP}OC.fits fi done + for CHIP in '$4' + test -f /home/albireo/14-17.02.2015/rgb//flat/flat_1.fits ls $1/$3/*_${CHIP}.fits ++ ls /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits + FILES='/home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits ']' + MAXX=3358 + MAXY=2536 + do_overscan=Y + '[' 0 = 0 ']' + '[' 0 = 0 ']' + do_overscan=N ls $1/$3/*_${CHIP}.fits | wc | ${P_GAWK} '{print $1}' ++ wc ++ /usr/bin/gawk '{print $1}' ++ ls /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits /home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits + numexp=24 + '[' 24 -lt 3 ']' + get_statminsection + local xmin + local ymin + local xmax + local ymax + '[' _A '!=' _A ']' + k=1 + statminsection= + xmax=3358 + ymax=2536 + '[' 1 -le 1 ']' + '[' 3358 -lt 3358 ']' + '[' 2536 -lt 2536 ']' + k=2 + '[' 2 -le 1 ']' + xmin=2 + ymin=2 + xmax=3357 + ymax=2535 + statminsection=2,3357,2,2535 + '[' 24 -gt 47 ']' + '[' N = Y ']' + '[' 24 -le 47 ']' + '[' N = N ']' ls /$1/$3/*_${CHIP}.fits ++ ls //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits + /home/albireo/THELI/theli//bin/Linux_64//preprocess //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits -c /home/albireo/THELI/theli//gui/reduction/imred.conf -MAXIMAGES 47 -OVERSCAN N -OVERSCAN_REGION 0,0 -OVERSCAN_REJECT 0,1 -BIAS N -BIAS_IMAGE //home/albireo/14-17.02.2015/rgb//BIAS/BIAS_1.fits -OUTPUT N -OUTPUT_DIR /home/albireo/14-17.02.2015/rgb//flat/ -TRIM Y -OUTPUT_SUFFIX OC.fits -TRIM_REGION 1,3358,1,2536 -COMBINE Y -COMBINE_RESCALE Y -COMBINE_OUTPUT /home/albireo/14-17.02.2015/rgb//flat/flat_1.fits -COMBINE_REJECT 0,1 -STATSSEC 2,3357,2,2535 -> Part 1: Preliminary work: -> Part 2: performing run processing -> 2.1: loading input images for run 1 loading framelist...: 1 out of 24 loading framelist...: 2 out of 24 loading framelist...: 3 out of 24 loading framelist...: 4 out of 24 loading framelist...: 5 out of 24 loading framelist...: 6 out of 24 loading framelist...: 7 out of 24 loading framelist...: 8 out of 24 loading framelist...: 9 out of 24 loading framelist...: 10 out of 24 loading framelist...: 11 out of 24 loading framelist...: 12 out of 24 loading framelist...: 13 out of 24 loading framelist...: 14 out of 24 loading framelist...: 15 out of 24 loading framelist...: 16 out of 24 loading framelist...: 17 out of 24 loading framelist...: 18 out of 24 loading framelist...: 19 out of 24 loading framelist...: 20 out of 24 loading framelist...: 21 out of 24 loading framelist...: 22 out of 24 loading framelist...: 23 out of 24 loading framelist...: 24 out of 24 -> 2.2: trimming images -> Part 3: stacking cube with rescaling determine median of input frames: Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-23_Red_0000_1.fits: 24411.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-28-45_Red_0001_1.fits: 24161.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-06_Red_0002_1.fits: 24070.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-28_Red_0003_1.fits: 23589.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-29-49_Red_0004_1.fits: 23506.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-10_Red_0005_1.fits: 23056.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-32_Red_0006_1.fits: 22713.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-30-53_Red_0007_1.fits: 22518.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-09_Green_0008_1.fits: 27264.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-29_Green_0009_1.fits: 26883.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-32-48_Green_0010_1.fits: 26450.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-07_Green_0011_1.fits: 26093.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-27_Green_0012_1.fits: 25810.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-33-46_Green_0013_1.fits: 25266.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-06_Green_0014_1.fits: 24857.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-26_Green_0015_1.fits: 24440.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-34-55_Blue_0016_1.fits: 19910.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-35-20_Blue_0017_1.fits: 19572.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-34_Blue_0018_1.fits: 18303.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-36-55_Blue_0019_1.fits: 17993.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-15_Blue_0020_1.fits: 17759.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-35_Blue_0021_1.fits: 17533.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-37-55_Blue_0022_1.fits: 17239.000000 Image //home/albireo/14-17.02.2015/rgb//flat/flat16.02.15_2015-02-16_21-38-15_Blue_0023_1.fits: 16943.000000 mean of medians; 22514.125000 + get_region + regionfile= + region= + test -f /home/albireo/.theli/tmp//regionfile_1.reg + '[' -f /home/albireo/.theli/instruments_user//Moravian10z.reg ']' + '[' -f /home/albireo/.theli/instruments_user//Moravian10z_1.reg ']' + '[' _A '!=' _A ']' + '[' '!' _A = _A ']' + '[' /home/albireo/14-17.02.2015/rgb//flat '!=' /home/albireo/14-17.02.2015/rgb//flat ']' + rm -f '//home/albireo/14-17.02.2015/rgb//flat/*_1OC.fits' + '[' /home/albireo/14-17.02.2015/rgb//flat '!=' /home/albireo/14-17.02.2015/rgb//flat ']' #!/bin/bash -xv # this script determines the flux ratios between chips in a # flat-field image, to be fed into preprocess. # $1: main directory # $2: Flat directory # preliminary work: INSTDIR=instruments_professional + INSTDIR=instruments_professional if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=instruments_commercial fi + '[' '!' -f instruments_professional/Moravian10z.ini ']' + INSTDIR=instruments_commercial if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=~/.theli/instruments_user/ fi + '[' '!' -f instruments_commercial/Moravian10z.ini ']' + INSTDIR=/home/albireo/.theli/instruments_user/ . ${INSTDIR}/${INSTRUMENT:?}.ini + . /home/albireo/.theli/instruments_user//Moravian10z.ini # ************************************************************** # Instrument config file created by THELI # ************************************************************** # Name of the instrument (must be unique!) INSTRUMENT=Moravian10z ++ INSTRUMENT=Moravian10z # Number of chips NCHIPS=1 ++ NCHIPS=1 # We need to export these export INSTRUMENT ++ export INSTRUMENT export NCHIPS ++ export NCHIPS # Load directories for programs . progs.ini ++ . progs.ini # If we include this script in bash verbose and trace mode we do not want # all the (verbose) output from this script: USE_X=0 +++ USE_X=0 if [[ ${-/x} != $- ]] ; then USE_X=1 set +xv fi +++ [[ hvB != hvxB ]] +++ USE_X=1 +++ set +xv # # this is needed for the GUI # P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh +++ P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh P_GETPMMOBJECTS=${BIN}/get_pmm_objects +++ P_GETPMMOBJECTS=/home/albireo/THELI/theli//bin/Linux_64//get_pmm_objects P_LDACDESC=${BIN}/ldacdesc +++ P_LDACDESC=/home/albireo/THELI/theli//bin/Linux_64//ldacdesc P_LDACDELTAB=${BIN}/ldacdeltab +++ P_LDACDELTAB=/home/albireo/THELI/theli//bin/Linux_64//ldacdeltab P_LDACDELKEY=${BIN}/ldacdelkey +++ P_LDACDELKEY=/home/albireo/THELI/theli//bin/Linux_64//ldacdelkey P_LDACPUTXY=${BIN}/ldacputxy +++ P_LDACPUTXY=/home/albireo/THELI/theli//bin/Linux_64//ldacputxy P_FITSFLIP=${BIN}/fitsflip +++ P_FITSFLIP=/home/albireo/THELI/theli//bin/Linux_64//fitsflip P_DEBLOOM=${BIN}/debloom +++ P_DEBLOOM=/home/albireo/THELI/theli//bin/Linux_64//debloom P_EQUALISEBAYERFLAT=${BIN}/equalisebayerflat +++ P_EQUALISEBAYERFLAT=/home/albireo/THELI/theli//bin/Linux_64//equalisebayerflat P_FITSCOLLAPSE=${BIN}/fitscollapse +++ P_FITSCOLLAPSE=/home/albireo/THELI/theli//bin/Linux_64//fitscollapse P_FITSWRAP=${BIN}/fitswrap +++ P_FITSWRAP=/home/albireo/THELI/theli//bin/Linux_64//fitswrap P_FITSCUBESLICE=${BIN}/fitscubeslice +++ P_FITSCUBESLICE=/home/albireo/THELI/theli//bin/Linux_64//fitscubeslice P_FITSADDKEY=${BIN}/fitsaddkey +++ P_FITSADDKEY=/home/albireo/THELI/theli//bin/Linux_64//fitsaddkey P_FITSCUT=${BIN}/fitscut +++ P_FITSCUT=/home/albireo/THELI/theli//bin/Linux_64//fitscut P_FITSBIN=${BIN}/fitsbin +++ P_FITSBIN=/home/albireo/THELI/theli//bin/Linux_64//fitsbin P_FITSBLOOMDETECT=${BIN}/fitsbloomdetect +++ P_FITSBLOOMDETECT=/home/albireo/THELI/theli//bin/Linux_64//fitsbloomdetect P_FITSFILL=${BIN}/fitsfill +++ P_FITSFILL=/home/albireo/THELI/theli//bin/Linux_64//fitsfill P_FITSGRAD=${BIN}/fitsgrad +++ P_FITSGRAD=/home/albireo/THELI/theli//bin/Linux_64//fitsgrad P_FITSDEMOSAICBAYER=${BIN}/fitsdemosaicbayer +++ P_FITSDEMOSAICBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitsdemosaicbayer P_FITS3MAX=${BIN}/fits3max +++ P_FITS3MAX=/home/albireo/THELI/theli//bin/Linux_64//fits3max P_FITSGAUSS=${BIN}/fitsgauss +++ P_FITSGAUSS=/home/albireo/THELI/theli//bin/Linux_64//fitsgauss P_FITS2TIFF=${BIN}/fits2tiff +++ P_FITS2TIFF=/home/albireo/THELI/theli//bin/Linux_64//fits2tiff P_TIFF2FITS=${BIN}/tiff2fits +++ P_TIFF2FITS=/home/albireo/THELI/theli//bin/Linux_64//tiff2fits P_CORRLIRIS=${BIN}/corrliris +++ P_CORRLIRIS=/home/albireo/THELI/theli//bin/Linux_64//corrliris P_CROSSTALK=${BIN}/crosstalk +++ P_CROSSTALK=/home/albireo/THELI/theli//bin/Linux_64//crosstalk P_FINDGSC=${BIN}/findgsc +++ P_FINDGSC=/home/albireo/THELI/theli//bin/Linux_64//findgsc P_FINDUSNOB1=${BIN}/findusnob1 +++ P_FINDUSNOB1=/home/albireo/THELI/theli//bin/Linux_64//findusnob1 P_FIND2MASS=${BIN}/find2mass +++ P_FIND2MASS=/home/albireo/THELI/theli//bin/Linux_64//find2mass P_FINDSDSS7=${BIN}/findsdss7 +++ P_FINDSDSS7=/home/albireo/THELI/theli//bin/Linux_64//findsdss7 P_FINDUCAC3=${BIN}/finducac3 +++ P_FINDUCAC3=/home/albireo/THELI/theli//bin/Linux_64//finducac3 P_FINDNOMAD1=${BIN}/findnomad1 +++ P_FINDNOMAD1=/home/albireo/THELI/theli//bin/Linux_64//findnomad1 P_VIZQUERY=${BIN}/vizquery +++ P_VIZQUERY=/home/albireo/THELI/theli//bin/Linux_64//vizquery P_SESAME=${BIN}/sesame +++ P_SESAME=/home/albireo/THELI/theli//bin/Linux_64//sesame P_ACLIENT=${BIN}/aclient +++ P_ACLIENT=/home/albireo/THELI/theli//bin/Linux_64//aclient P_SCAMP=${BIN}/scamp +++ P_SCAMP=/home/albireo/THELI/theli//bin/Linux_64//scamp P_RENAME=${BIN}/rename +++ P_RENAME=/home/albireo/THELI/theli//bin/Linux_64//rename P_MISSFITS=${BIN}/missfits +++ P_MISSFITS=/home/albireo/THELI/theli//bin/Linux_64//missfits P_SWARPFILTER=${BIN}/swarpfilter +++ P_SWARPFILTER=/home/albireo/THELI/theli//bin/Linux_64//swarpfilter P_FITSBLOCKEDIT=${BIN}/fitsblockedit +++ P_FITSBLOCKEDIT=/home/albireo/THELI/theli//bin/Linux_64//fitsblockedit P_FITSMASK=${BIN}/fitsmask +++ P_FITSMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsmask P_FITSMEDIAN=${BIN}/fitsmedian +++ P_FITSMEDIAN=/home/albireo/THELI/theli//bin/Linux_64//fitsmedian P_FITSNANMASK=${BIN}/fitsnanmask +++ P_FITSNANMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsnanmask P_FITSSTAT=${BIN}/fitsstat +++ P_FITSSTAT=/home/albireo/THELI/theli//bin/Linux_64//fitsstat P_FITSSPLITBAYER=${BIN}/fitssplitbayer +++ P_FITSSPLITBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitssplitbayer P_FITSGUESSGEOM=${BIN}/fitsguessgeom +++ P_FITSGUESSGEOM=/home/albireo/THELI/theli//bin/Linux_64//fitsguessgeom P_SPLITSUPRIMECAM=${BIN}/splitsuprimecam +++ P_SPLITSUPRIMECAM=/home/albireo/THELI/theli//bin/Linux_64//splitsuprimecam P_PASTEDUALCHANNEL=${BIN}/paste_dualchannel +++ P_PASTEDUALCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_dualchannel P_PASTEQUADCHANNEL=${BIN}/paste_quadchannel +++ P_PASTEQUADCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_quadchannel P_SPLITGMOSMULTIPORT=${BIN}/split_gmos_multiport +++ P_SPLITGMOSMULTIPORT=/home/albireo/THELI/theli//bin/Linux_64//split_gmos_multiport P_SPLITGROND=${BIN}/split_grond +++ P_SPLITGROND=/home/albireo/THELI/theli//bin/Linux_64//split_grond P_GETSHIFT=${BIN}/get_shift +++ P_GETSHIFT=/home/albireo/THELI/theli//bin/Linux_64//get_shift P_FITSNONLINEARITY=${BIN}/fitsnonlinearity +++ P_FITSNONLINEARITY=/home/albireo/THELI/theli//bin/Linux_64//fitsnonlinearity P_FITSPIC=${BIN}/fitspic +++ P_FITSPIC=/home/albireo/THELI/theli//bin/Linux_64//fitspic P_FITSPOLYGON=${BIN}/fitspolygon +++ P_FITSPOLYGON=/home/albireo/THELI/theli//bin/Linux_64//fitspolygon P_EXPANDSEXMASK=${BIN}/expand_sexmask +++ P_EXPANDSEXMASK=/home/albireo/THELI/theli//bin/Linux_64//expand_sexmask P_DCRAW=${BIN}/dcraw +++ P_DCRAW=/home/albireo/THELI/theli//bin/Linux_64//dcraw P_DECSEXCONV=${BIN}/decsexconv +++ P_DECSEXCONV=/home/albireo/THELI/theli//bin/Linux_64//decsexconv P_STIFF=${BIN}/stiff +++ P_STIFF=/home/albireo/THELI/theli//bin/Linux_64//stiff P_SUBSKY=${BIN}/subsky +++ P_SUBSKY=/home/albireo/THELI/theli//bin/Linux_64//subsky P_GIFSICLE=${BIN}/gifsicle +++ P_GIFSICLE=/home/albireo/THELI/theli//bin/Linux_64//gifsicle P_XY2SKY=${BIN}/xy2sky +++ P_XY2SKY=/home/albireo/THELI/theli//bin/Linux_64//xy2sky P_SKY2XY=${BIN}/sky2xy +++ P_SKY2XY=/home/albireo/THELI/theli//bin/Linux_64//sky2xy P_ISFITS=${BIN}/isfits +++ P_ISFITS=/home/albireo/THELI/theli//bin/Linux_64//isfits P_GETKEY=${BIN}/getkey +++ P_GETKEY=/home/albireo/THELI/theli//bin/Linux_64//getkey P_SKY2XY_HELPER=${BIN}/sky2xy_helper +++ P_SKY2XY_HELPER=/home/albireo/THELI/theli//bin/Linux_64//sky2xy_helper P_GET_ROTIMSIZE=${BIN}/get_rotimsize +++ P_GET_ROTIMSIZE=/home/albireo/THELI/theli//bin/Linux_64//get_rotimsize P_CLEANSTRING=${BIN}/cleanstring +++ P_CLEANSTRING=/home/albireo/THELI/theli//bin/Linux_64//cleanstring P_NUMCPU=${BIN}/numcpu +++ P_NUMCPU=/home/albireo/THELI/theli//bin/Linux_64//numcpu P_GETSTATS=${BIN}/getstats +++ P_GETSTATS=/home/albireo/THELI/theli//bin/Linux_64//getstats P_CDELT2CD=${BIN}/cdelt2cd +++ P_CDELT2CD=/home/albireo/THELI/theli//bin/Linux_64//cdelt2cd P_GETCDMATRIX=${BIN}/get_cdmatrix +++ P_GETCDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//get_cdmatrix P_GETPIXSCALE=${BIN}/get_pixscale +++ P_GETPIXSCALE=/home/albireo/THELI/theli//bin/Linux_64//get_pixscale P_GETPOSANGLE=${BIN}/get_posangle +++ P_GETPOSANGLE=/home/albireo/THELI/theli//bin/Linux_64//get_posangle P_GETREFCATRADIUS=${BIN}/get_refcat_radius +++ P_GETREFCATRADIUS=/home/albireo/THELI/theli//bin/Linux_64//get_refcat_radius P_MERGEDAT=${BIN}/mergedat +++ P_MERGEDAT=/home/albireo/THELI/theli//bin/Linux_64//mergedat P_ROTATECDMATRIX=${BIN}/rotate_cdmatrix +++ P_ROTATECDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//rotate_cdmatrix P_BACKGROUNDDYNAMIC=${BIN}/background_dynamic +++ P_BACKGROUNDDYNAMIC=/home/albireo/THELI/theli//bin/Linux_64//background_dynamic P_BACKGROUNDSTATIC=${BIN}/background_static +++ P_BACKGROUNDSTATIC=/home/albireo/THELI/theli//bin/Linux_64//background_static P_MJD=${BIN}/mjd +++ P_MJD=/home/albireo/THELI/theli//bin/Linux_64//mjd P_LST=${BIN}/lst +++ P_LST=/home/albireo/THELI/theli//bin/Linux_64//lst P_IMALYZER=${BIN}/imalyzer +++ P_IMALYZER=/home/albireo/THELI/theli//bin/Linux_64//imalyzer P_E2ELL=${BIN}/e2ell +++ P_E2ELL=/home/albireo/THELI/theli//bin/Linux_64//e2ell P_FITSSMOOTHEDGE=${BIN}/fitssmoothedge +++ P_FITSSMOOTHEDGE=/home/albireo/THELI/theli//bin/Linux_64//fitssmoothedge P_FITSEXTRACTEXT=${BIN}/fitsextractext +++ P_FITSEXTRACTEXT=/home/albireo/THELI/theli//bin/Linux_64//fitsextractext P_ERRTEST=${BIN}/errtest +++ P_ERRTEST=/home/albireo/THELI/theli//bin/Linux_64//errtest P_BUILDINDEX=${BIN_ANET}/build-astrometry-index +++ P_BUILDINDEX=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//build-astrometry-index P_SOLVEFIELD=${BIN_ANET}/solve-field +++ P_SOLVEFIELD=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//solve-field P_XCORR=${BIN}/xcorr +++ P_XCORR=/home/albireo/THELI/theli//bin/Linux_64//xcorr # # load the parameter files # . ~/.theli/param_set1.ini +++ . /home/albireo/.theli/param_set1.ini PROJECTNAME=Pferdekopf ++++ PROJECTNAME=Pferdekopf NPARA=1 ++++ NPARA=1 GUIVERSION=2.9.5 ++++ GUIVERSION=2.9.5 KERNEL=Linux ++++ KERNEL=Linux NFRAMES=47 ++++ NFRAMES=47 V_PRE_SPLITMIRCUBE=N ++++ V_PRE_SPLITMIRCUBE=N V_DO_BIAS=N ++++ V_DO_BIAS=N V_DO_FLAT=Y ++++ V_DO_FLAT=Y V_NONLINCORR=N ++++ V_NONLINCORR=N V_AP_SERVER=vizier.u-strasbg.fr ++++ V_AP_SERVER=vizier.u-strasbg.fr V_AP_MAGLIM=20 ++++ V_AP_MAGLIM=20 V_AP_RADIUS=30 ++++ V_AP_RADIUS=30 V_AP_REFRA=header ++++ V_AP_REFRA=header V_AP_REFDEC=header ++++ V_AP_REFDEC=header . ~/.theli/param_set2.ini +++ . /home/albireo/.theli/param_set2.ini V_PRE_RENAME_CHECKED=0 ++++ V_PRE_RENAME_CHECKED=0 V_RENAME_FITSKEY=ARCFILE ++++ V_RENAME_FITSKEY=ARCFILE V_SORT_FITSKEY=OBJECT ++++ V_SORT_FITSKEY=OBJECT V_SORT_BIAS= ++++ V_SORT_BIAS= V_SORT_DARK= ++++ V_SORT_DARK= V_SORT_DOMEFLAT= ++++ V_SORT_DOMEFLAT= V_SORT_SKYFLAT= ++++ V_SORT_SKYFLAT= V_SORT_STD= ++++ V_SORT_STD= V_PRE_XTALK_NOR_CHECKED=0 ++++ V_PRE_XTALK_NOR_CHECKED=0 V_PRE_XTALK_ROW_CHECKED=0 ++++ V_PRE_XTALK_ROW_CHECKED=0 V_PRE_XTALK_MULTI_CHECKED=0 ++++ V_PRE_XTALK_MULTI_CHECKED=0 V_PRE_XTALK_NOR_AMPLITUDE= ++++ V_PRE_XTALK_NOR_AMPLITUDE= V_PRE_XTALK_ROW_AMPLITUDE= ++++ V_PRE_XTALK_ROW_AMPLITUDE= V_PRE_XTALK_NOR_BUTTONID=-1 ++++ V_PRE_XTALK_NOR_BUTTONID=-1 V_PRE_XTALK_ROW_BUTTONID=-1 ++++ V_PRE_XTALK_ROW_BUTTONID=-1 V_PRE_XTALK_MULTI_BUTTONID=-1 ++++ V_PRE_XTALK_MULTI_BUTTONID=-1 V_PRE_XTALK_MULTI_NSECTION= ++++ V_PRE_XTALK_MULTI_NSECTION= V_CAL_OVSCANNLOW=0 ++++ V_CAL_OVSCANNLOW=0 V_CAL_OVSCANNHIGH=1 ++++ V_CAL_OVSCANNHIGH=1 V_CAL_BIASNLOW=0 ++++ V_CAL_BIASNLOW=0 V_CAL_BIASNHIGH=1 ++++ V_CAL_BIASNHIGH=1 V_CAL_DARKNLOW=0 ++++ V_CAL_DARKNLOW=0 V_CAL_DARKNHIGH=1 ++++ V_CAL_DARKNHIGH=1 V_CAL_FLATNLOW=0 ++++ V_CAL_FLATNLOW=0 V_CAL_FLATNHIGH=1 ++++ V_CAL_FLATNHIGH=1 V_CAL_USEDARK=FALSE ++++ V_CAL_USEDARK=FALSE V_BACK_NLOW=0 ++++ V_BACK_NLOW=0 V_BACK_NHIGH=1 ++++ V_BACK_NHIGH=1 V_BACK_DETECTTHRESH=1.3 ++++ V_BACK_DETECTTHRESH=1.3 V_BACK_DETECTMINAREA=5 ++++ V_BACK_DETECTMINAREA=5 V_BACK_MASKEXPAND= ++++ V_BACK_MASKEXPAND= V_BACK_ILLUMSMOOTH= ++++ V_BACK_ILLUMSMOOTH= V_BACK_FRINGESMOOTH= ++++ V_BACK_FRINGESMOOTH= V_BACK_WINDOWSIZE=0 ++++ V_BACK_WINDOWSIZE=0 V_BACK_GAPSIZE=1.0 ++++ V_BACK_GAPSIZE=1.0 V_BACK_APPLYMODE=0 ++++ V_BACK_APPLYMODE=0 V_BACK_COMBINEMETHOD=0 ++++ V_BACK_COMBINEMETHOD=0 V_BACK_SEXFILTER=Y ++++ V_BACK_SEXFILTER=Y V_BACK_ADJUSTGAINS=N ++++ V_BACK_ADJUSTGAINS=N V_BACK_FRINGESCALE=Y ++++ V_BACK_FRINGESCALE=Y V_BACK_SKY2PASS=Y ++++ V_BACK_SKY2PASS=Y V_COLLDETECTTHRESH= ++++ V_COLLDETECTTHRESH= V_COLLDETECTMINAREA= ++++ V_COLLDETECTMINAREA= V_COLLMASKEXPAND= ++++ V_COLLMASKEXPAND= V_COLLREJECTTHRESH=1.5 ++++ V_COLLREJECTTHRESH=1.5 V_COLLDIRECTION=x ++++ V_COLLDIRECTION=x V_COLLXMIN= ++++ V_COLLXMIN= V_COLLXMAX= ++++ V_COLLXMAX= V_COLLYMIN= ++++ V_COLLYMIN= V_COLLYMAX= ++++ V_COLLYMAX= V_COLLMASKACTION=1 ++++ V_COLLMASKACTION=1 V_COLLAUTOTHRESHOLD=0 ++++ V_COLLAUTOTHRESHOLD=0 V_WEIGHTBINMIN=-100 ++++ V_WEIGHTBINMIN=-100 V_WEIGHTBINMAX=500 ++++ V_WEIGHTBINMAX=500 V_WEIGHTBINSIZE=4 ++++ V_WEIGHTBINSIZE=4 V_WEIGHTLOWTHRESHOLD= ++++ V_WEIGHTLOWTHRESHOLD= V_WEIGHTHIGHTHRESHOLD= ++++ V_WEIGHTHIGHTHRESHOLD= V_GLOBWFLATMIN=0.5 ++++ V_GLOBWFLATMIN=0.5 V_GLOBWFLATMAX=1.6 ++++ V_GLOBWFLATMAX=1.6 V_GLOBWDARKMIN= ++++ V_GLOBWDARKMIN= V_GLOBWDARKMAX= ++++ V_GLOBWDARKMAX= V_GLOBWDARKDIR= ++++ V_GLOBWDARKDIR= V_COSMICSTHRESHOLD=0.1 ++++ V_COSMICSTHRESHOLD=0.1 V_COSMICDT=6 ++++ V_COSMICDT=6 V_COSMICDMIN=1 ++++ V_COSMICDMIN=1 V_DEFECT_KERNELSIZE= ++++ V_DEFECT_KERNELSIZE= V_DEFECT_ROWTOL= ++++ V_DEFECT_ROWTOL= V_DEFECT_COLTOL= ++++ V_DEFECT_COLTOL= V_DEFECT_CLUSTOL= ++++ V_DEFECT_CLUSTOL= V_DEFECT_KERNELSIZE_SF= ++++ V_DEFECT_KERNELSIZE_SF= V_DEFECT_ROWTOL_SF= ++++ V_DEFECT_ROWTOL_SF= V_DEFECT_COLTOL_SF= ++++ V_DEFECT_COLTOL_SF= V_DEFECT_CLUSTOL_SF= ++++ V_DEFECT_CLUSTOL_SF= V_BLOOMLOWLIMIT=20000 ++++ V_BLOOMLOWLIMIT=20000 V_BLOOMMINAREA=100 ++++ V_BLOOMMINAREA=100 V_BLOOMWIDTH=0 ++++ V_BLOOMWIDTH=0 V_WEIGHT_BINOUTLIER=FALSE ++++ V_WEIGHT_BINOUTLIER=FALSE V_GLOBW_UNIFORMWEIGHT=FALSE ++++ V_GLOBW_UNIFORMWEIGHT=FALSE V_MASKBLOOMSPIKE=0 ++++ V_MASKBLOOMSPIKE=0 V_AP_DETECTTHRESH=5.0 ++++ V_AP_DETECTTHRESH=5.0 V_AP_DETECTMINAREA=5 ++++ V_AP_DETECTMINAREA=5 V_DEBLENDMINCONT=0.0005 ++++ V_DEBLENDMINCONT=0.0005 V_AP_LOWNUM=1 ++++ V_AP_LOWNUM=1 V_SEXCATMINFWHM=1.5 ++++ V_SEXCATMINFWHM=1.5 V_SEXCATFLAG=0 ++++ V_SEXCATFLAG=0 V_SATURATION= ++++ V_SATURATION= V_SEXBACKLEVEL= ++++ V_SEXBACKLEVEL= V_AP_FILTER=N ++++ V_AP_FILTER=N V_SCAMP_MAXPOSANGLE=2.0 ++++ V_SCAMP_MAXPOSANGLE=2.0 V_SCAMP_MAXOFFSET=2.0 ++++ V_SCAMP_MAXOFFSET=2.0 V_SCAMP_MAXSCALE=1.05 ++++ V_SCAMP_MAXSCALE=1.05 V_SCAMP_SNLOW=5 ++++ V_SCAMP_SNLOW=5 V_SCAMP_SNHIGH=20 ++++ V_SCAMP_SNHIGH=20 V_SCAMP_FGROUPRADIUS=1.0 ++++ V_SCAMP_FGROUPRADIUS=1.0 V_SCAMP_CROSSIDRADIUS=0 ++++ V_SCAMP_CROSSIDRADIUS=0 V_SCAMP_ASTREFWEIGHT=1.0 ++++ V_SCAMP_ASTREFWEIGHT=1.0 V_SCAMP_ASTRINSTRUKEY=FILTER ++++ V_SCAMP_ASTRINSTRUKEY=FILTER V_SCAMP_PHOTINSTRUKEY=FILTER ++++ V_SCAMP_PHOTINSTRUKEY=FILTER V_SCAMP_POLYORDER=3 ++++ V_SCAMP_POLYORDER=3 V_SCAMP_POLYORDER2= ++++ V_SCAMP_POLYORDER2= V_SCAMP_DISTORTGROUPS= ++++ V_SCAMP_DISTORTGROUPS= V_SCAMP_DISTORTKEYS= ++++ V_SCAMP_DISTORTKEYS= V_SCAMP_STABILITY=INSTRUMENT ++++ V_SCAMP_STABILITY=INSTRUMENT V_SCAMP_MOSAICTYPE=UNCHANGED ++++ V_SCAMP_MOSAICTYPE=UNCHANGED V_SCAMP_FOCALPLANE=DEFAULT ++++ V_SCAMP_FOCALPLANE=DEFAULT V_SCAMP_MATCHFLIPPED=N ++++ V_SCAMP_MATCHFLIPPED=N V_SCAMP_MATCHWCS=Y ++++ V_SCAMP_MATCHWCS=Y V_ANET_MAXSCALE=1.05 ++++ V_ANET_MAXSCALE=1.05 V_ANET_POLYORDER=3 ++++ V_ANET_POLYORDER=3 V_ANET_SCAMPDISTORT=N ++++ V_ANET_SCAMPDISTORT=N V_ABSPHOT_FILTER=B ++++ V_ABSPHOT_FILTER=B V_ABSPHOT_COLORFIXED=0.1 ++++ V_ABSPHOT_COLORFIXED=0.1 V_ABSPHOT_EXTINCTION=0 ++++ V_ABSPHOT_EXTINCTION=0 V_ABSPHOT_ZPMIN=24 ++++ V_ABSPHOT_ZPMIN=24 V_ABSPHOT_ZPMAX=25 ++++ V_ABSPHOT_ZPMAX=25 V_ABSPHOT_STDMINMAG=0.0 ++++ V_ABSPHOT_STDMINMAG=0.0 V_ABSPHOT_MAXPHOTERR=0.05 ++++ V_ABSPHOT_MAXPHOTERR=0.05 V_ABSPHOT_WCSERR=10 ++++ V_ABSPHOT_WCSERR=10 V_ABSPHOT_WCSERRCHECKBOX=Y ++++ V_ABSPHOT_WCSERRCHECKBOX=Y V_ABSPHOT_STDCAT_DIRECT=SDSS ++++ V_ABSPHOT_STDCAT_DIRECT=SDSS V_ABSPHOT_STDFILTERM2=u ++++ V_ABSPHOT_STDFILTERM2=u V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI ++++ V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI V_ABSPHOT_STDFILTER=U ++++ V_ABSPHOT_STDFILTER=U V_ABSPHOT_STDCOLOR=UmB ++++ V_ABSPHOT_STDCOLOR=UmB V_ABSPHOT_CALMODE=RUNCALIB ++++ V_ABSPHOT_CALMODE=RUNCALIB V_ABSPHOT_FITTINGMETHODM2=CHIP ++++ V_ABSPHOT_FITTINGMETHODM2=CHIP V_COADD_REFRA= ++++ V_COADD_REFRA= V_COADD_REFDEC= ++++ V_COADD_REFDEC= V_COADD_IDENT=(null) ++++ V_COADD_IDENT=(null) V_COADD_SEEING= ++++ V_COADD_SEEING= V_COADD_RZP= ++++ V_COADD_RZP= V_COADD_PIXSCALE=0.89 ++++ V_COADD_PIXSCALE=0.89 V_COADD_SIZEX= ++++ V_COADD_SIZEX= V_COADD_SIZEY= ++++ V_COADD_SIZEY= V_COADD_SKYPOSANGLE= ++++ V_COADD_SKYPOSANGLE= V_COADD_CHIPS="" ++++ V_COADD_CHIPS= V_COADD_FILTERTHRESHOLD= ++++ V_COADD_FILTERTHRESHOLD= V_COADD_FILTERCLUSTERSIZE= ++++ V_COADD_FILTERCLUSTERSIZE= V_COADD_FILTERBORDERWIDTH= ++++ V_COADD_FILTERBORDERWIDTH= V_COADD_SMOOTHEDGE= ++++ V_COADD_SMOOTHEDGE= V_COADD_PROPMOTIONRA= ++++ V_COADD_PROPMOTIONRA= V_COADD_PROPMOTIONDEC= ++++ V_COADD_PROPMOTIONDEC= V_COADD_KERNEL=LANCZOS3 ++++ V_COADD_KERNEL=LANCZOS3 V_COADD_FILTER= ++++ V_COADD_FILTER= V_COADD_PROJECTION=TAN ++++ V_COADD_PROJECTION=TAN V_COADD_CELESTIALTYPE=EQUATORIAL ++++ V_COADD_CELESTIALTYPE=EQUATORIAL V_COADD_COMBINETYPE=WEIGHTED ++++ V_COADD_COMBINETYPE=WEIGHTED . ~/.theli/param_set3.ini +++ . /home/albireo/.theli/param_set3.ini # Geographical latitude of the observatory OBSLAT= ++ OBSLAT= # Pixel scale in x and y [degrees] (CDELT1/2 keywords) PIXSCX=-0 ++ PIXSCX=-0 PIXSCY=0 ++ PIXSCY=0 # Pixel scale in arcsec PIXSCALE= ++ PIXSCALE= # Detector gain GAIN=1.0 ++ GAIN=1.0 # Overscan regions. The values are the x-range. # The min is given by OVSCANX1, the max by OVSCANX2. # We need separate regions for each chip. OVSCANX1=([1]=0) ++ OVSCANX1=([1]=0) OVSCANX2=([1]=0) ++ OVSCANX2=([1]=0) # Sections for the cutting. Select those pixels that receive light. # Give min values for x and y, and the extent in x and y. CUTX=([1]=1) ++ CUTX=([1]=1) CUTY=([1]=1) ++ CUTY=([1]=1) SIZEX=([1]=3358) ++ SIZEX=([1]=3358) SIZEY=([1]=2536) ++ SIZEY=([1]=2536) # Reference pixels for each chip. REFPIXX=([1]=100) ++ REFPIXX=([1]=100) REFPIXY=([1]=100) ++ REFPIXY=([1]=100) # Camera type TYPE=OPT ++ TYPE=OPT . bash.include + . bash.include function get_wcs { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do if [ "${NCHIPS}" = "1" ]; then ext=0 else ext=${k} fi # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good" ] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_wcsext { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do ext=1 # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good" ] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_statminsection { local xmin local ymin local xmax local ymax if [ "${STATXMIN[${CHIP}]}"_A != "_A" ] && \ [ "${STATXMAX[${CHIP}]}"_A != "_A" ] && \ [ "${STATYMIN[${CHIP}]}"_A != "_A" ] && \ [ "${STATYMAX[${CHIP}]}"_A != "_A" ]; then xmin=${STATXMIN[${CHIP}]} xmax=${STATXMAX[${CHIP}]} ymin=${STATYMIN[${CHIP}]} ymax=${STATYMAX[${CHIP}]} else k=1 statminsection="" xmax=${SIZEX[1]} ymax=${SIZEY[1]} while [ "${k}" -le "${NCHIPS}" ] do if [ ${SIZEX[k]} -lt ${xmax} ]; then xmax=${SIZEX[k]} fi if [ ${SIZEY[k]} -lt ${ymax} ]; then ymax=${SIZEY[k]} fi k=$((${k}+1)) done xmin=2 ymin=2 xmax=$((${xmax}-1)) ymax=$((${ymax}-1)) fi statminsection="${xmin},${xmax},${ymin},${ymax}" } function p_ls { ls -l | grep ^- | ${P_GAWK} '{print $NF}' } # the following extracts keywords; # leading and trailing whitespace is removed (2nd call to gawk) # and blanks inside the keyword value are replaced with underscores (e.g. for OBJECT keyword) function get_key { ${P_DFITS} $1 | ${P_FITSORT} -d "$2" | ${P_GAWK} 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' | ${P_GAWK} '{$1=$1}1' | sed 's/ /_/g' } function get_keyext { ${P_DFITS} -x 1 $1 | ${P_FITSORT} -d "$2" | ${P_GAWK} '{if ($2!="KEY_N/A") print $2}' } function get_region { # merge global and local region files if they exist regionfile="" region="" test -f ${TEMPDIR}/regionfile_${CHIP}.reg && rm ${TEMPDIR}/regionfile_${CHIP}.reg if [ -f ${INSTDIR}/${INSTRUMENT}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}.reg >> \ ${TEMPDIR}/regionfile_${CHIP}.reg regionfile="${TEMPDIR}/regionfile_${CHIP}.reg" region="-d ${regionfile}" fi if [ -f ${INSTDIR}/${INSTRUMENT}_${CHIP}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}_${CHIP}.reg >> \ ${TEMPDIR}/regionfile_${CHIP}.reg regionfile="${TEMPDIR}/regionfile_${CHIP}.reg" region="-d ${regionfile}" fi } function get_region_chip { # merge global and local region files if they exist # for a different chip than the currently processed one # (constant sky subtraction) regionfile="" region="" test -f ${TEMPDIR}/regionfile_$1.reg_$$ && rm ${TEMPDIR}/regionfile_$1.reg_$$ if [ -f ${INSTDIR}/${INSTRUMENT}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}.reg >> \ ${TEMPDIR}/regionfile_$1.reg_$$ regionfile="${TEMPDIR}/regionfile_$1.reg_$$" region="-d ${regionfile}" fi if [ -f ${INSTDIR}/${INSTRUMENT}_$1.reg ]; then cat ${INSTDIR}/${INSTRUMENT}_$1.reg >> \ ${TEMPDIR}/regionfile_$1.reg_$$ regionfile="${TEMPDIR}/regionfile_$1.reg_$$" region="-d ${regionfile}" fi } function prepare_sex_for_region { backstring="" weightstring="" if [ -f ${regionfile} ] && [ "${regionfile}_A" != "_A" ]; then # estimate the background, ignoring blanked pixels background=`${P_FITSSTAT} -i ${FILE} -m median -q` backstring="-BACK_TYPE MANUAL -BACK_VALUE ${background}" weightstring="-WEIGHT_TYPE MAP_WEIGHT -WEIGHT_IMAGE /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits" # check if the dummy weight exists if [ ! -e /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits ]; then mkdir -p /$1/WEIGHTS/ ${P_IC} -p -32 -c ${SIZEX[${CHIP}]} ${SIZEY[${CHIP}]} '1' > ${TEMPDIR}/dummy_${CHIP}.fits ${P_FITSPOLYGON} -i ${TEMPDIR}/dummy_${CHIP}.fits -p ${regionfile} \ -o /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits rm ${TEMPDIR}/dummy_${CHIP}.fits fi fi } function prepare_imred_for_backgroundmodel { # The relevant settings depending on the apply mode of the background model flatstring="" fringestring="" #################################################### # subtract background model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 0 ]; then fringestring="-FLAT N -FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE ${FRINGEIMAGE}" fi #################################################### # divide background model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 1 ]; then flatstring="-FRINGE N -FLAT Y -FLAT_THRESHHOLD 0.1 -FLAT_SCALE ${V_BACK_ADJUSTGAINS}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FLATIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FLATIMAGE=/$1/$2/$2_${CHIP}.fits fi flatstring="${flatstring} -FLAT_IMAGE ${FLATIMAGE} -FLAT_SCALEIMAGE ${FLATIMAGE}" fi #################################################### # subtract fringing model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 2 ]; then fringestring="-FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_FRINGESMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe_smooth.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE /$1/$2/$2_${CHIP}_illum.fits" fi #################################################### # divide background model, subtract fringing model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 3 ]; then flatstring="-FLAT Y -FLAT_THRESHHOLD 0.1 -FLAT_SCALE ${V_BACK_ADJUSTGAINS}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FLATIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FLATIMAGE=/$1/$2/$2_${CHIP}.fits fi flatstring="${flatstring} -FLAT_IMAGE ${FLATIMAGE} -FLAT_SCALEIMAGE ${FLATIMAGE}" fringestring="-FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_FRINGESMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe_smooth.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE /$1/$2/$2_${CHIP}_illum.fits" fi } # setting up the (sorted) list of flats rm ${TEMPDIR}/flatlist_$$ + rm /home/albireo/.theli/tmp//flatlist_17200 rm: cannot remove '/home/albireo/.theli/tmp//flatlist_17200': No such file or directory k=1 + k=1 while [ "${k}" -le "${NCHIPS}" ] do echo /$1/$2/$2_${k}.fits >> ${TEMPDIR}/flatlist_$$ k=$((${k}+1)) done + '[' 1 -le 1 ']' + echo //home/albireo/14-17.02.2015/rgb//flat/flat_1.fits + k=2 + '[' 2 -le 1 ']' get_statminsection + get_statminsection + local xmin + local ymin + local xmax + local ymax + '[' _A '!=' _A ']' + k=1 + statminsection= + xmax=3358 + ymax=2536 + '[' 1 -le 1 ']' + '[' 3358 -lt 3358 ']' + '[' 2536 -lt 2536 ']' + k=2 + '[' 2 -le 1 ']' + xmin=2 + ymin=2 + xmax=3357 + ymax=2535 + statminsection=2,3357,2,2535 statminsection=`echo ${statminsection} | sed 's/,/ /g'` echo ${statminsection} | sed 's/,/ /g' ++ echo 2,3357,2,2535 ++ sed 's/,/ /g' + statminsection='2 3357 2 2535' # get the statistics of all flats ${P_IMSTATS} \ -i ${TEMPDIR}/flatlist_$$ \ -s ${statminsection} \ -t -1e8 1e8 \ -o ${TEMPDIR}/flatlist_imstats0_$$ + /home/albireo/THELI/theli//bin/Linux_64//imstats -i /home/albireo/.theli/tmp//flatlist_17200 -s 2 3357 2 2535 -t -1e8 1e8 -o /home/albireo/.theli/tmp//flatlist_imstats0_17200 # extract the mode ${P_GAWK} '($1 !~ /#/) {print $2}' ${TEMPDIR}/flatlist_imstats0_$$ > \ ${TEMPDIR}/flatlist_imstats1_$$ + /usr/bin/gawk '($1 !~ /#/) {print $2}' /home/albireo/.theli/tmp//flatlist_imstats0_17200 # extract the maximum maxlevel=`${P_GETSTATS} -i ${TEMPDIR}/flatlist_imstats1_$$ -m max` ${P_GETSTATS} -i ${TEMPDIR}/flatlist_imstats1_$$ -m max ++ /home/albireo/THELI/theli//bin/Linux_64//getstats -i /home/albireo/.theli/tmp//flatlist_imstats1_17200 -m max ./create_flat_ratio.sh: line 45: /home/albireo/THELI/theli//bin/Linux_64//getstats: No such file or directory + maxlevel= # normalise the fluxes with respect to the highest flux ${P_GAWK} '{printf "%.4f\n", $1/'"${maxlevel}"'}' ${TEMPDIR}/flatlist_imstats1_$$ > $1/$2//flatratios.dat + /usr/bin/gawk '{printf "%.4f\n", $1/}' /home/albireo/.theli/tmp//flatlist_imstats1_17200 gawk: cmd. line:1: {printf "%.4f\n", $1/} gawk: cmd. line:1: ^ syntax error echo ${maxlevel} > $1/$2/.maxlevel.dat + echo #!/bin/bash -xv # the script creates normalised images. # It puts the result into new directories ..._norm. # The script assumes a threshold. Pixels under # this threshhold will be assigned a value of threshold # in the resulting images. # # 30.05.04: # temporary files go to a TEMPDIR directory # # 12.04.2005: # - I rewrote the script so that it works with the # 'imstats' program instead with the FLIPS based # 'immode'. # $1: main directory # $2: directory from which normaliced images should be created # $3: chips to be processed INSTDIR=instruments_professional + INSTDIR=instruments_professional if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=instruments_commercial fi + '[' '!' -f instruments_professional/Moravian10z.ini ']' + INSTDIR=instruments_commercial if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=~/.theli/instruments_user/ fi + '[' '!' -f instruments_commercial/Moravian10z.ini ']' + INSTDIR=/home/albireo/.theli/instruments_user/ . ${INSTDIR}/${INSTRUMENT:?}.ini + . /home/albireo/.theli/instruments_user//Moravian10z.ini # ************************************************************** # Instrument config file created by THELI # ************************************************************** # Name of the instrument (must be unique!) INSTRUMENT=Moravian10z ++ INSTRUMENT=Moravian10z # Number of chips NCHIPS=1 ++ NCHIPS=1 # We need to export these export INSTRUMENT ++ export INSTRUMENT export NCHIPS ++ export NCHIPS # Load directories for programs . progs.ini ++ . progs.ini # If we include this script in bash verbose and trace mode we do not want # all the (verbose) output from this script: USE_X=0 +++ USE_X=0 if [[ ${-/x} != $- ]] ; then USE_X=1 set +xv fi +++ [[ hvB != hvxB ]] +++ USE_X=1 +++ set +xv # # this is needed for the GUI # P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh +++ P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh P_GETPMMOBJECTS=${BIN}/get_pmm_objects +++ P_GETPMMOBJECTS=/home/albireo/THELI/theli//bin/Linux_64//get_pmm_objects P_LDACDESC=${BIN}/ldacdesc +++ P_LDACDESC=/home/albireo/THELI/theli//bin/Linux_64//ldacdesc P_LDACDELTAB=${BIN}/ldacdeltab +++ P_LDACDELTAB=/home/albireo/THELI/theli//bin/Linux_64//ldacdeltab P_LDACDELKEY=${BIN}/ldacdelkey +++ P_LDACDELKEY=/home/albireo/THELI/theli//bin/Linux_64//ldacdelkey P_LDACPUTXY=${BIN}/ldacputxy +++ P_LDACPUTXY=/home/albireo/THELI/theli//bin/Linux_64//ldacputxy P_FITSFLIP=${BIN}/fitsflip +++ P_FITSFLIP=/home/albireo/THELI/theli//bin/Linux_64//fitsflip P_DEBLOOM=${BIN}/debloom +++ P_DEBLOOM=/home/albireo/THELI/theli//bin/Linux_64//debloom P_EQUALISEBAYERFLAT=${BIN}/equalisebayerflat +++ P_EQUALISEBAYERFLAT=/home/albireo/THELI/theli//bin/Linux_64//equalisebayerflat P_FITSCOLLAPSE=${BIN}/fitscollapse +++ P_FITSCOLLAPSE=/home/albireo/THELI/theli//bin/Linux_64//fitscollapse P_FITSWRAP=${BIN}/fitswrap +++ P_FITSWRAP=/home/albireo/THELI/theli//bin/Linux_64//fitswrap P_FITSCUBESLICE=${BIN}/fitscubeslice +++ P_FITSCUBESLICE=/home/albireo/THELI/theli//bin/Linux_64//fitscubeslice P_FITSADDKEY=${BIN}/fitsaddkey +++ P_FITSADDKEY=/home/albireo/THELI/theli//bin/Linux_64//fitsaddkey P_FITSCUT=${BIN}/fitscut +++ P_FITSCUT=/home/albireo/THELI/theli//bin/Linux_64//fitscut P_FITSBIN=${BIN}/fitsbin +++ P_FITSBIN=/home/albireo/THELI/theli//bin/Linux_64//fitsbin P_FITSBLOOMDETECT=${BIN}/fitsbloomdetect +++ P_FITSBLOOMDETECT=/home/albireo/THELI/theli//bin/Linux_64//fitsbloomdetect P_FITSFILL=${BIN}/fitsfill +++ P_FITSFILL=/home/albireo/THELI/theli//bin/Linux_64//fitsfill P_FITSGRAD=${BIN}/fitsgrad +++ P_FITSGRAD=/home/albireo/THELI/theli//bin/Linux_64//fitsgrad P_FITSDEMOSAICBAYER=${BIN}/fitsdemosaicbayer +++ P_FITSDEMOSAICBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitsdemosaicbayer P_FITS3MAX=${BIN}/fits3max +++ P_FITS3MAX=/home/albireo/THELI/theli//bin/Linux_64//fits3max P_FITSGAUSS=${BIN}/fitsgauss +++ P_FITSGAUSS=/home/albireo/THELI/theli//bin/Linux_64//fitsgauss P_FITS2TIFF=${BIN}/fits2tiff +++ P_FITS2TIFF=/home/albireo/THELI/theli//bin/Linux_64//fits2tiff P_TIFF2FITS=${BIN}/tiff2fits +++ P_TIFF2FITS=/home/albireo/THELI/theli//bin/Linux_64//tiff2fits P_CORRLIRIS=${BIN}/corrliris +++ P_CORRLIRIS=/home/albireo/THELI/theli//bin/Linux_64//corrliris P_CROSSTALK=${BIN}/crosstalk +++ P_CROSSTALK=/home/albireo/THELI/theli//bin/Linux_64//crosstalk P_FINDGSC=${BIN}/findgsc +++ P_FINDGSC=/home/albireo/THELI/theli//bin/Linux_64//findgsc P_FINDUSNOB1=${BIN}/findusnob1 +++ P_FINDUSNOB1=/home/albireo/THELI/theli//bin/Linux_64//findusnob1 P_FIND2MASS=${BIN}/find2mass +++ P_FIND2MASS=/home/albireo/THELI/theli//bin/Linux_64//find2mass P_FINDSDSS7=${BIN}/findsdss7 +++ P_FINDSDSS7=/home/albireo/THELI/theli//bin/Linux_64//findsdss7 P_FINDUCAC3=${BIN}/finducac3 +++ P_FINDUCAC3=/home/albireo/THELI/theli//bin/Linux_64//finducac3 P_FINDNOMAD1=${BIN}/findnomad1 +++ P_FINDNOMAD1=/home/albireo/THELI/theli//bin/Linux_64//findnomad1 P_VIZQUERY=${BIN}/vizquery +++ P_VIZQUERY=/home/albireo/THELI/theli//bin/Linux_64//vizquery P_SESAME=${BIN}/sesame +++ P_SESAME=/home/albireo/THELI/theli//bin/Linux_64//sesame P_ACLIENT=${BIN}/aclient +++ P_ACLIENT=/home/albireo/THELI/theli//bin/Linux_64//aclient P_SCAMP=${BIN}/scamp +++ P_SCAMP=/home/albireo/THELI/theli//bin/Linux_64//scamp P_RENAME=${BIN}/rename +++ P_RENAME=/home/albireo/THELI/theli//bin/Linux_64//rename P_MISSFITS=${BIN}/missfits +++ P_MISSFITS=/home/albireo/THELI/theli//bin/Linux_64//missfits P_SWARPFILTER=${BIN}/swarpfilter +++ P_SWARPFILTER=/home/albireo/THELI/theli//bin/Linux_64//swarpfilter P_FITSBLOCKEDIT=${BIN}/fitsblockedit +++ P_FITSBLOCKEDIT=/home/albireo/THELI/theli//bin/Linux_64//fitsblockedit P_FITSMASK=${BIN}/fitsmask +++ P_FITSMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsmask P_FITSMEDIAN=${BIN}/fitsmedian +++ P_FITSMEDIAN=/home/albireo/THELI/theli//bin/Linux_64//fitsmedian P_FITSNANMASK=${BIN}/fitsnanmask +++ P_FITSNANMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsnanmask P_FITSSTAT=${BIN}/fitsstat +++ P_FITSSTAT=/home/albireo/THELI/theli//bin/Linux_64//fitsstat P_FITSSPLITBAYER=${BIN}/fitssplitbayer +++ P_FITSSPLITBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitssplitbayer P_FITSGUESSGEOM=${BIN}/fitsguessgeom +++ P_FITSGUESSGEOM=/home/albireo/THELI/theli//bin/Linux_64//fitsguessgeom P_SPLITSUPRIMECAM=${BIN}/splitsuprimecam +++ P_SPLITSUPRIMECAM=/home/albireo/THELI/theli//bin/Linux_64//splitsuprimecam P_PASTEDUALCHANNEL=${BIN}/paste_dualchannel +++ P_PASTEDUALCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_dualchannel P_PASTEQUADCHANNEL=${BIN}/paste_quadchannel +++ P_PASTEQUADCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_quadchannel P_SPLITGMOSMULTIPORT=${BIN}/split_gmos_multiport +++ P_SPLITGMOSMULTIPORT=/home/albireo/THELI/theli//bin/Linux_64//split_gmos_multiport P_SPLITGROND=${BIN}/split_grond +++ P_SPLITGROND=/home/albireo/THELI/theli//bin/Linux_64//split_grond P_GETSHIFT=${BIN}/get_shift +++ P_GETSHIFT=/home/albireo/THELI/theli//bin/Linux_64//get_shift P_FITSNONLINEARITY=${BIN}/fitsnonlinearity +++ P_FITSNONLINEARITY=/home/albireo/THELI/theli//bin/Linux_64//fitsnonlinearity P_FITSPIC=${BIN}/fitspic +++ P_FITSPIC=/home/albireo/THELI/theli//bin/Linux_64//fitspic P_FITSPOLYGON=${BIN}/fitspolygon +++ P_FITSPOLYGON=/home/albireo/THELI/theli//bin/Linux_64//fitspolygon P_EXPANDSEXMASK=${BIN}/expand_sexmask +++ P_EXPANDSEXMASK=/home/albireo/THELI/theli//bin/Linux_64//expand_sexmask P_DCRAW=${BIN}/dcraw +++ P_DCRAW=/home/albireo/THELI/theli//bin/Linux_64//dcraw P_DECSEXCONV=${BIN}/decsexconv +++ P_DECSEXCONV=/home/albireo/THELI/theli//bin/Linux_64//decsexconv P_STIFF=${BIN}/stiff +++ P_STIFF=/home/albireo/THELI/theli//bin/Linux_64//stiff P_SUBSKY=${BIN}/subsky +++ P_SUBSKY=/home/albireo/THELI/theli//bin/Linux_64//subsky P_GIFSICLE=${BIN}/gifsicle +++ P_GIFSICLE=/home/albireo/THELI/theli//bin/Linux_64//gifsicle P_XY2SKY=${BIN}/xy2sky +++ P_XY2SKY=/home/albireo/THELI/theli//bin/Linux_64//xy2sky P_SKY2XY=${BIN}/sky2xy +++ P_SKY2XY=/home/albireo/THELI/theli//bin/Linux_64//sky2xy P_ISFITS=${BIN}/isfits +++ P_ISFITS=/home/albireo/THELI/theli//bin/Linux_64//isfits P_GETKEY=${BIN}/getkey +++ P_GETKEY=/home/albireo/THELI/theli//bin/Linux_64//getkey P_SKY2XY_HELPER=${BIN}/sky2xy_helper +++ P_SKY2XY_HELPER=/home/albireo/THELI/theli//bin/Linux_64//sky2xy_helper P_GET_ROTIMSIZE=${BIN}/get_rotimsize +++ P_GET_ROTIMSIZE=/home/albireo/THELI/theli//bin/Linux_64//get_rotimsize P_CLEANSTRING=${BIN}/cleanstring +++ P_CLEANSTRING=/home/albireo/THELI/theli//bin/Linux_64//cleanstring P_NUMCPU=${BIN}/numcpu +++ P_NUMCPU=/home/albireo/THELI/theli//bin/Linux_64//numcpu P_GETSTATS=${BIN}/getstats +++ P_GETSTATS=/home/albireo/THELI/theli//bin/Linux_64//getstats P_CDELT2CD=${BIN}/cdelt2cd +++ P_CDELT2CD=/home/albireo/THELI/theli//bin/Linux_64//cdelt2cd P_GETCDMATRIX=${BIN}/get_cdmatrix +++ P_GETCDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//get_cdmatrix P_GETPIXSCALE=${BIN}/get_pixscale +++ P_GETPIXSCALE=/home/albireo/THELI/theli//bin/Linux_64//get_pixscale P_GETPOSANGLE=${BIN}/get_posangle +++ P_GETPOSANGLE=/home/albireo/THELI/theli//bin/Linux_64//get_posangle P_GETREFCATRADIUS=${BIN}/get_refcat_radius +++ P_GETREFCATRADIUS=/home/albireo/THELI/theli//bin/Linux_64//get_refcat_radius P_MERGEDAT=${BIN}/mergedat +++ P_MERGEDAT=/home/albireo/THELI/theli//bin/Linux_64//mergedat P_ROTATECDMATRIX=${BIN}/rotate_cdmatrix +++ P_ROTATECDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//rotate_cdmatrix P_BACKGROUNDDYNAMIC=${BIN}/background_dynamic +++ P_BACKGROUNDDYNAMIC=/home/albireo/THELI/theli//bin/Linux_64//background_dynamic P_BACKGROUNDSTATIC=${BIN}/background_static +++ P_BACKGROUNDSTATIC=/home/albireo/THELI/theli//bin/Linux_64//background_static P_MJD=${BIN}/mjd +++ P_MJD=/home/albireo/THELI/theli//bin/Linux_64//mjd P_LST=${BIN}/lst +++ P_LST=/home/albireo/THELI/theli//bin/Linux_64//lst P_IMALYZER=${BIN}/imalyzer +++ P_IMALYZER=/home/albireo/THELI/theli//bin/Linux_64//imalyzer P_E2ELL=${BIN}/e2ell +++ P_E2ELL=/home/albireo/THELI/theli//bin/Linux_64//e2ell P_FITSSMOOTHEDGE=${BIN}/fitssmoothedge +++ P_FITSSMOOTHEDGE=/home/albireo/THELI/theli//bin/Linux_64//fitssmoothedge P_FITSEXTRACTEXT=${BIN}/fitsextractext +++ P_FITSEXTRACTEXT=/home/albireo/THELI/theli//bin/Linux_64//fitsextractext P_ERRTEST=${BIN}/errtest +++ P_ERRTEST=/home/albireo/THELI/theli//bin/Linux_64//errtest P_BUILDINDEX=${BIN_ANET}/build-astrometry-index +++ P_BUILDINDEX=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//build-astrometry-index P_SOLVEFIELD=${BIN_ANET}/solve-field +++ P_SOLVEFIELD=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//solve-field P_XCORR=${BIN}/xcorr +++ P_XCORR=/home/albireo/THELI/theli//bin/Linux_64//xcorr # # load the parameter files # . ~/.theli/param_set1.ini +++ . /home/albireo/.theli/param_set1.ini PROJECTNAME=Pferdekopf ++++ PROJECTNAME=Pferdekopf NPARA=1 ++++ NPARA=1 GUIVERSION=2.9.5 ++++ GUIVERSION=2.9.5 KERNEL=Linux ++++ KERNEL=Linux NFRAMES=47 ++++ NFRAMES=47 V_PRE_SPLITMIRCUBE=N ++++ V_PRE_SPLITMIRCUBE=N V_DO_BIAS=N ++++ V_DO_BIAS=N V_DO_FLAT=Y ++++ V_DO_FLAT=Y V_NONLINCORR=N ++++ V_NONLINCORR=N V_AP_SERVER=vizier.u-strasbg.fr ++++ V_AP_SERVER=vizier.u-strasbg.fr V_AP_MAGLIM=20 ++++ V_AP_MAGLIM=20 V_AP_RADIUS=30 ++++ V_AP_RADIUS=30 V_AP_REFRA=header ++++ V_AP_REFRA=header V_AP_REFDEC=header ++++ V_AP_REFDEC=header . ~/.theli/param_set2.ini +++ . /home/albireo/.theli/param_set2.ini V_PRE_RENAME_CHECKED=0 ++++ V_PRE_RENAME_CHECKED=0 V_RENAME_FITSKEY=ARCFILE ++++ V_RENAME_FITSKEY=ARCFILE V_SORT_FITSKEY=OBJECT ++++ V_SORT_FITSKEY=OBJECT V_SORT_BIAS= ++++ V_SORT_BIAS= V_SORT_DARK= ++++ V_SORT_DARK= V_SORT_DOMEFLAT= ++++ V_SORT_DOMEFLAT= V_SORT_SKYFLAT= ++++ V_SORT_SKYFLAT= V_SORT_STD= ++++ V_SORT_STD= V_PRE_XTALK_NOR_CHECKED=0 ++++ V_PRE_XTALK_NOR_CHECKED=0 V_PRE_XTALK_ROW_CHECKED=0 ++++ V_PRE_XTALK_ROW_CHECKED=0 V_PRE_XTALK_MULTI_CHECKED=0 ++++ V_PRE_XTALK_MULTI_CHECKED=0 V_PRE_XTALK_NOR_AMPLITUDE= ++++ V_PRE_XTALK_NOR_AMPLITUDE= V_PRE_XTALK_ROW_AMPLITUDE= ++++ V_PRE_XTALK_ROW_AMPLITUDE= V_PRE_XTALK_NOR_BUTTONID=-1 ++++ V_PRE_XTALK_NOR_BUTTONID=-1 V_PRE_XTALK_ROW_BUTTONID=-1 ++++ V_PRE_XTALK_ROW_BUTTONID=-1 V_PRE_XTALK_MULTI_BUTTONID=-1 ++++ V_PRE_XTALK_MULTI_BUTTONID=-1 V_PRE_XTALK_MULTI_NSECTION= ++++ V_PRE_XTALK_MULTI_NSECTION= V_CAL_OVSCANNLOW=0 ++++ V_CAL_OVSCANNLOW=0 V_CAL_OVSCANNHIGH=1 ++++ V_CAL_OVSCANNHIGH=1 V_CAL_BIASNLOW=0 ++++ V_CAL_BIASNLOW=0 V_CAL_BIASNHIGH=1 ++++ V_CAL_BIASNHIGH=1 V_CAL_DARKNLOW=0 ++++ V_CAL_DARKNLOW=0 V_CAL_DARKNHIGH=1 ++++ V_CAL_DARKNHIGH=1 V_CAL_FLATNLOW=0 ++++ V_CAL_FLATNLOW=0 V_CAL_FLATNHIGH=1 ++++ V_CAL_FLATNHIGH=1 V_CAL_USEDARK=FALSE ++++ V_CAL_USEDARK=FALSE V_BACK_NLOW=0 ++++ V_BACK_NLOW=0 V_BACK_NHIGH=1 ++++ V_BACK_NHIGH=1 V_BACK_DETECTTHRESH=1.3 ++++ V_BACK_DETECTTHRESH=1.3 V_BACK_DETECTMINAREA=5 ++++ V_BACK_DETECTMINAREA=5 V_BACK_MASKEXPAND= ++++ V_BACK_MASKEXPAND= V_BACK_ILLUMSMOOTH= ++++ V_BACK_ILLUMSMOOTH= V_BACK_FRINGESMOOTH= ++++ V_BACK_FRINGESMOOTH= V_BACK_WINDOWSIZE=0 ++++ V_BACK_WINDOWSIZE=0 V_BACK_GAPSIZE=1.0 ++++ V_BACK_GAPSIZE=1.0 V_BACK_APPLYMODE=0 ++++ V_BACK_APPLYMODE=0 V_BACK_COMBINEMETHOD=0 ++++ V_BACK_COMBINEMETHOD=0 V_BACK_SEXFILTER=Y ++++ V_BACK_SEXFILTER=Y V_BACK_ADJUSTGAINS=N ++++ V_BACK_ADJUSTGAINS=N V_BACK_FRINGESCALE=Y ++++ V_BACK_FRINGESCALE=Y V_BACK_SKY2PASS=Y ++++ V_BACK_SKY2PASS=Y V_COLLDETECTTHRESH= ++++ V_COLLDETECTTHRESH= V_COLLDETECTMINAREA= ++++ V_COLLDETECTMINAREA= V_COLLMASKEXPAND= ++++ V_COLLMASKEXPAND= V_COLLREJECTTHRESH=1.5 ++++ V_COLLREJECTTHRESH=1.5 V_COLLDIRECTION=x ++++ V_COLLDIRECTION=x V_COLLXMIN= ++++ V_COLLXMIN= V_COLLXMAX= ++++ V_COLLXMAX= V_COLLYMIN= ++++ V_COLLYMIN= V_COLLYMAX= ++++ V_COLLYMAX= V_COLLMASKACTION=1 ++++ V_COLLMASKACTION=1 V_COLLAUTOTHRESHOLD=0 ++++ V_COLLAUTOTHRESHOLD=0 V_WEIGHTBINMIN=-100 ++++ V_WEIGHTBINMIN=-100 V_WEIGHTBINMAX=500 ++++ V_WEIGHTBINMAX=500 V_WEIGHTBINSIZE=4 ++++ V_WEIGHTBINSIZE=4 V_WEIGHTLOWTHRESHOLD= ++++ V_WEIGHTLOWTHRESHOLD= V_WEIGHTHIGHTHRESHOLD= ++++ V_WEIGHTHIGHTHRESHOLD= V_GLOBWFLATMIN=0.5 ++++ V_GLOBWFLATMIN=0.5 V_GLOBWFLATMAX=1.6 ++++ V_GLOBWFLATMAX=1.6 V_GLOBWDARKMIN= ++++ V_GLOBWDARKMIN= V_GLOBWDARKMAX= ++++ V_GLOBWDARKMAX= V_GLOBWDARKDIR= ++++ V_GLOBWDARKDIR= V_COSMICSTHRESHOLD=0.1 ++++ V_COSMICSTHRESHOLD=0.1 V_COSMICDT=6 ++++ V_COSMICDT=6 V_COSMICDMIN=1 ++++ V_COSMICDMIN=1 V_DEFECT_KERNELSIZE= ++++ V_DEFECT_KERNELSIZE= V_DEFECT_ROWTOL= ++++ V_DEFECT_ROWTOL= V_DEFECT_COLTOL= ++++ V_DEFECT_COLTOL= V_DEFECT_CLUSTOL= ++++ V_DEFECT_CLUSTOL= V_DEFECT_KERNELSIZE_SF= ++++ V_DEFECT_KERNELSIZE_SF= V_DEFECT_ROWTOL_SF= ++++ V_DEFECT_ROWTOL_SF= V_DEFECT_COLTOL_SF= ++++ V_DEFECT_COLTOL_SF= V_DEFECT_CLUSTOL_SF= ++++ V_DEFECT_CLUSTOL_SF= V_BLOOMLOWLIMIT=20000 ++++ V_BLOOMLOWLIMIT=20000 V_BLOOMMINAREA=100 ++++ V_BLOOMMINAREA=100 V_BLOOMWIDTH=0 ++++ V_BLOOMWIDTH=0 V_WEIGHT_BINOUTLIER=FALSE ++++ V_WEIGHT_BINOUTLIER=FALSE V_GLOBW_UNIFORMWEIGHT=FALSE ++++ V_GLOBW_UNIFORMWEIGHT=FALSE V_MASKBLOOMSPIKE=0 ++++ V_MASKBLOOMSPIKE=0 V_AP_DETECTTHRESH=5.0 ++++ V_AP_DETECTTHRESH=5.0 V_AP_DETECTMINAREA=5 ++++ V_AP_DETECTMINAREA=5 V_DEBLENDMINCONT=0.0005 ++++ V_DEBLENDMINCONT=0.0005 V_AP_LOWNUM=1 ++++ V_AP_LOWNUM=1 V_SEXCATMINFWHM=1.5 ++++ V_SEXCATMINFWHM=1.5 V_SEXCATFLAG=0 ++++ V_SEXCATFLAG=0 V_SATURATION= ++++ V_SATURATION= V_SEXBACKLEVEL= ++++ V_SEXBACKLEVEL= V_AP_FILTER=N ++++ V_AP_FILTER=N V_SCAMP_MAXPOSANGLE=2.0 ++++ V_SCAMP_MAXPOSANGLE=2.0 V_SCAMP_MAXOFFSET=2.0 ++++ V_SCAMP_MAXOFFSET=2.0 V_SCAMP_MAXSCALE=1.05 ++++ V_SCAMP_MAXSCALE=1.05 V_SCAMP_SNLOW=5 ++++ V_SCAMP_SNLOW=5 V_SCAMP_SNHIGH=20 ++++ V_SCAMP_SNHIGH=20 V_SCAMP_FGROUPRADIUS=1.0 ++++ V_SCAMP_FGROUPRADIUS=1.0 V_SCAMP_CROSSIDRADIUS=0 ++++ V_SCAMP_CROSSIDRADIUS=0 V_SCAMP_ASTREFWEIGHT=1.0 ++++ V_SCAMP_ASTREFWEIGHT=1.0 V_SCAMP_ASTRINSTRUKEY=FILTER ++++ V_SCAMP_ASTRINSTRUKEY=FILTER V_SCAMP_PHOTINSTRUKEY=FILTER ++++ V_SCAMP_PHOTINSTRUKEY=FILTER V_SCAMP_POLYORDER=3 ++++ V_SCAMP_POLYORDER=3 V_SCAMP_POLYORDER2= ++++ V_SCAMP_POLYORDER2= V_SCAMP_DISTORTGROUPS= ++++ V_SCAMP_DISTORTGROUPS= V_SCAMP_DISTORTKEYS= ++++ V_SCAMP_DISTORTKEYS= V_SCAMP_STABILITY=INSTRUMENT ++++ V_SCAMP_STABILITY=INSTRUMENT V_SCAMP_MOSAICTYPE=UNCHANGED ++++ V_SCAMP_MOSAICTYPE=UNCHANGED V_SCAMP_FOCALPLANE=DEFAULT ++++ V_SCAMP_FOCALPLANE=DEFAULT V_SCAMP_MATCHFLIPPED=N ++++ V_SCAMP_MATCHFLIPPED=N V_SCAMP_MATCHWCS=Y ++++ V_SCAMP_MATCHWCS=Y V_ANET_MAXSCALE=1.05 ++++ V_ANET_MAXSCALE=1.05 V_ANET_POLYORDER=3 ++++ V_ANET_POLYORDER=3 V_ANET_SCAMPDISTORT=N ++++ V_ANET_SCAMPDISTORT=N V_ABSPHOT_FILTER=B ++++ V_ABSPHOT_FILTER=B V_ABSPHOT_COLORFIXED=0.1 ++++ V_ABSPHOT_COLORFIXED=0.1 V_ABSPHOT_EXTINCTION=0 ++++ V_ABSPHOT_EXTINCTION=0 V_ABSPHOT_ZPMIN=24 ++++ V_ABSPHOT_ZPMIN=24 V_ABSPHOT_ZPMAX=25 ++++ V_ABSPHOT_ZPMAX=25 V_ABSPHOT_STDMINMAG=0.0 ++++ V_ABSPHOT_STDMINMAG=0.0 V_ABSPHOT_MAXPHOTERR=0.05 ++++ V_ABSPHOT_MAXPHOTERR=0.05 V_ABSPHOT_WCSERR=10 ++++ V_ABSPHOT_WCSERR=10 V_ABSPHOT_WCSERRCHECKBOX=Y ++++ V_ABSPHOT_WCSERRCHECKBOX=Y V_ABSPHOT_STDCAT_DIRECT=SDSS ++++ V_ABSPHOT_STDCAT_DIRECT=SDSS V_ABSPHOT_STDFILTERM2=u ++++ V_ABSPHOT_STDFILTERM2=u V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI ++++ V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI V_ABSPHOT_STDFILTER=U ++++ V_ABSPHOT_STDFILTER=U V_ABSPHOT_STDCOLOR=UmB ++++ V_ABSPHOT_STDCOLOR=UmB V_ABSPHOT_CALMODE=RUNCALIB ++++ V_ABSPHOT_CALMODE=RUNCALIB V_ABSPHOT_FITTINGMETHODM2=CHIP ++++ V_ABSPHOT_FITTINGMETHODM2=CHIP V_COADD_REFRA= ++++ V_COADD_REFRA= V_COADD_REFDEC= ++++ V_COADD_REFDEC= V_COADD_IDENT=(null) ++++ V_COADD_IDENT=(null) V_COADD_SEEING= ++++ V_COADD_SEEING= V_COADD_RZP= ++++ V_COADD_RZP= V_COADD_PIXSCALE=0.89 ++++ V_COADD_PIXSCALE=0.89 V_COADD_SIZEX= ++++ V_COADD_SIZEX= V_COADD_SIZEY= ++++ V_COADD_SIZEY= V_COADD_SKYPOSANGLE= ++++ V_COADD_SKYPOSANGLE= V_COADD_CHIPS="" ++++ V_COADD_CHIPS= V_COADD_FILTERTHRESHOLD= ++++ V_COADD_FILTERTHRESHOLD= V_COADD_FILTERCLUSTERSIZE= ++++ V_COADD_FILTERCLUSTERSIZE= V_COADD_FILTERBORDERWIDTH= ++++ V_COADD_FILTERBORDERWIDTH= V_COADD_SMOOTHEDGE= ++++ V_COADD_SMOOTHEDGE= V_COADD_PROPMOTIONRA= ++++ V_COADD_PROPMOTIONRA= V_COADD_PROPMOTIONDEC= ++++ V_COADD_PROPMOTIONDEC= V_COADD_KERNEL=LANCZOS3 ++++ V_COADD_KERNEL=LANCZOS3 V_COADD_FILTER= ++++ V_COADD_FILTER= V_COADD_PROJECTION=TAN ++++ V_COADD_PROJECTION=TAN V_COADD_CELESTIALTYPE=EQUATORIAL ++++ V_COADD_CELESTIALTYPE=EQUATORIAL V_COADD_COMBINETYPE=WEIGHTED ++++ V_COADD_COMBINETYPE=WEIGHTED . ~/.theli/param_set3.ini +++ . /home/albireo/.theli/param_set3.ini # Geographical latitude of the observatory OBSLAT= ++ OBSLAT= # Pixel scale in x and y [degrees] (CDELT1/2 keywords) PIXSCX=-0 ++ PIXSCX=-0 PIXSCY=0 ++ PIXSCY=0 # Pixel scale in arcsec PIXSCALE= ++ PIXSCALE= # Detector gain GAIN=1.0 ++ GAIN=1.0 # Overscan regions. The values are the x-range. # The min is given by OVSCANX1, the max by OVSCANX2. # We need separate regions for each chip. OVSCANX1=([1]=0) ++ OVSCANX1=([1]=0) OVSCANX2=([1]=0) ++ OVSCANX2=([1]=0) # Sections for the cutting. Select those pixels that receive light. # Give min values for x and y, and the extent in x and y. CUTX=([1]=1) ++ CUTX=([1]=1) CUTY=([1]=1) ++ CUTY=([1]=1) SIZEX=([1]=3358) ++ SIZEX=([1]=3358) SIZEY=([1]=2536) ++ SIZEY=([1]=2536) # Reference pixels for each chip. REFPIXX=([1]=100) ++ REFPIXX=([1]=100) REFPIXY=([1]=100) ++ REFPIXY=([1]=100) # Camera type TYPE=OPT ++ TYPE=OPT # first calculate the modes: test -e ${TEMPDIR}/immode.dat_$$ && rm ${TEMPDIR}/immode.dat_$$ + test -e /home/albireo/.theli/tmp//immode.dat_17210 for CHIP in $3 do ${P_IMSTATS} -d 6 /$1/$2/$2_${CHIP}.fits -t -1e8 1e8 \ >> ${TEMPDIR}/immode.dat_$$ done + for CHIP in '$3' + /home/albireo/THELI/theli//bin/Linux_64//imstats -d 6 //home/albireo/14-17.02.2015/rgb//flat/flat_1.fits -t -1e8 1e8 # clean the modes (in case of strong vignetting they might be underestimated) # lower threshold: 25% of the initial mode estimate ${P_GAWK} '($1!="#") {print $2}' ${TEMPDIR}/immode.dat_$$ > \ ${TEMPDIR}/immode_clean.dat_$$ + /usr/bin/gawk '($1!="#") {print $2}' /home/albireo/.theli/tmp//immode.dat_17210 k=1 + k=1 for CHIP in $3 do thresh=`${P_GAWK} '(NR=='"${k}"') {print 0.25*$1}' ${TEMPDIR}/immode_clean.dat_$$` ${P_IMSTATS} -d 6 /$1/$2/$2_${CHIP}.fits -t ${thresh} 1e8 \ >> ${TEMPDIR}/immode_new.dat_$$ k=$((${k}+1)) done + for CHIP in '$3' ${P_GAWK} '(NR=='"${k}"') {print 0.25*$1}' ${TEMPDIR}/immode_clean.dat_$$ ++ /usr/bin/gawk '(NR==1) {print 0.25*$1}' /home/albireo/.theli/tmp//immode_clean.dat_17210 + thresh=5642.9 + /home/albireo/THELI/theli//bin/Linux_64//imstats -d 6 //home/albireo/14-17.02.2015/rgb//flat/flat_1.fits -t 5642.9 1e8 + k=2 # create the new directory for the normalized images if # it does not exist already if [ ! -d "/$1/$2_norm" ]; then mkdir "/$1/$2_norm" fi + '[' '!' -d //home/albireo/14-17.02.2015/rgb//flat_norm ']' + mkdir //home/albireo/14-17.02.2015/rgb//flat_norm MODES=`${P_GAWK} '($1!="#") {printf ("%f ", $2)}' ${TEMPDIR}/immode_new.dat_$$` ${P_GAWK} '($1!="#") {printf ("%f ", $2)}' ${TEMPDIR}/immode_new.dat_$$ ++ /usr/bin/gawk '($1!="#") {printf ("%f ", $2)}' /home/albireo/.theli/tmp//immode_new.dat_17210 + MODES='22571.605469 ' # Set the threshold to 20% of the smallest mode value THRESH=`${P_GAWK} '($1!="#") {print $2}' ${TEMPDIR}/immode_new.dat_$$ |\ ${P_GAWK} 'BEGIN{min=1000000} {if ($1 ?' /$1/$2/$2_${CHIP}.fits > \ ${RESULTDIR}/$2_norm_${CHIP}.fits i=$(( $i + 1 )) done + for CHIP in '$3' echo ${MODES} | ${P_GAWK} '{print $'${i}'}' ++ echo 22571.605469 ++ /usr/bin/gawk '{print $1}' + ACTUMODE=22571.605469 + '[' -L //home/albireo/14-17.02.2015/rgb//flat/flat_1.fits ']' + RESULTDIR=/home/albireo/14-17.02.2015/rgb//flat_norm/ + /home/albireo/THELI/theli//bin/Linux_64//ic -p -32 '%1 22571.605469 / 0.0 %1 4514.32 > ?' //home/albireo/14-17.02.2015/rgb//flat/flat_1.fits + i=2 #!/bin/bash -xv #$1: main directory (filter) #$2: Bias directory #$3: Flat directory #$4: Science directory #$5: chips to be processed # preliminary work: INSTDIR=instruments_professional + INSTDIR=instruments_professional if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=instruments_commercial fi + '[' '!' -f instruments_professional/Moravian10z.ini ']' + INSTDIR=instruments_commercial if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=~/.theli/instruments_user/ fi + '[' '!' -f instruments_commercial/Moravian10z.ini ']' + INSTDIR=/home/albireo/.theli/instruments_user/ . ${INSTDIR}/${INSTRUMENT:?}.ini + . /home/albireo/.theli/instruments_user//Moravian10z.ini # ************************************************************** # Instrument config file created by THELI # ************************************************************** # Name of the instrument (must be unique!) INSTRUMENT=Moravian10z ++ INSTRUMENT=Moravian10z # Number of chips NCHIPS=1 ++ NCHIPS=1 # We need to export these export INSTRUMENT ++ export INSTRUMENT export NCHIPS ++ export NCHIPS # Load directories for programs . progs.ini ++ . progs.ini # If we include this script in bash verbose and trace mode we do not want # all the (verbose) output from this script: USE_X=0 +++ USE_X=0 if [[ ${-/x} != $- ]] ; then USE_X=1 set +xv fi +++ [[ hvB != hvxB ]] +++ USE_X=1 +++ set +xv # # this is needed for the GUI # P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh +++ P_PIPELIBS=/home/albireo/THELI/PipeLibs//pipe.csh P_GETPMMOBJECTS=${BIN}/get_pmm_objects +++ P_GETPMMOBJECTS=/home/albireo/THELI/theli//bin/Linux_64//get_pmm_objects P_LDACDESC=${BIN}/ldacdesc +++ P_LDACDESC=/home/albireo/THELI/theli//bin/Linux_64//ldacdesc P_LDACDELTAB=${BIN}/ldacdeltab +++ P_LDACDELTAB=/home/albireo/THELI/theli//bin/Linux_64//ldacdeltab P_LDACDELKEY=${BIN}/ldacdelkey +++ P_LDACDELKEY=/home/albireo/THELI/theli//bin/Linux_64//ldacdelkey P_LDACPUTXY=${BIN}/ldacputxy +++ P_LDACPUTXY=/home/albireo/THELI/theli//bin/Linux_64//ldacputxy P_FITSFLIP=${BIN}/fitsflip +++ P_FITSFLIP=/home/albireo/THELI/theli//bin/Linux_64//fitsflip P_DEBLOOM=${BIN}/debloom +++ P_DEBLOOM=/home/albireo/THELI/theli//bin/Linux_64//debloom P_EQUALISEBAYERFLAT=${BIN}/equalisebayerflat +++ P_EQUALISEBAYERFLAT=/home/albireo/THELI/theli//bin/Linux_64//equalisebayerflat P_FITSCOLLAPSE=${BIN}/fitscollapse +++ P_FITSCOLLAPSE=/home/albireo/THELI/theli//bin/Linux_64//fitscollapse P_FITSWRAP=${BIN}/fitswrap +++ P_FITSWRAP=/home/albireo/THELI/theli//bin/Linux_64//fitswrap P_FITSCUBESLICE=${BIN}/fitscubeslice +++ P_FITSCUBESLICE=/home/albireo/THELI/theli//bin/Linux_64//fitscubeslice P_FITSADDKEY=${BIN}/fitsaddkey +++ P_FITSADDKEY=/home/albireo/THELI/theli//bin/Linux_64//fitsaddkey P_FITSCUT=${BIN}/fitscut +++ P_FITSCUT=/home/albireo/THELI/theli//bin/Linux_64//fitscut P_FITSBIN=${BIN}/fitsbin +++ P_FITSBIN=/home/albireo/THELI/theli//bin/Linux_64//fitsbin P_FITSBLOOMDETECT=${BIN}/fitsbloomdetect +++ P_FITSBLOOMDETECT=/home/albireo/THELI/theli//bin/Linux_64//fitsbloomdetect P_FITSFILL=${BIN}/fitsfill +++ P_FITSFILL=/home/albireo/THELI/theli//bin/Linux_64//fitsfill P_FITSGRAD=${BIN}/fitsgrad +++ P_FITSGRAD=/home/albireo/THELI/theli//bin/Linux_64//fitsgrad P_FITSDEMOSAICBAYER=${BIN}/fitsdemosaicbayer +++ P_FITSDEMOSAICBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitsdemosaicbayer P_FITS3MAX=${BIN}/fits3max +++ P_FITS3MAX=/home/albireo/THELI/theli//bin/Linux_64//fits3max P_FITSGAUSS=${BIN}/fitsgauss +++ P_FITSGAUSS=/home/albireo/THELI/theli//bin/Linux_64//fitsgauss P_FITS2TIFF=${BIN}/fits2tiff +++ P_FITS2TIFF=/home/albireo/THELI/theli//bin/Linux_64//fits2tiff P_TIFF2FITS=${BIN}/tiff2fits +++ P_TIFF2FITS=/home/albireo/THELI/theli//bin/Linux_64//tiff2fits P_CORRLIRIS=${BIN}/corrliris +++ P_CORRLIRIS=/home/albireo/THELI/theli//bin/Linux_64//corrliris P_CROSSTALK=${BIN}/crosstalk +++ P_CROSSTALK=/home/albireo/THELI/theli//bin/Linux_64//crosstalk P_FINDGSC=${BIN}/findgsc +++ P_FINDGSC=/home/albireo/THELI/theli//bin/Linux_64//findgsc P_FINDUSNOB1=${BIN}/findusnob1 +++ P_FINDUSNOB1=/home/albireo/THELI/theli//bin/Linux_64//findusnob1 P_FIND2MASS=${BIN}/find2mass +++ P_FIND2MASS=/home/albireo/THELI/theli//bin/Linux_64//find2mass P_FINDSDSS7=${BIN}/findsdss7 +++ P_FINDSDSS7=/home/albireo/THELI/theli//bin/Linux_64//findsdss7 P_FINDUCAC3=${BIN}/finducac3 +++ P_FINDUCAC3=/home/albireo/THELI/theli//bin/Linux_64//finducac3 P_FINDNOMAD1=${BIN}/findnomad1 +++ P_FINDNOMAD1=/home/albireo/THELI/theli//bin/Linux_64//findnomad1 P_VIZQUERY=${BIN}/vizquery +++ P_VIZQUERY=/home/albireo/THELI/theli//bin/Linux_64//vizquery P_SESAME=${BIN}/sesame +++ P_SESAME=/home/albireo/THELI/theli//bin/Linux_64//sesame P_ACLIENT=${BIN}/aclient +++ P_ACLIENT=/home/albireo/THELI/theli//bin/Linux_64//aclient P_SCAMP=${BIN}/scamp +++ P_SCAMP=/home/albireo/THELI/theli//bin/Linux_64//scamp P_RENAME=${BIN}/rename +++ P_RENAME=/home/albireo/THELI/theli//bin/Linux_64//rename P_MISSFITS=${BIN}/missfits +++ P_MISSFITS=/home/albireo/THELI/theli//bin/Linux_64//missfits P_SWARPFILTER=${BIN}/swarpfilter +++ P_SWARPFILTER=/home/albireo/THELI/theli//bin/Linux_64//swarpfilter P_FITSBLOCKEDIT=${BIN}/fitsblockedit +++ P_FITSBLOCKEDIT=/home/albireo/THELI/theli//bin/Linux_64//fitsblockedit P_FITSMASK=${BIN}/fitsmask +++ P_FITSMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsmask P_FITSMEDIAN=${BIN}/fitsmedian +++ P_FITSMEDIAN=/home/albireo/THELI/theli//bin/Linux_64//fitsmedian P_FITSNANMASK=${BIN}/fitsnanmask +++ P_FITSNANMASK=/home/albireo/THELI/theli//bin/Linux_64//fitsnanmask P_FITSSTAT=${BIN}/fitsstat +++ P_FITSSTAT=/home/albireo/THELI/theli//bin/Linux_64//fitsstat P_FITSSPLITBAYER=${BIN}/fitssplitbayer +++ P_FITSSPLITBAYER=/home/albireo/THELI/theli//bin/Linux_64//fitssplitbayer P_FITSGUESSGEOM=${BIN}/fitsguessgeom +++ P_FITSGUESSGEOM=/home/albireo/THELI/theli//bin/Linux_64//fitsguessgeom P_SPLITSUPRIMECAM=${BIN}/splitsuprimecam +++ P_SPLITSUPRIMECAM=/home/albireo/THELI/theli//bin/Linux_64//splitsuprimecam P_PASTEDUALCHANNEL=${BIN}/paste_dualchannel +++ P_PASTEDUALCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_dualchannel P_PASTEQUADCHANNEL=${BIN}/paste_quadchannel +++ P_PASTEQUADCHANNEL=/home/albireo/THELI/theli//bin/Linux_64//paste_quadchannel P_SPLITGMOSMULTIPORT=${BIN}/split_gmos_multiport +++ P_SPLITGMOSMULTIPORT=/home/albireo/THELI/theli//bin/Linux_64//split_gmos_multiport P_SPLITGROND=${BIN}/split_grond +++ P_SPLITGROND=/home/albireo/THELI/theli//bin/Linux_64//split_grond P_GETSHIFT=${BIN}/get_shift +++ P_GETSHIFT=/home/albireo/THELI/theli//bin/Linux_64//get_shift P_FITSNONLINEARITY=${BIN}/fitsnonlinearity +++ P_FITSNONLINEARITY=/home/albireo/THELI/theli//bin/Linux_64//fitsnonlinearity P_FITSPIC=${BIN}/fitspic +++ P_FITSPIC=/home/albireo/THELI/theli//bin/Linux_64//fitspic P_FITSPOLYGON=${BIN}/fitspolygon +++ P_FITSPOLYGON=/home/albireo/THELI/theli//bin/Linux_64//fitspolygon P_EXPANDSEXMASK=${BIN}/expand_sexmask +++ P_EXPANDSEXMASK=/home/albireo/THELI/theli//bin/Linux_64//expand_sexmask P_DCRAW=${BIN}/dcraw +++ P_DCRAW=/home/albireo/THELI/theli//bin/Linux_64//dcraw P_DECSEXCONV=${BIN}/decsexconv +++ P_DECSEXCONV=/home/albireo/THELI/theli//bin/Linux_64//decsexconv P_STIFF=${BIN}/stiff +++ P_STIFF=/home/albireo/THELI/theli//bin/Linux_64//stiff P_SUBSKY=${BIN}/subsky +++ P_SUBSKY=/home/albireo/THELI/theli//bin/Linux_64//subsky P_GIFSICLE=${BIN}/gifsicle +++ P_GIFSICLE=/home/albireo/THELI/theli//bin/Linux_64//gifsicle P_XY2SKY=${BIN}/xy2sky +++ P_XY2SKY=/home/albireo/THELI/theli//bin/Linux_64//xy2sky P_SKY2XY=${BIN}/sky2xy +++ P_SKY2XY=/home/albireo/THELI/theli//bin/Linux_64//sky2xy P_ISFITS=${BIN}/isfits +++ P_ISFITS=/home/albireo/THELI/theli//bin/Linux_64//isfits P_GETKEY=${BIN}/getkey +++ P_GETKEY=/home/albireo/THELI/theli//bin/Linux_64//getkey P_SKY2XY_HELPER=${BIN}/sky2xy_helper +++ P_SKY2XY_HELPER=/home/albireo/THELI/theli//bin/Linux_64//sky2xy_helper P_GET_ROTIMSIZE=${BIN}/get_rotimsize +++ P_GET_ROTIMSIZE=/home/albireo/THELI/theli//bin/Linux_64//get_rotimsize P_CLEANSTRING=${BIN}/cleanstring +++ P_CLEANSTRING=/home/albireo/THELI/theli//bin/Linux_64//cleanstring P_NUMCPU=${BIN}/numcpu +++ P_NUMCPU=/home/albireo/THELI/theli//bin/Linux_64//numcpu P_GETSTATS=${BIN}/getstats +++ P_GETSTATS=/home/albireo/THELI/theli//bin/Linux_64//getstats P_CDELT2CD=${BIN}/cdelt2cd +++ P_CDELT2CD=/home/albireo/THELI/theli//bin/Linux_64//cdelt2cd P_GETCDMATRIX=${BIN}/get_cdmatrix +++ P_GETCDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//get_cdmatrix P_GETPIXSCALE=${BIN}/get_pixscale +++ P_GETPIXSCALE=/home/albireo/THELI/theli//bin/Linux_64//get_pixscale P_GETPOSANGLE=${BIN}/get_posangle +++ P_GETPOSANGLE=/home/albireo/THELI/theli//bin/Linux_64//get_posangle P_GETREFCATRADIUS=${BIN}/get_refcat_radius +++ P_GETREFCATRADIUS=/home/albireo/THELI/theli//bin/Linux_64//get_refcat_radius P_MERGEDAT=${BIN}/mergedat +++ P_MERGEDAT=/home/albireo/THELI/theli//bin/Linux_64//mergedat P_ROTATECDMATRIX=${BIN}/rotate_cdmatrix +++ P_ROTATECDMATRIX=/home/albireo/THELI/theli//bin/Linux_64//rotate_cdmatrix P_BACKGROUNDDYNAMIC=${BIN}/background_dynamic +++ P_BACKGROUNDDYNAMIC=/home/albireo/THELI/theli//bin/Linux_64//background_dynamic P_BACKGROUNDSTATIC=${BIN}/background_static +++ P_BACKGROUNDSTATIC=/home/albireo/THELI/theli//bin/Linux_64//background_static P_MJD=${BIN}/mjd +++ P_MJD=/home/albireo/THELI/theli//bin/Linux_64//mjd P_LST=${BIN}/lst +++ P_LST=/home/albireo/THELI/theli//bin/Linux_64//lst P_IMALYZER=${BIN}/imalyzer +++ P_IMALYZER=/home/albireo/THELI/theli//bin/Linux_64//imalyzer P_E2ELL=${BIN}/e2ell +++ P_E2ELL=/home/albireo/THELI/theli//bin/Linux_64//e2ell P_FITSSMOOTHEDGE=${BIN}/fitssmoothedge +++ P_FITSSMOOTHEDGE=/home/albireo/THELI/theli//bin/Linux_64//fitssmoothedge P_FITSEXTRACTEXT=${BIN}/fitsextractext +++ P_FITSEXTRACTEXT=/home/albireo/THELI/theli//bin/Linux_64//fitsextractext P_ERRTEST=${BIN}/errtest +++ P_ERRTEST=/home/albireo/THELI/theli//bin/Linux_64//errtest P_BUILDINDEX=${BIN_ANET}/build-astrometry-index +++ P_BUILDINDEX=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//build-astrometry-index P_SOLVEFIELD=${BIN_ANET}/solve-field +++ P_SOLVEFIELD=/home/albireo/THELI/theli//gui/packages/astrometry/build/bin//solve-field P_XCORR=${BIN}/xcorr +++ P_XCORR=/home/albireo/THELI/theli//bin/Linux_64//xcorr # # load the parameter files # . ~/.theli/param_set1.ini +++ . /home/albireo/.theli/param_set1.ini PROJECTNAME=Pferdekopf ++++ PROJECTNAME=Pferdekopf NPARA=1 ++++ NPARA=1 GUIVERSION=2.9.5 ++++ GUIVERSION=2.9.5 KERNEL=Linux ++++ KERNEL=Linux NFRAMES=47 ++++ NFRAMES=47 V_PRE_SPLITMIRCUBE=N ++++ V_PRE_SPLITMIRCUBE=N V_DO_BIAS=N ++++ V_DO_BIAS=N V_DO_FLAT=Y ++++ V_DO_FLAT=Y V_NONLINCORR=N ++++ V_NONLINCORR=N V_AP_SERVER=vizier.u-strasbg.fr ++++ V_AP_SERVER=vizier.u-strasbg.fr V_AP_MAGLIM=20 ++++ V_AP_MAGLIM=20 V_AP_RADIUS=30 ++++ V_AP_RADIUS=30 V_AP_REFRA=header ++++ V_AP_REFRA=header V_AP_REFDEC=header ++++ V_AP_REFDEC=header . ~/.theli/param_set2.ini +++ . /home/albireo/.theli/param_set2.ini V_PRE_RENAME_CHECKED=0 ++++ V_PRE_RENAME_CHECKED=0 V_RENAME_FITSKEY=ARCFILE ++++ V_RENAME_FITSKEY=ARCFILE V_SORT_FITSKEY=OBJECT ++++ V_SORT_FITSKEY=OBJECT V_SORT_BIAS= ++++ V_SORT_BIAS= V_SORT_DARK= ++++ V_SORT_DARK= V_SORT_DOMEFLAT= ++++ V_SORT_DOMEFLAT= V_SORT_SKYFLAT= ++++ V_SORT_SKYFLAT= V_SORT_STD= ++++ V_SORT_STD= V_PRE_XTALK_NOR_CHECKED=0 ++++ V_PRE_XTALK_NOR_CHECKED=0 V_PRE_XTALK_ROW_CHECKED=0 ++++ V_PRE_XTALK_ROW_CHECKED=0 V_PRE_XTALK_MULTI_CHECKED=0 ++++ V_PRE_XTALK_MULTI_CHECKED=0 V_PRE_XTALK_NOR_AMPLITUDE= ++++ V_PRE_XTALK_NOR_AMPLITUDE= V_PRE_XTALK_ROW_AMPLITUDE= ++++ V_PRE_XTALK_ROW_AMPLITUDE= V_PRE_XTALK_NOR_BUTTONID=-1 ++++ V_PRE_XTALK_NOR_BUTTONID=-1 V_PRE_XTALK_ROW_BUTTONID=-1 ++++ V_PRE_XTALK_ROW_BUTTONID=-1 V_PRE_XTALK_MULTI_BUTTONID=-1 ++++ V_PRE_XTALK_MULTI_BUTTONID=-1 V_PRE_XTALK_MULTI_NSECTION= ++++ V_PRE_XTALK_MULTI_NSECTION= V_CAL_OVSCANNLOW=0 ++++ V_CAL_OVSCANNLOW=0 V_CAL_OVSCANNHIGH=1 ++++ V_CAL_OVSCANNHIGH=1 V_CAL_BIASNLOW=0 ++++ V_CAL_BIASNLOW=0 V_CAL_BIASNHIGH=1 ++++ V_CAL_BIASNHIGH=1 V_CAL_DARKNLOW=0 ++++ V_CAL_DARKNLOW=0 V_CAL_DARKNHIGH=1 ++++ V_CAL_DARKNHIGH=1 V_CAL_FLATNLOW=0 ++++ V_CAL_FLATNLOW=0 V_CAL_FLATNHIGH=1 ++++ V_CAL_FLATNHIGH=1 V_CAL_USEDARK=FALSE ++++ V_CAL_USEDARK=FALSE V_BACK_NLOW=0 ++++ V_BACK_NLOW=0 V_BACK_NHIGH=1 ++++ V_BACK_NHIGH=1 V_BACK_DETECTTHRESH=1.3 ++++ V_BACK_DETECTTHRESH=1.3 V_BACK_DETECTMINAREA=5 ++++ V_BACK_DETECTMINAREA=5 V_BACK_MASKEXPAND= ++++ V_BACK_MASKEXPAND= V_BACK_ILLUMSMOOTH= ++++ V_BACK_ILLUMSMOOTH= V_BACK_FRINGESMOOTH= ++++ V_BACK_FRINGESMOOTH= V_BACK_WINDOWSIZE=0 ++++ V_BACK_WINDOWSIZE=0 V_BACK_GAPSIZE=1.0 ++++ V_BACK_GAPSIZE=1.0 V_BACK_APPLYMODE=0 ++++ V_BACK_APPLYMODE=0 V_BACK_COMBINEMETHOD=0 ++++ V_BACK_COMBINEMETHOD=0 V_BACK_SEXFILTER=Y ++++ V_BACK_SEXFILTER=Y V_BACK_ADJUSTGAINS=N ++++ V_BACK_ADJUSTGAINS=N V_BACK_FRINGESCALE=Y ++++ V_BACK_FRINGESCALE=Y V_BACK_SKY2PASS=Y ++++ V_BACK_SKY2PASS=Y V_COLLDETECTTHRESH= ++++ V_COLLDETECTTHRESH= V_COLLDETECTMINAREA= ++++ V_COLLDETECTMINAREA= V_COLLMASKEXPAND= ++++ V_COLLMASKEXPAND= V_COLLREJECTTHRESH=1.5 ++++ V_COLLREJECTTHRESH=1.5 V_COLLDIRECTION=x ++++ V_COLLDIRECTION=x V_COLLXMIN= ++++ V_COLLXMIN= V_COLLXMAX= ++++ V_COLLXMAX= V_COLLYMIN= ++++ V_COLLYMIN= V_COLLYMAX= ++++ V_COLLYMAX= V_COLLMASKACTION=1 ++++ V_COLLMASKACTION=1 V_COLLAUTOTHRESHOLD=0 ++++ V_COLLAUTOTHRESHOLD=0 V_WEIGHTBINMIN=-100 ++++ V_WEIGHTBINMIN=-100 V_WEIGHTBINMAX=500 ++++ V_WEIGHTBINMAX=500 V_WEIGHTBINSIZE=4 ++++ V_WEIGHTBINSIZE=4 V_WEIGHTLOWTHRESHOLD= ++++ V_WEIGHTLOWTHRESHOLD= V_WEIGHTHIGHTHRESHOLD= ++++ V_WEIGHTHIGHTHRESHOLD= V_GLOBWFLATMIN=0.5 ++++ V_GLOBWFLATMIN=0.5 V_GLOBWFLATMAX=1.6 ++++ V_GLOBWFLATMAX=1.6 V_GLOBWDARKMIN= ++++ V_GLOBWDARKMIN= V_GLOBWDARKMAX= ++++ V_GLOBWDARKMAX= V_GLOBWDARKDIR= ++++ V_GLOBWDARKDIR= V_COSMICSTHRESHOLD=0.1 ++++ V_COSMICSTHRESHOLD=0.1 V_COSMICDT=6 ++++ V_COSMICDT=6 V_COSMICDMIN=1 ++++ V_COSMICDMIN=1 V_DEFECT_KERNELSIZE= ++++ V_DEFECT_KERNELSIZE= V_DEFECT_ROWTOL= ++++ V_DEFECT_ROWTOL= V_DEFECT_COLTOL= ++++ V_DEFECT_COLTOL= V_DEFECT_CLUSTOL= ++++ V_DEFECT_CLUSTOL= V_DEFECT_KERNELSIZE_SF= ++++ V_DEFECT_KERNELSIZE_SF= V_DEFECT_ROWTOL_SF= ++++ V_DEFECT_ROWTOL_SF= V_DEFECT_COLTOL_SF= ++++ V_DEFECT_COLTOL_SF= V_DEFECT_CLUSTOL_SF= ++++ V_DEFECT_CLUSTOL_SF= V_BLOOMLOWLIMIT=20000 ++++ V_BLOOMLOWLIMIT=20000 V_BLOOMMINAREA=100 ++++ V_BLOOMMINAREA=100 V_BLOOMWIDTH=0 ++++ V_BLOOMWIDTH=0 V_WEIGHT_BINOUTLIER=FALSE ++++ V_WEIGHT_BINOUTLIER=FALSE V_GLOBW_UNIFORMWEIGHT=FALSE ++++ V_GLOBW_UNIFORMWEIGHT=FALSE V_MASKBLOOMSPIKE=0 ++++ V_MASKBLOOMSPIKE=0 V_AP_DETECTTHRESH=5.0 ++++ V_AP_DETECTTHRESH=5.0 V_AP_DETECTMINAREA=5 ++++ V_AP_DETECTMINAREA=5 V_DEBLENDMINCONT=0.0005 ++++ V_DEBLENDMINCONT=0.0005 V_AP_LOWNUM=1 ++++ V_AP_LOWNUM=1 V_SEXCATMINFWHM=1.5 ++++ V_SEXCATMINFWHM=1.5 V_SEXCATFLAG=0 ++++ V_SEXCATFLAG=0 V_SATURATION= ++++ V_SATURATION= V_SEXBACKLEVEL= ++++ V_SEXBACKLEVEL= V_AP_FILTER=N ++++ V_AP_FILTER=N V_SCAMP_MAXPOSANGLE=2.0 ++++ V_SCAMP_MAXPOSANGLE=2.0 V_SCAMP_MAXOFFSET=2.0 ++++ V_SCAMP_MAXOFFSET=2.0 V_SCAMP_MAXSCALE=1.05 ++++ V_SCAMP_MAXSCALE=1.05 V_SCAMP_SNLOW=5 ++++ V_SCAMP_SNLOW=5 V_SCAMP_SNHIGH=20 ++++ V_SCAMP_SNHIGH=20 V_SCAMP_FGROUPRADIUS=1.0 ++++ V_SCAMP_FGROUPRADIUS=1.0 V_SCAMP_CROSSIDRADIUS=0 ++++ V_SCAMP_CROSSIDRADIUS=0 V_SCAMP_ASTREFWEIGHT=1.0 ++++ V_SCAMP_ASTREFWEIGHT=1.0 V_SCAMP_ASTRINSTRUKEY=FILTER ++++ V_SCAMP_ASTRINSTRUKEY=FILTER V_SCAMP_PHOTINSTRUKEY=FILTER ++++ V_SCAMP_PHOTINSTRUKEY=FILTER V_SCAMP_POLYORDER=3 ++++ V_SCAMP_POLYORDER=3 V_SCAMP_POLYORDER2= ++++ V_SCAMP_POLYORDER2= V_SCAMP_DISTORTGROUPS= ++++ V_SCAMP_DISTORTGROUPS= V_SCAMP_DISTORTKEYS= ++++ V_SCAMP_DISTORTKEYS= V_SCAMP_STABILITY=INSTRUMENT ++++ V_SCAMP_STABILITY=INSTRUMENT V_SCAMP_MOSAICTYPE=UNCHANGED ++++ V_SCAMP_MOSAICTYPE=UNCHANGED V_SCAMP_FOCALPLANE=DEFAULT ++++ V_SCAMP_FOCALPLANE=DEFAULT V_SCAMP_MATCHFLIPPED=N ++++ V_SCAMP_MATCHFLIPPED=N V_SCAMP_MATCHWCS=Y ++++ V_SCAMP_MATCHWCS=Y V_ANET_MAXSCALE=1.05 ++++ V_ANET_MAXSCALE=1.05 V_ANET_POLYORDER=3 ++++ V_ANET_POLYORDER=3 V_ANET_SCAMPDISTORT=N ++++ V_ANET_SCAMPDISTORT=N V_ABSPHOT_FILTER=B ++++ V_ABSPHOT_FILTER=B V_ABSPHOT_COLORFIXED=0.1 ++++ V_ABSPHOT_COLORFIXED=0.1 V_ABSPHOT_EXTINCTION=0 ++++ V_ABSPHOT_EXTINCTION=0 V_ABSPHOT_ZPMIN=24 ++++ V_ABSPHOT_ZPMIN=24 V_ABSPHOT_ZPMAX=25 ++++ V_ABSPHOT_ZPMAX=25 V_ABSPHOT_STDMINMAG=0.0 ++++ V_ABSPHOT_STDMINMAG=0.0 V_ABSPHOT_MAXPHOTERR=0.05 ++++ V_ABSPHOT_MAXPHOTERR=0.05 V_ABSPHOT_WCSERR=10 ++++ V_ABSPHOT_WCSERR=10 V_ABSPHOT_WCSERRCHECKBOX=Y ++++ V_ABSPHOT_WCSERRCHECKBOX=Y V_ABSPHOT_STDCAT_DIRECT=SDSS ++++ V_ABSPHOT_STDCAT_DIRECT=SDSS V_ABSPHOT_STDFILTERM2=u ++++ V_ABSPHOT_STDFILTERM2=u V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI ++++ V_ABSPHOT_STDCAT_INDIRECT=LANDOLT_UBVRI V_ABSPHOT_STDFILTER=U ++++ V_ABSPHOT_STDFILTER=U V_ABSPHOT_STDCOLOR=UmB ++++ V_ABSPHOT_STDCOLOR=UmB V_ABSPHOT_CALMODE=RUNCALIB ++++ V_ABSPHOT_CALMODE=RUNCALIB V_ABSPHOT_FITTINGMETHODM2=CHIP ++++ V_ABSPHOT_FITTINGMETHODM2=CHIP V_COADD_REFRA= ++++ V_COADD_REFRA= V_COADD_REFDEC= ++++ V_COADD_REFDEC= V_COADD_IDENT=(null) ++++ V_COADD_IDENT=(null) V_COADD_SEEING= ++++ V_COADD_SEEING= V_COADD_RZP= ++++ V_COADD_RZP= V_COADD_PIXSCALE=0.89 ++++ V_COADD_PIXSCALE=0.89 V_COADD_SIZEX= ++++ V_COADD_SIZEX= V_COADD_SIZEY= ++++ V_COADD_SIZEY= V_COADD_SKYPOSANGLE= ++++ V_COADD_SKYPOSANGLE= V_COADD_CHIPS="" ++++ V_COADD_CHIPS= V_COADD_FILTERTHRESHOLD= ++++ V_COADD_FILTERTHRESHOLD= V_COADD_FILTERCLUSTERSIZE= ++++ V_COADD_FILTERCLUSTERSIZE= V_COADD_FILTERBORDERWIDTH= ++++ V_COADD_FILTERBORDERWIDTH= V_COADD_SMOOTHEDGE= ++++ V_COADD_SMOOTHEDGE= V_COADD_PROPMOTIONRA= ++++ V_COADD_PROPMOTIONRA= V_COADD_PROPMOTIONDEC= ++++ V_COADD_PROPMOTIONDEC= V_COADD_KERNEL=LANCZOS3 ++++ V_COADD_KERNEL=LANCZOS3 V_COADD_FILTER= ++++ V_COADD_FILTER= V_COADD_PROJECTION=TAN ++++ V_COADD_PROJECTION=TAN V_COADD_CELESTIALTYPE=EQUATORIAL ++++ V_COADD_CELESTIALTYPE=EQUATORIAL V_COADD_COMBINETYPE=WEIGHTED ++++ V_COADD_COMBINETYPE=WEIGHTED . ~/.theli/param_set3.ini +++ . /home/albireo/.theli/param_set3.ini # Geographical latitude of the observatory OBSLAT= ++ OBSLAT= # Pixel scale in x and y [degrees] (CDELT1/2 keywords) PIXSCX=-0 ++ PIXSCX=-0 PIXSCY=0 ++ PIXSCY=0 # Pixel scale in arcsec PIXSCALE= ++ PIXSCALE= # Detector gain GAIN=1.0 ++ GAIN=1.0 # Overscan regions. The values are the x-range. # The min is given by OVSCANX1, the max by OVSCANX2. # We need separate regions for each chip. OVSCANX1=([1]=0) ++ OVSCANX1=([1]=0) OVSCANX2=([1]=0) ++ OVSCANX2=([1]=0) # Sections for the cutting. Select those pixels that receive light. # Give min values for x and y, and the extent in x and y. CUTX=([1]=1) ++ CUTX=([1]=1) CUTY=([1]=1) ++ CUTY=([1]=1) SIZEX=([1]=3358) ++ SIZEX=([1]=3358) SIZEY=([1]=2536) ++ SIZEY=([1]=2536) # Reference pixels for each chip. REFPIXX=([1]=100) ++ REFPIXX=([1]=100) REFPIXY=([1]=100) ++ REFPIXY=([1]=100) # Camera type TYPE=OPT ++ TYPE=OPT . bash.include + . bash.include function get_wcs { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do if [ "${NCHIPS}" = "1" ]; then ext=0 else ext=${k} fi # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good" ] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_wcsext { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do ext=1 # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good" ] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_statminsection { local xmin local ymin local xmax local ymax if [ "${STATXMIN[${CHIP}]}"_A != "_A" ] && \ [ "${STATXMAX[${CHIP}]}"_A != "_A" ] && \ [ "${STATYMIN[${CHIP}]}"_A != "_A" ] && \ [ "${STATYMAX[${CHIP}]}"_A != "_A" ]; then xmin=${STATXMIN[${CHIP}]} xmax=${STATXMAX[${CHIP}]} ymin=${STATYMIN[${CHIP}]} ymax=${STATYMAX[${CHIP}]} else k=1 statminsection="" xmax=${SIZEX[1]} ymax=${SIZEY[1]} while [ "${k}" -le "${NCHIPS}" ] do if [ ${SIZEX[k]} -lt ${xmax} ]; then xmax=${SIZEX[k]} fi if [ ${SIZEY[k]} -lt ${ymax} ]; then ymax=${SIZEY[k]} fi k=$((${k}+1)) done xmin=2 ymin=2 xmax=$((${xmax}-1)) ymax=$((${ymax}-1)) fi statminsection="${xmin},${xmax},${ymin},${ymax}" } function p_ls { ls -l | grep ^- | ${P_GAWK} '{print $NF}' } # the following extracts keywords; # leading and trailing whitespace is removed (2nd call to gawk) # and blanks inside the keyword value are replaced with underscores (e.g. for OBJECT keyword) function get_key { ${P_DFITS} $1 | ${P_FITSORT} -d "$2" | ${P_GAWK} 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' | ${P_GAWK} '{$1=$1}1' | sed 's/ /_/g' } function get_keyext { ${P_DFITS} -x 1 $1 | ${P_FITSORT} -d "$2" | ${P_GAWK} '{if ($2!="KEY_N/A") print $2}' } function get_region { # merge global and local region files if they exist regionfile="" region="" test -f ${TEMPDIR}/regionfile_${CHIP}.reg && rm ${TEMPDIR}/regionfile_${CHIP}.reg if [ -f ${INSTDIR}/${INSTRUMENT}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}.reg >> \ ${TEMPDIR}/regionfile_${CHIP}.reg regionfile="${TEMPDIR}/regionfile_${CHIP}.reg" region="-d ${regionfile}" fi if [ -f ${INSTDIR}/${INSTRUMENT}_${CHIP}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}_${CHIP}.reg >> \ ${TEMPDIR}/regionfile_${CHIP}.reg regionfile="${TEMPDIR}/regionfile_${CHIP}.reg" region="-d ${regionfile}" fi } function get_region_chip { # merge global and local region files if they exist # for a different chip than the currently processed one # (constant sky subtraction) regionfile="" region="" test -f ${TEMPDIR}/regionfile_$1.reg_$$ && rm ${TEMPDIR}/regionfile_$1.reg_$$ if [ -f ${INSTDIR}/${INSTRUMENT}.reg ]; then cat ${INSTDIR}/${INSTRUMENT}.reg >> \ ${TEMPDIR}/regionfile_$1.reg_$$ regionfile="${TEMPDIR}/regionfile_$1.reg_$$" region="-d ${regionfile}" fi if [ -f ${INSTDIR}/${INSTRUMENT}_$1.reg ]; then cat ${INSTDIR}/${INSTRUMENT}_$1.reg >> \ ${TEMPDIR}/regionfile_$1.reg_$$ regionfile="${TEMPDIR}/regionfile_$1.reg_$$" region="-d ${regionfile}" fi } function prepare_sex_for_region { backstring="" weightstring="" if [ -f ${regionfile} ] && [ "${regionfile}_A" != "_A" ]; then # estimate the background, ignoring blanked pixels background=`${P_FITSSTAT} -i ${FILE} -m median -q` backstring="-BACK_TYPE MANUAL -BACK_VALUE ${background}" weightstring="-WEIGHT_TYPE MAP_WEIGHT -WEIGHT_IMAGE /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits" # check if the dummy weight exists if [ ! -e /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits ]; then mkdir -p /$1/WEIGHTS/ ${P_IC} -p -32 -c ${SIZEX[${CHIP}]} ${SIZEY[${CHIP}]} '1' > ${TEMPDIR}/dummy_${CHIP}.fits ${P_FITSPOLYGON} -i ${TEMPDIR}/dummy_${CHIP}.fits -p ${regionfile} \ -o /$1/WEIGHTS/globalweight_dummy_${CHIP}.fits rm ${TEMPDIR}/dummy_${CHIP}.fits fi fi } function prepare_imred_for_backgroundmodel { # The relevant settings depending on the apply mode of the background model flatstring="" fringestring="" #################################################### # subtract background model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 0 ]; then fringestring="-FLAT N -FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE ${FRINGEIMAGE}" fi #################################################### # divide background model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 1 ]; then flatstring="-FRINGE N -FLAT Y -FLAT_THRESHHOLD 0.1 -FLAT_SCALE ${V_BACK_ADJUSTGAINS}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FLATIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FLATIMAGE=/$1/$2/$2_${CHIP}.fits fi flatstring="${flatstring} -FLAT_IMAGE ${FLATIMAGE} -FLAT_SCALEIMAGE ${FLATIMAGE}" fi #################################################### # subtract fringing model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 2 ]; then fringestring="-FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_FRINGESMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe_smooth.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE /$1/$2/$2_${CHIP}_illum.fits" fi #################################################### # divide background model, subtract fringing model #################################################### if [ "${V_BACK_APPLYMODE}" -eq 3 ]; then flatstring="-FLAT Y -FLAT_THRESHHOLD 0.1 -FLAT_SCALE ${V_BACK_ADJUSTGAINS}" if [ "${V_BACK_ILLUMSMOOTH}"_A != "_A" ]; then FLATIMAGE=/$1/$2/$2_${CHIP}_illum.fits else FLATIMAGE=/$1/$2/$2_${CHIP}.fits fi flatstring="${flatstring} -FLAT_IMAGE ${FLATIMAGE} -FLAT_SCALEIMAGE ${FLATIMAGE}" fringestring="-FRINGE Y -FRINGE_SCALE ${V_BACK_FRINGESCALE}" if [ "${V_BACK_FRINGESMOOTH}"_A != "_A" ]; then FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe_smooth.fits else FRINGEIMAGE=/$1/$2/$2_${CHIP}_fringe.fits fi fringestring="${fringestring} -FRINGE_IMAGE ${FRINGEIMAGE} -FRINGE_REFIMAGE /$1/$2/$2_${CHIP}_illum.fits" fi } # the resultdir is where the output coadded images # will go. If ONE image of the corresponding chip # is a link the image will go to THIS directory for CHIP in $5 do RESULTDIR[${CHIP}]="$1/$4" done + for CHIP in '$5' + RESULTDIR[${CHIP}]=/home/albireo/14-17.02.2015/rgb//science get_statminsection + get_statminsection + local xmin + local ymin + local xmax + local ymax + '[' _A '!=' _A ']' + k=1 + statminsection= + xmax=3358 + ymax=2536 + '[' 1 -le 1 ']' + '[' 3358 -lt 3358 ']' + '[' 2536 -lt 2536 ']' + k=2 + '[' 2 -le 1 ']' + xmin=2 + ymin=2 + xmax=3357 + ymax=2535 + statminsection=2,3357,2,2535 # perform preprocessing (overscan correction, # BIAS subtraction, first flatfield pass) for CHIP in $5 do # check whether the superflat exists. # if yes, then delete it. test -f $1/$4/$4_${CHIP}.fits && rm $1/$4/$4_${CHIP}.fits test -f $1/$4/$4_${CHIP}_illum.fits && rm $1/$4/$4_${CHIP}_illum.fits test -f $1/$4/$4_${CHIP}_fringe.fits && rm $1/$4/$4_${CHIP}_fringe.fits FILES=`ls $1/$4/*_${CHIP}.fits` for FILE in ${FILES} do BAYERPATTERN=`get_key ${FILE} BAYERPAT` if [ -L ${FILE} ]; then LINK=`${P_READLINK} ${FILE}` BASE=`basename ${LINK} .fits` DIR=`dirname ${LINK}` ln -s ${DIR}/${BASE}OFC.fits $1/$4/${BASE}OFC.fits RESULTDIR[${CHIP}]=`dirname ${LINK}` fi done MAXX=$(( ${CUTX[${CHIP}]} + ${SIZEX[${CHIP}]} - 1 )) MAXY=$(( ${CUTY[${CHIP}]} + ${SIZEY[${CHIP}]} - 1 )) do_overscan=Y if [ "${OVSCANX1[${CHIP}]}" = "0" ] && \ [ "${OVSCANX2[${CHIP}]}" = "0" ]; then do_overscan=N fi # flat scaling if [ "${V_DO_FLAT}" = "Y" ]; then FLATSCALEVALUE=`cat /$1/$3/.maxlevel.dat` FLATFLAG="-FLAT_SCALE ${V_DO_FLAT} -FLAT_SCALEVALUE ${FLATSCALEVALUE}" else FLATFLAG="-FLAT_SCALE ${V_DO_FLAT}" fi # check if a dark is subtracted, and if yes, # then if it needs to be rescaled # we subtract the mode of the dark to account for the bias level (on zero order) BIASDARKIMAGE=/$1/$2/$2_${CHIP}.fits if [ "${V_CAL_USEDARK}" = "TRUE" ] && [ "${V_DO_BIAS}" = "Y" ]; then DARKTIME=`get_key /$1/$2/$2_${CHIP}.fits EXPTIME` if [ "${DARKTIME}"_A != "_A" ]; then echo "ok, dark exposure time present; rescaling" testimage=`ls /$1/$4/*_${CHIP}.fits | ${P_GAWK} '(NR==1)'` EXPTIME=`get_key ${testimage} EXPTIME` MODE=`${P_IMSTATS} /$1/$2/$2_${CHIP}.fits | ${P_GAWK} '($0!~/#/) {print $2}'` RESCALE=`echo ${DARKTIME} ${EXPTIME} | \ ${P_GAWK} '{x=1.-$1/$2; if (x<0.01 && x> -0.01) print "NO"; else print $2/$1}'` if [ "${RESCALE}" != "NO" ]; then ${P_IC} '%1 '${MODE}' - '${RESCALE}' *' \ /$1/$2/$2_${CHIP}.fits > \ /$1/$2/$2_${CHIP}rescaled.fits BIASDARKIMAGE=/$1/$2/$2_${CHIP}rescaled.fits fi fi fi # health checks if [ "${V_DO_BIAS}" = "Y" ] && [ ! -f ${BIASDARKIMAGE} ]; then echo `${P_ERRTEST} 0`": Could not find master BIAS for chip ${CHIP}" continue fi if [ "${V_DO_FLAT}" = "Y" ] && [ ! -f /$1/$3/$3_${CHIP}.fits ]; then echo `${P_ERRTEST} 0`": Could not find master FLAT for chip ${CHIP}" continue fi # without nonlinearity correction if [ "${V_NONLINCORR}" = "N" ] || [ ! -f ${INSTDIR}/${INSTRUMENT}".nldat" ]; then # overscan correct, bias subtract and flatfield # science images: # enforce an upper limit to maximages (required by preprocess) if [ ${NFRAMES} -gt 1000 ]; then NFRAMES=1000 fi ${P_IMRED_ECL:?} `ls /$1/$4/*_${CHIP}.fits` \ -c ${CONF}/imred.conf \ -MAXIMAGES ${NFRAMES} \ -OVERSCAN ${do_overscan} \ -OVERSCAN_REGION ${OVSCANX1[${CHIP}]},${OVSCANX2[${CHIP}]} \ -OVERSCAN_REJECT ${V_CAL_OVSCANNLOW},${V_CAL_OVSCANNHIGH} \ -BIAS ${V_DO_BIAS} \ -BIAS_IMAGE ${BIASDARKIMAGE} \ -FLAT ${V_DO_FLAT} \ -FLAT_IMAGE /$1/$3/$3_${CHIP}.fits \ -COMBINE N \ -OUTPUT Y \ -OUTPUT_DIR /$1/$4/ \ -OUTPUT_SUFFIX OFC.fits \ -TRIM Y \ -FLAT_THRESHHOLD 0.1 \ -TRIM_REGION ${CUTX[${CHIP}]},${MAXX},${CUTY[${CHIP}]},${MAXY} \ -STATSSEC ${statminsection} ${FLATFLAG} fi # with nonlinearity correction if [ "${V_NONLINCORR}" = "Y" ] && [ -f ${INSTDIR}/${INSTRUMENT}".nldat" ]; then # overscan correct and bias subtract science images: ${P_IMRED_ECL:?} `ls /$1/$4/*_${CHIP}.fits` \ -c ${CONF}/imred.conf \ -MAXIMAGES ${NFRAMES} \ -OVERSCAN ${do_overscan} \ -OVERSCAN_REGION ${OVSCANX1[${CHIP}]},${OVSCANX2[${CHIP}]} \ -OVERSCAN_REJECT ${V_CAL_OVSCANNLOW},${V_CAL_OVSCANNHIGH} \ -BIAS ${V_DO_BIAS} \ -BIAS_IMAGE ${BIASDARKIMAGE} \ -COMBINE N \ -OUTPUT Y \ -OUTPUT_DIR /$1/$4/ \ -OUTPUT_SUFFIX O.fits \ -TRIM Y \ -TRIM_REGION ${CUTX[${CHIP}]},${MAXX},${CUTY[${CHIP}]},${MAXY} \ -STATSSEC ${statminsection} ${FLATFLAG} # do the nonlinearity correction ls /$1/$4/*_${CHIP}O.fits > ${TEMPDIR}/nonlinlist_$$ cat ${TEMPDIR}/nonlinlist_$$ |\ { while read file do BASE=`basename ${file} _${CHIP}O.fits` ${P_FITSNONLINEARITY} -i ${file} -o ${RESULTDIR[${CHIP}]}/${BASE}_${CHIP}ON.fits \ -c ${CHIP} -f ${INSTDIR}/${INSTRUMENT}".nldat" mv ${RESULTDIR[${CHIP}]}/${BASE}_${CHIP}ON.fits ${RESULTDIR[${CHIP}]}/${BASE}_${CHIP}O.fits done } # flatfield the images ${P_IMRED_ECL:?} `ls /$1/$4/*_${CHIP}O.fits` \ -c ${CONF}/imred.conf \ -MAXIMAGES ${NFRAMES} \ -FLAT ${V_DO_FLAT} \ -FLAT_IMAGE /$1/$3/$3_${CHIP}.fits \ -COMBINE N \ -OUTPUT Y \ -OUTPUT_DIR /$1/$4/ \ -OUTPUT_SUFFIX FC.fits \ -FLAT_THRESHHOLD 0.1 \ -STATSSEC ${statminsection} ${FLATFLAG} rm /$1/$4/*_${CHIP}O.fits fi # do the demosaicing if requested if [ ! "${BAYERPATTERN}_A" = "_A" ]; then BAYERFILES=`ls $1/$4/*_${CHIP}OFC.fits` for FILE in ${BAYERFILES} do BAYERPATTERN=`get_key ${FILE} BAYERPAT` BASE=`basename ${FILE} _${CHIP}OFC.fits` ${P_FITSDEMOSAICBAYER} -i ${FILE} -p ${BAYERPATTERN} -q 2 ${P_REPLACEKEY} $1/$4/${BASE}_${CHIP}OFC.R.fits "FILTER = 'Red'" FILTER ${P_REPLACEKEY} $1/$4/${BASE}_${CHIP}OFC.G.fits "FILTER = 'Green'" FILTER ${P_REPLACEKEY} $1/$4/${BASE}_${CHIP}OFC.B.fits "FILTER = 'Blue'" FILTER mv $1/$4/${BASE}_${CHIP}OFC.R.fits $1/$4/${BASE}_R_${CHIP}OFC.fits mv $1/$4/${BASE}_${CHIP}OFC.G.fits $1/$4/${BASE}_G_${CHIP}OFC.fits mv $1/$4/${BASE}_${CHIP}OFC.B.fits $1/$4/${BASE}_B_${CHIP}OFC.fits rm ${FILE} done fi if [ ! -d /$1/$4/SPLIT_IMAGES ]; then mkdir /$1/$4/SPLIT_IMAGES fi mv /$1/$4/*_${CHIP}.fits /$1/$4/SPLIT_IMAGES done + for CHIP in '$5' + test -f /home/albireo/14-17.02.2015/rgb//science/science_1.fits + test -f /home/albireo/14-17.02.2015/rgb//science/science_1_illum.fits + test -f /home/albireo/14-17.02.2015/rgb//science/science_1_fringe.fits ls $1/$4/*_${CHIP}.fits ++ ls /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-17-56_rot_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-29-10_Green_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-39-48_Red_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-50-06_Green_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-00-24_Blue_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-10-44_Red_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-21-02_Green_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-31-20_Blue_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-41-40_Red_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-51-58_Green_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-02-16_Blue_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-12-35_Red_1.fits + FILES='/home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-17-56_rot_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-29-10_Green_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-39-48_Red_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-50-06_Green_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-00-24_Blue_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-10-44_Red_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-21-02_Green_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-31-20_Blue_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-41-40_Red_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-51-58_Green_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-02-16_Blue_1.fits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-12-35_Red_1.fits' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-17-56_rot_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-17-56_rot_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-17-56_rot_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-29-10_Green_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-29-10_Green_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-29-10_Green_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-39-48_Red_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-39-48_Red_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-39-48_Red_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-50-06_Green_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-50-06_Green_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-50-06_Green_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-00-24_Blue_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-00-24_Blue_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-00-24_Blue_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-10-44_Red_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-10-44_Red_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-10-44_Red_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-21-02_Green_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-21-02_Green_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-21-02_Green_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-31-20_Blue_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-31-20_Blue_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-31-20_Blue_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-41-40_Red_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-41-40_Red_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /usr/bin/gawk '{$1=$1}1' ++ sed 's/ /_/g' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-41-40_Red_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-51-58_Green_1.fits BAYERPAT ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-51-58_Green_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-51-58_Green_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-02-16_Blue_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-02-16_Blue_1.fits + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-02-16_Blue_1.fits ']' + for FILE in '${FILES}' get_key ${FILE} BAYERPAT ++ get_key /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-12-35_Red_1.fits BAYERPAT ++ /home/albireo/THELI/theli//bin/Linux_64//dfits /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-12-35_Red_1.fits ++ /home/albireo/THELI/theli//bin/Linux_64//fitsort -d BAYERPAT ++ sed 's/ /_/g' ++ /usr/bin/gawk '{$1=$1}1' ++ /usr/bin/gawk 'BEGIN{FS="\t"} {if ($2!~/KEY_N\/A/) print $2}' + BAYERPATTERN= + '[' -L /home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-12-35_Red_1.fits ']' + MAXX=3358 + MAXY=2536 + do_overscan=Y + '[' 0 = 0 ']' + '[' 0 = 0 ']' + do_overscan=N + '[' Y = Y ']' cat /$1/$3/.maxlevel.dat ++ cat //home/albireo/14-17.02.2015/rgb//flat/.maxlevel.dat + FLATSCALEVALUE= + FLATFLAG='-FLAT_SCALE Y -FLAT_SCALEVALUE ' + BIASDARKIMAGE=//home/albireo/14-17.02.2015/rgb//BIAS/BIAS_1.fits + '[' FALSE = TRUE ']' + '[' N = Y ']' + '[' Y = Y ']' + '[' '!' -f //home/albireo/14-17.02.2015/rgb//flat/flat_1.fits ']' + '[' N = N ']' + '[' 47 -gt 1000 ']' ls /$1/$4/*_${CHIP}.fits ++ ls //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-17-56_rot_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-29-10_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-39-48_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-50-06_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-00-24_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-10-44_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-21-02_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-31-20_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-41-40_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-51-58_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-02-16_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-12-35_Red_1.fits + /home/albireo/THELI/theli//bin/Linux_64//preprocess //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-17-56_rot_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-29-10_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-39-48_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-50-06_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-00-24_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-10-44_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-21-02_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-31-20_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-41-40_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-51-58_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-02-16_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-12-35_Red_1.fits -c /home/albireo/THELI/theli//gui/reduction/imred.conf -MAXIMAGES 47 -OVERSCAN N -OVERSCAN_REGION 0,0 -OVERSCAN_REJECT 0,1 -BIAS N -BIAS_IMAGE //home/albireo/14-17.02.2015/rgb//BIAS/BIAS_1.fits -FLAT Y -FLAT_IMAGE //home/albireo/14-17.02.2015/rgb//flat/flat_1.fits -COMBINE N -OUTPUT Y -OUTPUT_DIR //home/albireo/14-17.02.2015/rgb//science/ -OUTPUT_SUFFIX OFC.fits -TRIM Y -FLAT_THRESHHOLD 0.1 -TRIM_REGION 1,3358,1,2536 -STATSSEC 2,3357,2,2535 -FLAT_SCALE Y -FLAT_SCALEVALUE \_/ **** fatal error in pid 17302 '-0-' **** segmentation fault or bus error --0-- .-0-. + '[' N = Y ']' + '[' '!' _A = _A ']' + '[' '!' -d //home/albireo/14-17.02.2015/rgb//science/SPLIT_IMAGES ']' + mkdir //home/albireo/14-17.02.2015/rgb//science/SPLIT_IMAGES + mv //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-17-56_rot_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-29-10_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-39-48_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_19-50-06_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-00-24_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-10-44_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-21-02_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-31-20_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-41-40_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_20-51-58_Green_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-02-16_Blue_1.fits //home/albireo/14-17.02.2015/rgb//science/pferdekopf_2015-02-16_21-12-35_Red_1.fits //home/albireo/14-17.02.2015/rgb//science/SPLIT_IMAGES