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.
______________________________________________________
¶m
<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 ¶m 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://barre.oso.chalmers.se/hgs/ttimm/ttimm4kl.ins
another one:
TK>
¶m
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.
¶m
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>
¶m
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 ¶m, 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