## Comparison of Swan results with measurements

Two Swan computations were made for the Haringvliet area, one without current, the other with zero current. It turned out that the two sets of results were not identical, unfortunately. See the Swan control file for the case with zero current.
With the OPGraph control file shown below, both computational results are compared with observations also available for this area.

```project 'SWAN conf. test' 'f01'

! this procedure generates a plot where results of 2 Swan runs
! are compared with observations
! the upper panel is for Hs, the lower panel for period

outfile 'f01_parms.eps'
option  frame=0  plot EPS
page panels 1, 2

transf 'sta' lin 0, 9
transf 'Hs'   lin 0, 4.2
transf 'Tm'   lin 0, 8
CHAR .27
appearance 33 pen=2 thickn=0.05 syms=0.32

for     '[var]' is 'Hs' 'Tm'   &
coupled '[num]' is '1' '2'     &
coupled '[Scol]' is '4' '5'    &
coupled '[Vcol]' is '2' '3'    &
do
panel  num [num] x 'sta' y '[var]'
axis x '   stations ' LABEL TICKS AT 1 2 3 4 5 6 7 8
axis y '[var]' auto
if '[num]' is '1'
curve 1 panel 0.36, 1, 0.36, 0.88, 1, 0.88
text panel 0.4, 0.9, 0. 'F01, Haringvliet: 23:00 UTC, Oct. 14, 1982'
endif

file 'verif.dat\hvl2300a.dat' 1 simple
obtain 'st' 1 1
obtain '[var]' 1 [Vcol]
graph 33 x 'st'  y '[var]'
if '[num]' is '2'
legend 'Observations around 23:00 hrs' panel 0.1, 0.21, 1
curve app=1 panel 0, 0.27, 0.6, 0.27, 0.6, 0.
endif
file 'verif.dat\hvl2300b.dat' 1 simple
obtain 'st' 1 1
obtain '[var]' 1 [Vcol]
graph 33 x 'st'  y '[var]'

appearance 22 pen=3 thickn=0.05 syms=0.32
appearance 24 pen=3 thickn=0.05 syms=0.32
for '[case]' is 'a', 'b'       &
coupled '[ltyp]' is '22','24'  &
do
file 'f01[case].tab' 2 simple
obtain '[var]' 2 [Scol]
assignto 'stat' series 1, 'numcomp([var])', 'numcomp([var])'
graph [ltyp] x 'stat'  y '[var]'
if '[num]' is '2'
if '[case]' is 'a'
legend 'SWAN, no current' panel 0.1, 0.09, 1
endif
if '[case]' is 'b'
legend 'SWAN, current=0' panel 0.1, 0.03, 1
endif
endif
endfor
endfor

INCLUDE '..\incl\copyrgt.opg'
test 0

! the procedure below generates a plot file showing 1d spectra
! computed by runs F01a and F01b
! with observations of spectra at 7 buoys and one wave staff
! there is a panel for each buoy

outfile 'f01_spec.eps'
option  frame=0  plot EPS
page panels= 2, 4
char 0.2
transf 'f' lin 0, 0.6

for     '[pan]' is '1' '2' '3' '4' '5' '6' '7' '8'        &
coupled '[sta]' is 'Wav' 'Wr1' 'Wr2' 'Wr3' 'Wr4' 'Wr5' 'Wr6' 'E75'   &
coupled '[vrt]' is '20' '15' '10' '10' '1' '2' '6' '1'        &
coupled '[hr1]' is '23' '21' '22' '22' '22' '22' '23' '23'    &
coupled '[mi1]' is '00' '00' '30' '50' '50' '50' '10' '00'    &
coupled '[hr2]' is '  ' '01' '  ' '  ' '23' '23' '  ' '  '    &
coupled '[mi2]' is '  ' '30' '  ' '  ' '10' '10' '  ' '  '   do
transf 'e' lin 0, [vrt]
panel num [pan] x 'f' y 'e'
axis x 'f' auto
axis y 'E' auto
text pan 0.5, 0.87, 0, 'station [sta]'

file 'f01B.spc' 1 swan
obtain 'AFREQ' From 1
obtain 'VaDens' from 1 for 'LOCATIONS' [pan]
graph 2 x  'AFREQ'  Y  'VaDens'
if '[pan]' is '5'
legend 'SWAN, current=0' panel 0.15, 1.05, 1
endif

file 'f01A.spc' 1 swan
obtain 'AFREQ' From 1
obtain 'VaDens' from 1 for 'LOCATIONS' [pan]
graph 3 x  'AFREQ'  Y  'VaDens'
if '[pan]' is '1'
legend 'SWAN, no current' panel 0.15, 1.05, 1
endif

file 'verif.dat\[sta]_[hr1][mi1].dat'  1 simple
obtain 'AFREQ' From 1 1
graph 4 x  'AFREQ'  Y  'VaDens'
legend '[hr1]:[mi1]' panel 0.5, 0.77, 1.
if '[hr2]' isnot '    ' then
file 'verif.dat\[sta]_[hr2][mi2].dat'  1 simple
obtain 'AFREQ' From 1 1
graph 4 x  'AFREQ'  Y  'VaDens'
legend '[hr2]:[mi2]' panel 0.5, 0.67, 1.
endif
endfor

INCLUDE '..\incl\copyrgt.opg'
stop
```
Two pictures result from this OPGraph job; one is called 'F01_parms.ps'; it shows the deviation between the computed wave parameters Hs and Tm01 and the measurements made by 7 wave buoys and a wave staff. One buoy provided directional information; this was used as boundary condition. It is noted that the computations are made for time 23.00; observations for this time were not present at all buoys; therefore in some locations two observations at instants around 23.00 were used.

The second picture is called 'F01_spec.ps'; it shows a comparison between computed 1-d spectra and observed spectra. These observed spectra were read from files of which the name consisted of the station name and the time of observation.

go to top