1 Prepare

   cd ~/gps/Upto_01/blt
   setenv TAPDIR ~/tap/Upto_01/blt
   setenv DDIR ~/gps/Upto_01/blt
   setenv SDIR /net_mnt/Bifrost/gps/RESULTS
   setenv SD /net_mnt/Bifrost/gps/RESULTS/ITRF_YY
   setenv REFF /net_mnt/Bifrost/gps/RESULTS/all_itrf00

2 Running blanalysis
(2.1) Making the blt.info table from the stacov file
   cd $DDIR
   blanalysis.out blanalysis.ins '>BI>'

(2.2) Doing all baseline collection
   blanalysis.out

3 Running all tap
   cd ~/tap
   setenv SITEEVENTS ~/tap/Upto_01/siteevents.dat

   foreach bl ( `the-sites-in -a $DDIR _bl.mc` )
   set sita=`echo $bl | sed 's/^\(....\).*/\1/'`
   set sitb=`echo $bl | sed 's/^....\(....\).*/\1/'`
   urtap-bl L $sita $sitb -poi >! urtap-bl.prt
   urtap-bl T $sita $sitb -poi >! urtap-bt.prt
   echo $sita $sitb
   end

4 Extracting the rate information
   cd $TAPDIR
   rm blt.rsl
   touch blt.rsl
   foreach bl ( `the-sites-in -a ./ .bl.prl` )
   set sita=`echo $bl | sed 's/\(....\).*/\1/'`
   set sitb=`echo $bl | sed 's/....\(....\).*/\1/'`
   echo $sita $sitb
   fgrep '/[year]' ${sita}${sitb}.b?.prl >> blt.rsl
   end

5 Getting the observed rates, converting to mm/yr, and scaling errors up with 2

   awk '/\.bl\./{sub("WETB","WETT",$1);sub(".*/","",$1);sub("\.b.\.prl:","",$1); \
       printf "%10.3f%10.3f %s\n", 1000*$6,2000*$8,$1}' blt.rsl >! bl_obs_sca2.rsl

   awk '/\.bt\./{sub("WETB","WETT",$1);sub(".*/","",$1);sub("\.b.\.prl:","",$1); \
       printf "%10.3f%10.3f %s\n", 1000*$6,2000*$8,$1}' blt.rsl >! bt_obs_sca2.rsl

Tasks 5 and 6 are combined in the tcsh-source script
blt_collect_rates.env

6 Computing model rates
   cd ~/gps
   setenv GAMMODEL 120p810
   mblanalysis.out

Plotting the results. These plots need csh and awk scripts in ~/gps/bin
   cd $DDIR
   site-clr-key-plot `the-sites-of SWEPOS FINNNET` WTZR
   bt-rr-plot
   btr-rr-plot