Skip to content
Snippets Groups Projects
Commit 8f3ee773 authored by Andreas Korpi-Lagg's avatar Andreas Korpi-Lagg
Browse files

added fe6173_synth.ipt

parent ed5e0ce5
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 193 additions and 0 deletions
input/fe617_synth.ipt 0 → 100644
+ 193
− 0
View file @ 8f3ee773
;--------- COMMENTS ---------------------------------------------------
;--------- DIRECTORIES ------------------------------------------------
PS ./ps/ ;directory for postscript output
SAV ./sav/ ;directory for storing results (sav-file)
PROFILE_ARCHIVE ./profile_archive/ ;input directory for profiles / observations
ATM_ARCHIVE ./atm_archive/ ;output directory for results (atmospheres)
ATM_SUFFIX mo-0 ;add a suffix to the atm-directory to identify this
; run
WGT ./wgt/ ;directory for wgt-files
ATOM ./atom/ ;directory for atomic data files
;--------- CONTROL ----------------------------------------------------
DISPLAY_MAP 1 ;if 1 then display maps of parameters after
; successful run (multiple pixels only!)
DISPLAY_PROFILE 1 ;if 1 then display fitted & observed profiles
DISPLAY_COMP 0 ;if 1 then display individual atm. components of
; fitted profiles
VERBOSE 2 ;verbosity of output: 0=quiet, 1=normal, 2=women
SAVE_FITPROF 1 ;if 1 then save fitted profile as sav-file (can be
; used as input profile)
FITSOUT 1 ;if 1 then write out map results in FITS format
OUTPUT X ;set to 'PS' for postscript output
OLD_NORM 0 ;use code with old normalization of the various
; atmospheric components
OLD_VOIGT 0 ;if 1 use old normalization of Faraday/Voigt function
;--------- DATA SET ---------------------------------------------------
OBSERVATION ;observation in one of the formats:
; 1.) sav file of the format
; 'obsname'.profiles.sav, 2. spinor-profile, 3.
; sav-file created with helix 'SAVE_FITPROF'
; keyword, 4. reduced TIP data file (extension
; '.cc') and auxliary file ('.ccx')
;WL_RANGE 6169.600000 6177.000000 ;WL-range to be used for analysis (may be
WL_RANGE 6172.600000 6174.000000 ;WL-range to be used for analysis (may be
; only a part of the observation WL-range)
WL_NUM 400 ;# of WL-points (for synthesis only)
WL_DISP 0.000000 ;WL-calibration: dispersion per WL-bin number
WL_OFF 0.000000 ;WL-calibration: offset (used if != 0)
WL_BIN 1 ;wavelength binning
XPOS 0 0 ;two-elements vector containing xmin,xmax of the
; observation map to be analyzed
YPOS 0 0 ;two-elements vector containing ymin,ymax of the
; observation map to be analyzed
PROFILE_LIST ;filename of profile-list (2 columns: xxx yyy). If
; present, XPOS and YPOS are ignored. File must be
; located in ./ or PROFILE_ARCHIVE.
STEPX 5 ;step size for going from xmin to xmax
STEPY 5 ;step size for going from ymin to ymax
AVERAGE 1 ;if 1 then average observation over the stepx/stepy
; size
SCANSIZE 0 ;stepsize of multiple scans within one observation
SYNTH 1 ;if 1 then create synthetic profile
NOISE 0.00000 ;noise level for adding random noise (IQUV)
INOISE 0.0 ;noise level for adding random noise to Stokes I only
SMOOTH 0 0 ;smooth-value for profiles and smooth-method:
; (0=IDL-smooth function,1=FFT Low-Pass)
MEDIAN 0 ;median filter for observed profiles (0=off, >=2
; median filter width)
SPIKE 0 ;if >=1 then remove electronic spikes. Higher values
; remove broader spikes.
MIN_QUV 0.0 ;min. mag. signal. If sqrt(Q^2+U^2+V^2)<MIN_QUV then
; the fit will be done with B=0.
STRAYPOL_AMP 0.00000 ;amplitude for stray-polarization (only used for
; synthesis)
STRAYPOL_CORR 0 B ;iteration steps and orientation of
; scattering-polarization correction. Orientation:
; 'B' = along B-field, 'X' = a number defining an
; angle manually
HANLE_AZI 0. ;Offset for Hanle calc. (experimental)
SLIT_ORIENTATION 0.00 ;slit-orientation (used for straypol-corr)
M1ANGLE 0 ;angle of coelostat mirror M1
SOLAR_POS 0.0 0.0 ;pos. of observation (x and y in arcseconds)
SOLAR_RADIUS 950. ;radius of sun in arcsec for time of obs.
HELIO_ANGLE 0. ;heliocentric angle of observation in degrees (for
; full Hanle). For TIP: theta=acos(sqrt(x^2+y^2)/r,
; (x,y) = position of observation in terrestrial
; coordinates.
+Q2LIMB 0. ;pos. Q to limb direction in degrees (for full
; Hanle). For TIP: Q2L=atan(y,x), emission vector
; angle gamma=360.-Q2L+90 .
HIN_SCANNR 1 ;Hinode Scan# (repetetive scans)
NORM_STOKES IC ;normalization of Stokes profiles: I or IC
NORM_CONT LOCAL ;type of cont. normalization: local, slit or image
IC_LEVEL 0. ;set level for Ic. A positive number overwrites the
; Ic value contained in the observational data.
LOCALSTRAY_RAD 0. ;add a local straylight component: take the
; average Stokes profile of the surrounding
; profiles within the radius specified here and mix
; it via alpha to the other components. The radius
; is in units of the pixel size defined with the
; STEPX / STEPY keywords.
LOCALSTRAY_FWHM 0. ;FWHM of the Gaussian weighting used for the
; localstraylight profiles
LOCALSTRAY_CORE 0. ;exclude core profiles within the radius
; specified here
LOCALSTRAY_POL 1 ;1=polarized / 0=unpolarized local straylight
;--------- POST PROCESSING --------------------------------------------
CONV_FUNC ./convolve/fp_617.3.hw100.dat ;convolution function of instrument
CONV_NWL 256 ;# of bins for convolution (used if number of
; WL_bins in data is small
CONV_MODE FFT ;convolution method: FFT or MUL
CONV_WLJITTER 0.0 ;uncertainty in \AA in knowledge of exact
; wavelength position for every filter position
; (e.g. caused by error in etalon voltage). Make
; sure that the wavelength resolution is high
; enough (i.e. increase CONV\_NWL)! Must be set to
; ZERO for real data!
CONV_OUTPUT 0 ;flag to control the display/output of the fitted
; profile. If 1 then the profile is displayed with
; CONV_NWL WL pixels, if 0 the WL-pixels of the
; observation are used. NOTE: for local straylight
; correction this flag is set to 0.
PREFILTER ;prefilter curve
PREFILTER_WLERR 0.0 ;error in \AA in knowledge of prefilter curve
; (e.g. caused by unknown temperature). Must be set
; to ZERO for real data!
;--------- ATMOSPHERES ------------------------------------------------
NCOMP 1 ;number of components
; --- atmospheric component 1 ---
; NAME Value SCL_MIN SCL_MAX FIT %RG FIT
BFIEL 2000.000 0.000 3000.000 4 ;magnetic field value in Gauss
AZIMU 00.000 -90.000 90.000 3 ;azimut of B-vector [deg]
GAMMA 60.000 0.000 180.000 2 ;inclination of B-vector [deg]
VELOS 0.000 -5000.000 5000.000 1 ;line-of-sight velocity in m/s
VDAMP 0.050 0.000 0.990 1 ;damping constant (Voigt +
; hanle-slab)
VDOPP 0.050 0.000 0.990 1 ;doppler broadening (Voigt +
; hanle-slab)
EZERO 5.000 0.000 110.000 0 ;amplitude of components of
; propagation matrix (Voigt profile only!)
SGRAD 1.000 0.000 9.000 6 ;gradient of source function
ALPHA 1.000 0.010 0.990 0 ;Filling factor for this component
USE_ATOM fe617.3.dat ;atomic data file(s) for this component
;--------- LINE PARAMETERS --------------------------------------------
; --- fit parameters for line 1 ---
; NAME Value SCL_MIN SCL_MAX FIT
LINE_ID 6173.3412 ;wavelength to identify the line
LINE_STRENGTH 1.00 0.00 0.00 0 ;factor for line strength
LINE_WLSHIFT 0.00 0.00 0.00 0 ;adjust central WL of line
;--------- BLENDS -----------------------------------------------------
NBLEND 1 ;number of telluric blends
; --- telluric blend 1 ---
; NAME Value SCL_MIN SCL_MAX FIT
BLEND_WL 0.000 0.000 0.000 0 ;WL and WL-range for blend
BLEND_WIDTH 0.000 0.000 0.000 0 ;width of voigt-profile
BLEND_DAMP 0.000 0.000 0.000 0 ;damping of voigt-profile
BLEND_A0 0.000 0.000 0.000 0 ;amplitude of voigt-profile
; (blend is not used if BLEND_A0=0
;--------- GENERAL FIT PARAMETERS -------------------------------------
; NAME Value SCL_MIN SCL_MAX FIT
CCORR 1.00 1.00 1.00 0 ;factor for cont. correction
STRAYLIGHT 0.00 0.00 1.00 0 ;spectral straylight
; correction (non-dispersive, non-polarized)
RADCORRSOLAR 1.00 0.10 10.00 0 ;Correction for the intensity
; of the solar radiation field (for Hanle mode only)
;--------- ANALYSIS METHOD --------------------------------------------
IQUV_WEIGHT 1.00000 1.00000 1.00000 1.00000 ;4-element vector defining
; relative weighting of IQUV (in that order).
; Additionally the code does an automatic weighting
; according to the strength of the I signal
; compared to the QUV signals. This weighting
; scheme is used in the PIKAIA fit routine.
WGT_FILE fe_6173.wgt ;file with WL-dependent weighting function for IQUV
PROFILE voigt ;functional form for pi- and sigma components of
; spectral line. Available: gauss, voigt,
; voigt_phys, voigt_szero
MAGOPT 0 ;include magneto-optical effects (dispersion
; coefficients, (Voigt profile only!))
USE_GEFF 0 ;use effective Lande factor (=1) or real Zeeman
; pattern (=0)
USE_PB 0 ;if set, the zeeman-splitting and strength includes
; the Paschen-Back effect (from the table by
; Socas-Navarro)
PB_METHOD poly ;use polynomials (=poly) or table interpolations
; (=table) to calculate the PB-effect
;--------- PIKAIA PARAMATERS ------------------------------------------
CODE FORTRAN ;PIKAIA code to use. Available: FORTRAN (=fast) or
; IDL (=platform independent).
METHOD PIKAIA ;minimization method: PIKAIA, POWELL (fast), LMDIF
; (fast) or PIK_LM (combined, for maps)
NCALLS 0 ;number of iterations in PIKAIA routine / max.
; number of calls for POWELL or LMDIF
CHI2MODE 0 ;method to calc. chi2 (JM = Borrero method, PLAIN =
; weight only depends on IQUV_WEIGHT values).
; Default: Weight depends on 1/signal strength
PIXELREP 1 ;number of repetitions per pixel (to perform
; statistical analysis)
KEEPBEST 1 ;store result only if fitness is better than
; existing result (FITSOUT only)
PIKAIA_STAT NONE ;switch on and control statistic module based on
; PIKAIA populations. Requires FITSOUT=1 and
; PIXELREP=1. Example: POP75
PIKAIA_POP 0 ;number of PIKAIA populations (0=automatic)
;--------- END OF INPUT FILE ------------------------------------------
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment