HOW TO  scale GPS rate errors after the fact:

look at ~/gps/errsca.env

First you need to copy sasm03.csh into $TAPDIR and eventually adjust. Here is a template
#!/bin/csh
set site=$1
set comp=$2
unsetenv term
unsetenv tty
sasm03 << END
TER:
21 ^ ./${site}.${comp}.wr.ts
   Q
 &param
  lvllen=-1, qshlvl=.false., pef_order=8
  rec_mrs=-99999.9, dat_mrs=-99999.9
  l_section=128, logxax=.true.
  target='BIN'
  units='[mm]\$'
  mx_miss=0, dff1=0.d0, ibeg_burg=24, qrfsp=.false.
  qrepair=.true., qgra=.false.
 &end
END

cd $TAPDIR
source ~/gps/errsca.env

Three files are generated for each component:
acv.$co.dat, chi.$co.dat and sca.$co.dat
The latter is read in by averotm.out and contains the scaling factors to obtain a-posteriori sigmas from the a-priori sigmas.

To create horizontal motion vectors for map drawing:

cd ~/gps/Upto_98/std
ss_rates_data .ra > rates.ra.dat
ss_rates_data .ea > rates.ea.dat
ss_rates_data .no > rates.no.dat
 

cat rates.ra.dat rates.ea.dat rates.no.dat >! actual_rates.dat

setenv WETT_alias WTZR

averotm.out '>STD>'
$

Here is a template for the averotm.ins STD section:
STD>
21 O actual_rates.dat
22 O ~/Maxwell/glenn_rates.dat
25 O $TAPDIR/sca.ea.dat
26 O $TAPDIR/sca.no.dat
31 B ~/gps/maps/bstdhor-vecs.dat
32 B ~/gps/maps/glennhor-vecs.dat
33 B bstd_adj_glenn.dat
42 B ~/gps/maps/bstdhor-vels.dat
43 B bstd_adj_glenn.rsl
   Q
 &param
 stacov_file='$REFF.stacov'
 legend_obs  ='12 0.0 0 1 Obs. rot''d # mm/yr '
 legend_model='12 0.0 0 1 Model # mm/yr '
 explain_obs  ='Rate # mm/yr ','Observed Bifrost'
 explain_model='Rate # mm/yr ','Model Milne (1999)'
 size=300., hordir=180., nmsg=0
 xlatl=40.0, xlonl=22., dlatl=.2, dlonl=+.5, sample_obs=2.
 xlatm=38.0, xlonm=22., dlatm=.2, dlonm=+.5, sample_model=2.
 q_incl_sites=.false.
 drop_sites='WETB WTZR SAAR MADR'
 scale_sigma_ea=2., scale_sigma_no=2.
 &end
 

To create vertical motion data,

fira -r -s > vertical_rates.dat

.bye