;;;;;;;;;;;;;;
 ;;;
 ;;;  This file is read by the routine INIT_VARIABLES.
 ;;;  It defines all the variables necessary for the speckle 
 ;;;  routines to work correctly.  It is the only file you have to
 ;;;  change for different setups!
 ;;;  The syntax has to be valid IDL, as the file gets executed.  So
 ;;;  remember placing comments only after a ';'
 ;;;
 ;;;;;;;;;;;;;;

RawDir    = './raw/'                     ;;; Location of raw data
HelpDir   = './support/'                 ;;; Dir for additional files
ResDir	  = './results/'                 ;;; Dir for corrected results
f_raw	  = 'rawflat_g'			 ;;; filename raw flats
z_raw	  = 'rawdark_g'			 ;;;    "      "  darks
d_raw	  = 'rawdata'			 ;;;    "      "  data
FileOffset= 0			 ;;; Byte offset for files with
DarkOffset= 0                            ;;; headers, 0 for raw, 2880
FlatOffset= 0			 ;;; for FITS
FileEndian= 1                            ;;; Format of raw data,
                                         ;;; 1=PC, 0=Sun,HP
Bitpix    = 16				 ;;; byte(8) or integer(16) raw data
Rot_Img	  = 2				 ;;; Rotation of image for sim. series
FitsAdd   = ''                           ;;; File with aditional FITS entries
Mask_Depth= 9				 ;;; 
I_Limit   = 1023			 ;;; Max. Intensity level
N_Image   =  100			 ;;; # of images per burst
N_Flat	  =  100			 ;;; # of flats in burst
N_Dark	  =  100			 ;;; # of darks in burst
X_Size_0  = 1296                         ;;; Raw data X axis
Y_Size_0  = 1028		         ;;;  "    "  Y axis
X_Size_S  =    0                          ;;; used data X axis start
X_Size_E  = 1295                         ;;; used data X axis end
Y_Size_S  =    0		         ;;; used   "  Y axis start
Y_Size_E  = 1027		         ;;; used   "  Y axis end
X_Width	  = X_Size_E - X_Size_S + 1
Y_Width	  = Y_Size_E - Y_Size_S + 1
S_Size	  = 64 				 ;;; Size of subimages
  ;;;  To compute subimage numbers automatically, for an overlap of
  ;;;  approx 33% use the following formula:
;S_nx      = fix(X_Width * 1.5 / S_Size - 0.5)
;S_ny      = fix(Y_Width * 1.5 / S_Size - 0.5)
  ;;; More general: Overlap of xx (here: xx=0.4)
;S_nx      = fix((X_Width-S_Size)/(S_Size*(1.-0.4))+1.5)
;S_ny      = fix((Y_Width-S_Size)/(S_Size*(1.-0.4))+1.5)
S_nx	  = 35				 ;;; # of sub-images in X
S_ny	  = 28				 ;;; # of sub-images in Y
S_n	  = S_nx * S_ny			 ;;; total # of sub-images
A_Wid	  = 7				 ;;; rimwidth of apodisation window
Tel_dia   = 43.8                         ;;; Telescope aperture in cm
Lambda    = 4305.                        ;;; Observing wavelength in A
Arcs_pp   = 0.071                        ;;; Arcsec per pixel
Quiet	  = 1				 ;;; Be verbose
ScaleNoise= 0				 ;;; Scale noise power to data level
UseOptFilt= 1                            ;;; Use optimum filter?
UseXDR    = 0
TmpEndian = 1

 ;;; The following variables normally don't need to be changed

StfDir	  = '~/idl/stf/'		 ;;; The precomputed STF's
STF_2D    = 'DOT/'			 ;;; For which aperture?
flatname  = '../flat+noise/flat_g'
darkname  = '../flat+noise/dark_g'
noisename = '../flat+noise/rauschpower_g'
dimname   = 'basisdimensionen'
longname  = 'langzeitbild'
refname   = 'referenzbild'
rmsname   = 'kontrast'
scalename = 'flatnorm'
partsclnam= 'partnorm'
flatokname= 'flatbilder'
splitname = 'teilbildserie'
resname   = 'teilbild'
res1name  = 'allbild'
regname	  = 'flaeche'
motionname= 'imagemotion'
restdim	  = 'ruhdimensionen'
trackname = 'ruhbilder'
stfname	  = 'stf'
srname	  = 'ratio'
sqname	  = 'quotient'
filtername= 'rauschfilter'
powername = 'powerspectrum'
stepname  = 'schrittweite'
rawsumname= 'rawsum'