SCG Data - How to get it, how to compile it, where to put it
We have scripts for these tasks: HOW.TO files that can be sourced
- rsync memo
- GWR-DDAS
- Monument
temperatures and Invar rod
- Onsala meteo
station
- Onsala
mareograph
- Onsala rain data
- SMHI rain
data
- SMHI
tide
gauge data
- USGS earthquake events
- IERS Polar
motion
- BKG Atmacs
rsync:
You are running cygwin on ore
#
# From froste:
#
cd
/cygdrive/f/tidedata/GWR_o054/2009o054/RAW_o054
rsync -avp
hgs@holt.oso.chalmers.se:/home/hgs/tidedata/GWR_o054/2009o054/RAW_o054/
.
cd
/cygdrive/f/tidedata/GWR_o054/2009o054/MON_o054
rsync -avp
hgs@holt.oso.chalmers.se:/home/hgs/tidedata/GWR_o054/2009o054/MON_o054/
.
#
# To froste:
cd
/cygdrive/f/tidedata/GWR_o054/2009o054/RAW_o054
rsync -avp *
hgs@holt.oso.chalmers.se:/home/hgs/tidedata/GWR_o054/2009o054/RAW_o054/
cd
/cygdrive/f/tidedata/GWR_o054/2009o054/MON_o054
rsync -avp *
hgs@holt.oso.chalmers.se:/home/hgs/tidedata/GWR_o054/2009o054/MON_o054/
#
# for transferring gravity
data only; there are files and directories for processing.
#
rsync -lptgoDvp
hgs@holt.oso.chalmers.se:/home/hgs/tidedata/GWR_o054/2009o054/PM/
PM/
rsync -rlptgoDvp
hgs@holt.oso.chalmers.se:/home/hgs/tidedata/GWR_o054/2009o054/tgg/
tgg/
#
# first without subdirectory
recursion, specific directory, then with recursion
# from froste to local
#
GWR-DDAS Gravity files
HOW TO Copy the DDAS files to the directory Tidedata and to the
Tidearchive on machine Elder.
You'll be on Elder Cygwin.
Login to Elder by ssh.
cd Tidedata
cp-from-sg
-X #MM ['day-expression']
for month #MM, example
cp-from-sg
-X 02 '0[1-3]'
for the first three days of February the current year.
..
Collects data and compresses files in ~/Tidedata and ~/Tidearchive
(the latter is on the external harddisk).
You have the same directory structure as on the DDAS PC, namely
yyyyo054/DAY_o054 ../MON_o054 ../RAW_o054
On Holt there is a remote-copy script in ~/TD. To
copy "Raw" data from Elder, use
scpe -R 'filename-expression'
example
scpe -R '??1302??.054*'
to copy the files for February 2013. The trailing `*´ globs
the name extension of compressed files.
Monument temperatures and invar rod
Get an update from OSO's weather server, wx:
cd ~/wx/TEMP
get-monutemps [#year]
Collects T_2_$year.ts ... T_5_$year.ts
and T_11_$year.ts, four
temperature sensors and the Invar rod.
The cron job:
crontab -l
# m
h dom mon dow command
40 01 * * *
/home/hgs/wx/TEMP/get-monutemps >
/home/hgs/wx/TEMP/get-monutemps.log 2>&1
Onsala meteo
OSO Rain data
Here is how to get annual ts-files for
rain rate (5min averages)
rain-1h (1h accumulation)
24h-rain (24h accumulation)
NOTE that from 2012 on the sampling
interval in the weather logs is 60s and there is an offset of
30s.
Thus, except for the daily rain ORD${y}.ts,
60s is the sampling rate of the *.ts files
cd ~/wx/ARC
get-oso-weather [COLLECT] [#year]
Without the COLLECT option
the script will wget new files from the weather server.
With or without the option, the script collects
ORD$y.ts Daily rain [mm]
ORH$y-1h.ts
Hourly rain [mm]
ORTH$y.ts TEMP_________VAL [deg C]
HUMI_________VAL [%]
OWSD$y.ts
WSPEED_______VAL [m/s]
WDIR_________VAL [deg]
ORR$y-5m.ts 5-min rain rate [mm/h]
for the year
specified, or the current year as the default.
We run a cron job:
crontab -l
# m
h dom mon dow command
30 01 * * *
/home/hgs/wx/ARC/get-oso-weather >
/home/hgs/wx/ARC/get-oso-weather.log 2>&1
OSO weather
ReadMe
OSO mareograph - home page
URGENTLY INCLUDE THIS IN THE DAILY CRON
Data is available from Nov 23, 2011.
SMHI rain data
weekly ftp, see local README file at
file:///c:/cygwin/home/hgs/SMHI/RAIN/HOW.TO
froste: ~/SMHI/RAIN/HOW.TO
The data part looks like this
rm -f
ROX-1d.txt
touch
ROX-1d.txt
foreach
file ( onsala*.txt )
awk -F,
'(NR>1){printf "%s %9.2f\n",$1,$2}' $file >> ROX-1d.txt
end
SMHI tide gauge data
Here, a cron job automatically downloads and decodes SMHI html files
crontab -l
# m
h dom mon dow command
45 07 *
* * /home/hgs/SMHI/getsealevel -g r3g.txt >>
/home/hgs/SMHI/crontouch 2>&1
45 19 *
* * /home/hgs/SMHI/getsealevel -g r3g.txt >>
/home/hgs/SMHI/crontouch 2>&1
getsealevel is in ~/SMHI/AUTO
The data is stored in ~/SMHI/AUTO
The collecting ts-files are in ~/SMHI
m
RINGHALS-2010-04-01-06.00.dat
2010 04
01 06 0.0
2010 04
01 07 0.0
2010 04
01 08 -7.0
...
cd ~/SMHI
cat AUTO/RING* | tslist -
-gI4,3i3,f7.0]W- -k1 -I -o tgg-rin-100301.ts
cat AUTO/GOT* |
tslist - -gI4,3i3,f7.0]W- -k1 -I -o
tgg-got-100301.ts
Copies for Johan Löfgren
tslist tgg-rin-100301.ts -C1 -Ff10.1 -qqq > !
~/www/4JL/tgg-rin-100901.tsf
tslist tgg-got-100301.ts -C1 -Ff10.1 -qqq > !
~/www/4JL/tgg-got-100901.tsf
USGS Earthquake events
cd ~/Seism/USGS
cp usgs-JD-M.dat usgs-JD-M.bup
tail usgs-JD-M.bup
JDC -j -m 55245
2010
02 18
set bdate = ( 2010 02 19 )
set sdate = ( 2010 04 27 )
source wget-epic
cp epic.dat epic-100427.bup
cat usgs-JD-M.dat
>> usgs-JD-M.bup
cp usgs-JD-M.bup usgs-JD-M.dat
This gives you only the MJD and the moment magnitude. Anything else
must process the epic-files.
A csh-script can be found in ~/Seism/USGS/EXTENDREC
EXTENDREC with no
arguments appends at the last record on usgs-JD-M.dat and
through the current date.
Absolute gravity
This is a little obsolete. If you only want to
process sets, have a look at /TD/a/Allsets
The remainder of this § considers drops.
Put the set.txt files into a subdir of ~/TD, e.g. ~/TD/cal1004
hinting at the date of the experiment
Find out the first day, find out the instrument number
Hannover: 220 "Olga"
Lantmateriet: 233 "Greta"
cd ~/TD
ls cal1004/*.set.txt
locate 'cal1004/*.project.txt'
-d ~/CYGHGS/grav + 'Meter S/N'
Do the conversion
cd ~/TD
FG5-set2ts -m -i FG220 -d a
100419 cal1004/AA_*.set.txt
This needs still improvement. If we relate everything to platform AA
we need the gravity ties.
Provide SCG data at AG times: G1_1_yymm_resmp4AG.ts
Use TD/cal1106/synch-with-drops.tse
tslist <appended 1-s SCG>
-Ecal1106/synch-with-drops.tse,DROPS -I -o G1_1_1106_resmp4olga.ts
TD/cal was a bit special: the set file was not reduced for
predictable effects.
The following is to no avail, since there will be an arbitrary
offset.
setenv FILESUB cal/G1_1_0907_resmp4greta.ts
tslist
cal/Onsala_AC_20090706b.set.ts -Ecal/ag-scg.tse,DIFF -L'A|V'
-L'A|S' -C3 -Y-9817159000 -SO10,10 -F1p,2d12.4 -qqq -w
cal/ag090706-scg.tsf
FG5-set2ts -d a -u nm/s2 -f
i4,2i3,1x,3i3,t37 -i FG233 -t cal/ag090706-scg.tsf 090706
but it shows the use of other options to FG5-set2ts
Next attempt, compute own corrections. Possibly pertaining problems,
e.g. whether the M0S0 tide is subtracted. Anyway:
get baro data, theo-tide and polmot, and subtract
tsf2ts-decim -r 6 10 -o cal -c
3 -d 600,300 -a RAW_o054/ G1 0907
tslist cal/G1_3_090706.dat
-A'2009 07 06 13 25]' -g'(i4,2i3,i4,2i3,t36,e19.0)' -k3 -rs600 -C3
-F1p,e14.5 -I -Y-1012.35 -Un563 -o cal/G1_3_090706-600s.ts
(the 563 is the number of 10-min samples through the following date
range; do tslist
cal/Onsala_AC_20090706b.set.ts -I -C3 -qbe )
run_urtip cal/ag090706.mc ; cp o/pt.ts cal/THG090706-600s.ts
polmotm -d 1.16 -l 11.9,57.6
cal/G1_3_090706-600s.ts | fgrep -v '<' | tslist -
-gI4,4i3,t66,f10.0 -k2 -C1 -I -o cal/PMG090706-600s.ts
tsl
cal/Onsala_AC_20090706b.set.ts -Ecal/ag-scg.tse,THT -L'A|V'
-C3 -Y-9817159000 -I -o cal/ag090706-tht.ts # for plotting with
tsd
tslist
cal/Onsala_AC_20090706b.set.ts -Ecal/ag-scg.tse,THT -L'A|V'
-L'A|S' -C3 -Y-9817159000 -F2f10.3 -SO1,10 -qqq -w
cal/ag090706-tht.tsf
look at
/home/hgs/tslist-logs/tidedata%GWR_o054%2009o054%cal%ag090706-tht.ts.log
IERS Polar motion
After a major strategy change, I now recommend to renew the IERS
latest version EOPC IAU2000 file,
long-term (since 1962 and daily). It is called 214 (? always?) and
available with
wget
http://datacenter.iers.org/eop/-/somos/5Rgv/latest/214
RENEW -W includes the wget step and a backup to ~/TD/PM/old/
cd ~/TD/PM
RENEW -W 2009 06 01
To calculate a series parallel with a measurement file, do
polmotm -l 11.9,57.6 d/g090701-140131-1h.mc |\
fgrep -v '<' |\
tslist - -gi4,5i3,t66,f10.0 -k3 -I -o
d/PMG090701-140131-1h.ts
DONE!
Table header EOP Finals, annual:
FORMAT(3(I4),I7,2(F11.6),2(F12.7),2(F11.6),2(F11.6),2(F11.7),2F12.6)
**********************************************************************************
Date
MJD
x
y
UT1-UTC
LOD
dX
dY x
Err y Err UT1-UTC Err
LOD Err dX
Err dY Err
"
"
s
s
"
"
"
"
s
s
" "
(0h UTC)
2012 1 1 55927
0.118700 0.263209 -0.4190307
0.0012643 0.000520 -0.000064
0.000019 0.000021 0.0000070
0.0000055 0.000134 0.000215
2012 1 2 55928
0.117478 0.262407 -0.4202660
0.0011668 0.000506 -0.000044
0.000018 0.000020 0.0000043
0.0000055 0.000107 0.000169
2012 1 3 55929
0.115947 0.261649 -0.4213914
0.0010732 0.000487 -0.000035
0.000018 0.000020 0.0000046
0.0000055 0.000080 0.000127
Table header Bulletin B, monthly:
1 - DAILY FINAL VALUES OF x, y, UT1-UTC, dX, dY
Angular unit is milliarcsecond (mas), time unit is
millisecond (ms).
Upgraded solution from February 1 2011 - consistent
with ITRF 2008.
DATE
MJD
x
y
UT1-UTC dX
dY x err y
err UT1 err X err Y err
(0 h
UTC)
mas
mas
ms
mas mas
mas
mas ms
mas mas
Final values
Mean formal error
0.028 0.027
0.0098 0.028 0.036
2012 8 2 56141
147.315 399.345 409.5109 -0.259
-0.018 0.020
0.023 0.0248 0.018 0.023
2012 8 3 56142
148.473 398.802 409.0457 -0.267
0.054 0.020
0.023 0.0056 0.024 0.030
2012 8 4 56143
149.249 398.140 408.4970 -0.253
0.091 0.020
0.023 0.0047 0.027 0.034
BKG Atmacs - This is basically
~/TD/Atmacs/HOW.TO
# We have a cron job script ~/TD/get-atmacs doing the
following
#
#
# To get the wget part done, you must specify
#
# get-atmacs -w
#
# It will update the MC file Atmacs/atmac.mc [sic!]
#
# Then you can play with
# atmacs2ts [-r 6]
# to get MC file Atmacs/atmacs-3h.mc
[default] or,
# with -r 6, MC file Atmacs/atmacs-3h.mc
#
#########################################################################
#
# http://atmacs.bkg.bund.de/docs/computation.php
says a.o.:
#
# Extract the model air pressure (col 2) from the total
atmospheric effect
# using a constant admittance factor being typical for the
local zone
# (e.g. 2...2.5 nm/s2/hPa). Interpolate and and
re-intruduce the air pressure
# basing on your high resolution data using the same
admittance factor. It
# could be necessary to vary the admittance factor in order
to obtain the
# best result.
#
# Data is explained here:
http://atmacs.bkg.bund.de/docs/data.php
# "gm" - glocal model
# "lm" - local model
#
# ( Newton ) + ( Load ) = (lm#3 + lm#4 + gm#3) + (gm$4)
#
# To get data sets, issue
wget
http://atmacs.bkg.bund.de/data/results/lm/os_lm2_12km_13deg.grav
#
1)
wget
http://atmacs.bkg.bund.de/data/results/gm/os_gm384_13deg.grav
#
2)
wget
http://atmacs.bkg.bund.de/data/results/gm/os_gm256_13deg.grav
# 3)
# 1) is from Feb 2., 2010 --
# 2) is from March 2012 --
# 3) is from Feb 2., 2010 through Feb. 2012
#
# If you would like to start all over:
set atmac_new =
if ( $?atmac_new == 1 ) then
rm -f atmac.mc
tslist
os_gm256_13deg.grav -g'i4,3i2,f13.0)' -k1 -C1 -S1d9 -I \
-O:`label ATMG,GMN` atmac.mc
tslist
os_gm256_13deg.grav -g'i4,3i2,t24,f13.0)' -k1 -C1 -S1d9 -I \
-O:`label ATMG,GML` atmac.mc
awk
'\!/^201001/||/^20100201/{print}' os_lm2_12km_13deg.grav |\
tslist -
-g'i4,3i2,t25,f13.0)' -k1 -C1 -S1d9 -I \
-O:`label ATMG,LMN` atmac.mc
awk
'\!/^201001/||/^20100201/{print}' os_lm2_12km_13deg.grav |\
tslist -
-g'i4,3i2,t38,f13.0)' -k1 -C1 -S1d9 -I \
-O:`label ATMG,RMN` atmac.mc
awk
'\!/^201001/||/^20100201/{print}' os_lm2_12km_13deg.grav |\
tslist -
-g'i4,3i2,f14.0)' -k1 -C1 -S0.01 -I \
-O:`label BARO,VAL` atmac.mc
# Check with tslist atmac.mc `tslql -L
atmac.mc` -C1 -F5f10.2 | m
# Labels:
# _GMN - Global Newton
# _GML - Global Load
# _LMN - Local Newton
# _RMN - Region Newton
# _VAL - Local model pressure
#
# Adding together and output the sum:
tslist atmac.mc
`tslql -L atmac.mc` -C1 -F10f10.2 -RMA:1,2,3,4=1 -I \
-o atmac-sum.ts
endif # atmac_new
#
Atmacs interpolation with the gain method
#
# Only the regional and local effects - determine daily gain
factors
tslist atmac.mc -L'A|L' -L'A|R' -L'B|V' \
-RMLe:1.,1,1.,2=1
-RMA8,0,1:1,3,W=HANN -C3 \
-o atmacs-LoRe-adm-3h.ts
# Explanations:
# -L - import local (L) and regional (R)
loading and their barometer
# -RMLe - add the first two columns together and store in
column 1.
# Do "early"
(before tsfedit)
# -RMA - time-dependent admittance,
# half-width
8, start 0, step 1 : (c1 * c3)/|c3|2 with a HANN
window
# Since there are 8 ordinates per day, the output rate is 3h.
# Removing DC-levels is intrinsic to the function
short_term_admit, sas/p/shrtadms.f
#
#
# Apply the gain factors on our barometer series, bar =
baro-atmacs-regional
setenv PRESSRC -LoRe-adm-3h
tslist d/g090701-130515-1h.mc -L'B|V' -BHc2010,2,2,11 -I
-Eatmacs.tse,GAIN -OEA \
-o d/bar100202-130515-1h.ts
# Explanations:
# The -OEA option keeps the epoch and fills with MRS's up
to the first sample.
# -Eatmacs.tse,GAIN:
TSF EDIT GAIN
REMDC
OPEN 22 < Atmacs/atmacs${PRESSRC}-adm-1d.ts
22, 'BIN]',-99999.0,'U', 0, 0, 'K*',
1.0
TDGAIN D
END
#
# A later tsfedit stage may replace the leader with zeroes.
# REMDC
# POKE 0.0 To 2010 02 02 00 05 00
#
# If our barometer is to supplement Atmacs,
# TRUNC At 2010
02 02 00 10 00
# REMDC
# EXTEND V=0.0 To yyyy mm dd hh mm ss fff
# ; i.e. to end-of-series
#
#
# Interpolate the far-field effects, bag = baro-atmacs-global
tslist Atmacs/atmac.mc -L'A|GMN' -L'A|GML'
-RMLe:1.,1,1.,2=1 -Eatmacs.tse,R1H -I \
-BHc2010,2,2,11
-U2013,05,16,19 -o d/bag100202-130515-1h.ts
# with
TSF EDIT R
RESAMPLE 3 0. ${RSMPDT:[600.]} [s]
END
# For comparison, interpolate the near-field effects
tslist Atmacs/atmac.mc -L'A|LMN' -L'A|RMN'
-RMLe:1.,1,1.,2=1 -Eatmacs.tse,R1H -I \
-BHc2010,2,2,11
-U2013,05,16,19 -o d/bari100202-130515-1h.ts
# In urtap-bigtides*.ins the Atmacs series have abbreviated file
names.
# Take gravity residual, no air pressure correction,
subtract polar motion:
tslist ~/TD/d/garb-600s.mc -L'G|R' -C3 -Bc2010,2,1 -H24
-Edeci.tse,DRD -I -o gr-6h.ts
tslist ~/TD/d/garb-600s.mc -L'B|V' -C3 -Bc2010,2,1 -H24
-Edeci.tse,DD -I -o b-6h.ts
tslist ~/TD/d/garb-600s.mc -L'B|V' -C3 -Bc2010,1,1 -H24
-Edeci.tse,D6 -I -o b-1h.ts
# Continue with polar motion reduction:
polmotm -l 11.9,57.6 -d 1.16 gr-6h.ts |\
tslist - -gI4,5i3,t66,f11.0 -k3 -C3 -I
-o pm-6h.ts
setenv FNSUB pm-6h.ts
tslist gr-6h.ts -E~/TD/subts.tse,SUBU -I -o grp-6h.ts
tslist ~/TD/d/garbp-600s.ts -C3 -Bc2010,2,1 -H24
-Edeci.tse,DD -I -o gbr-6h.ts
#
# Subtract atmac-sum.ts from grp:
setenv FNSUB atmac-sum.ts
tslist grp-6h.ts -E~/TD/subts.tse,ADDV -I -o gar-6h.ts
#
# Continue with reducing sea level
effects
# Get the Ringhals tide gauge file
if ( ! -r ~/SMHI/tgg-rin-tmp.ts ) then
cd ~/SMHI
tslist-app -I -o tgg-rin-tmp.ts +
SMHI_2010_JAN_NOV/slrin100101-1h.ts tgg-rin-100301.ts
cd ~/TD/Atmacs
endif
touch rin-1h.tse
urtap-outlit -pc urtap-rin.ins
#
# decimate and date-cut the residual
tslist rin-1h.ra.ts -C3 -Bc2010,2,1 -H24 -Edeci.tse,D6 -I
-o rin-6h.ts
# Lsq analysis (with outlier iteration):
urtap-outlit urtap.ins
urtap-outlit urtap-with-rin.ins
urtap-outlit urtap-baro.ins
# and plot
plot-atmac -w
#
# That's for manual import of SMHI-data:
#awk '/2010-02/||/2010-03-01/{sub(/,/,".",$0); print
substr($0,9,38)}' RINGHALS.RW.Hour.txt | tslist -
-gi4,1x,i2,1x,i2,2x,i2,t30,f9.0 -k1 -C1 -o ../tgg-rin-100201.ts
#
.bye
PS
OLD!
DDAS-DATA old methods
To copy:
RAW_o054 MON_o054
DAY_o054 to the respective directories on auda
(...\hgs\RAW_o054 etc.)
We don't distinguish years on auda.
On froste, copy the data on auda 1:1 to ~/TD/
.../hgs/RAW_o054
to ~/TD/RAW_o054
using scp ( ~/bin/scpa )
e.g.
cd ~/TD/RAW_o054
scpa 'RAW_o054/??1004??.054' .
bzip2 *.054
cd ZIP
scpa 'RAW_o054/??1004??.ZIP' .
cd ~/TD/MON_o054
scpa 'MON_o054/??1004??.TSF' .
bzip2 *.TSF
For a local backup, start cycgwin and connect the 'My Book'
(or another USB disk device) and find out the drive letter.
Assume E:
cd
/cygdrive/e/tidedata/GWR_o054/2009o054
copy the new files from froste to the USB device like above.
Using the DOS-Window
In order to prevent Windows Explorer to allocate but not free
workspace, try using copy commands in a DOS Window
cd
d:\GRW_o054\2011o054\RAW_o054
copy ??11MM??.054
\\Hgs-dator\Users\HGS\Tidedata\2011o054\RAW_o054
cd
..\MON_o054
copy GW11MM*.*
\\Hgs-dator\Users\HGS\Tidedata\2011o054\MON_o054
# Wind speed:
cd ~/wx/ws
wget
ftp://wx.oso.chalmers.se/pub/weather/2009/wind_speed/'ws20091[6-9]?.log'
#
# concatenate and save all as
binary
awk '{printf "%s %2.2d %8.4f
%6.2f\n", substr(FILENAME,3,4)" 0 ",substr(FILENAME,7,3),$1,$2}'
ws2009*.log |\
tslist - -a
-g'(i4,i2,i5,f9.4,f7.2),DH+' -k1 -C3 -rs600 -I -o ws0906.ts
# Weather, general, old:
wget
http://wx.oso.chalmers.se/wview/Archive/ARC-2010-03-27.txt
Result:
--Timestamp--- Temp Chill HIndex Humid Dewpt Wind HiWind WindDir Rain Barom Solar ET UV
-------------- ---- ----- ------ ----- ----- ------ ------ ------- ----- ----- ----- ----- ----
20100327 00:00 43.2 43.2 43.2 98 42.7 2 4 342 0.0 29.5 -100000 -100000.000 -100000.0
20100327 00:05 43.1 43.1 43.1 98 42.5 3 5 324 0.0 29.5 -100000 -100000.000 -100000.0
20100327 00:10 42.3 40.1 42.3 97 41.6 4 5 309 0.0 29.5 -100000 -100000.000 -100000.0
...
# Weather, general, new:
This is essentially part of OSO Rain data
cd ~/wx/ARC
http://wx.oso.chalmers.se/data/xweather/2010/
foreach d ( `fromto -f '%4i' 1 117` )
wget
http://wx.oso.chalmers.se/data/xweather/2010/weather_2010_$d.log
end
Check the README
file (http://wx.oso.chalmers.se/data/xweather/README)
Super-script:
source HOW.TO [COLLECT] 2010
(without COLLECT the
script executes a loop with wget for the new available days)
If you need wind speed in X and Y, follow these suggestions (yields
60s-data)
tslist OWSD2010.ts -L'WS|V'
-Ere.tse,R -C3 -F1p,e12.4 -O:`label WSPEED,VAL` OWSD2010-60s.mc -I
tslist OWSD2010.ts -L'WD|V'
-Ere.tse,R -C3 -F1p,e12.4 -O:`label WDIR,VAL` OWSD2010-60s.mc -I
tslist OWSD2010-60s.mc
-L'WS|V' -L'WD|V' -TCD:1,2]180. -F1p,2e12.4 -C3 -O1:`label
WSX,VAL`'|'2:`label WSY,VAL` OWSXY2010-60s.mc -I
or all in one:
tslist OWSD2010.ts -L'WS|V'
-L'WD|V' -E1,2:re.tse,R -TCD180 -O1:`label WSX,VAL`'|'2:`label
WSY,VAL` OWSXY2010-60s.mc -I
.bye