HMI Ring Diagrams|
Pipeline Module Specifications
1.0 19 VIII 2010
pspec3 is designed to calculate the full 3-dimensional power spectra of
a set of time series of 2-dimensional images. The input data are optionally
apodized by a (1-x^2)^2 taper operating between specified inner and
outer radii of the images and a similar taper operating on the time axis
(image number). The input data series must consist of 3-dimensional
data cubes (such as those produced by mtrack), each treated as a
separate data record.
Included in JSOC release 5.10
The source code and Makefile are in the JSOC CVS repository
pspec3_v10 is a compiled standalone DRMS binary built with NetDRMS 2.4
to run in the JSOC configuration. Versions are in:
The source code and Makefile is in ~rick/src/pspec
pspec3 [-lvx] in= InputDescriptor pspec= OutputDescriptor
[param= val ...]
-l output direct power spectrum rather than scaled log
-n do not save output record (diagnostics only)
-v run verbose
-x use double-precision calculation internally
in The name of an input dataset (record selection); No default.
segment The name of the segment in the input data set to be used;
ignored if there is only one 3-dimensional segment in the
data set. Default: not specified
pspec The name of the output data series, which must already
exist; record specification is based on the input data record(s)
and the series prime keywords. No default.
mask_in Inner radial edge of spatial taper; Default: 0.9375
mask_ex Outer radial edge of spatial taper; Default: 1.0
apodize Inner edge of temporal taper; Default: 0.96875
fbin Output frequency binning; 0 -> no binning; Default: 0
copy A comma-separated list of keywords to be propagated as-is from
input to output; Default: "+"
Input Data series class:
tracked mapped data cubes from solar image data
Sample input series:
Input keys used or inspected:
CDELT3, CUNIT3, MapScale, T_STOP, T_START
Either CDELT3 or both T_STOP and T_START must be present. If CDELT3 is
present it is used; if CUNIT3 is missing and needed to interpret CDELT3,
the value of CUNIT3 is assumed to be "s"
MapScale is required for the WCS keywords in the output spectrum to be
set; its values are expected to be in deg/pixel.
Output Data series class:
power spectra of data cubes
Output keys set (if possible):
CRPIX1, CRPIX2, CRPIX3, CRVAL1, CRVAL2, CRVAL3, CTYPE1, CTYPE2, CTYPE3,
CDELT1, CDELT2, Delta_k, Delta_kx, Delta_ky, CDELT3, Delta_nu, D_OMEGA,
Cadence, Apode_k_min, APOD_MIN, Apode_k_max, APOD_MAX, apode_f, LOG_BASE,
DataMin, DataMax, Module, BLD_VERS, Source, Input, Created
Sample output series:
Input keys propagated to output (if possible) if "+" in copy list:
CarrTime, CarrRot, CMLon, LonHG, LatHG, LonCM, MidTime, Duration, LonSpan,
T_START, T_STOP, Coverage, ZonalTrk, ZonalVel, MeridTrk, MeridVel,
MapScale, Size, Width, Height, MapProj, Map_PA, PosAng, RSunRef, MAI,
Significant changes from Previous Version
— Added Ident to list of keys propagated by default; added option for
providing an alternate or additional list of keywords for propagation
Bugs in current version
For large input data cubes (~1 Gpxl) the results may be garbage,
for larger ones (~2 Gpxl) the input may fail altogether. The exact
causes of these problems and consequently the exact limits of
validity are not known, but the code has been verified to work with
"normal" input cubes (Dopplergrams as compressed scaled shorts)
of up to 800 Mpxl.
There is no checking that the input data series corresponds with a link
in the output series (if present)
The value of CDELT3 in the input set overrides the value of Cadence if
propagated, but is not checked against the units field of the Cadence
keyword, hence it will be incorrect if the Cadence units field differs
from the value of CUNIT3.
WCS keywords should be an acceptable (even preferred) alternative to
MapScale as input keywords.