\hypertarget{namespacesloshing__initialization}{}\doxysection{sloshing\+\_\+initialization Module Reference}
\label{namespacesloshing__initialization}\index{sloshing\_initialization@{sloshing\_initialization}}


\doxysubsection{Detailed Description}
Initialization for the \char`\"{}sloshing\char`\"{} internal waves configuration. 

The module configures the model for the non-\/rotating sloshing test case. \doxysubsection*{Functions/\+Subroutines}
\begin{DoxyCompactItemize}
\item 
subroutine, public \mbox{\hyperlink{namespacesloshing__initialization_a0ecf019d4e63415fc16ce5ca4158f7f3}{sloshing\+\_\+initialize\+\_\+topography}} (D, G, param\+\_\+file, max\+\_\+depth)
\begin{DoxyCompactList}\small\item\em Initialization of topography. \end{DoxyCompactList}\item 
subroutine, public \mbox{\hyperlink{namespacesloshing__initialization_a5f999277fd41bc2730d6ef1dffffdb59}{sloshing\+\_\+initialize\+\_\+thickness}} (h, G, GV, US, param\+\_\+file, just\+\_\+read\+\_\+params)
\begin{DoxyCompactList}\small\item\em Initialization of thicknesses This routine is called when T\+H\+I\+C\+K\+N\+E\+S\+S\+\_\+\+C\+O\+N\+F\+IG is set to \textquotesingle{}sloshing\textquotesingle{}. \end{DoxyCompactList}\item 
subroutine, public \mbox{\hyperlink{namespacesloshing__initialization_a7806e9a52f01ac49906e7285e2c2bb3b}{sloshing\+\_\+initialize\+\_\+temperature\+\_\+salinity}} (T, S, h, G, GV, param\+\_\+file, eqn\+\_\+of\+\_\+state, just\+\_\+read\+\_\+params)
\begin{DoxyCompactList}\small\item\em Initialization of temperature and salinity. \end{DoxyCompactList}\end{DoxyCompactItemize}


\doxysubsection{Function/\+Subroutine Documentation}
\mbox{\Hypertarget{namespacesloshing__initialization_a7806e9a52f01ac49906e7285e2c2bb3b}\label{namespacesloshing__initialization_a7806e9a52f01ac49906e7285e2c2bb3b}} 
\index{sloshing\_initialization@{sloshing\_initialization}!sloshing\_initialize\_temperature\_salinity@{sloshing\_initialize\_temperature\_salinity}}
\index{sloshing\_initialize\_temperature\_salinity@{sloshing\_initialize\_temperature\_salinity}!sloshing\_initialization@{sloshing\_initialization}}
\doxysubsubsection{\texorpdfstring{sloshing\_initialize\_temperature\_salinity()}{sloshing\_initialize\_temperature\_salinity()}}
{\footnotesize\ttfamily subroutine, public sloshing\+\_\+initialization\+::sloshing\+\_\+initialize\+\_\+temperature\+\_\+salinity (\begin{DoxyParamCaption}\item[{real, dimension( g \%isd\+: g \%ied, g \%jsd\+: g \%jed,  g \%ke), intent(out)}]{T,  }\item[{real, dimension( g \%isd\+: g \%ied, g \%jsd\+: g \%jed,  g \%ke), intent(out)}]{S,  }\item[{real, dimension( g \%isd\+: g \%ied, g \%jsd\+: g \%jed,  g \%ke), intent(in)}]{h,  }\item[{type(ocean\+\_\+grid\+\_\+type), intent(in)}]{G,  }\item[{type(verticalgrid\+\_\+type), intent(in)}]{GV,  }\item[{type(param\+\_\+file\+\_\+type), intent(in)}]{param\+\_\+file,  }\item[{type(eos\+\_\+type), pointer}]{eqn\+\_\+of\+\_\+state,  }\item[{logical, intent(in), optional}]{just\+\_\+read\+\_\+params }\end{DoxyParamCaption})}



Initialization of temperature and salinity. 

This subroutine initializes linear profiles for T and S according to reference surface layer salinity and temperature and a specified range. Note that the linear distribution is set up with respect to the layer number, not the physical position). 
\begin{DoxyParams}[1]{Parameters}
\mbox{\texttt{ in}}  & {\em g} & Ocean grid structure. \\
\hline
\mbox{\texttt{ in}}  & {\em gv} & The ocean\textquotesingle{}s vertical grid structure. \\
\hline
\mbox{\texttt{ out}}  & {\em t} & Potential temperature \mbox{[}degC\mbox{]}. \\
\hline
\mbox{\texttt{ out}}  & {\em s} & Salinity \mbox{[}ppt\mbox{]}. \\
\hline
\mbox{\texttt{ in}}  & {\em h} & Layer thickness \mbox{[}H $\sim$$>$ m or kg m-\/2\mbox{]}. \\
\hline
\mbox{\texttt{ in}}  & {\em param\+\_\+file} & A structure indicating the open file to parse for model parameter values. \\
\hline
 & {\em eqn\+\_\+of\+\_\+state} & Equation of state structure. \\
\hline
\mbox{\texttt{ in}}  & {\em just\+\_\+read\+\_\+params} & If present and true, this call will only read parameters without changing h. \\
\hline
\end{DoxyParams}


Definition at line 179 of file sloshing\+\_\+initialization.\+F90.


\begin{DoxyCode}{0}
\DoxyCodeLine{181   \textcolor{keywordtype}{type}(ocean\_grid\_type),                     \textcolor{keywordtype}{intent(in)}  :: G\textcolor{comment}{ !< Ocean grid structure.}}
\DoxyCodeLine{182   \textcolor{keywordtype}{type}(verticalGrid\_type),                   \textcolor{keywordtype}{intent(in)}  :: GV\textcolor{comment}{ !< The ocean's vertical grid structure.}}
\DoxyCodeLine{183 \textcolor{keywordtype}{  real}, \textcolor{keywordtype}{dimension(SZI\_(G),SZJ\_(G), SZK\_(G))}, \textcolor{keywordtype}{intent(out)} :: T\textcolor{comment}{ !< Potential temperature [degC].}}
\DoxyCodeLine{184 \textcolor{keywordtype}{  real}, \textcolor{keywordtype}{dimension(SZI\_(G),SZJ\_(G), SZK\_(G))}, \textcolor{keywordtype}{intent(out)} :: S\textcolor{comment}{ !< Salinity [ppt].}}
\DoxyCodeLine{185 \textcolor{keywordtype}{  real}, \textcolor{keywordtype}{dimension(SZI\_(G),SZJ\_(G), SZK\_(G))}, \textcolor{keywordtype}{intent(in)}  :: h\textcolor{comment}{ !< Layer thickness [H \string~> m or kg m-\/2].}}
\DoxyCodeLine{186   \textcolor{keywordtype}{type}(param\_file\_type),                     \textcolor{keywordtype}{intent(in)}  :: param\_file\textcolor{comment}{ !< A structure indicating the}}
\DoxyCodeLine{187 \textcolor{comment}{                                                            !! open file to parse for model}}
\DoxyCodeLine{188 \textcolor{comment}{                                                            !! parameter values.}}
\DoxyCodeLine{189   \textcolor{keywordtype}{type}(EOS\_type),                            \textcolor{keywordtype}{pointer}     :: eqn\_of\_state\textcolor{comment}{ !< Equation of state structure.}}
\DoxyCodeLine{190   \textcolor{keywordtype}{logical},       \textcolor{keywordtype}{optional}, \textcolor{keywordtype}{intent(in)}  :: just\_read\_params\textcolor{comment}{ !< If present and true, this call will}}
\DoxyCodeLine{191 \textcolor{comment}{                                                      !! only read parameters without changing h.}}
\DoxyCodeLine{192 }
\DoxyCodeLine{193   \textcolor{keywordtype}{integer} :: i, j, k, is, ie, js, je, nz}
\DoxyCodeLine{194 \textcolor{keywordtype}{  real}    :: delta\_S, delta\_T}
\DoxyCodeLine{195 \textcolor{keywordtype}{  real}    :: S\_ref, T\_ref;      \textcolor{comment}{! Reference salinity and temerature within}}
\DoxyCodeLine{196                                 \textcolor{comment}{! surface layer}}
\DoxyCodeLine{197 \textcolor{keywordtype}{  real}    :: S\_range, T\_range;  \textcolor{comment}{! Range of salinities and temperatures over the}}
\DoxyCodeLine{198                                 \textcolor{comment}{! vertical}}
\DoxyCodeLine{199   \textcolor{keywordtype}{integer} :: kdelta}
\DoxyCodeLine{200 \textcolor{keywordtype}{  real}    :: deltah}
\DoxyCodeLine{201 \textcolor{keywordtype}{  real}    :: xi0, xi1}
\DoxyCodeLine{202   \textcolor{keywordtype}{logical} :: just\_read    \textcolor{comment}{! If true, just read parameters but set nothing.}}
\DoxyCodeLine{203   \textcolor{keywordtype}{character(len=40)}  :: mdl = \textcolor{stringliteral}{"initialize\_temp\_salt\_linear"} \textcolor{comment}{! This subroutine's}}
\DoxyCodeLine{204                                                             \textcolor{comment}{! name.}}
\DoxyCodeLine{205 }
\DoxyCodeLine{206   is = g\%isc ; ie = g\%iec ; js = g\%jsc ; je = g\%jec ; nz = g\%ke}
\DoxyCodeLine{207 }
\DoxyCodeLine{208   just\_read = .false. ; \textcolor{keywordflow}{if} (\textcolor{keyword}{present}(just\_read\_params)) just\_read = just\_read\_params}
\DoxyCodeLine{209 }
\DoxyCodeLine{210   \textcolor{keyword}{call }get\_param(param\_file, mdl, \textcolor{stringliteral}{"S\_REF"}, s\_ref, \textcolor{stringliteral}{'Reference value for salinity'}, \&}
\DoxyCodeLine{211                  default=35.0, units=\textcolor{stringliteral}{'1e-\/3'}, do\_not\_log=just\_read)}
\DoxyCodeLine{212   \textcolor{keyword}{call }get\_param(param\_file, mdl, \textcolor{stringliteral}{"T\_REF"}, t\_ref, \textcolor{stringliteral}{'Reference value for temperature'}, \&}
\DoxyCodeLine{213                  units=\textcolor{stringliteral}{'degC'}, fail\_if\_missing=.not.just\_read, do\_not\_log=just\_read)}
\DoxyCodeLine{214 }
\DoxyCodeLine{215   \textcolor{comment}{! The default is to assume an increase by 2 ppt for the salinity and a uniform temperature.}}
\DoxyCodeLine{216   \textcolor{keyword}{call }get\_param(param\_file, mdl,\textcolor{stringliteral}{"S\_RANGE"},s\_range,\textcolor{stringliteral}{'Initial salinity range.'}, \&}
\DoxyCodeLine{217                  units=\textcolor{stringliteral}{'1e-\/3'}, default=2.0, do\_not\_log=just\_read)}
\DoxyCodeLine{218   \textcolor{keyword}{call }get\_param(param\_file, mdl,\textcolor{stringliteral}{"T\_RANGE"},t\_range,\textcolor{stringliteral}{'Initial temperature range'}, \&}
\DoxyCodeLine{219                  units=\textcolor{stringliteral}{'degC'}, default=0.0, do\_not\_log=just\_read)}
\DoxyCodeLine{220 }
\DoxyCodeLine{221   \textcolor{keywordflow}{if} (just\_read) \textcolor{keywordflow}{return} \textcolor{comment}{! All run-\/time parameters have been read, so return.}}
\DoxyCodeLine{222 }
\DoxyCodeLine{223   \textcolor{comment}{! Prescribe salinity}}
\DoxyCodeLine{224   \textcolor{comment}{!delta\_S = S\_range / ( G\%ke -\/ 1.0 )}}
\DoxyCodeLine{225 }
\DoxyCodeLine{226   \textcolor{comment}{!S(:,:,1) = S\_ref}}
\DoxyCodeLine{227   \textcolor{comment}{!do k = 2,G\%ke}}
\DoxyCodeLine{228   \textcolor{comment}{!  S(:,:,k) = S(:,:,k-\/1) + delta\_S}}
\DoxyCodeLine{229   \textcolor{comment}{!enddo}}
\DoxyCodeLine{230 }
\DoxyCodeLine{231   deltah = g\%max\_depth / nz}
\DoxyCodeLine{232   \textcolor{keywordflow}{do} j=js,je ; \textcolor{keywordflow}{do} i=is,ie}
\DoxyCodeLine{233     xi0 = 0.0}
\DoxyCodeLine{234     \textcolor{keywordflow}{do} k = 1,nz}
\DoxyCodeLine{235       xi1 = xi0 + deltah / g\%max\_depth \textcolor{comment}{! =  xi0 + 1.0 / real(nz)}}
\DoxyCodeLine{236       s(i,j,k) = 34.0 + 0.5 * s\_range * (xi0 + xi1)}
\DoxyCodeLine{237       xi0 = xi1}
\DoxyCodeLine{238 \textcolor{keywordflow}{    enddo}}
\DoxyCodeLine{239 \textcolor{keywordflow}{  enddo} ;\textcolor{keywordflow}{ enddo}}
\DoxyCodeLine{240 }
\DoxyCodeLine{241   \textcolor{comment}{! Prescribe temperature}}
\DoxyCodeLine{242   delta\_t = t\_range / ( g\%ke -\/ 1.0 )}
\DoxyCodeLine{243 }
\DoxyCodeLine{244   t(:,:,1) = t\_ref}
\DoxyCodeLine{245   \textcolor{keywordflow}{do} k = 2,g\%ke}
\DoxyCodeLine{246     t(:,:,k) = t(:,:,k-\/1) + delta\_t}
\DoxyCodeLine{247 \textcolor{keywordflow}{  enddo}}
\DoxyCodeLine{248   kdelta = 2}
\DoxyCodeLine{249   t(:,:,g\%ke/2 -\/ (kdelta-\/1):g\%ke/2 + kdelta) = 1.0}
\DoxyCodeLine{250 }

\end{DoxyCode}
\mbox{\Hypertarget{namespacesloshing__initialization_a5f999277fd41bc2730d6ef1dffffdb59}\label{namespacesloshing__initialization_a5f999277fd41bc2730d6ef1dffffdb59}} 
\index{sloshing\_initialization@{sloshing\_initialization}!sloshing\_initialize\_thickness@{sloshing\_initialize\_thickness}}
\index{sloshing\_initialize\_thickness@{sloshing\_initialize\_thickness}!sloshing\_initialization@{sloshing\_initialization}}
\doxysubsubsection{\texorpdfstring{sloshing\_initialize\_thickness()}{sloshing\_initialize\_thickness()}}
{\footnotesize\ttfamily subroutine, public sloshing\+\_\+initialization\+::sloshing\+\_\+initialize\+\_\+thickness (\begin{DoxyParamCaption}\item[{real, dimension( g \%isd\+: g \%ied, g \%jsd\+: g \%jed, gv \%ke), intent(out)}]{h,  }\item[{type(ocean\+\_\+grid\+\_\+type), intent(in)}]{G,  }\item[{type(verticalgrid\+\_\+type), intent(in)}]{GV,  }\item[{type(unit\+\_\+scale\+\_\+type), intent(in)}]{US,  }\item[{type(param\+\_\+file\+\_\+type), intent(in)}]{param\+\_\+file,  }\item[{logical, intent(in), optional}]{just\+\_\+read\+\_\+params }\end{DoxyParamCaption})}



Initialization of thicknesses This routine is called when T\+H\+I\+C\+K\+N\+E\+S\+S\+\_\+\+C\+O\+N\+F\+IG is set to \textquotesingle{}sloshing\textquotesingle{}. 

This routine initializes layer positions to set off a sloshing motion in the zonal direction in a rectangular basin. All layers have initially the same thickness but all interfaces (except bottom and sea surface) are displaced according to a half-\/period cosine, with maximum value on the left and minimum value on the right. This sets off a regular sloshing motion. 
\begin{DoxyParams}[1]{Parameters}
\mbox{\texttt{ in}}  & {\em g} & The ocean\textquotesingle{}s grid structure. \\
\hline
\mbox{\texttt{ in}}  & {\em gv} & The ocean\textquotesingle{}s vertical grid structure. \\
\hline
\mbox{\texttt{ in}}  & {\em us} & A dimensional unit scaling type \\
\hline
\mbox{\texttt{ out}}  & {\em h} & The thickness that is being initialized \mbox{[}H $\sim$$>$ m or kg m-\/2\mbox{]}. \\
\hline
\mbox{\texttt{ in}}  & {\em param\+\_\+file} & A structure indicating the open file to parse for model parameter values. \\
\hline
\mbox{\texttt{ in}}  & {\em just\+\_\+read\+\_\+params} & If present and true, this call will only read parameters without changing h. \\
\hline
\end{DoxyParams}


Definition at line 56 of file sloshing\+\_\+initialization.\+F90.


\begin{DoxyCode}{0}
\DoxyCodeLine{57   \textcolor{keywordtype}{type}(ocean\_grid\_type),   \textcolor{keywordtype}{intent(in)}  :: G\textcolor{comment}{           !< The ocean's grid structure.}}
\DoxyCodeLine{58   \textcolor{keywordtype}{type}(verticalGrid\_type), \textcolor{keywordtype}{intent(in)}  :: GV\textcolor{comment}{          !< The ocean's vertical grid structure.}}
\DoxyCodeLine{59   \textcolor{keywordtype}{type}(unit\_scale\_type),   \textcolor{keywordtype}{intent(in)}  :: US\textcolor{comment}{          !< A dimensional unit scaling type}}
\DoxyCodeLine{60 \textcolor{keywordtype}{  real}, \textcolor{keywordtype}{dimension(SZI\_(G),SZJ\_(G),SZK\_(GV))}, \&}
\DoxyCodeLine{61                            \textcolor{keywordtype}{intent(out)} :: h\textcolor{comment}{           !< The thickness that is being initialized [H \string~> m or kg m-\/2].}}
\DoxyCodeLine{62   \textcolor{keywordtype}{type}(param\_file\_type),   \textcolor{keywordtype}{intent(in)}  :: param\_file\textcolor{comment}{  !< A structure indicating the open file}}
\DoxyCodeLine{63 \textcolor{comment}{                                                      !! to parse for model parameter values.}}
\DoxyCodeLine{64   \textcolor{keywordtype}{logical},       \textcolor{keywordtype}{optional}, \textcolor{keywordtype}{intent(in)}  :: just\_read\_params\textcolor{comment}{ !< If present and true, this call will}}
\DoxyCodeLine{65 \textcolor{comment}{                                                      !! only read parameters without changing h.}}
\DoxyCodeLine{66 }
\DoxyCodeLine{67 \textcolor{keywordtype}{  real}    :: displ(SZK\_(G)+1)   \textcolor{comment}{! The interface displacement in depth units.}}
\DoxyCodeLine{68 \textcolor{keywordtype}{  real}    :: z\_unif(SZK\_(G)+1)  \textcolor{comment}{! Fractional uniform interface heights [nondim].}}
\DoxyCodeLine{69 \textcolor{keywordtype}{  real}    :: z\_inter(SZK\_(G)+1) \textcolor{comment}{! Interface heights, in depth units.}}
\DoxyCodeLine{70 \textcolor{keywordtype}{  real}    :: a0                 \textcolor{comment}{! The displacement amplitude in depth units.}}
\DoxyCodeLine{71 \textcolor{keywordtype}{  real}    :: weight\_z           \textcolor{comment}{! A (misused?) depth-\/space weighting, in inconsistent units.}}
\DoxyCodeLine{72 \textcolor{keywordtype}{  real}    :: x1, y1, x2, y2     \textcolor{comment}{! Dimensonless parameters.}}
\DoxyCodeLine{73 \textcolor{keywordtype}{  real}    :: x, t               \textcolor{comment}{! Dimensionless depth coordinates?}}
\DoxyCodeLine{74   \textcolor{keywordtype}{logical} :: use\_IC\_bug         \textcolor{comment}{! If true, set the initial conditions retaining an old bug.}}
\DoxyCodeLine{75   \textcolor{keywordtype}{logical} :: just\_read          \textcolor{comment}{! If true, just read parameters but set nothing.}}
\DoxyCodeLine{76   \textcolor{comment}{! This include declares and sets the variable "version".}}
\DoxyCodeLine{77 \textcolor{preprocessor}{\# include "version\_variable.h"}}
\DoxyCodeLine{78 \textcolor{preprocessor}{}  \textcolor{keywordtype}{character(len=40)}  :: mdl = \textcolor{stringliteral}{"sloshing\_initialization"}\textcolor{comment}{ !< This module's name.}}
\DoxyCodeLine{79 }
\DoxyCodeLine{80   \textcolor{keywordtype}{integer} :: i, j, k, is, ie, js, je, nx, nz}
\DoxyCodeLine{81 }
\DoxyCodeLine{82   is = g\%isc ; ie = g\%iec ; js = g\%jsc ; je = g\%jec ; nz = g\%ke}
\DoxyCodeLine{83 }
\DoxyCodeLine{84   just\_read = .false. ; \textcolor{keywordflow}{if} (\textcolor{keyword}{present}(just\_read\_params)) just\_read = just\_read\_params}
\DoxyCodeLine{85   \textcolor{keywordflow}{if} (.not.just\_read) \textcolor{keyword}{call }log\_version(param\_file, mdl, version, \textcolor{stringliteral}{""})}
\DoxyCodeLine{86   \textcolor{keyword}{call }get\_param(param\_file, mdl, \textcolor{stringliteral}{"SLOSHING\_IC\_AMPLITUDE"}, a0, \&}
\DoxyCodeLine{87                  \textcolor{stringliteral}{"Initial amplitude of sloshing internal interface height "}//\&}
\DoxyCodeLine{88                  \textcolor{stringliteral}{"displacements it the sloshing test case."}, \&}
\DoxyCodeLine{89                  units=\textcolor{stringliteral}{'m'}, default=75.0, scale=us\%m\_to\_Z, do\_not\_log=just\_read)}
\DoxyCodeLine{90   \textcolor{keyword}{call }get\_param(param\_file, mdl, \textcolor{stringliteral}{"SLOSHING\_IC\_BUG"}, use\_ic\_bug, \&}
\DoxyCodeLine{91                  \textcolor{stringliteral}{"If true, use code with a bug to set the sloshing initial conditions."}, \&}
\DoxyCodeLine{92                  default=.false., do\_not\_log=just\_read)}
\DoxyCodeLine{93 }
\DoxyCodeLine{94   \textcolor{keywordflow}{if} (just\_read) \textcolor{keywordflow}{return} \textcolor{comment}{! All run-\/time parameters have been read, so return.}}
\DoxyCodeLine{95 }
\DoxyCodeLine{96   \textcolor{comment}{! Define thicknesses}}
\DoxyCodeLine{97   \textcolor{keywordflow}{do} j=g\%jsc,g\%jec ; \textcolor{keywordflow}{do} i=g\%isc,g\%iec}
\DoxyCodeLine{98 }
\DoxyCodeLine{99     \textcolor{comment}{! Define uniform interfaces}}
\DoxyCodeLine{100     \textcolor{keywordflow}{do} k = 0,nz}
\DoxyCodeLine{101       z\_unif(k+1) = -\/real(k)/real(nz)}
\DoxyCodeLine{102 \textcolor{keywordflow}{    enddo}}
\DoxyCodeLine{103 }
\DoxyCodeLine{104     \textcolor{comment}{! 1. Define stratification}}
\DoxyCodeLine{105     \textcolor{keywordflow}{do} k = 1,nz+1}
\DoxyCodeLine{106 }
\DoxyCodeLine{107       \textcolor{comment}{! Thin pycnocline in the middle}}
\DoxyCodeLine{108       \textcolor{comment}{!z\_inter(k) = (2.0**(n-\/1)) * (z\_unif(k) + 0.5)**n -\/ 0.5}}
\DoxyCodeLine{109 }
\DoxyCodeLine{110       \textcolor{comment}{! Thin pycnocline in the middle (piecewise linear profile)}}
\DoxyCodeLine{111       x1 = 0.30; y1 = 0.48; x2 = 0.70; y2 = 0.52}
\DoxyCodeLine{112 }
\DoxyCodeLine{113       x = -\/z\_unif(k)}
\DoxyCodeLine{114 }
\DoxyCodeLine{115       \textcolor{keywordflow}{if} ( x <= x1 ) \textcolor{keywordflow}{then}}
\DoxyCodeLine{116         t = y1*x/x1}
\DoxyCodeLine{117       \textcolor{keywordflow}{elseif} ( (x > x1 ) .and. ( x < x2 )) \textcolor{keywordflow}{then}}
\DoxyCodeLine{118         t = y1 + (y2-\/y1) * (x-\/x1) / (x2-\/x1)}
\DoxyCodeLine{119       \textcolor{keywordflow}{else}}
\DoxyCodeLine{120         t = y2 + (1.0-\/y2) * (x-\/x2) / (1.0-\/x2)}
\DoxyCodeLine{121 \textcolor{keywordflow}{      endif}}
\DoxyCodeLine{122 }
\DoxyCodeLine{123       t = -\/ z\_unif(k)}
\DoxyCodeLine{124 }
\DoxyCodeLine{125       z\_inter(k) = -\/t * g\%max\_depth}
\DoxyCodeLine{126 }
\DoxyCodeLine{127 \textcolor{keywordflow}{    enddo}}
\DoxyCodeLine{128 }
\DoxyCodeLine{129     \textcolor{comment}{! 2. Define displacement}}
\DoxyCodeLine{130     \textcolor{comment}{! a0 is set via get\_param; by default a0 is a 75m Displacement amplitude in depth units.}}
\DoxyCodeLine{131     \textcolor{keywordflow}{do} k = 1,nz+1}
\DoxyCodeLine{132 }
\DoxyCodeLine{133       weight\_z = -\/ 4.0 * ( z\_unif(k) + 0.5 )**2 + 1.0}
\DoxyCodeLine{134 }
\DoxyCodeLine{135       x = g\%geoLonT(i,j) / g\%len\_lon}
\DoxyCodeLine{136       \textcolor{keywordflow}{if} (use\_ic\_bug) \textcolor{keywordflow}{then}}
\DoxyCodeLine{137         displ(k) = a0 * cos(acos(-\/1.0)*x) + weight\_z * us\%m\_to\_Z}
\DoxyCodeLine{138       \textcolor{keywordflow}{else}}
\DoxyCodeLine{139         displ(k) = a0 * cos(acos(-\/1.0)*x) * weight\_z}
\DoxyCodeLine{140 \textcolor{keywordflow}{      endif}}
\DoxyCodeLine{141 }
\DoxyCodeLine{142       \textcolor{keywordflow}{if} ( k == 1 ) \textcolor{keywordflow}{then}}
\DoxyCodeLine{143         displ(k) = 0.0}
\DoxyCodeLine{144 \textcolor{keywordflow}{      endif}}
\DoxyCodeLine{145 }
\DoxyCodeLine{146       \textcolor{keywordflow}{if} ( k == nz+1 ) \textcolor{keywordflow}{then}}
\DoxyCodeLine{147         displ(k) = 0.0}
\DoxyCodeLine{148 \textcolor{keywordflow}{      endif}}
\DoxyCodeLine{149 }
\DoxyCodeLine{150       z\_inter(k) = z\_inter(k) + displ(k)}
\DoxyCodeLine{151 }
\DoxyCodeLine{152 \textcolor{keywordflow}{    enddo}}
\DoxyCodeLine{153 }
\DoxyCodeLine{154     \textcolor{comment}{! 3. The last interface must coincide with the seabed}}
\DoxyCodeLine{155     z\_inter(nz+1) = -\/g\%bathyT(i,j)}
\DoxyCodeLine{156     \textcolor{comment}{! Modify interface heights to make sure all thicknesses are strictly positive}}
\DoxyCodeLine{157     \textcolor{keywordflow}{do} k = nz,1,-\/1}
\DoxyCodeLine{158       \textcolor{keywordflow}{if} ( z\_inter(k) < (z\_inter(k+1) + gv\%Angstrom\_Z) ) \textcolor{keywordflow}{then}}
\DoxyCodeLine{159         z\_inter(k) = z\_inter(k+1) + gv\%Angstrom\_Z}
\DoxyCodeLine{160 \textcolor{keywordflow}{      endif}}
\DoxyCodeLine{161 \textcolor{keywordflow}{    enddo}}
\DoxyCodeLine{162 }
\DoxyCodeLine{163     \textcolor{comment}{! 4. Define layers}}
\DoxyCodeLine{164     \textcolor{keywordflow}{do} k = 1,nz}
\DoxyCodeLine{165       h(i,j,k) = gv\%Z\_to\_H * (z\_inter(k) -\/ z\_inter(k+1))}
\DoxyCodeLine{166 \textcolor{keywordflow}{    enddo}}
\DoxyCodeLine{167 }
\DoxyCodeLine{168 \textcolor{keywordflow}{  enddo} ;\textcolor{keywordflow}{ enddo}}
\DoxyCodeLine{169 }

\end{DoxyCode}
\mbox{\Hypertarget{namespacesloshing__initialization_a0ecf019d4e63415fc16ce5ca4158f7f3}\label{namespacesloshing__initialization_a0ecf019d4e63415fc16ce5ca4158f7f3}} 
\index{sloshing\_initialization@{sloshing\_initialization}!sloshing\_initialize\_topography@{sloshing\_initialize\_topography}}
\index{sloshing\_initialize\_topography@{sloshing\_initialize\_topography}!sloshing\_initialization@{sloshing\_initialization}}
\doxysubsubsection{\texorpdfstring{sloshing\_initialize\_topography()}{sloshing\_initialize\_topography()}}
{\footnotesize\ttfamily subroutine, public sloshing\+\_\+initialization\+::sloshing\+\_\+initialize\+\_\+topography (\begin{DoxyParamCaption}\item[{real, dimension(g\%isd\+:g\%ied,g\%jsd\+:g\%jed), intent(out)}]{D,  }\item[{type(dyn\+\_\+horgrid\+\_\+type), intent(in)}]{G,  }\item[{type(param\+\_\+file\+\_\+type), intent(in)}]{param\+\_\+file,  }\item[{real, intent(in)}]{max\+\_\+depth }\end{DoxyParamCaption})}



Initialization of topography. 


\begin{DoxyParams}[1]{Parameters}
\mbox{\texttt{ in}}  & {\em g} & The dynamic horizontal grid type \\
\hline
\mbox{\texttt{ out}}  & {\em d} & Ocean bottom depth in the units of depth\+\_\+max \\
\hline
\mbox{\texttt{ in}}  & {\em param\+\_\+file} & Parameter file structure \\
\hline
\mbox{\texttt{ in}}  & {\em max\+\_\+depth} & Maximum ocean depth in arbitrary units \\
\hline
\end{DoxyParams}


Definition at line 31 of file sloshing\+\_\+initialization.\+F90.


\begin{DoxyCode}{0}
\DoxyCodeLine{32   \textcolor{keywordtype}{type}(dyn\_horgrid\_type),  \textcolor{keywordtype}{intent(in)}  :: G\textcolor{comment}{ !< The dynamic horizontal grid type}}
\DoxyCodeLine{33 \textcolor{keywordtype}{  real}, \textcolor{keywordtype}{dimension(G\%isd:G\%ied,G\%jsd:G\%jed)}, \&}
\DoxyCodeLine{34                            \textcolor{keywordtype}{intent(out)} :: D\textcolor{comment}{ !< Ocean bottom depth in the units of depth\_max}}
\DoxyCodeLine{35   \textcolor{keywordtype}{type}(param\_file\_type),   \textcolor{keywordtype}{intent(in)}  :: param\_file\textcolor{comment}{ !< Parameter file structure}}
\DoxyCodeLine{36 \textcolor{keywordtype}{  real},                    \textcolor{keywordtype}{intent(in)}  :: max\_depth\textcolor{comment}{ !< Maximum ocean depth in arbitrary units}}
\DoxyCodeLine{37 }
\DoxyCodeLine{38   \textcolor{comment}{! Local variables}}
\DoxyCodeLine{39   \textcolor{keywordtype}{integer}   :: i, j}
\DoxyCodeLine{40 }
\DoxyCodeLine{41   \textcolor{keywordflow}{do} i=g\%isc,g\%iec ; \textcolor{keywordflow}{do} j=g\%jsc,g\%jec}
\DoxyCodeLine{42     d(i,j) = max\_depth}
\DoxyCodeLine{43 \textcolor{keywordflow}{  enddo} ;\textcolor{keywordflow}{ enddo}}
\DoxyCodeLine{44 }

\end{DoxyCode}
