\hypertarget{namespacew__module}{}\section{w\+\_\+module Module Reference}
\label{namespacew__module}\index{w\+\_\+module@{w\+\_\+module}}


2次元球面\+: 球面調和関数展開  


\subsection*{Functions/\+Subroutines}
\begin{DoxyCompactItemize}
\item 
subroutine, public \hyperlink{namespacew__module_a34cf65ed64bfeec316d31e8764e8ffc0}{w\+\_\+initial} (n\+\_\+in, i\+\_\+in, j\+\_\+in, np\+\_\+in)
\begin{DoxyCompactList}\small\item\em スペクトル変換の格子点数, 波数および O\+P\+E\+N\+MP 使用時の 最大スレッド数を設定する. \end{DoxyCompactList}\item 
subroutine, public \hyperlink{namespacew__module_a8fa63b83c6aee6f23c4560de239475e6}{w\+\_\+finalize}
\begin{DoxyCompactList}\small\item\em モジュールの終了処理(割り付け配列の解放)をおこなう. \end{DoxyCompactList}\end{DoxyCompactItemize}


\subsection{Detailed Description}
2次元球面\+: 球面調和関数展開 

\begin{DoxyAuthor}{Author}
Shin-\/ichi Takehiro, Youhei S\+A\+S\+A\+KI 
\end{DoxyAuthor}
\begin{DoxyCopyright}{Copyright}
Copyright (C) S\+P\+M\+O\+D\+EL Development Group, 2019. License is M\+I\+T/\+X11. see \hyperlink{COPYRIGHT}{C\+O\+P\+Y\+R\+I\+G\+HT} in detail
\end{DoxyCopyright}
spml/w\+\_\+module モジュールは球面上での 2 次元流体運動を 球面調和函数を用いたスペクトル法によって数値計算するための Fortran90 関数を提供する. \hyperlink{namespacew__module}{w\+\_\+module} は実際には以下の下部モジュールからなる．
\begin{DoxyItemize}
\item \hyperlink{namespacew__module__base}{w\+\_\+module\+\_\+base}\+: 基本変換
\item \hyperlink{namespacew__module__deriv}{w\+\_\+module\+\_\+deriv}\+: 微分計算
\item \hyperlink{namespacew__module__integral}{w\+\_\+module\+\_\+integral}\+: 積分・平均計算
\item \hyperlink{namespacew__module__spectrum}{w\+\_\+module\+\_\+spectrum}\+: スペクトル解析
\item \hyperlink{namespacew__module__interpolate}{w\+\_\+module\+\_\+interpolate}\+: 補間計算
\end{DoxyItemize}

内部で I\+S\+P\+A\+C\+K3 の S\+X\+P\+A\+C\+Kサブルーチンを呼んでいる. スペクトルデータお よび格子点データの格納方法や変換の詳しい計算法については I\+S\+P\+A\+C\+K3/\+S\+X\+P\+A\+CK のマニュアルを参照されたい.

\begin{DoxyWarning}{Warning}
後方非互換である点に注意されたい
\begin{DoxyItemize}
\item 旧版と異なり (I\+S\+P\+A\+C\+K/\+S\+X\+P\+A\+C\+Kのレベルで) 球面調和関数の正規化が異なる．
\item 旧版と異なり, l\+\_\+nm, nm\+\_\+l は初期化したのちにしか使うことができない.
\end{DoxyItemize}
\end{DoxyWarning}


\paragraph*{命名法}


\begin{DoxyItemize}
\item 関数名の先頭 ({\ttfamily w\+\_\+}, {\ttfamily nm\+\_\+}, {\ttfamily n\+\_\+}, {\ttfamily xy\+\_\+}, {\ttfamily x\+\_\+}, {\ttfamily y\+\_\+}) は, 返す値の形を示している.
\begin{DoxyItemize}
\item {\ttfamily w\+\_\+} \+: スペクトルデータ
\item {\ttfamily xy\+\_\+} \+: 2 次元格子点データ
\item {\ttfamily nm\+\_\+} \+: スペクトルデータの並んだ 3 次元配列
\item {\ttfamily n\+\_\+} \+: スペクトルデータの並んだ 2 次元配列
\item {\ttfamily x\+\_\+} \+: 経度方向 1 次元格子点データ
\item {\ttfamily y\+\_\+} \+: 緯度方向 1 次元格子点データ
\end{DoxyItemize}
\item 関数名の間の文字列({\ttfamily D\+Lon}, {\ttfamily Grad\+Lat}, {\ttfamily Grad\+Lat}, {\ttfamily Div\+Lon}, {\ttfamily Div\+Lat}, {\ttfamily Lapla}, {\ttfamily Lapla\+Inv}, {\ttfamily Jacobian})は, その関数の作用を表している.
\item 関数名の最後 ({\ttfamily \+\_\+w\+\_\+w}, {\ttfamily \+\_\+w}, {\ttfamily \+\_\+xy}, {\ttfamily \+\_\+x}, {\ttfamily \+\_\+y}) は, 入力変数の形スペクトルデータおよび格子点データであることを示している.
\begin{DoxyItemize}
\item {\ttfamily \+\_\+w} \+:\+: スペクトルデータ
\item {\ttfamily \+\_\+w\+\_\+w} \+:\+: 2 つのスペクトルデータ
\item {\ttfamily \+\_\+xy} \+:\+: 2 次元格子点データ
\item {\ttfamily \+\_\+x} \+:\+: 経度方向 1 次元格子点データ
\item {\ttfamily \+\_\+y} \+:\+: 緯度方向 1 次元格子点データ
\end{DoxyItemize}
\end{DoxyItemize}

\paragraph*{各データの種類の説明}


\begin{DoxyItemize}
\item 以下, {\ttfamily DP = kind(1.\+0\+D0)}
\item {\ttfamily xy} \+: 2 次元格子点データ.
\begin{DoxyItemize}
\item 変数の種類と次元は {\ttfamily real(\+D\+P), dimension(0\+:im-\/1,1\+:jm)}.
\item {\ttfamily im}, {\ttfamily jm} はそれぞれ経度, 緯度座標の格子点数であり, サブルーチン \hyperlink{namespacew__module_a34cf65ed64bfeec316d31e8764e8ffc0}{w\+\_\+\+Initial} にてあらかじめ設定しておく.
\end{DoxyItemize}
\item {\ttfamily w} \+: スペクトルデータ.
\begin{DoxyItemize}
\item 変数の種類と次元は {\ttfamily real(\+D\+P), dimension((nm+1)$\ast$(nm+1))}.
\item nm は球面調和函数の最大全波数であり, サブルーチン w\+\_\+\+Initial にて あらかじめ設定しておく.
\item スペクトルデータの格納のされ方は関数 l\+\_\+nm, nm\+\_\+l によって 調べることができる.
\end{DoxyItemize}
\item {\ttfamily nm} \+: スペクトルデータの並んだ 2 次元配列.
\begin{DoxyItemize}
\item 変数の種類と次元は {\ttfamily real(\+D\+P), dimension(0\+:nm,-\/nm\+:nm)}. 第 1 次元が水平全波数, 第 2 次元が帯状波数を表す.
\item {\ttfamily nm} は球面調和函数の最大全波数であり, サブルーチン \hyperlink{namespacew__module_a34cf65ed64bfeec316d31e8764e8ffc0}{w\+\_\+\+Initial} にてあらかじめ設定しておく.
\end{DoxyItemize}
\item {\ttfamily n} \+: スペクトルデータの並んだ 1 次元配列.
\begin{DoxyItemize}
\item 変数の種類と次元は {\ttfamily real(\+D\+P), dimension(0\+:nm)}.
\item 第 1 次元が水平全波数を表す. nm は球面調和函数の最大全波数であり, サブルーチン \hyperlink{namespacew__module_a34cf65ed64bfeec316d31e8764e8ffc0}{w\+\_\+\+Initial} にてあらかじめ設定しておく.
\end{DoxyItemize}
\item {\ttfamily x}, {\ttfamily y} \+: 経度, 緯度方向 1 次元格子点データ.
\begin{DoxyItemize}
\item 変数の種類と次元はそれぞれ {\ttfamily real(\+D\+P), dimension(0\+:im-\/1)} および {\ttfamily real(\+D\+P), dimension(1\+:jm)}.
\end{DoxyItemize}
\item {\ttfamily w\+\_\+} で始まる関数が返す値はスペクトルデータに同じ.
\item {\ttfamily xy\+\_\+} で始まる関数が返す値は 2 次元格子点データに同じ.
\item {\ttfamily x\+\_\+}, {\ttfamily y\+\_\+} で始まる関数が返す値は 1 次元格子点データに同じ.
\item スペクトルデータに対する微分等の作用とは, 対応する格子点データに 微分などを作用させたデータをスペクトル変換したものことである.
\end{DoxyItemize}

\paragraph*{初期化, 終了処理}

\hyperlink{namespacew__module__base}{w\+\_\+module\+\_\+base} にて定義


\begin{DoxyItemize}
\item \hyperlink{namespacew__module_a34cf65ed64bfeec316d31e8764e8ffc0}{w\+\_\+\+Initial} \+: スペクトル変換の格子点数, 波数, 領域の大きさの設定
\item \hyperlink{namespacew__module_a8fa63b83c6aee6f23c4560de239475e6}{w\+\_\+\+Finalize} \+: モジュールの終了処理(割り付け配列の解放)をおこなう.
\end{DoxyItemize}

\paragraph*{座標変数}


\begin{DoxyItemize}
\item \hyperlink{namespacew__module__base_ae875323468dc412db617d265fcf88b2a}{x\+\_\+\+Lon} \+: 経度座標格子点を格納した 1 次元配列
\item \hyperlink{namespacew__module__base_a4c34a1556ac18599ef73000b6b583579}{y\+\_\+\+Lat} \+: 緯度座標格子点を格納した 1 次元配列
\item \hyperlink{namespacew__module__base_a58d8184c4227f4fa64acabca0ac5c977}{x\+\_\+\+Lon\+\_\+\+Weight} \+: 重み座標を格納した 1 次元配列
\item \hyperlink{namespacew__module__base_a52e4db808630002b9e3e75b55f5efcb2}{y\+\_\+\+Lat\+\_\+\+Weight} \+: 重み座標を格納した 1 次元配列
\item \hyperlink{namespacew__module__base_a916da409070cdd92318081f7e13c1861}{xy\+\_\+\+Lon} \+: 経度座標(\+X,\+Y)の格子点データ(2 次元配列)
\item \hyperlink{namespacew__module__base_a263eec8100074f932519e03dcdc7fe70}{xy\+\_\+\+Lat} \+: 緯度座標(\+X,\+Y)の格子点データ(2 次元配列)
\end{DoxyItemize}

\paragraph*{基本変換}


\begin{DoxyItemize}
\item \hyperlink{namespacew__module__base_a8c988f23cc273731d93e9d3e249b64e5}{xy\+\_\+w} \+: スペクトルデータから格子データへの変換
\item \hyperlink{namespacew__module__base_ac2c97a298ba072c23d3a148018998ac5}{w\+\_\+xy} \+: 格子データからスペクトルデータへの変換
\item l\+\_\+nm, nm\+\_\+l \+: スペクトルデータの格納位置と全波数・帯状波数の変換
\item \hyperlink{namespacew__module__base_a496a5089e89bf920c12c56f29089d51b}{w\+\_\+\+Stream\+Potential2\+Vector} \+: 流線ポテンシャルから速度場計算
\item \hyperlink{namespacew__module__base_af7cb0ea7552281e416b5045af7ab3837}{w\+\_\+\+Vector2\+Vor\+Div} \+: 速度場から渦度発散を計算
\item \hyperlink{namespacew__module__base_a972eb288433b29554df9c889f8ef54e6}{w\+\_\+\+Vector\+Cos\+Lat2\+Vor\+Div} \+: 速度場から渦度発散を計算
\end{DoxyItemize}

\paragraph*{微分}

\hyperlink{namespacew__module__deriv}{w\+\_\+module\+\_\+deriv} にて定義


\begin{DoxyItemize}
\item \hyperlink{namespacew__module__deriv_a0bf8cc78dc8dde2691985c05044d9c4b}{w\+\_\+\+Lapla\+\_\+w} \+: スペクトルデータにラプラシアンを作用させる
\item \hyperlink{namespacew__module__deriv_aa475b204ff39d1a9d988d1893de349fd}{rn} \+: スペクトルデータのラプラシアンを計算するための係数.
\item \hyperlink{namespacew__module__deriv_aa247a4a281588e92be2ec6a11da04a5b}{w\+\_\+\+Lapla\+Inv\+\_\+w} \+: スペクトルデータにラプラシアンの逆変換を作用させる
\item \hyperlink{namespacew__module__deriv_a1676aabd52064c64b765e84b140b8a82}{w\+\_\+\+D\+Lon\+\_\+w} \+: スペクトルデータに経度微分 $∂/∂λ$を作用させる
\item \hyperlink{namespacew__module__deriv_a88efad865d2eb89d3a069ac9c45b3f9b}{xy\+\_\+\+Grad\+Lon\+\_\+w} \+: スペクトルデータに勾配型経度微分 $1/\cosφ・∂/∂λ$を作用させる
\item \hyperlink{namespacew__module__deriv_a79fbe256e0d5145fc01ca38435860741}{xy\+\_\+\+Grad\+Lat\+\_\+w} \+: スペクトルデータに勾配型緯度微分 $∂/∂φ$を作用させる
\item \hyperlink{namespacew__module__deriv_a4ad86b3db507fb0bff454e007e515ba0}{w\+\_\+\+Div\+Lon\+\_\+xy} \+: 格子データに発散型経度微分 $1/\cosφ・∂/∂λ$を作用させる
\item \hyperlink{namespacew__module__deriv_ab13eefdf5b81060ce83e3065b1fea953}{w\+\_\+\+Div\+Lat\+\_\+xy} \+: 格子データに発散型緯度微分 $1/\cosφ・∂(g \cosφ)/∂φ$を作用させる
\item \hyperlink{namespacew__module__deriv_a4ac727e6389247d61f99e5358aebf069}{w\+\_\+\+Div\+\_\+xy\+\_\+xy} \+: ベクトル成分である 2 つの格子データに発散を作用させる
\item \hyperlink{namespacew__module__deriv_a7a61436e763799d6fb0b980e84e1b232}{w\+\_\+\+Jacobian\+\_\+w\+\_\+w} \+: 2 つのスペクトルデータからヤコビアンを計算する
\end{DoxyItemize}

$(λ,μ)$ 座標 $(μ=\sinφ)$ での微分


\begin{DoxyItemize}
\item \hyperlink{namespacew__module__deriv_aafdccb5625e83fc74e07f28c29c663ae}{xy\+\_\+\+Grad\+Lambda\+\_\+w} \+: スペクトルデータに勾配型経度微分 $∂/∂λ$を作用させる
\item \hyperlink{namespacew__module__deriv_a602369ddd0f5e6b89163ac558e0d71b1}{xy\+\_\+\+Grad\+Mu\+\_\+w} \+: スペクトルデータに勾配型緯度微分 $(1-μ^2)∂/∂μ$を作用させる
\item \hyperlink{namespacew__module__deriv_a6dc1edaff2ff00ce89f907370a0b0292}{w\+\_\+\+Div\+Lambda\+\_\+xy} \+: 格子データに発散型経度微分 $1/(1-μ^2)・∂/∂λ$を作用させる
\item \hyperlink{namespacew__module__deriv_a984edeb7ba46736df32910567e30ed8e}{w\+\_\+\+Div\+Mu\+\_\+xy} \+: 格子データに発散型緯度微分 $∂/∂μ$を作用させる
\end{DoxyItemize}

\paragraph*{補間}

\hyperlink{namespacew__module__interpolate}{w\+\_\+module\+\_\+interpolate} にて定義


\begin{DoxyItemize}
\item Interpolate\+\_\+w \+: スペクトルデータから任意の点での値を求める.
\end{DoxyItemize}

\paragraph*{積分・平均}

\hyperlink{namespacew__module__integral}{w\+\_\+module\+\_\+integral} にて定義


\begin{DoxyItemize}
\item \hyperlink{namespacew__module__integral_a681cc02a96b569a528fd0e9cff0a878b}{Int\+Lon\+Lat\+\_\+xy} \+: 2 次元格子点データの全領域積分
\item \hyperlink{namespacew__module__integral_ae81714fa06a9b31798ef10f8ab12a414}{Avr\+Lon\+Lat\+\_\+xy} \+: 2 次元格子点データの全領域平均
\item \hyperlink{namespacew__module__integral_a38a305bdc56997297d5d6f2e46ab38b3}{y\+\_\+\+Int\+Lon\+\_\+xy} \+: 2 次元格子点データの経度方向積分
\item \hyperlink{namespacew__module__integral_ada18b7030e3bf35699bd0af8611a183d}{y\+\_\+\+Avr\+Lon\+\_\+xy} \+: 2 次元格子点データの経度方向平均
\item \hyperlink{namespacew__module__integral_ac7ba110474a4ab162f9b4e7397b2dc91}{Int\+Lon\+\_\+x} \+: 1 次元(\+X)格子点データの経度方向積分
\item \hyperlink{namespacew__module__integral_a6f252f1df33c794abdf86114346bf1c5}{Avr\+Lon\+\_\+x} \+: 1 次元(\+X)格子点データの経度方向平均
\item \hyperlink{namespacew__module__integral_a63f03bca68fcc1d72b665e5a1ef12b68}{x\+\_\+\+Int\+Lat\+\_\+xy} \+: 2 次元格子点データの緯度方向積分
\item \hyperlink{namespacew__module__integral_aec162c0bf27dd164bc686a38e9c084cb}{x\+\_\+\+Avr\+Lat\+\_\+xy} \+: 2 次元格子点データの緯度方向平均
\item \hyperlink{namespacew__module__integral_aee775f5db056a1cbe1e0d12459c9f1ea}{Int\+Lat\+\_\+y} \+: 1 次元(\+Y)格子点データの緯度方向積分
\item \hyperlink{namespacew__module__integral_a6d7a3294fdda72a4cc589a6630414d74}{Avr\+Lat\+\_\+y} \+: 1 次元(\+Y)格子点データの緯度方向平均
\end{DoxyItemize}

\paragraph*{スペクトル解析}

\hyperlink{namespacew__module__spectrum}{w\+\_\+module\+\_\+spectrum} にて定義


\begin{DoxyItemize}
\item \hyperlink{namespacew__module__spectrum_a549029558183c587d43bdb0a3592df12}{nm\+\_\+\+Energy\+From\+Streamfunc\+\_\+w} \+: 流線関数からエネルギースペクトルを計算する(水平全波数 n, 帯状波数 m 空間)
\item \hyperlink{namespacew__module__spectrum_a7619d87ba2a8eb0ba4c9cb1bfe51ace8}{n\+\_\+\+Energy\+From\+Streamfunc\+\_\+w} \+: 流線関数からエネルギースペクトルを計算する(水平全波数 n 空間)
\item \hyperlink{namespacew__module__spectrum_abaac50d008bc5f3a7e31ac39a8955929}{nm\+\_\+\+Enstrophy\+From\+Streamfunc\+\_\+w} \+: 流線関数からエンストロフィースペクトルを計算する (水平全波数 n, 帯状波数 m 空間)
\item \hyperlink{namespacew__module__spectrum_a2e2a246f41310d2bd1ae0db84d73b499}{n\+\_\+\+Enstrophy\+From\+Streamfunc\+\_\+w} \+: 流線関数からエンストロフィースペクトルを計算する (水平全波数 n 空間)
\item \hyperlink{namespacew__module__spectrum_a29b88d42fa75011df629c02dcd7422d4}{w\+\_\+spectrum\+\_\+\+V\+Miss} \+: 欠損値 
\end{DoxyItemize}

\subsection{Function/\+Subroutine Documentation}
\mbox{\Hypertarget{namespacew__module_a8fa63b83c6aee6f23c4560de239475e6}\label{namespacew__module_a8fa63b83c6aee6f23c4560de239475e6}} 
\index{w\+\_\+module@{w\+\_\+module}!w\+\_\+finalize@{w\+\_\+finalize}}
\index{w\+\_\+finalize@{w\+\_\+finalize}!w\+\_\+module@{w\+\_\+module}}
\subsubsection{\texorpdfstring{w\+\_\+finalize()}{w\_finalize()}}
{\footnotesize\ttfamily subroutine, public w\+\_\+module\+::w\+\_\+finalize (\begin{DoxyParamCaption}{ }\end{DoxyParamCaption})}



モジュールの終了処理(割り付け配列の解放)をおこなう. 

\begin{DoxyWarning}{Warning}
解像度を変更する際にはこのサブルーチンを呼んで終了処理を おこなったのちに再度 w\+\_\+\+Initial で初期設定しなければ ならない. 
\end{DoxyWarning}


Definition at line 419 of file w\+\_\+module.\+f90.



References w\+\_\+module\+\_\+base\+::w\+\_\+base\+\_\+finalize(), and w\+\_\+module\+\_\+deriv\+::w\+\_\+deriv\+\_\+finalize().


\begin{DoxyCode}
419     \textcolor{keyword}{call }w\_base\_finalize
420     \textcolor{keyword}{call }w\_deriv\_finalize
421 
422     \textcolor{keyword}{call }messagenotify(\textcolor{stringliteral}{'M'},\textcolor{stringliteral}{'w\_Finalize'},&
423       \textcolor{stringliteral}{'w\_module (2019/10/02) is finalized'})
424 
\end{DoxyCode}
Here is the call graph for this function\+:
\nopagebreak
\begin{figure}[H]
\begin{center}
\leavevmode
\includegraphics[width=350pt]{namespacew__module_a8fa63b83c6aee6f23c4560de239475e6_cgraph}
\end{center}
\end{figure}
\mbox{\Hypertarget{namespacew__module_a34cf65ed64bfeec316d31e8764e8ffc0}\label{namespacew__module_a34cf65ed64bfeec316d31e8764e8ffc0}} 
\index{w\+\_\+module@{w\+\_\+module}!w\+\_\+initial@{w\+\_\+initial}}
\index{w\+\_\+initial@{w\+\_\+initial}!w\+\_\+module@{w\+\_\+module}}
\subsubsection{\texorpdfstring{w\+\_\+initial()}{w\_initial()}}
{\footnotesize\ttfamily subroutine, public w\+\_\+module\+::w\+\_\+initial (\begin{DoxyParamCaption}\item[{integer, intent(in)}]{n\+\_\+in,  }\item[{integer, intent(in)}]{i\+\_\+in,  }\item[{integer, intent(in)}]{j\+\_\+in,  }\item[{integer, intent(in), optional}]{np\+\_\+in }\end{DoxyParamCaption})}



スペクトル変換の格子点数, 波数および O\+P\+E\+N\+MP 使用時の 最大スレッド数を設定する. 

他の関数を呼ぶ前に, 最初にこのサブルーチンを呼んで初期設定を しなければならない.

np\+\_\+in に 1 より大きな値を指定すれば I\+S\+P\+A\+CK の球面調和函数変換 O\+P\+E\+N\+MP 並列計算ルーチンが用いられる. 並列計算を実行するには, 実行時に環境変数 O\+M\+P\+\_\+\+N\+U\+M\+\_\+\+T\+H\+R\+E\+A\+DS を np\+\_\+in 以下の数字に設定する等の システムに応じた準備が必要となる.

np\+\_\+in に 1 より大きな値を指定しなければ並列計算ルーチンは呼ばれない.


\begin{DoxyParams}[1]{Parameters}
\mbox{\tt in}  & {\em i\+\_\+in} & 格子点数(東西)\\
\hline
\mbox{\tt in}  & {\em j\+\_\+in} & 格子点数(南北)\\
\hline
\mbox{\tt in}  & {\em n\+\_\+in} & 切断波数の設定\\
\hline
\mbox{\tt in}  & {\em np\+\_\+in} & O\+P\+E\+N\+MP での最大スレッド数 \\
\hline
\end{DoxyParams}


Definition at line 387 of file w\+\_\+module.\+f90.



References w\+\_\+module\+\_\+base\+::w\+\_\+base\+\_\+initial(), and w\+\_\+module\+\_\+deriv\+::w\+\_\+deriv\+\_\+initial().


\begin{DoxyCode}
387 
388     \textcolor{keywordtype}{integer},\textcolor{keywordtype}{intent(in)} :: i\_in
390     \textcolor{keywordtype}{integer},\textcolor{keywordtype}{intent(in)} :: j\_in
392     \textcolor{keywordtype}{integer},\textcolor{keywordtype}{intent(in)} :: n\_in
394     \textcolor{keywordtype}{integer},\textcolor{keywordtype}{intent(in)}, \textcolor{keywordtype}{optional} :: np\_in
395 
396     \textcolor{keywordflow}{if} ( \textcolor{keyword}{present} (np\_in) ) \textcolor{keywordflow}{then}
397       \textcolor{keyword}{call }messagenotify(\textcolor{stringliteral}{'M'},\textcolor{stringliteral}{'w\_Initial'}, &
398         & \textcolor{stringliteral}{'Optional argnument np is dummy.'})
399 \textcolor{keywordflow}{    endif}
400 
401     \textcolor{keyword}{call }w\_base\_initial(n\_in,i\_in,j\_in)
402     \textcolor{keyword}{call }w\_deriv\_initial
403 
404     \textcolor{keyword}{call }messagenotify(\textcolor{stringliteral}{'M'},\textcolor{stringliteral}{'w\_initial'},&
405       \textcolor{stringliteral}{'w\_module (2019/10/02) is initialized'})
406 
\end{DoxyCode}
Here is the call graph for this function\+:
\nopagebreak
\begin{figure}[H]
\begin{center}
\leavevmode
\includegraphics[width=345pt]{namespacew__module_a34cf65ed64bfeec316d31e8764e8ffc0_cgraph}
\end{center}
\end{figure}
