run_urtip

USAGE:
          run_urtip [options] [ts-modelfile]
 
PURPOSE:

 
          OBS!!!   Presently this script runs urtipgt by default.
                   Use  -C urtipg to change to the Bonn-version.
 
          Compute tide predictions from a table, produce
          a time-series file exactly synchronous with
          the specified ts-model file (t0, dt, epoch, length)
 
          The table is prepared by urtap (.trs, unit 13),
          default file name t/urtap.trs
 
          The output is found in t/pt.tsf by default.
          For converting ASCII -> BIN, use
             tslist o/pt.tsf -gi4,5i3,t33,d18.0 -k3 -r#rate
          (the -r#rate option might be needed if the
          auto-detect rate is not precise enough).

          The timing parameters are derived from the model file.
          If no model file is specified, a source-script file
          must be prepared. See the section on FILES below.

          If you modify a trs-table (see -t ), observe the note.
 
OPTIONS:

 -e          - test run, the instruction input to urtip is shown.
 
 -C prog     - the program name for predicting                         [urtipgt]
 
 -E src      - If no ts-file is specified, source src file to
               set the runtime parameters.                         [./urtip.env]
               A drift signal can be produced with prepared code,
               by adding a line
               set addrift=-Erun_urtip,EXP
               set addy=1.0
              
(if the drift signal is to be added with factor 1.0   [0.0])
 
 -t table    - the table file                                      [t/urtap.trs]
 
 -o output   - the output file (ASCII)                                [o/pt.tsf]
                "    "      "  (BIN) accordingly:              
                                     extension .ts                     [o/pt.ts]
 
 -PM         - add polar motion                                          [don't]
 -LL         - longitude,latitude for PM calculation                 [11.9,57.6]
 
 -b [output] - produce ASCII and BIN file                           [ASCII only]
 -B [output] - Retain only BIN file                                    [o/pt.ts]              
 -bm         - like b, produce GRAV|PRED column
               in MC-file, however
 -BM         - like -bm, retain only BIN file, however
 
-rmc        - remove the MC-file first                                  [don't]

 -s scale    - rescale the output                                        [1.0d0]
 -s u        - unscale the output (1/calfactor from the table file)
 -u          - unscale the output (1/calfactor ...))
 -dudt       - compute the time-derivative of the effect

 -prf file,dt[,trg]                                              [no pre-filter]
             - apply a pre-filter (the filtering that was involved in
               the model file, for instance).
               file - produced with a FILTER command and OUTFS in tsfedit
               dt   - the input sampling interval [h]
               trg  - the label of the filter section if the file
                      contains more than one. [' ']

 -besselprf n[,f0]
                                            [no Bessel filter]
             - apply the Bessel-filter that was used on the analogue
               side of the data acquisition, GWR: 8,0.0174
               n    - Filter order
               f0   - Scaling frequency in Hertz.


ENVIRONMENT:
 
   RESIDEDIT - an option (or a series of options) for tslist when inspecting
               the model file. Typically:  setenv RESIDEDIT "-Efile.tse,target".
               run_urtip spawns for epoch, time step, and length

FILES:

  ./urtip.ins - output, a copy of the run-time instructions read by urtipgt. 

  -E option, source file coded in csh, example
set date=2012,10,25  
set hms=0,0,0  
set dt=1.     
set tscale=60.
set nout=1440
   The time step in hours used by urtipg[t] is:  dt/tscale; thus, the example
   generates 1-minute data.

   Here is the list of namelist parameters that can be (must) set with the source file:
 &param
 date=$date, hms=$hms
 dt=$dt, tscale=$tscale
 nout=$nout
 iun=31

 iun_prf=$iunprf, dt_prf=$dtprf, trg_prf='$trgprf'
 n_bessel_prf=$nbessel, f_bessel_prf=$fbessel
 qdudt=$dudt
 $locstr
 &end
 
  -PM option, Polar Motion:
    This script will run the polar motion program and add the time series

A NOTE ON trs-FILES:
Modifying a tide table (trs-file), it's important that the numbers right next left of the `;ยด-sign point to the tide harmonic table's leading wave of the band.
The latter is obtained after the removal of the static tide (a.k.a. M0S0 and its degree-4 companion). Thus, the lunar nodal tide refers to number 1. Here is
an example that generates the tide band with frequencies below 1 cyc/yr:

(2,0,0,0,*)                      =     Nod2    1 ;  -3.56995845D+03  0.0
(3,0,0,0,*)                      =     Nod3    5 ;  -5.08756D+03     0.0
END WAVEGROUPS
   11.926000   57.396400    7.500000    0 = site long/lat/alt & IALTG
    9.798529    6378.137    6363.006      = g-used [m/s^2] r-used r-geocent [km]
   72     = n_eigenvalues
    0     = whitening filter length
S_GRAVIT  = core model
  -774.421 [nm/s^2]$        = CALFACTOR
 


EXAMPLES:
   Producing using nodal tides ("gnt...") for subtracting later from observations  d/g${begd}-OPNEND-1h.mc 
    set begd=100202
    run_urtip -t t/nodal-only.trs -rmc -BM d/gnt${begd}-OPNEND-1h.ts d/g${begd}-OPNEND-1h.mc
    setenv FNSUB
d/gnt${begd}-OPNEND-1h.ts
    setenv FACSUB 774.47
    tslist 
d/g${begd}-OPNEND-1h.mc -LG+V -I -Esubts.tse,ADD \
           -O:`label GRAV,NTSUB`
d/g${begd}-OPNEND-1h.mc

  Producing tides at 1-minute rate for a day:
    prtide4day -D 2013,01,02 -rm 1. -n 1440 > file.env
    run_urtip -A PT130102-60s.tsf -E file.env -u 

 
A rather standard type of task: 1-s
    run_urtip -t tt/g090615-OPNEND-1h-nnt-ochy.trs -b tmp/G1_p_150204-1s.ts lpfu/G1_g_150204-1s.ts