HMI_Logo HMI Ring Diagrams
Pipeline Module Specifications

rdfitc

Module Version
	1.3	18 IV 2017

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
	In use for ring-diagram pipelines from 2017.03.04 – (CR 2170–2171, 2187 –)
	Included in JSOC release 9.1
	The source code and Makefile are in the JSOC CVS repository
	This version supersedes Version 1.2

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 spherical-harmonic 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 kx in input spectrum. Default: no minimum
	kxmax		Maximum accepted value of kx in input spectrum. Default: no maximum
	kymin		Minimum accepted value of ky in input spectrum. Default: no minimum
	kymax		Maximum accepted value of ky 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 comma-separated list of keywords to be propagated as-is 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
	— Use template record for output series checking
	— Log host/timing info even when not in drms


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
	power-law 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 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.


Valid HTML 4.01 Transitional HMI Ring Diagrams 7 May 2019, 11:22-0700