Skip to content

Commit

Permalink
Modified manual Sec. "3.4.9 Rotated grids"
Browse files Browse the repository at this point in the history
  • Loading branch information
CarstenHansen committed Oct 23, 2019
1 parent 35ea698 commit 6a0ffac
Show file tree
Hide file tree
Showing 5 changed files with 74 additions and 84 deletions.
42 changes: 29 additions & 13 deletions manual/num/rotagrid.tex
Original file line number Diff line number Diff line change
Expand Up @@ -4,8 +4,8 @@ \subsubsection{~Rotated grids} \label{sub:num_space_rotagrid}

\noindent
The rotated grid is a latitude-longitude (lat-lon) grid and is obtained by
rotating the North Pole to a new position at latitude $\phi_{p}$ and
longitude $\lambda_{p}$ in the standard latitude-longitude system. The new
rotating the North Pole along a longitude $\lambda_{p}$ to a new position at
latitude $\phi_{p}$ in the standard latitude-longitude system. The new
pole position is chosen so that the model domain of interest may be placed
around the rotated equatorial area for a evenly-spaced lat-lon mesh. For this
reason the rotated grid is also known as \emph{Equatorial grid}. For instance,
Expand All @@ -15,22 +15,38 @@ \subsubsection{~Rotated grids} \label{sub:num_space_rotagrid}
grid allows a much more evenly spaced lat-lon mesh in the NAEW domain than the
standard lat-lon grid in the same area.

In \ws\, the rotated grid is implemented with minimum changes to the original
In \ws, the rotated grid is implemented with minimum changes to the original
lat-lon grid. In fact, the rotated grid is treated just like the standard
lat-lon grid inside the model. To set up and run a rotated grid model configuration,
users should choose the regular lat-lon grid along with the {\code RTD} switch.
The rotated pole position is set using the {\code PLAT} and {\code PLON} variables in the
{\bf ww3\_grid.inp} namelist {\code ROTD}. Model input files, like wind, current and ice files
should be mapped on to the rotated grid. For convenience of nesting in standard
lat-lon grid frameworks, boundary conditions provided to and output from the
The rotated pole position is set using the {\code PLAT} and {\code PLON} variables
under the namelist {\code ROTD} in the input file {\bf ww3\_grid.inp}
(see Sect.~\ref{sec:config011}). If the pole is set as {\code PLAT = 90.0},
{\code PLON = -180.0}, the grid is treated as a standard lat-lon system.

Model input files, like wind, current and ice files
should be mapped on to the rotated grid. For convenience of nesting in standard
lat-lon grid frameworks, boundary conditions data provided to and output from the
rotated grid use spectra referenced to a standard grid north and standard lat-lon
grid points values, which are converted into rotated grid lat-lon inside \ws\.
The list of 2D spectral output locations in {\bf ww3\_shel.inp} are also specified in
standard lat-lon.
grid points values, which are converted into rotated grid lat-lon inside \ws.
The list of 2D spectral output locations in {\bf ww3\_shel.inp} or
{\bf ww3\_shel.nml} are also specified in standard lat-lon. When the boundary
conditions are applied to a model on a standard lat-lon grid, the {\code RTD}
switch is not required in compiling the executables for that model.

Output of spectra at boundary points to one-way nested inner grids are transferred
as described in Appendix\ \ref{app:nest}. Output b.c.\ may be defined in the input
file (see {\bf ww3\_grid.inp}, Sect.~\ref{sec:config012}) as a sequence of straight
lines given in coordinates of the inner grid. If the inner grid is rotated, it's
pole position must be defined as the values of the array elements
{\code BPLAT(\sl{n})} and {\code BPLON(\sl{n})} under namelist {\code ROTB}. The
array index {\sl{n}} is the index ({\sl{1:9}}) of the boundary conditions file
{\file nest{\sl{n}}.ww3}.

Model directional and x-y vector outputs can be converted to a standard grid
north reference by setting the UNROT variable in the {\bf ww3\_grid.inp} namelist
ROTD to True. With this set, for point outputs lat-lon locations all directional
Model directional and x-y vector outputs from a rotated grid can be converted
to a standard grid north reference by setting the UNROT variable in the
{\bf ww3\_grid.inp} namelist ROTD to True. With this set, for point outputs
lat-lon locations all directional
values such as wind direction, current direction and 2D spectra are converted
into standard lat-lon orientation. Functions to de-rotate gridded
fields are applied in {\bf ww3\_ounf}, {\bf ww3\_outf} and {\bf ww3\_grib}.
Expand Down
38 changes: 20 additions & 18 deletions model/inp/ww3_grid.inp
Original file line number Diff line number Diff line change
Expand Up @@ -345,19 +345,20 @@ $ SEAWND : Use sea-point only wind input. .FALSE.
$ &PSMC DTIME = 39600.0, LATMIN=85.0, RFMAXD = 36.0, LvSMC=3, JEQT=1344 /
$
$ Rotated pole ------------------------------------------------------ $
$ Pole parameters : Namelist ROTD
$ PLAT : Rotated pole latitude
$ PLON : Rotated pole longitude
$ UNROT : Logical, un-rotate directions to true north
$ Pole parameters : Namelist ROTD
$ PLAT : Rotated pole latitude
$ PLON : Rotated pole longitude
$ UNROT : Logical, un-rotate directions to
$ true north
$
$ Compile switch /RTD required.
$ Compile switch /RTD required.
$
$ These will be used to derive rotation angle corrections in the
$ model. The corrections are used for rotation of boundary spectra
$ and for restoring conventional lat/lon orientation of the
$ output spectra, mean direction or any related variables.
$ The PLAT/LON example below is a standard setting for Met
$ Office UK regional models.
$ These will be used to derive rotation angle corrections in the
$ model. The corrections are used for rotation of boundary spectra
$ and for restoring conventional lat/lon orientation of the
$ output spectra, mean direction or any related variables.
$ The PLAT/LON example below is a standard setting for Met
$ Office UK regional models.
$
$ &ROTD PLAT = 37.5, PLON = 177.5, UNROT = .TRUE. /
$
Expand All @@ -367,17 +368,18 @@ $ &ROTD PLAT = 90.0, PLON = -180.0, UNROT = .FALSE. /
$
$ Output boundary conditions to rotated pole grids ------------------ $
$
$ Pole parameters : Namelist ROTB
$ BPLAT(1:9) : Dest. grid rotated pole latitude
$ BPLON(1:9) : Rotated pole longitude
$ Pole parameters : Namelist ROTB
$ BPLAT(1:9) : Dest. grid rotated pole latitude
$ BPLON(1:9) : Rotated pole longitude
$
$ Compile switch /RTD required.
$
$ In the section 'Output boundary points' below, the points must be given
$ as straight lines in the reference system of each destination grid. The
$ poles of each destination grids have default values of standard lat/lon.
$ In the section 'Output boundary points' below, the points must be
$ given as straight lines in the reference system of each destination
$ grid. The poles of each destination grids have default values of
$ standard lat/lon.
$
$ &ROTB BPLAT(1)=90., BPLON(1)=-180., BPLAT(2)=90., BPLON(2)=-180. / $ (etc.)
$ &ROTB BPLAT(1)=90., BPLON(1)=-180., BPLAT(2)=90. / $ (etc.)
$
$ Output of 3D arrays------------------------------------------------- $
$ In order to limit the use of memory, arrays for 3D output fiels (i.e.
Expand Down
24 changes: 24 additions & 0 deletions model/nml/namelists.nml
Original file line number Diff line number Diff line change
@@ -0,0 +1,24 @@
$ -------------------------------------------------------------------- $
$ WAVEWATCH III - model parameters namelist
$ -------------------------------------------------------------------- $
$ This file is to be used with ww3_grid.nml where its name is the value
$ of GRID%NML
$
$ The contents is equal to the namelist input of ww3_grid.inp
$ (traditional form) ending with a mandatory string 'END OF NAMELISTS'
$
$ Start of namelist input section ------------------------------------ $
$
$ [ ... ]
$
$ In the 'Out of the box' test setup we run with sub-grid obstacles
$ and with continuous ice treatment.
$
&MISC CICE0 = 0.25, CICEN = 0.75, FLAGTR = 4 /
&FLX3 CDMAX = 3.5E-3 , CTYPE = 0 /
$
END OF NAMELISTS
$
$ -------------------------------------------------------------------- $
$ End of model parameters namelist file
$ -------------------------------------------------------------------- $
52 changes: 0 additions & 52 deletions model/nml/namelists_GULF.nml

This file was deleted.

2 changes: 1 addition & 1 deletion model/nml/ww3_grid.nml
Original file line number Diff line number Diff line change
Expand Up @@ -140,7 +140,7 @@
! -------------------------------------------------------------------- !
&GRID_NML
GRID%NAME = 'GULF OF NOWHERE'
GRID%NML = 'namelists_GULF.nml'
GRID%NML = 'namelists.nml'
GRID%TYPE = 'RECT'
GRID%COORD = 'SPHE'
GRID%CLOS = 'SMPL'
Expand Down

0 comments on commit 6a0ffac

Please sign in to comment.