Program:
~hgs/gps/p/m/apl_reduce_stacov.f           (source code)
/Bifrost/hgs/gps/bin/apl_reduce_stacov.out  (executable)

Instruction file:
~hgs/Oload/smhi/aplj/apl_reduce_stacov.ins

This instruction file is only for a test; it references mapped stacov files. The real job should be done on unmapped stacov-files.

1994>
 &param
 rename='/netsol/','/netsol-apl/'
 &end
31 ^ &/all_apl_24h.xyz
   q
/Bifrost/gps/netsol/ITRF96_1994/94jan17_itrf_vf.stacov
/Bifrost/gps/netsol/ITRF96_1994/94jan19_itrf_vf.stacov
...etc....
$$$$

1995>
...etc....
 

Run:
set path = ( $path /Bifrost/hgs/gps/bin/ )
setenv path ( $path )
set PATH = $PATH:/Bifrost/hgs/gps/bin/
setenv PATH $PATH

# For Wettzell we need an alias:
setenv WETB_alias WTZR

cd somewhere
mkdir /Bifrost/gps/netsol-apl
mkdir /Bifrost/gps/netsol-apl/ITRF96_1994
uncompressdir /Bifrost/gps/netsol/ITRF96_1994/

apl_reduce_stacov.out ~hgs/Oload/smhi/aplj/ '>1994>'
$

That's what you see on the stdout:

G[hgs]~/Oload/smhi> apl_reduce_stacov.out aplj/
 <OpenIn>>> Your terminal: /dev/pty/ttyqe (I suppose)
 <OpenIn>>> Open Terminal = /dev/pty/ttyqe, unit was   2
 > *.ins file (apl_reduce_stacov.ins)> $

 <OpenIn>>>  4:  aplj/apl_reduce_stacov.ins           - on 00-Mar-10 16:47:51
 <OpenIn>>>  4:  aplj/apl_reduce_stacov.ins           - INSTRUCTION FILE
 <OPENF->>> 31:^ &/all_apl_24h.xyz                    -
 <SetEpo>>> 1993 08 01, JUL=49200
 <Main-->>>  IN: /Bifrost/gps/netsol/ITRF96_1994/94jan17_itrf_vf.stacov
 <Main-->>> OUT: /Bifrost/gps/netsol-apl/ITRF96_1994/94jan17_itrf_vf.stacov-apl
 <Main-->>> 1994 01 17 12 00 00 00, id=   170
 <Main-->>> site    1 ARJE
 <Main-->>> site    4 HASS
 <Main-->>> site    5 HERS
etc.
--

... id=   170   tells it's the 170'th day after the epoch of   all_apl_24h.xyz
If a site in the stacov file is not in the xyz-file, you have the opportunity to specify an alias, else the stacov record will not be changed. The alias specification is done via the shell environment. See above.

 <Main-->>> OUT:  Notice that we add  '-apl'  to the output file name.

The file   all_apl_24h.xyz  was made the following way:

cd ~hgs/Oload/smhi

(1) extend loading time series, appending  apls/SITE_apl.tsf

(2) make binary files

   foreach site ( `the-sites-in apls/ _apl.tsf` )
   tsf2ts apls/ $site
   end

(3) join time series in  apl/ and  apls/ :

   foreach site ( `the-sites-in apl/ _apl.ts` )
   ts-join97 apl/ apls/ aplj/ $site
   end

(4) make 24h-decimated binary files:

   foreach site ( `the-sites-in aplj/ _apl.ts` )
   24h-decimate -B aplj/ $site
   end

(5) make the big file:

   apl_ren2xyz.out aplj/apl_ren2xyz.ins
   $

 &param
 dir='aplj'
 &end
31 < &/all_apl_24h.xyz
   q
ARJE
BORA
BRUS
... (all the sites)

The file contains currently all air pressure loading displacements from 1993AUG03 until 2000JAN29. It's binary to keep its size low.