Exciting basins with time- and space-dependent potential - keyword SLOPE

/*
 &param
 idttep=1                                     Compute TGP at every time step
 dt_fixed=90.0                                90 s fixed time step
 etd_slope_opt(1)='N:1.0,84.0,400.0,2.0'      A very wide slope, wavelength 400 grid units, cutoff radius twide as wide
 etd_slope_opt(2)='E:1.0,95.0,400.0,2.0'       ... the amplitude could by increased (don't forget higher watchdog thresholds).
 etd_slope_opt(3)='T:3.,33.,0.25'              ... start at hour 3.0, continue for 33 hours, sine period 33/0.25 hours long
                                              (i.e. the slope is suddenly removed)
More examples

 Greensf='mc00egbc'
 Gpath='/home/hgs/Oload/Greenstb/] '
 ctide='M2 '                                  S2 would be o.k. too; however, it's not important since we prescribe the time step.
 tide='M2'
 tideh=' '
 qtake_all=.true., qtake(1)=.true.            All tides are included (we'll switch all off, doesn't work with qtake)
 nrcyc=2                                      Two major cycles (implies two more namelist blocks will be input). 
 cycwb=2                                      Array output will start from cycle 2                                
 nsol=-1                                      Auto number the solutions
 iusap=11                                     File unit for TGP input
 tstep=0.6                                    Sub-criticality of time step (not in effect here)
 nvu=46                                       Max number of tides to input for TGP and AB
 gef=1.00                                     Effective gravity value: full force
 slp=0.02                                     SAL parameter
 nonlin=0                                     Completely linear equations.
 path='./BALTIC/]'
 imon=30
 jmon=21
 itgra=10                                     Graphic display interval for Elevation
 range=1.6                                    ... data range
 hmin=1.0                                     Wet nodes with shallower depth are clipped at 1 m.
 Fricp=0.0, 0.0, 0.0, 0.0                     No friction at all.
 n_ramp=100,n_hold_damp=0,n_close_damp=101    
 geocxyz=.true.                               Eval-loading site coordinate type
 qfb=.false.                                  Free active boundary if .true. (remove currents, apply elevation)
 qspeca=.false.                                
 ispeca=23,44,26,47
 pspeca(1)=.1, pspeca(2)=.5, pspeca(3)=200.
 qbt=.false.                                  No tides on active boundary 
 qlt=.false.                                  No excitation by body tides
 start_date=2010,10,1, start_time=0,0,0,0
 bark_wd=3., bite_wd=10., qpoint_wd=.true.    These figures are rather low
 samp_int_tgg=1800.0                          not too sparse output, Nyquist should'nt be near suspected mode frequencies
 samp_int_evl=1800.0
 &end
31 < BALTIC/FLZUb.DAT
32 < BALTIC/FLMH.DAT
33 < BALTIC/ABOUND.DAT
03 * BALTIC/d1.dmp
   Q
8
Ringhals          30  47
StPeters         156  97
Haparand         100 171
Gavle             69 100
Stockhlm          82  88
Kiel              11  17
Gdansk            77  13
Mariehmn          89  91
8, -1, -1
Ringhals          30  47
StPeters         156  97
Haparand         100 171
Gavle             69 100
Stockhlm          82  88
Kiel              11  17
Gdansk            77  13
Mariehmn          89  91
 1, -1, -1
ONSALA      7205      3370658.676     711877.029    5349786.868
02 < BALTIC/dmp1.dat
41 < BALTIC/M2v4xt.dat
42 < BALTIC/M4v4xt.dat
75 B BALTIC/tgg-sAll.daf                       No other output at present but tide gauges and "potential sensors".
76 B BALTIC/tgp-sAll.daf                       (later we'll use the _EL_OUT feature)
65 B BALTIC/evlv-sAll.daf
66 B BALTIC/evlv-sAll.daf
67 B BALTIC/evlv-sAll.daf
   Q
Namelist for the start-up cycle:
 &param
 idttep=1
 iutgg=75
 iuevl=65
 iutgp=76
 nsver=500
 jsver=0
 Fricp=0.0, 0.0, 0.0, 0.e0
 nonlin=0
 Gra_Opt='Y'
 phacut=0.
 qspeca=.false.
 qmyprc=.false.
 kcyc=200                                      Integration time, in periods of the M2 here. 
 qadj_tend=.false.                             ... Don't adjust kcyc; the M2 is not interesting here.
 &end

 &param
 idttep=1
 iutgg=75
 iuevl=65
 nsver=0
 jsver=0
 Fricp=0.0, 0.0, 0.0, 0.0
 nonlin=0
 Gra_Opt='Y'
 kcyc=500
 qadj_tend=.false.
 &end

etd_slope_opt(1)='N:1.0,146.0,50.0,1.0'
etd_slope_opt(2)='E:0.1,103.0,30.0,1.0'
etd_slope_opt(3)='T:3.,33.,0.25'

etd_slope_opt(1)='N:1.0,10.0,30.0,1.0'
etd_slope_opt(2)='E:1.0,10.0,30.0,1.0'
etd_slope_opt(3)='T:3.,33.,0.25'

etd_slope_opt(1)='N:1.0,92.0,30.0,1.0'
etd_slope_opt(2)='E:1.0,156.0,50.0,1.0'
etd_slope_opt(3)='T:3.,33.,0.25'