Scribd
Upload
89 views33 pages

Wavewatch Lecture Ibnu

This document discusses using wave models like Wavewatch III and SWAN to model ocean waves. It provides information on: 1. Wavewatch III and SWAN are commonly used wave models, with Wavewatch III being better for large global and regional domains while SWAN is more advanced for shallow waters. 2. Creating model grids involves nesting smaller high resolution regional domains within larger lower resolution global domains to save on computation. 3. Running Wavewatch III involves preparing inputs like wind data, generating grids, initializing the model run, and post-processing outputs.

Uploaded by

lhaftc2
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
89 views33 pages

Wavewatch Lecture Ibnu

This document discusses using wave models like Wavewatch III and SWAN to model ocean waves. It provides information on: 1. Wavewatch III and SWAN are commonly used wave models, with Wavewatch III being better for large global and regional domains while SWAN is more advanced for shallow waters. 2. Creating model grids involves nesting smaller high resolution regional domains within larger lower resolution global domains to save on computation. 3. Running Wavewatch III involves preparing inputs like wind data, generating grids, initializing the model run, and post-processing outputs.

Uploaded by

lhaftc2
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 33

IbnuSofian

-The National Coordinating Agency for Surveys and Mapping (Bakosurtanal)


-Coastal and marine research center, Bandung Institute of Technology (ITB)
TMI/TRMM
TMI/TRMM
danQuickscat
danQuickscat

Altimeter
Altimeter (Jason-
(Jason-
1&2,
1&2, Envisat,
Envisat, GFO)
GFO)

-Wave
-Wave models
models
(WW
(WW IIIIII and
and
SWAN)
SWAN)
-Tidal
-Tidal Forcing
Forcing
Wave Model
 Wavewatch III
 Global to Regional Wave Model
 Uniform Grid

 User friendly

 SWAN
 Regional to Local Wave Model
 Advanced model for shallow water

 Uniform and unstructured grid

 Not user friendly


Sample of wave model (1)
Sample of wave model (2)
1. High spatial resolutions
using nesting method
2. Need computer resources
and calculation time

Model

1. Spatial resolution 1 deg


2. Only available from late
December 2005

Altimeter SWH
System requirement
 Operating System (OS)
 Linux OS
 Unix OS  Mac OS X, Sun OS, etc

 Fortran Compiler
 Portland Group Fortran Compiler
 Trial version can be downloaded from www.pgroup.com
 Intel Fortran Compiler
 Non-commercial license for Linux from www.intel.com
System requirement (2)
 Matlab both for Unix and Windows
 Mexnc Library from
http://mexcdf.sourceforge.net/downloads/R2008a.ht
ml for ver. 2008a or below
 http://mexcdf.sourceforge.net/downloads/ for ver.
2008b or higher
 Grads from http://www.iges.org/grads/ for post-
processing
Step by step WWIII preparation
 WWIII code request
 Uncompress and unzip the package
 WWIII installation
 Run the model
WWIII code
 Request code from
http://polar.ncep.noaa.gov/waves/wavewatch/lice
nse.shtml
 Request username and password
 Download all package
 Grid generator:
gridgen.bin.tar.gzgridgen.ref.tar.gzgridgen.exp.tar.gz
 WWII code: install_wwatch3.sh wwatch3.ftn.tar
wwatch3.aux.tar wwatch3.inp.tar wwatch3.tst.tar
WWIII installation (1)
 Execute the installer install_wwatch3.sh
 And compile the wavewatch code by exec.
/WWII/bin/w3_make
WWIII installation (2)

ibnusofian$ sh install_wwatch3.sh
How to run WWIII?
 Installation
 Edit link and comp in /WWIII LOCATION/bin/link and comp
 Compiling the WWIII w3_make
 Generating the domain grid
 Create the nested region (small domain)
 Define the boundary buffer zone with no redundant point
 Note: the boundary started from grid no. 2 until n-1, the outer most
boundaries are read as land
 Create the global domain
 Define the boundaries of nested region using the long – lat format
 Preparation
 Wind input using GFS Final analysis or GFS forecasting 1º or 0.5º
rotating
 [Air Temp – Skin Temp] input (optional)
 Cold restartww3_strt
 Model Exec. ww3_shel
How to run WWIII? (2)
 Temporal forcing
 Hindcast 3 days  stable model
 Forecast 7 days
Grid generation
NxM grid matrix

The outer most grid is read as land and not used


Grid generation
 Global to regional wave model using the nesting
method
 Create the regional domain using
create_grid_regional.m
 Define the boundary input from the regional domain
with the I, j format
 Create the global domain using create_grid_global.m

 Define the boundary output for regional model with


long lat format
Grid generation (1)

 Expand the gridgen file


gridgen.bin.tar.gzgridgen.ref.tar.gzgridgen.exp.tar.gz
 Tar zxffile_name.tar.gz

 Open matlab in example directory /YOUR_WWIII_DIR/examples


Grid generation (2)

 Copy the original create_grid_global.m and


regional.m to another filename to keep the orig. file
 Edit the matlabm file to generate the domain grid
Grid generation (3)
Grid generation (3)..Results
$Boundry input
$South
Grid generation (4) 2 2 F
140 2 T
 Ww3_grid.inp 187 2 F
207 2 T
 Boundary Input 225 2 F
242 2 T
$North
2 122 F
32 122 T
81 122 F
120 122 T
130 122 F
242 122 T
$West
2 3 F
2 121 T
$East
242 3 F
242 18 T
242 27 F
242 35 T
242 41 F
242 121 T
Grid generation (5) Global..Results
Grid generation (6)
 Ww3_grid.inp
 Boundary Output
$Boundry output
$North
90 15 0.25 0 31
109.75 15 0.25 0 40
12215 0.25 0 113
$South
90 -15 0.25 0 139
Starting point  136.25 -15 0.25 0 21
Total points
145.75 -15 0.25 0 18
$West
90 -14.75 0 0.25 119
$East
150 -14.75 0 0.25 16
150 -8.75 0 0.25 9
150 -5.25 0 0.25 81
increments
Wavewatch preprocessing (ww3_grid)
Wavewatch preprocessing
 Run ww3_strt

 Run ww3_prep
Ww3_prep
 Wind input
 Download the GFS final (1 deg) or GFS forecasting (1 deg)
 Convert grib to netcdf using NCL (ncar common
language) nclconvert2nc <file.grb>
 Using panoply for reading the netcdf of GFS file
Wind input
Wind input
 Wind format
Copy and paste wind data array
 YYYYMMDD HHMMSS
from panoply, and change [tab]
 [360x181]  u wind data to [space] using hidemaru editor
 [360x181] v wind data (windows), gedit(Centos/RHEL
5.0 or greater), and nedit
 YYYYMMDD HHMMSS (centos/RHEL 4.x)
 360x181]  u wind data
Note: gedit in RHEL 4.x is little bit
 [360x181] v wind data
“buggy”
…
For mac OSX user can use Editra
 ..
(from http://editra.org/download)
 .. rather than “stupid” Textedit.app
Sample of ww3_prep.inp

'WND' 'LL' T
0 359 360 -90 90 181
'NAME' 3 1 '(..T..)' '(..F..)'
20 'wind_2000.txt’
Ww3_shel
FF Water levels
FF Currents
TF Winds
F Ice concentrations
F Assimilation data : Mean parameters
F Assimilation data : 1-D spectra
F Assimilation data : 2-D spectra.
20000123 070000
20000131 010000
$ Type 1 : Fields of mean wave parameters
20000123 070000 10800 20000131 010000
T TTTTTTTTTTTTTTTTT
$ Type 2 : Point output
20000123 070000 900 20000131 010000
0.0 0.0 'STOPSTRING'
$ Type 3 : Output along track.
20000123 070000 900 20000131 010000
T
$ Type 4 : Restart files (no additional data required).
20000123 070000 86400 20000131 010000
$ Type 5 : Boundary data (no additional data required).
20000123 070000 3600 20000131 010000
'STP‘
$ End of input file
Happy modeling

You might also like