| 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
 tide='..', nver=0, qpacked=.true.
31 ^ &/FLZU.DAT
32 ^ &/${TIDE}${VCODE}.dat
33 ^ &/sg00egbc.rad
41 < &/${TIDE}${VCODE}_rad.dat

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        "      "    "     "

otem64 - instruction file
C Instruction file:
C (1) Namelist &param                              [default] or  * = required

C COMP         - char*4               Component                          *
C TIDE         - char*4               Tide symbol                        *
C NVER         - integer              Version number                     *
C NEXPGR       - integer              Expansion length Green's         [2000]

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 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 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 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                                     land/sea mask before output   [.false.]

C OPTOZ          - char*1             Option for retrieving the

C                                     Outer-Zone data (call GETOZ).
C                                     'F' - full                        ['F']
C TARGET         - char*8             Replacement OZ tides             ['O]']
C TARGET_XC      - char*8             For explicit convolution         ['S]']
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 (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 (*) I have no actual example where this input is used.

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 (1) Namelist &param:                                  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 (2) File open block for output, log.unit 21: Greens WVN spectrum
c ____________________________________________________________________________
Example  MDL22h/otem67.ins
 nf=512, scale=8.
21 < MDL22h/sg00egbc.rad

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


nes MDL22/prflm.prt  # alias nes nedit -fn spc08x14h \!*
                     # unsolved: print on paper (except
                     # using old MS-DOS)


  (MDL22/otem16.ins) -> MDL22/ETD.DMP

  (MDL22/otemt01.ins) -> MDL22/evlr.daf evle.daf evln.daf tgg.daf

EVL -o MDL22/WESTFORD_evl.mc -I MDL22/evl 1

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

  (MDL22/urtapmco.ins) -> MDL22/urtapmco.prl

prl2blq.pl -p MDL22/ MDL22/urtapmco.prl >> MDL22/ih.blq

  (MDL22/olstm.ins) >LEPR> -> MDL22/GWA4LEPR.DAT
  (MDL22/olstm.ins) >SCHW> -> MDL22/GWA4SCHW.DAT

  (MDL22/olfg.ins) -> MDL22/lepr_i.blq

  (MDL22/cnfload.ins) -> MDL22/lepr_c.blq

  (program got lost!)

| Back to HOW_TO  |  Back to General Information  |