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


Run BALTIC - Free Oscillations, 46 tides possible but switched off.
/*
&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


EXAMPLES:

1. Forcing the Bay of Bothnia, files named -BBaySlope
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'
Plot the results (time series and periodograms) using
    plot-tgg -t 336- BALTIC/tgg-BBaySlope.daf 


2. Forcing the Bay of Kiel, files named -KielSlope
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'

3. Forcing the Bay of Finland, files named -GoFSlope
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'
Gulf of Riga, GoRSplash ?  Slope?