--- manual/s_phys_pkgs/text/obcs.tex 2011/10/26 08:29:48 1.14 +++ 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:} ~ \\ @@ -356,9 +374,9 @@ in \code{data}, see \code{verification/exp4} for an example. \paragraph{OBCS\_CALC\_STEVENS:} ~ \\ -(THE IMPLEMENTATION OF THESE BOUNDARY CONDITIONS IS NOT COMPLETE. SO -FAR ONLY EASTERN AND WESTERN BOUNDARIES ARE SUPPORTED. PASSIVE TRACERS -AND NON-LINEAR FREE SURFACE ARE NOT SUPPORTED.) \\ +(THE IMPLEMENTATION OF THESE BOUNDARY CONDITIONS IS NOT +COMPLETE. PASSIVE TRACERS, SEA ICE AND NON-LINEAR FREE SURFACE ARE NOT +SUPPORTED PROPERLY.) \\ The boundary conditions following \citet{stevens:90} require the vertically averaged normal velocity (originally specified as a stream function along the open boundary) $\bar{u}_{ob}$ and the tracer fields