=begin JA = Fortran 90 ライブラリ STPK デモ ここでは, STPK を用いたデモファイル群, demo の結果を示します. 順次追加します. ((:

:)) デモの一部を動作させるためには, (()) および, (()) がインストールされている必要があります. ((:
:)) また, 計算結果は netcdf 形式で出力されますので, この形式をサポートしている可視化ツールを用いれば, 結果の図化が可能です. ((:
:)) 本サンプルプログラムは, (()) がインストールされていれば, 可視化ツールもサンプルと同じディレクトリに置いてありますので, ご利用下さい. ((:
:)) 使い方は同梱の PDF マニュアルを見ていただければ幸いです. == poisson ポアソン方程式をヤコビ法を用いて解く. 強制項は領域中心に有限領域をもつテーブル関数型. * ((<結果|URL:./poisson.png>)) * 使用サブルーチン : Ellip_Jacobi (use Ellip_Slv) == cov とある分布データを用いて, 相関係数を計算し, 回帰直線をグラフ化する * ((<結果|URL:./cov.png>)) * 使用サブルーチン : Reg_Line_1d, Cor_Coe_1d (use Statistics) == wind 解析的な気圧場を与えて, 地衡風, 渦度, 流線を計算します. * ((<気圧場|URL:./wind_pres.png>)) * ((<地衡風場|URL:./wind_velo.png>)) * 使用サブルーチン : grad_1d (use analy) * ((<渦度場|URL:./wind_rot.png>)) * 使用サブルーチン : curl (use analy) * ((<流線場|URL:./wind_SL.png>)) * 使用サブルーチン : Forward_Trajectory (use Trajectory) == matrix_test no_draw ディレクトリ内に同梱. data.dat に格納されているデータを用いて, * 連立方程式の解 * 使用サブルーチン (すべて use Algebra) : * gausss, LU_devs, Gau_Sei, Jacobi_algebra, SOR_Gau_Sei, SOR_Jacobi_algebra * 逆行列 * 使用サブルーチン : invert_mat (use Algebra) * 固有値 * 使用サブルーチン : eigenvalue_power (use Algebra) * 固有値については, べき乗法を用いて最大固有値と固有ベクトルを計算する. を求める. 実行すると, 1 桁の整数を入力するよう求められる. それぞれの数字が何に対応するかは, 実行時に表示されるので, それにしたがって入力. == adjust 1 次元線形浅水方程式系モデルから地衡流調節と重力波放射をシミュレーションするモデル. * 使用サブルーチン (use Derivation) * grad_1d * 初期値 * 静止状態で流体の高さに不連続を与えておく. * ((<初期の流体の高さ|URL:./adjust_1.png>)) * 設定パラメータ * ロスビー変形半径 : 300 m. * 重力波の伝播速度 : 30 m/s. * 計算結果 * ((<流体の高さの時系列図|URL:./adjust_2.png>)) * 黒実線がロスビー変形半径. * ((<地衡流の時系列図|URL:./adjust_3.png>)) * 定常状態での((<流体の高さ|URL:./adjust_4.png>)), 調節で引き起こされた((<地衡流|URL:./adjust_5.png>)) * 計算初期での顕著な重力波の((<外向き放射|URL:./adjust_6.png>)) * ロスビー変形半径と重力波の伝播速度は理論値と整合的. == Karman 2 次元流内におけるカルマン渦シミュレーションモデル. * ((<詳細|URL:./Karman/index.html>)) * 使用サブルーチン : Ellip_Jacobi (use Ellip_Slv) == WK04 水平鉛直 2 次元流内における冷気外出流シミュレーションモデル. * ((<詳細|URL:./WK04/index.html>)) * 使用サブルーチン : Ellip_Jacobi (use Ellip_Slv) == SEQ Pendergrass and Willoughby (2009) の 2 次元軸対称台風モデル. * ((<詳細|URL:./SEQ/index.html>)) * 使用サブルーチン : Ellip_Gausei (use Ellip_Slv) == NM01 浅水系における渦ロスビー波伝播シミュレーション用数値モデル. 軸対称流に埋め込まれた非軸対称成分の時間発展を計算する. 接線方向はフーリエスペクトル, 動径方向は中心差分のハイブリッドスペクトルモデルである. その様から, ((::))洗濯機モデル((::))と作者は命名. * ((<詳細|URL:./NM01/index.html>)) * 使用サブルーチン : ffttp_1d (use ffts), grad_* (use derivation) == BAROVOR 順圧非発散渦度モデルを空間方向にスペクトル法を用いて数値積分する. * ((<詳細|URL:./BAROVOR/index.html>)) * 使用サブルーチン : ffttp_2d (use ffts) ((:

:)) ((<1 つ上に戻る|URL:../index.html>)) =end JA =begin EN = Demo Some results, which are got by using programs in "demo" directory, are shown in this page. ((:

:)) To use these programs, you must install some libraries : (()), (()). ((:
:)) Some results produced by programs are output as NetCDF format. You can visualized these results to use some visualized software. ((:
:)) Tools for DCL are included in "demo" directory. Then, if you install (()), you can use the tools. Please see the document how to use these tools. == poisson The program to solve a two-dimensional Poisson equation by the Jacobi iteration method. A forcing in the Poisson equation is given by a shape of a circle. * (()) * using subroutine : Ellip_Jacobi (use Ellip_Slv) == cov The program to calculate correlation coefficient in a distribution and draw the result. * (()) * using subroutine : Reg_Line_1d, Cor_Coe_1d (use Statistics) == wind The program to calculate a geostrophic wind, vorticity and stream line for an analytic pressure distribution. * (()) * (()) * using subroutine : grad_1d (use analy) * (()) * using subroutine : curl (use analy) * (()) * using subroutine : Forward_Trajectory (use Trajectory) == matrix_test By using the data in the "data.dat", the program to * solve a simultaneous equation * using subroutine (use Algebra) : * gausss, LU_devs, Gau_Sei, Jacobi_algebra, SOR_Gau_Sei, SOR_Jacobi_algebra * invert a matrix * using subroutine : invert_mat (use Algebra) * solve an eigenvalue problem * using subroutine : eigenvalue_power (use Algebra) * the program calculates the maximum eigenvalue and the eigenvector correponding to the eigenvalue by using the power method. When you run this program, select one from some single digits. The digits correspond to some calculating methods. == adjust A one-dimensional, linear and shallow-water model to simulate processes of geostropic adjustment and radiation of gravity wave. * using subroutine (use Derivation) * grad_1d * Initial condition * At rest, a discontinuous distribution for the height of the fluid, whose state is unstable. * (()) * some parameters * Rossby deformation radius : 300 m. * propagating speed of gravity wave : 30 m/s. * Results * ((

:)) (()) =end EN