| Back to HOW_TO
|
Back
to General
Information | In OTEQ/TTEQ | Post-proc of tide harmonic field
|
Post
processing of a tide harmonic solution (plane grid,
result
of OTEQ and TTEQ)
A
new project (Aug. 2014): Using GoogleEarth to create a local
hi-resolution array for the computation of Newtonian attraction at a
station at altitude.
Intro remarks:
Main program otemt1 , e.g., may have produced
a tide solution, e.g.
M2v5.DAT. This is usually a packed CMPX data set.
Use otem64 to compute loading effects.
Example MDL22h/otem64.ins
¶m
target='S'
comp='RADI'
nf=512
tide='..',
nver=0,
qpacked=.true.
&end
31 ^
&/FLZU.DAT
32 ^
&/${TIDE}${VCODE}.dat
33 ^
&/sg00egbc.rad
34 ^
&/OOTIDES.DAT
41 <
&/${TIDE}${VCODE}_rad.dat
q
Obs! Need to define environment parameters TIDE and VCODE!
Components (taken from
PROG/s/greensf.f):
C
GRAVity
C
TILT of
crust w r t
plumb line
C
ASTRonomical deflection of plumb line w r t star ...
C
ASTN ...
Newtonian
normalized (not recommended)
C
RADIal displacement
C
TANGential "
C
PTAN
horizontal displacement quasi-potential
C
POTEntial in free space
C
OTEP (ocean
tide
effective) potential
C
on
deformed surface
C
NN-Strain along azimuth
C
XX-Strain cross "
C
AREAl strain
C
ARES AREAl
strain,
normalized by sin(dis/2)
C
NGRA
Newton+GRAvity,
usually trunc.Love series
C
NOTP
Newton+OTeP
" "
" "
C
otem64 - instruction file:
C____________________________________________________________________________
C Instruction file:
C
C (1) Namelist
¶m
[default] or * = required
C----------------------------------------------------------------------------
C
COMP -
char*4
Component
*
C
TIDE -
char*4
Tide
symbol
*
C
NVER -
integer
Version
number
*
C
C NEXPGR -
integer
Expansion
length Green's
[2000]
C
C NF -
integer
The
side lengh of the entire
C
space
domain, incl. padding space.
C
Must
not exceed the limit
C
allowed
by the Green's spectrum
C
else
it will be reduced. (#)
C
Must
be a power of 2.
C
C NADDW, NADDE -
integer
On
output add version numbers ...
C
to
wavenumber
transform
[100]
C
to
explicit
convolution
[200]
C
COMP='PTAN':
1000 is added for 'E'ast
C
2000
is added for 'N'orth
C
C
ZVALUE,KFL,LFL,QCHFLG
Used
if no file is opened on unit 34:
C
Complex
Zvalue is submitted where
C
KFL
<= Flagvalue <= LFL
C
QCHFLG=.true.:
The flag array
C
is
updated to contain 1 at nodes
C
that
pass the criterion
C
and
0 elsewhere.
C
[(1.0,.0), 1, 1, .FALSE.]
C
C QRead_Spectrum - logical -
.false.: Transf.Spectrum = const. [.true.]
C
QWave_Trans -
logical - .false.:
Skip
[.true.]
C QExplicit_Convolve - logical
-
.f.:
Skip
[.false.]
C
QSalMask -
logical
For
SAL iteration, impose
C
C
land/sea
mask before output [.false.]
C
C OPTOZ -
char*1
Option
for retrieving the
C
Outer-Zone
data (call GETOZ).
C
'F'
-
full
['F']
C
C
TARGET -
char*8
Replacement
OZ tides
['O]']
C
TARGET_XC -
char*8
For
explicit
convolution ['S]']
C
C
C (#) If NF is less than
2*max(N,M) the padding spaces are less than
C
quadrant-size. When NF approaches max(N,M) the sacrifice is an
C
increasing degree of circular correlation.
C
C
C (2) File opening block:
C Unit 21 - Input, global
loading
effect (output from otem92, optional (*))
C
31
- Z-flags (resolved) array and area
C
32
- Tide solution array
C
33
- Green's fct. spectrum (output from otem67)
C
34
- OZ = Outer zone tides, optional (output from otem92)
C
41
- Output (1 or 2 OTSAL-potential arrays)
C
or
E and N (case PTAN)
C
42
- If E and N go to 41, the displacement potential goes here
C
44
- optional: tides after resolving OZ
C
45
- optional: flags after resolving OZ
C
C (*) I have no actual example where this input is used.
C____________________________________________________________________________
otem67 produces a 2-D spectral loading
transfer function, an array that
is read in by otem64 on unit 33.
The transfer spectrum depends on the specific grid spacing of the
model
(parameter SCALE); therefore,
otem64 should find the spectrum in a thematic
subdirectory.
otem67 -
instruction file
c
____________________________________________________________________________
c Instruction file:
c
c (1) Namelist
¶m:
Required
* | [default]
c
----------------------------------------------------------------------------
c NF, SCALE - domain
full
width, grid spacing
(km)
*
c
GFN - Green's table name
(excluding file suffix)
*
c
COMP - Component, usually
OTEP
*
c
QREAL - take the real part of the
Green's
function
[.true.]
c
EPSIT - accuracy limit for iteration
(integral)
[0.001]
c
MAXIT - maximum number of
iterations
[128]
c
c (2) File open block for
output,
log.unit 21: Greens WVN spectrum
c
____________________________________________________________________________
c
Example MDL22h/otem67.ins
R>
¶m
nf=512,
scale=8.
gfn='mc00fpt'
comp='RADI'
qreal=.true.
maxit=256
&end
21 <
MDL22h/sg00egbc.rad
q
OTEQ - the route to ocean
loading coefficients via
a regional ocean basin.
Create a tide model, solve explicitly time-dependent, compute
loading effects on the fly. Fit constituents solution and
combine with global loading effects.
Programs in /home/hgs/bin or
/home/hgs/PC/OTEQ/bin
Most instruction and data files in a model-specific catalog.
cd /home/hgs/OTEQ
mkdir MDL22 # a model-specific catalog
create
(MDL22/create.ins)
nes MDL22/prflm.prt # alias nes
nedit -fn spc08x14h \!*
#
unsolved:
print on paper (except
#
using
old MS-DOS)
otem92
(MDL22/otem92.ins)
otem16
(MDL22/otem16.ins) -> MDL22/ETD.DMP
otemt1
(MDL22/otemt0.ins)
otemt1
(MDL22/otemt01.ins) -> MDL22/evlr.daf evle.daf
evln.daf
tgg.daf
EVL -o MDL22/WESTFORD_evl.mc -I MDL22/evl 1
TGL
After air pressure modelling, run aploading and...
EVL -I -o GOMSL/WESTFORD_evlap.mc -c 'RAD EAS NOR' -r 6. -lcc
3
-x - GOMSL/evlap2. 1
~/tap/m/urtap
(MDL22/urtapmco.ins) -> MDL22/urtapmco.prl
prl2blq.pl -p MDL22/ MDL22/urtapmco.prl >> MDL22/ih.blq
olstm
(MDL22/olstm.ins) >LEPR> -> MDL22/GWA4LEPR.DAT
(MDL22/olstm.ins) >SCHW> -> MDL22/GWA4SCHW.DAT
olfg
(MDL22/olfg.ins) -> MDL22/lepr_i.blq
cnfload.out
(MDL22/cnfload.ins) -> MDL22/lepr_c.blq
(program got lost!)
| Back to HOW_TO
|
Back
to General
Information |
.bye