Running ttimm  |  Namelist runtime parameters   Opening files - coding instructions  | Update notes  |  Examples  |  DUT - The UTC file

     Program  ttimm


 computes time series of Body tides (BT) and/or Ocean Loading tides (OLT).
 Can compute
             displacement tides
             strain tides
             tilt tides
             gravity tides (rigid earth, elastic earth with liquid core)
             inertial accelerations due to displacements
             Coriolis       "               "    "           "

Call this program with command line

     ttimm @ <ins-file> [:LOC]

 where LOC designates a location target in the instruction file referring to a section:
 LOC>
 
 Upon starting ttimm without any argument will cause the program to prompt for an instruction file name.
 Entering a `$ยด at this prompter means the default file  ./ttimm.ins  and the first section of instructions.

 If  nsites=1,  ttimm computes tide displacement time series  
 (alt. other components) at a location along components Radial East North (REN)
 (also UEN, means the same as REN).

 If nsites=2, ttimm computes baseline changes along geocentric XYZ components.
 or, alt., baseline LTU components

 ______________________________________________________

Instruction file
['ttimms.ins']:

 Consists of one NAMELIST and one OPEN FILE BLOCK.
  See examples.
 ______________________________________________________
  &param
  <namelist parameter=value expressions>
  &end
    * <Text>                                            
Openf-block, ignored line
    C <Text>                                            
Line printed to STDOUT
 iu o <file-name> [<comment(iu)>]    <- 21 B result.tsf 
Output time series
    Q                                                   
End block
 ______________________________________________________

 
About files to open

Namelist &param  parameters
 Not all parameters are documented here. See the complete list.
 ___________________________________________________________________________

 Name      - type    - meaning                                   [defaults]

 ___________________________________________________________________________

 tpath     - char64  - Path/filename of the harmonic tide development        !
                                                      ['etgtab/etcpot.dat '] !

 utc_file  - char64  - Ascii-file with dates and leap-seconds
                                                     ['./etgtab/ttutcf.dat']

 idate(3)  - integer - year, month, day  (e.g. idate=1993,12,6 
                       for Dec 6, 1993)                         [1994,01,01]

 Defining the site(s):
 nsites    - integer  - 1 or 2 for single site or baseline mode,
                        respectively.                                    [1]
 xsite(2)
 ysite(2)
 zsite(2)  - real*8  - coordinates for sites 1 and 2
 q_xyz     - logical - .true.  -  coordinates are geocentric XYZ in meters.
                                                                   [.false.]
 q_lola    - logical - .true.  -  coordinates are longitude latitude height
                                  in degrees                       [.false.]

 Defining the component:
 q_se_tide - logical - compute solid earth tide.                    [.true.]
 q_ol_tide - logical - compute ocean loading tide.                  [.true.]
 q_gravity - logical - gravity tides elast earth w core            [.false.]
 q_re_grav - logical - gravity tides rigid earth                   [.false.]
 q_re_pot  - logical - rigid-earth tide potential/g                [.false.]
 q_strain  - logical - strain components                           [.false.]
 q_tilt    - logical - tilt components                             [.false.]
 q_permanent_tide    - logical -                                   [.false.]
 q_acc     - logical - compute linear accelerations                [.false.]
 q_corio   - logical - compute coriolis accelerations              [.false.]

 q_acc and q_corio cannot be combined with q_gravity, q_re_grav or q_strain
 

 The sites can be taken from the ocean loading table:
 snm_opt(2)- char*4  - option for finding site by name (call ocean_load_site):
                       'M' - monument code
                       'S' - site name
                       'N' - site name AND monu code
                       'U' - name comparison in upper case chars
                       'R' - one rewind of file per call.
                       ' ' - don't call ocean_load_site
                             (one of q_lola and q_xyz must be set .true.
                             and the coordinates must be given explicitly.)
 site_name(2) - char*16 -
 monu_code(2) - char*6  -
                      
 pmt_opt   - char*16 - options for print_major_tides
                      '+H' - show relative phase w.r.t. astronomical potential
                      '+P' - like +H, but evaluate in extra column
                      '-L' - with +P, don't care about longitude
                       See p/prtides.f for more (bizarre) options.

 Ocean loading options:
 ol_opt    - char*4  - option parameter for subroutine Ocean_Load_tide,
                       composed of
                       'I' - interpolate spectrum
                       'P' - extended print:  loading parameters
                       'R' - even more print: response parameters
                       'Q' - quiet, 'q' - nearly quiet
                       '.' - use default.
                       Opposed to 'I'-option, the sparse loading spectrum
                       is returned, consisting only of the tabelled data.
                                                                    ['I...']

 The output time series:
 nt        - integer  - requested length of time series output.        [300]
 dt        - real*8   - time interval in seconds.                  [3600.d0]
 t0        - real*8   - time from midnight in seconds                [0.0d0]

 ncomp_out - integer  - 1 or 3
                        For BLQ gravity with tilt, if you don't
                        want tilt, specify 1                   [BLQ decides]

 out       - char*3   -          default =                           ['UEN']
                        'LTU'  - baseline tides:  Length Transverse Up
                        'XYZ'  - baseline or single site, geocentric tides.

                        For gravity or potential 'UEN' must be used. 
                        'UEN'  - single site tides:  Up East North
                                 'REN' can be used as a synonym.

                        'UEN'  - baseline tides: Up East North;
                                 you can select the reference site:
 n_refsite - integer  - 1 or 2 - the site where the unit vectors UEN
                                 are computed
                        other  - the point above the baseline center
                                 is used.                                [0]

 fmtout    - char*64  - fortran format code  
                                      ['(i6,f6.5,1x,i4,2i3.2,#i3.2,3f10.4)']
                        Inspect the default:
                          i6,f6.5 is for a fractional Julian date
                          i4,2i3.2 is for  year, month, day
                          # will become the number khms
                          the remainder is for the ordinate values.

 khms      - integer  - indicates whether the hour(1) minute(2) second(3)
                        part is to be written                            [1]

 sco(3)    - real*8   - scale output                                 [1.0d0]
____________________________________________________________________________

 Files to open (subroutine openf)

 11 - input, ocean loading coefficients, BLQ-format

 21 - output, ASCII, time series, optional:
      Time [h]  radial east north.

 22 - output, ASCII, optional
      horizontal vs vertical.

 31 - input, optional
      A previous output file, in order to produce differences.

 32 - input, optional
      First component (radial), added to this procedure's output
 33,34 second and third component


The UTC-file with DUT's:
...
1999 01 01 +1
2006 01 01 +1
2009 01 01 +1
2012 06 30 +1
2999 12 31 00
-1
must be maintained manually.


Update notes:

 2014-06-18  Updated: Tilt added. Not tested yet.
             If ocean loading is added, it's BLQ rows 2 & 3
             that are added to BT columns 2 & 3, and 1 remains empty.
             OBS! signs in BT, signs in OLT not checked yet.
             OBS! Are the functions PNM DPNM o.k. for n=1..6 ???


EXAMPLES:
1. Loading effects, displacement - a collection of segments
 http://holt.oso.chalmers.se/hgs/ttimm/ttimm4kl.ins

another one:

TK>
 &param
 idate=2009,06,17     
 nsites=1
 dt=300.d0, t0=0.d0                        
300 s interval, start at 00:00:00 on idate
 nt=600             
                        600 samples
 khms=3              
                       Output time tag: hours, minutes and seconds
 site_name(1)='3JADE]'     
                 Locate this site in the BLQ-file
 snm_opt(1)='SRF'          
                 Locate by site name, one rewind of BLQ-file,
                           
                   fail if the site cannot be found.
 ol_opt='IPR.'        
                      Interpolate the spectrum to (almost) all
                       
                       degree-2 tides, diagnostic Print,
                         
                     include the Response parameters
 Q_ol_tide=.true.           
                Ocean loading tide: yes
 Q_se_tide=.false.           
               Solid earth tide: no
 out='REN'                    
              Output in Radial (=vertical), East, North
 tpath='/home/hgs/Ttide/etgtab/etcpot.dat '
 fmtout='(i6,f6.5,1x,i4,2i3.2,#i3.2,3f10.4)'
 sco=1000.,1000.,1000.            
          Output scaling to millimetres
 &end
11 R tk.blq                       
          Input
21 B 3JADE.ren                    
          Output
   Q


... with the BLQ-file tk.blq:
$$ Ocean loading displacement
...
$$ END HEADER
$$
  3JADE             
$$ GOT00.2_PP ID: 2010-09-30 20:02:37
$$ Computed by OLCMC/OLMPP by H G Scherneck, Onsala Space Observatory, 2010
$$ 3JADE,                               lon/lat:    7.7170   53.8790    0.00
  .00423 .00256 .00051 .00064 .00181 .00160 .00059 .00041 .00051 .00021 .00035
  .00299 .00065 .00061 .00016 .00213 .00126 .00069 .00027 .00006 .00005 .00001
  .00290 .00051 .00055 .00019 .00358 .00210 .00117 .00052 .00040 .00016 .00029
   157.7 -104.4 -167.1  -89.6  155.7  106.1  151.8  118.3   45.0   34.3    5.4
   148.0 -156.1  113.4 -163.3   72.0   47.5   70.6   19.1  -79.6  178.1   14.9
   108.7  160.0   83.8  105.1 -102.2 -112.0 -102.9 -107.3   59.3   31.8    5.6
$$

2. Body and ocean loading tide, gravity:
  In this example, the station coordinates are retrieved from a BLQ-file.
 &param
 idate=2017,12,1
 nsites=1
 dt=3600.
 nt=2000
 khms=3
 q_lola=.true.
 Q_ol_tide=.true.
 Q_se_tide=.true.
 Q_re_grav=.false.
 Q_gravity=.true.
 q_permanent_tide=.false.
 site_name(1)='TSUK]'
 snm_opt(1)='SRF'
 ol_opt='IPR.'
 out='REN'
 sco=10.d0,1.d0,1.d0
 myunits='[nm/s^2] [nrad] [nrad]'
 fmtout='(i6,f7.6,1x,i4,2i3.2,#i3.2,3f10.1)'
 tpath='/home/hgs/Ttide/etgtab/etcpot.dat '
 &end
   C This example computes solid earth tide displacements for Onsala
   C 1000 hourly values starting on May 19, 2011.
11 R Tsukuba-g.blq
21 B Tsukuba_g_171201.tsf                              The output time series
   Q


3. Loading tide, gravity: ttimm-ohig.ins

ttimm @
ttimm-ohig.ins
ttimm @ ttimm-ohig.ins :OHIG_TPXO72N  # is not necessary since the section is first.

OHIG_TPXO72N>
 &param
 msg_level=0
 idate=2015,06,28
 nsites=1
 dt=3600.d0, t0=0.d0
 nt=144
 site_name(1)='OHIGGINS_SHIFT', q_lola=.true.,q_xyz=.false.
 snm_opt(1)='S'
 ol_opt='IPR.'
 Q_ol_tide=.true.
 Q_se_tide=.false.
 q_re_grav=.false.
 q_gravity=.true.
 q_inject_leaps=.true.
 out='UEN'
 tpath='etgtab/etcpot.dat '
 fmtout='(i6,f9.8,1x,i4,2i3.2,#i3.2,3f10.4)', sco=1.d0,1.d0,1.d0
 myunits='[nm/s**2]'
 khms=3
 ncomp_out=1
 &end
11 O ~/Oload/BLQ/ohig-holt-20121029.blq   GRAV ocean loading coefficients
21 B ohig_tpxo72n_150628_i_1h.tsf         The output time series
   Q


with BLQ-file
ohig-holt-20121029.blq
$$ Ocean loading gravity+tilt
$$
$$ Calculated on holt using olfg/olmpp of H.-G. Scherneck
$$
$$ COLUMN ORDER:  M2  S2  N2  K2  K1  O1  P1  Q1  MF  MM SSA
$$
$$ ROW ORDER:
$$ AMPLITUDES (um/s**2 / 0.1 urad / 0.1 urad )
$$   RADIAL
$$   TILT E
$$   TILT N
$$  PHASES (degrees)
$$   RADIAL
$$   TILT E
$$   TILT N
$$
$$ Gravity is defined positive upwards.
$$ The phase lag is relative to Greenwich and lags positive. The
$$ Gutenberg-Bullen Greens function is used. In the ocean tide model the
$$ deficit of tidal water mass has been corrected by subtracting a uniform
$$ layer of water with a certain phase lag globally.
$$
$$ Complete <model name> : No interpolation of ocean model was necessary
$$ <model name>_PP       : Ocean model has been interpolated near the station
$$                         (PP = Post-Processing)
$$
$$ Ocean tide model: TPXO.7.2, Ssa tide from TPXO.6.2
$$
$$ END HEADER
$$
  OHIGGINS_SHIFT
$$ TPXO.7.2_PP ID: 2012-10-29 16:09:46 SCF=1.0E+01
$$ Computed by OLMPP by H G Scherneck, Onsala Space Observatory, 2012
$$ OHIGGINS_SHIFT,            GRAV TILT  lon/lat:  302.0992  -63.3230    8.000
  .06152 .02700 .00940 .00798 .03946 .04193 .01291 .00957 .00488 .00272 .00218
  .29225 .27132 .05909 .07781 .11999 .10878 .02909 .02613 .00336 .00159 .00151
  .22680 .20905 .02911 .06397 .15256 .10874 .03003 .01924 .00684 .00353 .00379
    84.1  121.4   56.6  122.6 -116.3 -134.8 -114.0 -147.2   14.3   11.0    2.0
    85.2  101.3   93.6  104.7  159.7  175.9  179.9  170.4   77.5   46.0    9.8
  -134.7  -94.2 -151.8  -88.3    9.2   -6.9   19.8  -10.9 -139.0 -160.4 -175.4
$$
  OHIGGINS_VLBI
$$ TPXO.7.2_PP ID: 2012-10-29 16:10:17 SCF=1.0E+01
$$ Computed by OLMPP by H G Scherneck, Onsala Space Observatory, 2012
$$ OHIGGINS_VLBI,             GRAV TILT  lon/lat:  302.0992  -63.3211    8.000
  .06082 .02661 .00932 .00786 .03898 .04143 .01275 .00946 .00484 .00270 .00217
  .29211 .27121 .05907 .07778 .11996 .10871 .02907 .02612 .00336 .00158 .00151
  .22639 .20876 .02904 .06387 .15238 .10858 .02999 .01920 .00683 .00352 .00379
    83.9  120.9   56.4  122.0 -116.3 -134.8 -113.9 -147.2   14.2   11.0    2.0
    85.2  101.3   93.6  104.7  159.7  175.9  179.8  170.4   77.6   46.1    9.8
  -134.7  -94.2 -151.8  -88.3    9.1   -6.9   19.8  -10.9 -138.9 -160.4 -175.4
$$ END TABLE



Namelist &param, complete set of parameters:

idate
xsite
ysite
zsite
xsiteg
ysiteg

q_geogr
q_xyz
q_lola
myunits

nsites
nt
crtdis
out
fmt32
a1
a2

pmt_opt
ol_opt
n_refsite
ncomp_out

t0
dt
q_inject_leaps

site_name
monu_code
snm_opt
tpath
eff_ave

delete
select
ndel
nsel
q_abs_tide
utc_file

q_ol_tide
q_se_tide
q_daily_ave
q_gravity
q_re_grav
q_strain
q_acc
q_corio
q_tilt
q_re_pot
q_permanent_tide

fmtout
khms
sco
n_print
msg_level

iun_spectrum
qraw_amp
amp_zp21
pha_zp21
btsndr

nstop


.bye