
HMI Ring Diagrams
Pipeline Module Specifications

rdfitc
Module Version
1.2 9 VI 2011
General Description
rdfitc implements the procedure of Basu and Antia (ApJ 512, 458;
Solar Phys 192, 469) for fitting 3d solar acoustic power spectra to
a multiparameter model for line and background characteristics
Status
This version has been superseded by Version 1.3
In use for ringdiagram pipelines 2011.06.20 – 2017.02.20
(CR 2096–2097, 2109–2186)
Included in JSOC release 6.0
The source code and Makefile are in the JSOC CVS repository
This version supersedes Version 1.1
rdfitc_v12 is a compiled standalone DRMS binary built with NetDRMS 9.0
to run in the JSOC configuration. Versions are in:
~rick/bin/linux_avx/rdfitc_v12
~rick/bin/rick/rdfitc_v12
Usage
rdfitc [nvx2] in= InputDescriptor out= OutputDescriptor
[param= val ...]
Flags
n no fitting (diagnostics only)
v run verbose
x extra reporting of number of iterations
2 calculate second derivatives explicitly
Parameters
in the name of an input data set (record selection); No default
guessfile Name of a file containing the initial frequency guesses;
should be a rooted pathname. No default:
out Name of the output data series (or file)
Default: fit.out
The value of the out parameter is interpreted as the name of
of a data series, and if such a series exists with insert
permission, the output fit files and logs will be inserted as
records in that series; otherwise, the value will be interpreted
as the name of a file to which the fit parameters are written,
and the log information goes to stderr.
nmin Minimum radial (vertical) order to fit. Default: 0
nmax Maximum radial order to fit. Default: 10
lmin Minimum equivalent sphericalharmonic degree to fit. Default: 80
lmax Maximum degree to fit. Default: 3000
fmin Minimum frequency to fit [μHz]. Default: 900
fmax Maximum frequency to fit [μHz]. Default: 5500
kxmin Minimum accepted value of k_{x} in input spectrum. Default: no minimum
kxmax Maximum accepted value of k_{x} in input spectrum. Default: no maximum
kymin Minimum accepted value of k_{y} in input spectrum. Default: no minimum
kymax Maximum accepted value of k_{y} in input spectrum. Default: no maximum
bfgsct Number of iterations for bfgs function.. Default: 125
linminct Number of iterations for linmim function.. Default: 15
ux Initial guess for zonal flow speed [m/s]. Default: 0.0
uy Initial guess for meridional flow speed [m/s]. Default: 0.0
A1_guess Initial guess for amplitude term A1. Default: 0.0
A1 Fixed value for amplitude term A1 in lieu of fit. Default: "fit"
S_guess Initial guess for asymmetry term S. Default: 200.0
S Fixed value for asymmetry term S in lieu of fit. Default: "fit"
copy A commaseparated list of keywords to be propagated asis from
input to output; Default: "+"
Input Data series class:
power spectra of data cubes
Input keys used or inspected:
WCSNAME, CTYPE1, CTYPE2, CTYPE3, CUNIT1, CUNIT2, CUNIT3, CDELT1, CDELT2, CDELT3
DELTA_K*, RSunRef, DELTA_NU*, D_OMEGA*, LOG_BASE
* only if preferred keys not found
Sample input series:
hmi.rdVpspec_fd05
hmi.rdVpspec_fd15
hmi.rdVpspec_fd30
mdi.rdVpspec_dp
Output Data series class:
ring diagram mode fits
Output keys set (if possible):
CarrTime*, Module, Source, Created, BLD_VERS, n_x_fits, n_x_fmin, n_x_fmax
* if not propagated
Input keys propagated to output (if possible) if "+" in copy list:
CarrTime, CarrRot, CMLon, LonHG, LatHG, LonCM, MAI, Ident, MidTime,
Duration, Cadence, LonSpan, T_START, T_STOP, Coverage, Size, Width, Height,
MapProj, MapScale, Map_PA, PosAng, ZonalTrk, ZonalVel, MeridTrk, MeridVel,
Apode_f, Apode_k_min, Apode_k_max, APODIZNG, APOD_MIN, APOD_MAX
Sample output series:
hmi.rdVfitsc_fd05
hmi.rdVfitsc_fd15
hmi.rdVfitsc_fd30
mdi.rdVfitsc_dp
Changes from Previous Version
changed contents of column 5 from w0 to D_w0
set values in qual column based on autoweeding used in rdvinv
fixed bug in determination of default input segment when input series
contains multiple segments
added optional rectangular filtering on input spectrum
Bugs in current version
The output column for δν (5) is labeled D_w_0 (and identical with
column 20), since there is no formal uncertainty on the frequency.
The quality parameters sn (amplitude/background) and pk (the fit
powerlaw exponent for ν(l) are not printed (simply because they
were not included in the standard).
If either of the flags 2x are given the output will not conform to
the standard; the flag 2 is not correctly interpreted in any case.
It is implicitly assumed (not checked) that the input power spectrum
is centred, with CRPIX1 = cols / 2 + 1, CRPIX2 = rows / 2 + 1, CRPIX3 = 1,
CRVALi = 0.0, and that the units are in Mm1 and sec1; only the units
are verified. The code will not work properly if the number of rows or
columns is odd.