|
HMI Ring Diagrams
Pipeline Module Specifications
|
rdfitc
Module Version
0.8 25 VIII 2009
General Description
rdfitc implements the procedure of Basu and Antia (ApJ 512, 458;
Solar Phys 192, 469) for fitting 3-d solar acoustic power spectra to
a multi-parameter model for line and background characteristics
Status
This version has been superseded by Version 0.9
Never used for ring-diagram pipelines
Not included in any JSOC release
The source code and Makefile is in ~rick/src/ringfit/rdfitc
This version supersedes Version 0.7
Usage
rdfitc_v08 [-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)
guessfile Name of a file containing the initial frequency guesses;
should be a rooted pathname.
Default: /home/rick/rings/fit/sb_ringfit/freq.guess
out Name of the output file; should be a rooted pathname
Default: fit.out
nmin Minimum radial (vertical) order to fit. Default: 0
nmax Maximum radial order to fit. Default: 12
lmin Minimum equivalent spherical-harmonic degree to fit. Default: 100
lmax Maximum degree to fit. Default: 1400
fmin Minimum frequency to fit [μHz]. Default: 900
fmax Maximum frequency to fit [μHz]. Default: 7000
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
Input Data series class:
power spectra of data cubes
Sample input series:
mdi.rdVpspec_dp
Output Data series class:
ring diagram mode fits
Sample output series:
mdi.rdVfitsc_dp
Testing
Changes from Previous Version
— fixed scaling of ux_guess and uy_guess and test on par[2]
— trimmed output format slightly
— fixed bug in setting of Dnu from Domega
Bugs in current version
Certain keywords expected to be present in the input are not (MAPSCALE
and MAP_PROJ); they are ignored, with default values set. WCS keywords
should be used instead.
The output is only to a named file, not a DRMS record segment.
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 Mm-1 and sec-1; only the units
are verified. The code will not work properly if the number of rows or
columns is odd.