Commented instructions

 Run test   MDL22 - B a y  o f  F u n d y   - 39 tides    Comment
                                                               "
 Start                                                         "
                                                               "
/*                                                       Comment end

 &param                                                  Namelist block #1
 Greensf='mc00egbc'
 Gpath='/home/hgs/Oload/Greenstb/] '
 ctide='M2'                                      Harmonic analysis tide
 tide='S2'                                       Time-step determination tide
 tideh='M4'                                Upper-harmonic analysis (nonlinear product) tide
 nrcyc=2                                         Two namelist blocks = two iteration stages
 cycwb=2                                         Start to output in the second stage
 nsol=-1                                         Don't read a dump file = start from scratch
 tstep=0.9729                          Subcriticality factor for time step (a qualified guess)
 nvu=40                                          Number of tides in the spectrum of forcing fields
 gef=1.00                                        Factor on gravity
 slp=0.02                                        Self-attraction and -loading parameter
 nonlin=3                                        Both kinds of nonlinear terms in Diff.eq.
 path='./]'                                      Path to ETD.DMP file
 imon=12, jmon=21                                Monitor point
 range=2.0                                               Graphic display: amplitude range
 Fricp=0.1,0.0,0.18,1.3e4                                Damping and Friction parameters
 n_ramp=200,n_hold_damp=100,n_close_damp=500             Start-up softly
 geocxyz=.true.                                  Geodynamic site record (below) in geocentric coordinates
 qfb=.true.                                      Free the inner (double) row of active boundary on Z-grid
 qspeca=.false.                                  No "Special Area" exceptions
 ispeca=23,44,26,47                              Has no effect because of qspeca=.false.
 pspeca(1)=.1, pspeca(2)=.5, pspeca(3)=200.      Has no effect because of qspeca=.false.
 qbt=.true.                                      Excite with body tide
 qlt=.true.                                        "         loading "
 start_date=91,3,9, start_time=0,0,0,0           Date and time for start (should be near otem16 date)
 &end
31 < MDL22h/FLZU.DAT                             Open-file block: input
32 < MDL22h/FLMH.DAT
33 < MDL22h/ABOUND.DAT
   Q
10                                               10 tide gauge records follow: result table
Boston            12  21
Sandwich CapeCod  20  13
Portland Maine    19  37
Yarmouth NScotia  59  42
St.John NBrunsw1  57  58
St.John NBrunsw2  57  59
Wolfville         80  61
Truro             87  63
Sackville         78  69
Port Greville     77  63
10, -1                                           10 tide gauge records follow: time-series
Boston            12  21
Sandwich CapeCod  20  13
Portland Maine    19  37
Yarmouth NScotia  59  42
St.John NBrunsw1  57  58
St.John NBrunsw2  57  59
Wolfville         80  61
Truro             87  63
Sackville         78  69
Port Greville     77  63
 1, -1                                            1 geodynamic site follows
WESTFORD    7205      1492404.7400  -4457266.5170   4296881.7740
02 < MDL22h/dmp1.dat                              Open-file block, output
41 < MDL22h/M2v1.dat
42 < MDL22h/M4v1.dat
   Q
 &param                                           Second namelist block
 kcyc=2                                           2 tidal cycles to iterate through
 jsver=0                                          Extra increment for version number counting
 Fricp=0.1,0.0,0.18,1.3e4
 nonlin=-1                                        Don't change
 Gra_Opt='n'                                      No graphics
 n_steps=20                                       Add 20 steps to cycle length before rounding
 qadj_tend=.false.                                Round to get full tidal cycles? Don't
 &end

 &param
 nsver=0
 jsver=0
 Fricp=0.1,0.0,0.18,1.3e4
 nonlin=-1
 Gra_Opt='n'
 kcyc=14
 qadj_tend=.true.
 &end

With the above choice of  ctide, kcyc, n_steps  and  qadj_tend   the program will go first through two
complete cycles of S2 tide + 20 steps, then 14 cycles of S2.

How to continue that job? Click here.

.bye