| Back to HOW_TO  |  Back to General Information  |  otemt1 |
 

otemt01.ins - Continuing a model using otemt1


To see the startup set of instructions, click here.

Commented instructions
 

 Run   MDL22h - B a y  o f  F u n d y   - 39 tides
 Version qf: less turbulent and less quadratic friction, Hmin = 30 m
 Continue
 Harm.sol.: M2
                                             "
/*
 &param
 Greensf='mc00egbc'
 Gpath='/home/hgs/Oload/Greenstb/] '
 ctide='M2'
 tide='S2'
 tideh='NO'                                      Don't compute harmonic result for nonlinear product tide
 nrcyc=2
 cycwb=2
 nsol=999                                        Number of dump to retrieve from the file on unit 3 (*)
 nvu=40
 gef=1.00
 slp=0.02
 nonlin=3
 path='./]'
 imon=12
 jmon=21
 range=5.0
 Fricp=0.0, 0.0, 0.12, 0.7e4
 Hmin=30.
 n_ramp=200,n_hold_damp=100,n_close_damp=500
 geocxyz=.true.
 qfb=.true.
 qbt=.true.
 qlt=.true.
 qload_and_stop=.false.
 &end
31 < MDL22h/FLZU.DAT
32 < MDL22h/FLMH.DAT
33 < MDL22h/ABOUND.DAT
03 < MDL22h/dmp1.dat                             (*) Dump file, result from otemt0.ins
   Q
10
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
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
WESTFORD    7205      1492404.7400  -4457266.5170   4296881.7740
02 > MDL22h/dmpqf.dat                               New dump file; qf stands for quadratic friction
41 > MDL22h/M2qf.dat
42 * > MDL22h/M4qf.dat
75 B MDL22h/tggqf.daf                               time-series output: tide gauges
65 B MDL22h/evlqfr.daf                                   "         "  : crustal dynamics, vertical
66 B MDL22h/evlqfe.daf                                   "         "       "        "     east
67 B MDL22h/evlqfn.daf                                   "         "       "        "     north displacement
   C                                                evl stands for "evaluate loading"
   Q
 &param
 iutgg=-75                                          prepare tide gauge output, but don't write yet
 iuevl=-65
 nsver=3
 jsver=0
 Fricp=0.0, 0.0, 0.12, 0.7e4
 nonlin=3
 Gra_Opt='e'
 kcyc=7                                             7 cycles
 qadj_tend=.false.
 &end

 &param
 iutgg=75                                            activate output to tide gauge file
 iuevl=65                                            activate output to crustal dynamics file
 Fricp=0.0, 0.0, 0.12, 0.7e4
 nonlin=3
 Gra_Opt='e'
 kcyc=30                                             30 cycles
 qadj_tend=.false.
 &end

We have two iteration phases. The first asks for 7 cycles of S2 (really S2? Not harmonic solution tide M2? Well, I'm uncertain at the moment). The first asks for 7 cycles to allow a transition of the solution from one set of parameters to another. For instance we might change a friction parameter. The solution will generally produce larger tides, but it will take a series of time-steps until the new solution has established itself.
   The second block asks for 30 cycles. This might still be too little to obtain good frequency resolution, particularly if we want to solve for small harmonic tides. For obtaining an approximate M2 solution in a bay where M2 is the dominating tide the 30 cycles should be
enough.
   The time-series output to the tide-gauge file and the crustal-dynamics files  is active during the second block of iterations (cycwb=2)
The nature of your expiriments will be reflected in the naming of the files. This job highlights the use of quadratic friction. We might have a series of results to compare where we varied the eddy viscosity coefficient. Those jobs could result in files with names containing edv and
maybe a number or character to desginate the set of actual parameters used. There will be some need for  bookkeeping on paper or in  private documents.

To see the startup set of instructions, click here.

.bye