USAGE:

        polmotm  [-l lon,lat] [-f file] [-o ascii_output]
                [{-p lon,lat | -p= }] [ -s scale ]
                 [-a adde ] [-b addb] [ -d delta ] [ -g gammma ]
                 {-t date(3),hms(4) -r[s|d] n,rate | model-file.ts }

 Creates a time series of polar motion effects, gravity and/or potential.

 OUTPUT with gravity:

    polmotm -l 11.9,57.6 d/g090615-OPNEND-1h.mc
      YYYY MM DD HH MM SS MJD Dx Dy Dlat[deg] Dg[nm/s2]

   For conversion to BIN, pipe into fgrep -v '<' and
      | tslist - -gi4,5i3,t66,f10.0 -k3 -I -o out.ts

 OUTPUT with gravity and potential to ascii file:

    polmotm -l 11.9,57.6 -p= -a 144 -b 144 -d 1.160 -g 0.695 -o tmp/pmtest.dat d/g090615-OPNEND-1h.mc
<               date              >   DX    [mas]    DY   angle[o]  Dg nm/s2 local N [m]   N cos [m]   N sin [m]
2009 06 09 12 00 00 00 54991.500000   0.06306   0.54028   0.04971  -8.56905  1.6671E-03 -2.3373E-03 -2.0026E-02
2009 06 09 13 00 00 00 54991.541667   0.06316   0.54030   0.04961  -8.55198  1.6638E-03 -2.3412E-03 -2.0027E-02
2009 06 09 14 00 00 00 54991.583333   0.06327   0.54032   0.04951  -8.53498  1.6605E-03 -2.3450E-03 -2.0028E-02
...

N sin  and  N cos specify factors with which pole tide loading effects can be composed.
The factors can be extracted from ~/Oload/HAMTIDE-PREMY15-s3deg/olfg-eqlt-ugal.results

dg(t) = -125.07*(Ncos) - 3.04*(Nsin)
calc '-2.3373E-03*(-125.7)+(-2.0026E-02)*(-3.04)'
      = 0.35 nm/s^2
      4% of the direct effect. Ha ha

There's a maximum Dg:

2009 09 19 06 00 00 00 55093.250000   0.27544   0.36513  -0.19423  33.48352 -6.5143E-03 -1.0210E-02 -1.3534E-02
calc '-1.0210E-02*(-125.7) + (-1.3534E-02)*(-3.04)'
      = 1.32 nm/s^2
      4% again

 OPTIONS:
  -f  file name for input of a Julian-X-Y polar motion series.
      Default /home/hgs/TD/PM/pmjxy.dat (assigned in
      Ttide/p/polmots.f , generated by TD/PM/RENEW )

  -a  extend duration with adde samples after the end
  -b  extend duration with addb samples before the beginning

  -l  site location longitude, latitude
      If given, latitude change and gravity nm/s^2 are printed
      else only pole X Y is printed. Latitude change and X Y
      are in units of seconds of arc

  -p  Like -l, outputs equipotential surface however

  -p= will add the equipotential effect to the output, site location is copied from -l

  -d  delta    gravimetric factor, default = 1.0

  -g  gamma    Love number combination for potential,
                default=1.0; a resonable value is 0.695

  -s  scale    scale factor, default = 1.0

  -s/ scale    reciprocal... For output in volt, specify
                -s/ -777.4

  -t  yyyy-mm-dd_hh:mm:ss:fff  set epoch and start time
      OBS! no blanks.

  -r  n,sampling interval [h]. Default = 10,1.0
  -rs n,sampling interval [s]
  -rd n,sampling interval [d]
      n specifies amount of samples

  -o  ascii_output file. Default = STDOUT

  -q  Test mode, only protocol is produced

      file.ts:  a time series file for which a parallel polar motion series
      is computed. -a and -b can be used to extend output.