export FITS_KEY_TEMPERATURE="CCD-TEMP"

mkdir -p 20170227/
cd 20170227/
rm -f TYC*.fits
rm -f TYC*.fits~

wget ftp://anonymous@integral.physics.muni.cz/pub/praktikum/20170227/*
munipack dark -o d45.fits --verbose M67_*_Zadny_022?.fits

rename 's/ /_/' TYC*.fits

# Geograficke souradnice Vyskova 49.2836214N, 17.0224619E
for A in TYC_*.fits; do
    munipack fits --update-key LATITUDE=49.2836214 $A
    munipack fits --update-key LONGITUD=-17.0224619 $A
done

munipack phcorr -dark d45.fits TYC_*.fits

munipack find -f 7 TYC_*.fits
munipack aphot TYC_*.fits

munipack cone -r 0.2 -- 188.266288 +44.915328
munipack astrometry --sig 5 --disable-flux-check TYC_*.fits

python ucac_jmuc.py cone.fits cone_John.fits

REF=TYC_3020-2221-1_2017-02-27_22-24-47_R_0149.fits
echo REF: $REF
munipack phcal --photsys-ref Johnson --area 0.2 \
	 -f R --col-ra RAJ2000 --col-dec DEJ2000 --col-mag Rmag --col-magerr e_Rmag \
	 -c cone_John.fits -O --mask '\1_catcal.\2' $REF

for A in TYC_3020-2221-1_*_R_????.fits; do
    munipack phcal --photsys-ref Johnson --area 0.2 $SAPER -f R \
	     -c ${REF%.fits}_catcal.fits\[PHOTOMETRY\] \
       --col-ra RA --col-dec DEC --col-mag MAG --col-magerr MAGERR \
       -O --mask '\1_refcal.\2' $A;
done

munipack timeseries --stdout 188.266288,44.915328 188.281,44.937  -c MAG,MAGERR TYC_*_refcal.fits > lc.dat