cd ~/gps
source gdh2000.env

For yy = 93, 94, ..., 00
copy  SITE.rad SITE.lon SITE.lat to `SD yy`  by ftp

cd `SD yy`/..
collect_many_NEV `ls -d HEFLIN_9?` HEFLIN_00

cd ~/IGS
posvel2ssc.pl heflin_00_00.posvel > heflin_00_00.ssc

(2) Some gipsy nev files contained overrange values (***)

(3) New procedure ~/gps/bin/gd2ren-0mh

(4)   cat $SDIR/WETB.nev $SDIR/WTZR.nev > $SDIR/WETT.nev

      Some more stations that can be concatenated:

cd $DDIR
mv TROM.gd.ep TROO.gd.ep
cat TROO.gd.ep TRO1.gd.ep > TROM.gd.ep

 
cd ~/gps
many-gd2ren -0mh `the-sites-of SWEPOS FINNNET IGS` WETT

cd ~/tap
many-urtap-ss `the-sites-of SWEPOS FINNNET IGS` WETT

cd $TAPDIR
source ~/gps/errsca.env

cd $RDIR
ss_rates_data .ra -v \'HOFN REYK MATE\' > ! rates.rsl
ss_rates_data .ea -v \'HOFN REYK MATE\' >> rates.rsl
ss_rates_data .no -v \'HOFN REYK MATE\' >> rates.rsl

cd $RDIR
averotm.out

ss_rates_data .ra > ! ra_rates.rsl
awk '{printf "%.4s %s\n", $1,$2}' ra_rates.rsl | xargs -l1 lola-of-site > vertical_rates.xyz

cd maps
mindist.out -X0,40,1 -Y50,75,1 -F../vertical_rates.xyz | xyz2grd -R4/36/52/72 -I1/1 -G../distance.grd

~/gps/maps/clr_vrates ../vertical_rates.grd ../distance.grd bifrost

source bstdghes-vels-1.env

~/gps/define-apl.env
~/gps/define-tgg.env

Model versus Bifrost plots for the three components:

~/gps/bin/gps_vs_model_plot -o gps_vs_model.ps -r bstd_adj_glenn.rsl
~/gps/bin/gps_vs_model_plot -o gps_vs_model.ps -e bstd_adj_glenn.rsl
~/gps/bin/gps_vs_model_plot -o gps_vs_model.ps -n bstd_adj_glenn.rsl

EOF
  prepare (jobs without clim, but with air press load)
      source ~/gps/define-apl.env
      setenv TAPVARIANT /NoClim
      setenv KINDEVENTS all
      urtap-ss LEKS -poi R
      setenv BEGIN_DATE 1993,10,10
      collect-resid `cat the_EOF_sites_north `
      unseat BEGIN_DATE
      collect-resid `cat the_EOF_sites_south `
      cd $TAPDIR$TAPVARIANT
      source errsca.env
      setenv RDIR ~/gps/Upto_00/NoClim/with_apl
      cd $RDIR
      ss_rates_data -d $TAPDIR$TAPVARIANT .ra > ! rates.rsl
      rates_unc_sca -H -r rates.rsl > ! ra.rates_unc_sca.tbl
      setenv TDIR $TAPDIR/NCA
      setenv PSDIR /Bifrost/hgs/gps/Upto_00/vfh00-93/NoClim/with_apl
      setenv GDPLOT_DO R
      cd ~/gps
      gdplot-ts LEKS

 EOF proper
      setenv TAPVARIANT /NCA
      setenv TAPPPDIR /EOF
      setenv RDIR ~/gps/Upto_00/EOF
      setenv PSDIR $RDIR/PS
      setenv MARGIN 0.4
      setenv APLDIR /home/hgs/Oload/smhi/aplj
      setenv TAPAPL
      urtap_eof_north LEKS -pu
      cd $TAPDIR$TAPPPDIR
      rm -f acv.*.dat sca.*.dat chi.*.dat
      source errsca.env
      cd $RDIR
      setenv TDIR $TAPDIR$TAPPPDIR
      ss_rates_data -d $TDIR .apl.ra >! ra.rates.dat
      rates_unc_sca -H -r ra.rates.dat > ! ra.rates_unc_sca.tbl
      setenv TADDR .apl
      gdplot-ts LEKS
 

 prep horizontal (jobs without clim, no air press load)
      unseat TAPAPL APLDIR
      setenv TAPVARIANT /NoClim
      urtap_ss LEKS -pu E N
 

      setenv TAPPPJUMP
      setenv TAPPPNET N
      urtap-eof-horiz -pu N LEKS
      urtap-eof-horiz -pu N LEKS

      setenv TADD
      setenv TADDR .apl
      setenv TDIR $TAPDIR$TAPPPDIR
      cd gps
      gdplot-ts LEKS
 
 

BOBW (Best of both worlds (Lulu))
      cd /Bifrost/hgs/gps/Upto_00/bobw/
      source HOW.TO
 
 

PLOT RATES vs RATES
      setenv GEOIDRATES ~/Maxwell/prof/Bifrost-geoid.xyz
      setenv RATESTBL ra.rates_vs_drtgg.dat
      rates_unc_sca -r rates.rsl > ra.rates_unc_sca.tbl
      awk -f ~/gps/bin/scarates.awk ra.rates_unc_sca.tbl > sca_ra_rates.dat
      rates_vs_rates_data -sf sca_ra_rates.dat
      rates_vs_rates_plot -s $RATESTBL

      averotm.out
      sshrrv.out
      awk -v l=0.375 -v s=5 -f glennver.awk ~/Maxwell/GAM/rates_120120.dat > maps/glennver.vecs.dat
 

.bye

.bye