|
HMI Ring Diagrams
Pipeline Script Specifications
|
avgpsbycr
Script Version
1.6 11 VIII 2015
General Description
Produces time averaged power spectra over a Carrington rotation
from the power spectra of tracked ring-diagram tiles
Status
This version has been superseded by Version 1.7
Used for pipeline production from Aug 11, 2015 through Sep 25, 2017, including all of
hmi.rdVavgpspec_fd05[2161-2177, 2192-2193], hmi.rdVavgpspec_fd15[2168-2194], and hmi.rdVavgpspec_fd30[2168-2194]
Included in JSOC release 8.11
The script is in the JSOC CVS repository
This version supersedes Version 1.5
Usage
avgpsbycr [-test] Size CarRot
Arguments
-test If specified, the output series will be in namespace hmi_test;
otherwise it is in hmi
Size The tile size of spectra to be averaged: must be 5, 15, or 30
CarRot The Carrington rotation number of the averaging interval
Input Data Series:
hmi.rdVpspec_fd*
Output Data Series:
hmi.rdVavgpspec_fd*
hmi_test.rdVavgpspec_fd*
Modules invoked:
datavg
Auxiliary Files used:
~rick/hmi/notify (for addresses of email notifications)
~rick/hmi/rings/testoptions (for test options if -test specified)
~rick/hmi/tmpfs (optional locations for roots of temp and scratch directories)
if not present, both default to /tmp
Fixed Parameter Values for module datavg:
flags -v
mean logP
pkey CMLon
mscale 0.0005
mzero -11.5
copy "CarrRot,LatHG,LonCM,DELTA_K,DELTA_NU,D_OMEGA,Apode_f,APOD_MIN,APOD_MAX"
| param | Size 5 | Size 15 | Size 30 |
| out |
hmi[_test].rdVavgpspec_fd05 |
hmi[_test].rdVavgpspec_fd15 |
hmi[_test].rdVavgpspec_fd30 |
| in |
hmi.rdVpspec_fd05[CarRot] |
hmi.rdVpspec_fd15[CarRot] |
hmi.rdVpspec_fd30[CarRot] |
Sample commands:
avgpsbycr 5 2161
avgpsbycr 15 2096
avgpsbycr 30 2096
Changes from Previous Version
— add version number and report in completion message and runlog
— use cluster 1 for 5-deg tiles
— abort if fewer than 2 available spectra to average
Bugs in current version
Script fails if there are no power spectra for CM Longitude 180.
The tile target list is based on the union of all available tiles for the rotation,
so some tiles may be unfilled if they are not set for CMLon = 180, and some may be
based on a small number of spectra if the tile set changed near the beginning or end
of the rotation.