|
HMI Ring Diagrams
Pipeline Module Specifications
|
rdvinv
Module Version
0.9 24 IV 2012
General Description
rdvinv performs an OLA inversion on a set of ring-fit parameters Ux and Uy
for transverse velocity functions of depth.
If output is to a named DRMS series, the module combines the output files
corresponding to the input records selected into a single SUMS directory
or a pair of directories, depending whether the output data series has
one or two segments defined. If the output series has a single segment,
all output files will be placed in its SUMS directory; otherwise, separate
subdirectories for the uxseg and uyseg parameters will be created.
In either case, the ouput inversion files will be named with extensions
.Ux and .Uy. The filenames are constructed from the Carrington time,
latitude and (Stonyhurst) longitude information for the records, e.g.
2099:060.0_15.0E45.0S.Ux
If values are provided for any combination of parameters cr, clon, lon, and
lat, the module will select all records from the input record set or series
matching the selected values.
Status
This version has been superseded by Version 0.9
Included in JSOC release 6.2
In use for ring-diagram pipelines 2012.05.12 – 2012.08.27 (CR 2096, 2122–2126)
This version supersedes Version 0.8
The source code and Makefile are in the JSOC CVS repository
Usage
rdvinv [-v] [param= val ...]
Flags
-v run verbose
Parameters
in Name of an input data series or dataset (record selection);
No Default
seg Name of the input segment; Default: fit.out
out Output filename or directory; Default: fort.10.hmixy
uxseg Output data segment name for Ux inversions; Default: Ux
uyseg Output data segment name for Uy inversions; Default: Uy
cr Carrington rotation number common to all regions;
Default: Not Specified
clon Central meridian longitude common to all regions;
Default: Not Specified
lon Carrington longitude common to all regions;
Default: Not Specified
lat Heliographic latitude common to all regions;
Default: Not Specified
amu Error trade-off parameter; Default: 0.005
ob Lower frequency limit (mHz); Default: 1.0
oe Upper frequency limit (mHz); Default: 5.2
rb Lower radius limit; Default: 0.97
re Upper radius limit; Default: 1.00
num Number of target inversion points; Default: 40
kernel Kernel dataset record or filename; No Default
ave Output filename for averaging kernels; if unspecified,
kernels are not written. Default: Not Specified
coef Output filename for inversion coefficients; if unspecified,
coefficients are not written. Default: Not Specified
Output Data series class:
ring diagram flow inversions
Sample input series:
hmi.rdVfitsf_fd15
hmi.rdVfitsc_fd30
mdi.rdVfitsc_dp
Output Data series class:
ring diagram mode fits
Sample output series:
hmi.rdVflows_fd15_frame
hmi.rdVflows_fd30_frame
Changes from Previous Version
— Added (preferred) option for specifying kernel as DRMS data record
rather than filename
Bugs in current version:
Certain arrays in the FORTRAN are not reinitialized between calls,
so the results depend to a (small) extent on the order of processing
when multiple input records are processed.
The mode weeding for abnormal values of Ui is only applied to modes of order
n <= 8; all higher-order modes are rejected. The rejection criteria in the mode
weeding algorithm are hard-coded.
If any of the parameters cr, clon, lon, and lat are specified, the input data
must be in a DRMS series with the appropriate keywords CarrRot, CMLon, LonHG,
and LatHG.
It is not possible to restrict the DRMS input record set based on Stonyhurst
longitude via a calling parameter, only via the recordset descriptor.
When writing output to a DRMS series, there is no verification that the
requested output segments have the correct protocol (generic); but it probably
doesn't matter in any case.
When writing output to a file, rather than a drms record set, only one
set of files can be specified — if the input data set contains more than
one record, each inversion will overwrite the last. Likewise, specification
of the parameters ave and coef for the optional output of the averaging kernels
and the inversion coefficients is only applicable when there is a single input
record; otherwise, the result for each record will overwrite the last.
The current version includes code from ola_xy.c that generates compiler
warnings, but this code is not actually used.