--- manual/s_phys_pkgs/text/obcs.tex 2011/11/09 11:48:40 1.15 +++ manual/s_phys_pkgs/text/obcs.tex 2016/06/15 13:40:53 1.18 @@ -233,9 +233,9 @@ and for each meridional position $j=1,\ldots,N_y$, a zonal index $i$ specifies the Eastern/Western OB position. For Northern/Southern OB this defines an $N_x$-dimensional -``row'' array $\tt OB\_Jnorth(Ny)$ / $\tt OB\_Jsouth(Ny)$, +``row'' array $\tt OB\_Jnorth(Nx)$ / $\tt OB\_Jsouth(Nx)$, and an $N_y$-dimenisonal -``column'' array $\tt OB\_Ieast(Nx)$ / $\tt OB\_Iwest(Nx)$. +``column'' array $\tt OB\_Ieast(Ny)$ / $\tt OB\_Iwest(Ny)$. Positions determined in this way allows Northern/Southern OBs to be at variable $j$ (or $y$) positions, and Eastern/Western OBs at variable $i$ (or $x$) positions. @@ -268,6 +268,24 @@ eg. $\tt OB\_Jnorth(3)=-1$ means that the point $\tt (3,Ny)$ is a northern OB. +\noindent\textbf{Simple examples:} For a model grid with $ N_{x}\times +N_{y} = 120\times144$ horizontal grid points with four open boundaries +along the four egdes of the domain, the simplest way of specifying the +boundary points in \code{data.obcs} is: +\begin{verbatim} + OB_Ieast = 144*-1, +# or OB_Ieast = 144*120, + OB_Iwest = 144*1, + OB_Jnorth = 120*-1, +# or OB_Jnorth = 120*144, + OB_Jsouth = 120*1, +\end{verbatim} +If only the first $50$ grid points of the southern boundary are +boundary points: +\begin{verbatim} + OB_Jsouth(1:50) = 50*1, +\end{verbatim} + \noindent \textsf{Add special comments for case \#define NONLIN\_FRSURF, see obcs\_ini\_fixed.F} @@ -279,7 +297,7 @@ \paragraph{OBCS\_READPARMS:} ~ \\ Set OB positions through arrays -{\tt OB\_Jnorth(Ny), OB\_Jsouth(Ny), OB\_Ieast(Nx), OB\_Iwest(Nx)}, +{\tt OB\_Jnorth(Nx), OB\_Jsouth(Nx), OB\_Ieast(Ny), OB\_Iwest(Ny)}, and runtime flags (see Table \ref{tab:pkg:obcs:runtime_flags}). \paragraph{OBCS\_CALC:} ~ \\