peace moments andrew lahiff 29 november 2004 10 january 2005 24 january 2005

PEACE momentsPEACE moments

Andrew LahiffAndrew Lahiff

29 November 200429 November 2004

10 January 200510 January 2005

24 January 200524 January 2005

peacemoments - peacemoments - summarysummary

Cluster & Double Star PEACE moments: Cluster & Double Star PEACE moments: number number

density, velocity, pressure, heat flow vector, density, velocity, pressure, heat flow vector,

temperaturetemperature

GSE or B-field aligned coordinate systemsGSE or B-field aligned coordinate systems any distribution, any sensor combinationany distribution, any sensor combination S/C potential correctionsS/C potential corrections plain ASCII, CEF & CDF output filesplain ASCII, CEF & CDF output files

if you run it without thinking, i.e. looking at spectrograms, your if you run it without thinking, i.e. looking at spectrograms, your

results will probably be wrong - more on this later…results will probably be wrong - more on this later…

peacemoments - GUIpeacemoments - GUI

There are 3 pages:There are 3 pages: basicbasic S/C potentialS/C potential integrationintegration

Note: there are helpNote: there are help

bubblesbubbles

peacemoments - GUIpeacemoments - GUI

There are 3 pages:There are 3 pages: basicbasic S/C potentialS/C potential integrationintegration

peacemoments - recent peacemoments - recent changeschanges

v1.2 & below:v1.2 & below:

- called exportIDFS, then Qtran:- called exportIDFS, then Qtran:

cdf cdf qft qft read into peacemoments read into peacemoments

- temporary files were ASCII- temporary files were ASCII

- this was very slow- this was very slow now peacemoments can read data directly now peacemoments can read data directly

from the IDFS archive; temporary files are from the IDFS archive; temporary files are binarybinary

- this is - this is muchmuch faster faster

peacemomentspeacemoments written in C/C++; GUI in Perl/Tk. Needs:written in C/C++; GUI in Perl/Tk. Needs:

- SDDAS- SDDAS- Qtran for CDF output files- Qtran for CDF output files- gnuplot for generating plots- gnuplot for generating plots

the latest versions use the latest versions use IDFS libraries (not in standard SDDAS installation)IDFS libraries (not in standard SDDAS installation) ““cIDFSSource” code from Joey Mukherjee (SwRI)cIDFSSource” code from Joey Mukherjee (SwRI)

to promote & read data from the databaseto promote & read data from the database

available as an executable for Linux, Solaris & Mac OS Xavailable as an executable for Linux, Solaris & Mac OS X

peacemomentspeacemomentsLEEA 3D LEEA 3D velocity velocity distribution functiondistribution function

HEEA 3D HEEA 3D velocity velocity distribution functiondistribution function

sweep mode, sweep mode, sun pulse timessun pulse times

EFW PP potentialEFW PP potential

GSE spin-axis GSE spin-axis anglesangles

FGM PP B-field FGM PP B-field vectorvector

.mom file .mom file moments moments

.mom files gzipped and stored - used next time.mom files gzipped and stored - used next timeyou calculate moments for the same time periodyou calculate moments for the same time period(faster than getting data from database).(faster than getting data from database).

IDFS routines - promotion & reading data from the IDFS archiveIDFS routines - promotion & reading data from the IDFS archive

PEACE - distributionsPEACE - distributions 3D distributions: 3DR, 3DXP, 3DXE, 3DF; 2D: SPINPAD3D distributions: 3DR, 3DXP, 3DXE, 3DF; 2D: SPINPAD

LAR: Low Angular ResoutionLAR: Low Angular ResoutionMAR: Medium Angular ResolutionMAR: Medium Angular ResolutionHAR: High Angular ResolutionHAR: High Angular Resolution

3DF - highest resolution3DF - highest resolution• 15 degree polar res.15 degree polar res.• 22.5 degree azimuth res.22.5 degree azimuth res.

Double Star only has 3DFDouble Star only has 3DF

SPINPAD determined onSPINPAD determined onthe ground from PAD datathe ground from PAD data

peacemoments - distributionspeacemoments - distributions

SPINPADSPINPAD 3DR3DR 3DF3DF 3DX1P, 3DX1E, 3DX13DX1P, 3DX1E, 3DX1 3DX1P + 3DX2P3DX1P + 3DX2P 3DX1P + 3DR3DX1P + 3DR 3DX1E + 3DR3DX1E + 3DR also PEACE onboard (MOMS_L2), HIA PP, also PEACE onboard (MOMS_L2), HIA PP,

CODIF PP, WHISPER PPCODIF PP, WHISPER PP

peacemoments - sensor peacemoments - sensor

combinationscombinations 1-sensor:1-sensor: L1L2B, TH1H2, H1H2, L1L2L1L2B, TH1H2, H1H2, L1L2 2-sensor:2-sensor: TL1L2B, TH1H2B, TL1L2, H1H2BTL1L2B, TH1H2B, TL1L2, H1H2B normally use either normally use either TL1L2B or L1L2BTL1L2B or L1L2B

(HEEA more likely to be (HEEA more likely to be saturated than LEEA)saturated than LEEA)

usually like this (but not always)usually like this (but not always)

Double Star - sweep option 1Double Star - sweep option 1 26 keV

1 keV

34 eV

1 eV Spin 1,3,5,7,… Spin 2,4,6,8,…

Double Star - alternating sweepsDouble Star - alternating sweeps 26 keV

1 keV

34 eV

1 eV Spin 1,3,5,7,… Spin 2,4,6,8,…

Double Star - alternating sweepsDouble Star - alternating sweeps

TopTop

OverlapOverlap

BottomBottom

HH

LL

moments coveringmoments coveringfull energy range - full energy range - time resolutiontime resolution8-seconds8-seconds

Can have similarCan have similar““sensor combinations”sensor combinations”to Cluster, e.g.to Cluster, e.g.TL1L2B, L1L2, …TL1L2B, L1L2, …

Double Star - sensor Double Star - sensor combinationscombinations

use TL1L2B for alternating sweep modeuse TL1L2B for alternating sweep mode

don’t use TH1H2B (peacemoments don’t use TH1H2B (peacemoments won’t actually let you)won’t actually let you)

- frequently the overlap part of the ”high - frequently the overlap part of the ”high sweep preset” sweeps aren’t sweep preset” sweeps aren’t transmittedtransmitted “ “H1H2” doesn’t existH1H2” doesn’t exist

timetagstimetags 3 options:3 options:

PEACE interval averagedPEACE interval averaged PEACE sunpulsePEACE sunpulse CSDS interval averaged (same as PP’s)CSDS interval averaged (same as PP’s)

ttPPPP = t = tPEACE; sunPEACE; sun + t + tdurationduration (0.5 + (SPOS - 75/16)/1024) (0.5 + (SPOS - 75/16)/1024)

For CL, SPOS = 176 or 1200; for DS SPOS = 1104.For CL, SPOS = 176 or 1200; for DS SPOS = 1104.

Note: from v1.41, these timetags also apply to SPINPAD,Note: from v1.41, these timetags also apply to SPINPAD,PEACE onboard, HIA, CODIF & WHISPERPEACE onboard, HIA, CODIF & WHISPER

Available Data Listing - ClusterAvailable Data Listing - Cluster peacemoments has an “Available Data peacemoments has an “Available Data

Listing” command Listing” command (this uses CPMODES)(this uses CPMODES)

Example Output: 17th November 2002, CL-2Example Output: 17th November 2002, CL-2

SPINPAD (LEEA)START: 2002-11-17 00:00:02.549 STOP: 2002-11-17 03:09:03.581 MAR 79START: 2002-11-17 03:10:03.692 STOP: 2002-11-17 05:11:00.002 MAR 71

SPINPAD (HEEA)START: 2002-11-17 00:00:02.549 STOP: 2002-11-17 03:09:03.581 MAR 91START: 2002-11-17 03:10:03.692 STOP: 2002-11-17 05:11:00.002 MAR 89

3DR (LEEA)START: 2002-11-17 00:00:02.549 STOP: 2002-11-17 03:09:03.581

3DR (HEEA)START: 2002-11-17 00:00:02.549 STOP: 2002-11-17 03:09:03.581

3DX1 (LEEA)START: 2002-11-17 03:10:03.692 STOP: 2002-11-17 05:09:51.950 [N]

Note that each time interval listed has a single sweep mode and sweepNote that each time interval listed has a single sweep mode and sweeppreset.preset.

Available Data Listing - Double Available Data Listing - Double StarStar

Example Output: 13:20 to 14:20 11th March 2004, TC-1Example Output: 13:20 to 14:20 11th March 2004, TC-1

3DF (HEEA)3DF (HEEA)--------------------START: 2004-03-11 13:20:00.156 STOP: 2004-03-11 14:19:56.630 MAR 91 63START: 2004-03-11 13:20:00.156 STOP: 2004-03-11 14:19:56.630 MAR 91 63

Output: GSEOutput: GSE can select GSE or B-field aligned outputcan select GSE or B-field aligned output running the program with GSE output running the program with GSE output

gives all of the following moments:gives all of the following moments:* n* n* * vv* * PP* * QQ* T* T

plain ASCII, cef & cdf output filesplain ASCII, cef & cdf output files

Output: B-field alignedOutput: B-field aligned running the program gives all of the following running the program gives all of the following

moments: moments: * n* n

* * vvpara para vvperpperp

* * QQpara para QQperpperp

* P* Pparapara P Pperp perp PPperp1 perp1 P Pperp2perp2

* T* Tparapara T Tperpperp T Tperp1perp1 T Tperp2perp2

P rotated so that the z-axis is parallel to the B-fieldP rotated so that the z-axis is parallel to the B-field PPpara para = P= Pzzzz P Pperp1perp1 = P = Pxxxx P Pperp2perp2 = P = Pyy yy PPperp perp = (P= (Pperp1 perp1 + P+ Pperp2perp2)/2)/2

Also,Also, T = P/(nk)T = P/(nk)

where, e.g.,where, e.g.,

vvparapara = = BB ( (vv..BB)/()/(BB..BB))

vvperpperp = = vv - - vvparapara

BB is the FGM PP B-field vector. is the FGM PP B-field vector.

Configuration filesConfiguration files can save/load configuration filescan save/load configuration files

- stores start & stop times & dates, distribution, sensor - stores start & stop times & dates, distribution, sensor combination, S/C number, coordinate system, S/C potential combination, S/C number, coordinate system, S/C potential offsets, etc.offsets, etc.

- preferences are not saved in these files (temperature units, - preferences are not saved in these files (temperature units, timetag convention, etc)timetag convention, etc)

preferences can be savedpreferences can be saved ‘‘Calculate All in Current Directory’ commandCalculate All in Current Directory’ command

- calculates moments using all configuration files in the - calculates moments using all configuration files in the current directorycurrent directory

running from the command linerunning from the command line

peacemoment can be run from the peacemoment can be run from the command line - example:command line - example:

peacemoments cl1_25Jan2005.cpmpeacemoments cl1_25Jan2005.cpm

peacemoments_cmd 2002 025 0 0 2002 025 0 0 CL1 3DR peacemoments_cmd 2002 025 0 0 2002 025 0 0 CL1 3DR TL1L2B ON 1 2 bins NONE CSDS_AVERAGED eVTL1L2B ON 1 2 bins NONE CSDS_AVERAGED eV

Restricting the integration Restricting the integration rangesranges

the energy, polar & azimuth integration the energy, polar & azimuth integration

ranges can be restrictedranges can be restricted the lower & upper limits can be different the lower & upper limits can be different

for each spin if necessary for each spin if necessary

- useful if you want to isolate beams- useful if you want to isolate beams peacemoments reads in the limits from peacemoments reads in the limits from

text filestext files

Calculation of momentsCalculation of moments moments (relative to the bulk flow) of the moments (relative to the bulk flow) of the

velocity distribution function, e.g.:velocity distribution function, e.g.:

n = n = ∫∫ f( f(vv).d).d33vv

nnuu = = ∫∫ f( f(vv).).vv.d.d33vv

ppijij = m = m ∫∫ (v (vii - u - uii).(v).(vjj - u - ujj).f().f(vv).d).d33vv

qq = (m/2) = (m/2) ∫∫ ( (vv - - uu))22.(.(vv - - uu).f().f(vv).d).d33vv

Integration is over the whole of velocity space.Integration is over the whole of velocity space.

Calculation of momentsCalculation of moments

standard moments (these do standard moments (these do not contain the drift velocity):not contain the drift velocity):

n = n = ∫∫ f( f(vv).d).d33vv

nnuu = = ∫∫ f( f(vv).).vv.d.d33vv

PPijij = m = m ∫∫ v vii .v .vjj . f( . f(vv).d).d33vv

QQ = (m/2) = (m/2) ∫∫ vv22. . vv . . f(f(vv).d).d33vv

relative moments:relative moments:

n = nn = n

uu = = uu

pij = Pij - m n upij = Pij - m n uii u ujj

qq = = QQ - - uu Tr(P)/2 - Tr(P)/2 - uu.P.P + m n u+ m n u22 u u

peacemoments calculates moments the same way peacemoments calculates moments the same way as the on-board PEACE moments are calculatedas the on-board PEACE moments are calculated

Calculation of momentsCalculation of moments

we only have 4 integrals to calculate:we only have 4 integrals to calculate:

MOMENTMOMENTii = = ∫∫ f( f(vv).M).Mii((vv)d)d33vv

wherewhere

MM11((vv) = 1) = 1 MM22((vv) = ) = vv

MM3ij3ij((vv) = mv) = mviivvjj MM44((vv) = (m/2) v) = (m/2) v22 vv

Note:Note:

Calculation of momentsCalculation of moments we know we know f(f(vv)) for a finite number of energy, polar & azimuth for a finite number of energy, polar & azimuth

binsbins

e.g. for 3DRe.g. for 3DR15 energy bins covering E15 energy bins covering Eminmin to E to Emaxmax

6 polar bins covering 0 - 180 degrees6 polar bins covering 0 - 180 degrees16 azimuth bins covering 0 - 360 degrees16 azimuth bins covering 0 - 360 degrees

rewrite each integral as a sum over integrals for each bin, e.g.rewrite each integral as a sum over integrals for each bin, e.g.

Calculation of momentsCalculation of moments Therefore:Therefore:

becomesbecomes

wherewhere

Calculation of momentsCalculation of moments Assume Assume f(f(, , , v), v) doesn’t vary much inside each bin, therefore doesn’t vary much inside each bin, therefore

becomesbecomes

These integrals are carried out analytically.These integrals are carried out analytically.

For more information see For more information see Extracting the Bulk Parameters of a Particle Extracting the Bulk Parameters of a Particle DistributionDistribution by R.A. Gowen & M.A. Birdseye by R.A. Gowen & M.A. Birdseye

Spacecraft potential Spacecraft potential corrections corrections (1)(1)

when carrying out the energy integration:when carrying out the energy integration:

- ignore any energy bins if- ignore any energy bins if

VVSCSC > E > Elower lower (here V(here VSCSC is the S/C potential) is the S/C potential)

This removes the photoelectrons. This removes the photoelectrons.

- remaining energy bins are rescaled- remaining energy bins are rescaled

EElowerlower E Elowerlower - V - VSCSC

EEupperupper E Eupperupper - V - VSCSC

This removes the acceleration of the electrons by the S/C This removes the acceleration of the electrons by the S/C

potential.potential.

Spacecraft potential Spacecraft potential corrections corrections (2)(2)

can use EFW PP probe potential, or an ASCII filecan use EFW PP probe potential, or an ASCII file

two S/C potential offsetstwo S/C potential offsets- - offset 1offset 1: difference between EFW probe: difference between EFW probeand spacecraft potential (value from EFW and spacecraft potential (value from EFW team, typically 1 eV)team, typically 1 eV)

VVSCSC = V = VEFWEFW + V + Voffset 1offset 1

- - offset 2offset 2: for removing photoelectrons : for removing photoelectrons above the EFW probe potential (either in above the EFW probe potential (either in bins or eV). Cuts off the integration at a bins or eV). Cuts off the integration at a higher energy.higher energy.

Important to set offset 2 correctly. Usually leave offset1 at default Important to set offset 2 correctly. Usually leave offset1 at default value.value.

Increasing offset 1:Increasing offset 1:Increases VIncreases VSCSC

Increasing offset 2:Increasing offset 2:- removes more removes more energy bins above Venergy bins above VSCSC

User-defined S/C potential fileUser-defined S/C potential file

if there’s no EFW data:if there’s no EFW data:

ask someone from Cassini for helpask someone from Cassini for help

look at spectrogram & determine the potential by eyelook at spectrogram & determine the potential by eye

if the potential for another S/C exists & is similar:if the potential for another S/C exists & is similar:

efwgenerate 2002-025 10:00:00 2002-025 12:00:00 CL1efwgenerate 2002-025 10:00:00 2002-025 12:00:00 CL1

generates an ASCII file containing EFW datagenerates an ASCII file containing EFW data

User-defined S/C potential fileUser-defined S/C potential file making your own S/C potential ASCII file:making your own S/C potential ASCII file:

2002-01-25 10:00:00 10.02002-01-25 10:00:00 10.0

2002-01-25 12:00:00 15.02002-01-25 12:00:00 15.0

2002-01-25 10:00:00 10.02002-01-25 10:00:00 10.0

2002-01-25 10:59:59 10.02002-01-25 10:59:59 10.0

2002-01-25 11:00:00 15.02002-01-25 11:00:00 15.0

2002-01-25 12:00:00 15.02002-01-25 12:00:00 15.0

VV

12:00:0012:00:0010:00:0010:00:00

VV

12:00:0012:00:0010:00:0010:00:00 11:00:0011:00:00

Examples & removing photoelectronsExamples & removing photoelectrons

SPINPAD 2-sec resolution SPINPAD 2-sec resolution momentsmoments

• when HEEA and LEEA cover the same energy rangewhen HEEA and LEEA cover the same energy range

• this has slightly better time resolution than 3DR!this has slightly better time resolution than 3DR!

It’s sometimes easy… It’s sometimes easy… (1)(1)

• in this example (offset1 = 1 & offset2 = 0), i.e. Vin this example (offset1 = 1 & offset2 = 0), i.e. VSCSC = V = VEFWEFW + 1 + 1

PEACEPEACEHIAHIA

It’s sometimes easy… It’s sometimes easy… (3)(3)

PEACE (TH1H2B)PEACE (TH1H2B)PEACE (TL1L2B)PEACE (TL1L2B)HIAHIA

• in this example (offset1 = 1 & offset2 = 0), i.e. Vin this example (offset1 = 1 & offset2 = 0), i.e. VSCSC = V = VEFWEFW + 1 + 1

……but this is more commonbut this is more common

PEACEPEACEHIAHIA

Comparisons with CISComparisons with CIS

If PEACE and CIS disagree:If PEACE and CIS disagree: CIS might be in the wrong mode (look at the status CIS might be in the wrong mode (look at the status

bytes)bytes)

you need to adjust the S/C potential offset to ensureyou need to adjust the S/C potential offset to ensure photoelectrons have been removed completelyphotoelectrons have been removed completely ……

Importance of S/C potential Importance of S/C potential correctionscorrections

Removing photoelectrons Removing photoelectrons (1)(1)

removal of all photoelectrons (this can’t removal of all photoelectrons (this can’t easily be automated)easily be automated)

an example spectrogram for < 100 eVan example spectrogram for < 100 eV

• using just EFW probe potential won’t remove all photoelectronsusing just EFW probe potential won’t remove all photoelectrons• offset 2 in peacemoments needs to be used to remove photoelectrons aboveoffset 2 in peacemoments needs to be used to remove photoelectrons above the EFW probe potentialthe EFW probe potential• need to be careful if there are cold electronsneed to be careful if there are cold electrons


1.1.

2.2.

3. 3.

• may need to use Steve Schwartz’ QTMC software to fill in the gap between the S/C potential and the lowest energymay need to use Steve Schwartz’ QTMC software to fill in the gap between the S/C potential and the lowest energy measured by peacemeasured by peace

• need to be careful if there are cold electronsneed to be careful if there are cold electrons


Example 1:Example 1:

Importance of setting offset 2 correctlyImportance of setting offset 2 correctly

green & bluegreen & bluecuvers overlapcuvers overlap


Example 2:Example 2:

Importance of setting offset 2 correctlyImportance of setting offset 2 correctly


peacemoments can generate peacemoments can generate SpectroScalar layouts showing PEACE SpectroScalar layouts showing PEACE data & EFW potentialdata & EFW potential- meant to encourage users to look at - meant to encourage users to look at spectrograms!spectrograms!

- use “lower band edge” not “centre” energies!- use “lower band edge” not “centre” energies!

Problems & interference from other Problems & interference from other instrumentsinstruments

Vz problemVz problemPEACEPEACEHIAHIA

CL-1CL-1 CL-3CL-3

in 3D & on-board momentsin 3D & on-board moments

HEEAHEEAusuallyusuallyworse thanworse thanLEEALEEA

HEEA-LEEA velocity mirror HEEA-LEEA velocity mirror effecteffect

PEACE (TH1H2B)PEACE (TH1H2B)PEACE (TL1L2B)PEACE (TL1L2B)HIAHIA

Vz problem alsoVz problem alsocan be seen here!can be seen here!

HEEA saturationHEEA saturation sometimes HEEA can be saturatedsometimes HEEA can be saturated

- use TL1L2B instead of TH1H2B- use TL1L2B instead of TH1H2B

- use L1L2B instead of TH1H2- use L1L2B instead of TH1H2 (if LEEA & HEEA cover same (if LEEA & HEEA cover same energy range)energy range)

HEEAHEEALEEALEEA

PEACE secondary PEACE secondary electronselectrons

• need to work out what to do about them…need to work out what to do about them…

PEACE internal electronsPEACE internal electrons

When the sun is within the field of view of the analyser, the influx of solar ultraviolet radiation produces low-energy photoelectrons.

WHISPER interference WHISPER interference (1)(1)

sometimes moments have spikes at sometimes moments have spikes at regular intervals:regular intervals:

Problems caused by EDIProblems caused by EDI

• in very low density plasmas there can sometimes be severe EDI-inducedin very low density plasmas there can sometimes be severe EDI-induced spacecraft chargingspacecraft charging

additional moments-related softwareadditional moments-related software

peaceplotpeaceplot Command-line program for instantly Command-line program for instantly

plotting moments comparisons (no plotting moments comparisons (no endless menus and mouse-clicking)endless menus and mouse-clicking)

plots: number density, speed, velocity, plots: number density, speed, velocity, temperature, anisotropy, current densitytemperature, anisotropy, current density

PS or PDF output filesPS or PDF output files

All the moments plots in this talk were generated using peaceplot.All the moments plots in this talk were generated using peaceplot.

peaceplot - examplepeaceplot - examplepeaceplot -date 2003-007 -cl 4 -npeaceplot -date 2003-007 -cl 4 -n

peacecurrentpeacecurrent

calculates current densities from calculates current densities from peacemoments output filespeacemoments output files

current density = electron charge x number density x velocitycurrent density = electron charge x number density x velocity

has ASCII and CDF output fileshas ASCII and CDF output files

gives parallel & perpendicular currents from gives parallel & perpendicular currents from peacemoments B-field aligned output filespeacemoments B-field aligned output files

peacetppeacetp Calculates transition parameters as described in Calculates transition parameters as described in Re-Re-

ordered Electron Data in the Low-Latitude Boundary Layerordered Electron Data in the Low-Latitude Boundary Layer by M.A. Hapgood and D.A. Bryantby M.A. Hapgood and D.A. Bryant

The transition parameter (TP) indicates the range of The transition parameter (TP) indicates the range of progression between magnetosheath (TP = 0) and the progression between magnetosheath (TP = 0) and the magnetosphere (TP = 100)magnetosphere (TP = 100)

•fit polynomial to log[Tfit polynomial to log[Tee] versus log[N] versus log[Nee]]•project each data point onto the curveproject each data point onto the curve•calculate distance from each projectedcalculate distance from each projected data point along the curve to an arbitrarydata point along the curve to an arbitrary origin beyond the magnetosheath endorigin beyond the magnetosheath end of the plot of the plot raw TP raw TP•raw TP normalized to the range 0-100raw TP normalized to the range 0-100

peacetppeacetp needs a peacemoments ASCII output file needs a peacemoments ASCII output file

or any ASCII file containing date, time, or any ASCII file containing date, time, number density & temperaturenumber density & temperature

generates an ASCII file containing generates an ASCII file containing timetags and transition parameterstimetags and transition parameters

peacecdfcombinepeacecdfcombine

combines any number of cdf files, e.g.combines any number of cdf files, e.g.

peacecdfcombine file1.cdf file2.cdf file3.cdf …peacecdfcombine file1.cdf file2.cdf file3.cdf …

uses the header from the first file and the uses the header from the first file and the data from all filesdata from all files

Other miscellaneous utilitiesOther miscellaneous utilities

day-of-year day-of-year month & month & dayday

to convert YYYY-MM-DD to YYYY-to convert YYYY-MM-DD to YYYY-DOYDOY

daytodoy 2004-05-03daytodoy 2004-05-03

to convert YYYY-DOY to YYYY-MM-to convert YYYY-DOY to YYYY-MM-DDDD

doytoday 2003-127doytoday 2003-127

ssleditssledit command-line program for instantly making command-line program for instantly making

global changes to SpectroScalar layoutsglobal changes to SpectroScalar layouts

Very useful whenVery useful whenyou have manyyou have manyspectrograms!spectrograms!

(previously called peacess)

ssledit - examplesssledit - examples change start & stop timeschange start & stop times

ssledit layout.lay -st 12:00:00 -sp 15:00:00ssledit layout.lay -st 12:00:00 -sp 15:00:00

change start time, date & S/C numberchange start time, date & S/C numberssledit layout.lay -st 10:00:00 2002-025 -sc 3ssledit layout.lay -st 10:00:00 2002-025 -sc 3

change y-axis limitschange y-axis limitsssledit layout.lay -ymin 10 -ymax 1000ssledit layout.lay -ymin 10 -ymax 1000

change colour scalechange colour scalessledit layout.lay -cmin 1e-06 -cmax 32000ssledit layout.lay -cmin 1e-06 -cmax 32000

also can change units, logalso can change units, loglinear scales, etc.linear scales, etc.

sslbatchsslbatch

for generating many spectrograms using the same layout, for generating many spectrograms using the same layout, but differentbut different- times & dates- times & dates- S/C numbers- S/C numbers

gif & ps output optionsgif & ps output options

Example:Example:

sslbatch file.lay listfilesslbatch file.lay listfile

where listfile is of the form:where listfile is of the form:

2 2002-055 14:00:00 2002-055 14:55:002 2002-055 14:00:00 2002-055 14:55:004 2003-05-04 18:30:00 2003-05-04 19:40:004 2003-05-04 18:30:00 2003-05-04 19:40:001 2003-05-04 22:30:15 2003-05-04 23:55:001 2003-05-04 22:30:15 2003-05-04 23:55:00

a batch interface to SpectroScalar

peacepromotepeacepromote

For very easily promoting all data required to calculate moments onto For very easily promoting all data required to calculate moments onto e.g. a laptope.g. a laptop

peacepromote file.txtpeacepromote file.txt

An example file.txt:An example file.txt:

2004 39 TC1 FP MOM 2004 39 TC1 FP MOM <- gets TC-1, 3DF + pitch angles, + everything for moments<- gets TC-1, 3DF + pitch angles, + everything for moments

2004 39 CL1 S 2004 39 CL1 S <- gets CL-1, SPINPAD<- gets CL-1, SPINPAD

2004 39 CL2 R MOM2004 39 CL2 R MOM <- gets CL-2, 3DR, + everything for moments<- gets CL-2, 3DR, + everything for moments

For Cluster, will promote : CPCGP, CPMODES, AUX, EFW & FGMFor Cluster, will promote : CPCGP, CPMODES, AUX, EFW & FGMFor Double Star, will promote: DPGEN, AUX, & FGMFor Double Star, will promote: DPGEN, AUX, & FGM

batch promotion toolbatch promotion tool

moments in the IDFS archivemoments in the IDFS archive

moments in the databasemoments in the database

MOMENTSMOMENTS

- on-board moments with problems- on-board moments with problems MOMS_L2MOMS_L2

- on-board moments with corrections- on-board moments with corrections MOMS_L3MOMS_L3

- current version of the “best” - current version of the “best” momentsmoments

status of MOMS_L2status of MOMS_L2 do not use MOMENTSdo not use MOMENTS

MOMS_L2 : onboard moments but withMOMS_L2 : onboard moments but with- corrected geometric factors- corrected geometric factors- corrected factor of 2- corrected factor of 2- available in GSE coordinates- available in GSE coordinates- (now) has correct pressure tensor & temperature- (now) has correct pressure tensor & temperature

MOMS_L2 currently available for 2001 to MOMS_L2 currently available for 2001 to 2004-02-132004-02-13

status of MOMS_L2status of MOMS_L2 problems still in MOMS_L2:problems still in MOMS_L2:

based on on-board moments which use on-based on on-board moments which use on-board energy tables, efficiencies & response board energy tables, efficiencies & response surfaces. Any errors in these parameters surfaces. Any errors in these parameters cannot be corrected in ground processing.cannot be corrected in ground processing.

status of MOMS_L3status of MOMS_L3

currently corresponds to version 3 PEACE currently corresponds to version 3 PEACE PPPP

- generated using MOMS_L2 moments- generated using MOMS_L2 moments

- rules to decide which sensor - rules to decide which sensor combination to use depending on the S/C combination to use depending on the S/C potentialpotential

e.g. if B contains mostly photoelectrons, e.g. if B contains mostly photoelectrons, use TL1L2 (not TL1L2B)use TL1L2 (not TL1L2B)

• Example number density comparions with HIAExample number density comparions with HIA

the endthe end

peace moments andrew lahiff 29 november 2004 10 january 2005 24 january 2005

Documents