1. Selecting files by weekday
2. RAW data |
   
3. 1-min gravity, baro, predicted tides -> MC |
4. 1-s gravity -> BIN ts
5. 30-s MC ASCII from G1 |
6. 10-min data from G1  |  10-min data from e.g. garb MC |
   
7. Empirical residual
   
8. Polar motion data:  |wget IERS data; extend eopc04_IAU2000.orig;  compute gravity
9. Subtract polar motion: | see item 7.
   
10. Drift determination: |  tide analysis, air pressure, polar motion reduction
11. Tide analysis |
12. Gravimeter: sliding residual RMS |
13. Tilt controllers: | sliding RMS and long-term drift; and a power spectrum.
• Tilt controllers and the Maule earthquake
  
14. Si Diode thermometers | at the coldhead ("neck")
15. Plot | Plot many curves |
16. GGS | 1-s data for mode detection and ftp to GGP
17. Free Oscillations |
18. Daily residuals, micro-seismic spectrum, shocks, and cleaned low-rate data |
19. Daily RMS |
20. Daily PSP |
    
21. A sequence of periodograms | (every hour for 1 day, 1-s data) is an example in the tslist man page
• ... see also below on this page
22. Determine the daily "human" signature (working days versus holidays)
23. A diagnostic plot with GMT |

# ... get 1-min gravity, baro, and predicted tides (from GGP-files) into one MC-file
#

cd ~/TD
rm -f d/MC0906-1m.ts
cat MON_o054/GW09??00.GGP |\
 awk '/^2009/{printf "%s %s %11.7f\n",$1,$2,$3}' |\
   tslist - -gi4,2i2,i3,i2,t16,f12.0 -k2 -Bc2009,06,13 -I -O:`label GRAV,VAL` d/MC0906-1m.ts
cat MON_o054/GW09??00.GGP |\
 awk '/^2009/{printf "%s %s %11.4f\n",$1,$2,$5}' |\
   tslist - -gi4,2i2,i3,i2,t16,f12.0 -k2 -Bc2009,06,13 -I -O:`label BARO,VAL` d/MC0906-1m.ts
run_urtip d/MC0906-1m.ts
tslist o/pt.ts -I -O:`label GRAV,PRED` d/MC0906-1m.ts

#
# See script  myfo  for a recipe how to generate de-tided and de-baroed gravity
#
# See the new script resids e.g.
#

resids -xp -d 10 2009 06 15 2009 10 31

# creates 10-min data d/MC0906-10m.ts and a plot in plot/MC0906-10m.ps

foreach day ( 8 9 10 )
tsf2ts-decim -r $day,$day -W -c 5 -d 1 -O GRAV,VAL d RAW_o054/ A2 0909
end

#
#or

sasm03-tilt.ins
21 ^ ${sasm03_input}
41 < sasm03.pef
42 B ${sasm03_output}
   Q
 &param
  lvllen=-1, qshlvl=.false.
  rec_mrs=-99999.9, dat_mrs=-99999.9
  l_section=8192, logxax=.true.
  target='BIN', fmt='U'
  units='[mV]$'
  mx_miss=10000, dff1=-.999d0, ibeg_burg=24, qrfsp=.false.
  qrepair=.true.
  lpef=21,pef_order=19, lpef_wh=-1, iun_wh=-1, qrpefsp=.false.
  pef_show=11, pef_order_show=20,15,12,9,7,6,5,4,3,2,1
  t_window='KAIS', p_window=2.2
  qgra=.true.
  cutlvl=1.d-8
  scale=1.d3
  q_tsf_edit=.true., tsf_edit_name='SHOCKS'
  iuspout=42
 &end
TSF EDIT SHOCKS
CONT 43 R tilt/shocks.tse
END
#
#  Result:



The tilt controllers went over-range  in the surface waave action after Maule, Feb. 27, 2010. The following recipe ~accomplished a restoration of the signal:

setenv RESIDEDIT -Etilt-unclip.tse,U
daily-resid -x 11,19,12,20 -L TX-PWR,VAL -L TY-PWR,VAL -L TX-BAL,VAL -L TY-BAL,VAL -NT -T tbxy RAW_o054/A2100227.054

tilt-unclip.tse:
TSF EDIT U
DEL > 9.9
UNCLIP T=1.d-5 I=1000 S=4 O=P+
END

Result:


The tick interval on the time axis is 1 minute. The "After" and the "Interpolated" curves have been moved by -2 and +2 Volt, respectively.
For more on the subject of unclipping, see "Displacements from gravity"  and "Signal Restoration".

For the plot, use the following two lines in tilt-unclip.tse to prepare output:
OUTTS   U=51 < o/Chile-TXPWR-unclip-1s.ts
OUTTS W U=51 < o/Chile-TXPWR-uncpre-1s.ts
and the plot-init file plot/tilt-unclip.init

Knowing the X- and Y-signals, we would like to obtain the East and North motion/acceleration or the motion in the direction of wave propagation and at right-angle.  For this purpose there is a routine that write a tse-file with the appropriate parameters
~/TD/tiltrot
tiltrot -a 65 -c PWR -o d/A2_tiq_

GGS data

Daily residuals, micro-seismic spectrum, shocks, and cleaned low-rate data

has been moved