(57)【要約】【課題】定常状態の達成を海域平均運動エネルギーの時系列変化を表示させることによって客観的に判定する、海洋モデル定常状態判定方法、海洋モデル定常状態判定装置及びそのプログラムを提供する。【解決手段】風応力データ生成部１が生成した風応力を入力として海洋環境算出部２で算出した流速、密度に基づき運動エネルギー算出部３で各点の運動エネルギーを算出する。海域平均運動エネルギー算出部４はこれら運動エネルギーを海域全体での平均値を算出する。傾圧ロスビー波伝搬時間算出部５では海域の一方の岸で発生した擾乱が、他方の岸に達するまでの時間を算出する。海域平均運動エネルギー表示部６では海域平均運動エネルギーの時系列変化を、傾圧ロスビー波伝搬時間表示部７では傾圧ロスビー波伝搬時間を表示し、定常状態達成の有無の判断を支援する。【選択図】図１
2005321269.tif 000002
【特許請求の範囲】
【請求項１】
定常状態の風応力データを生成し、
前記定常状態の風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出し、
算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出し、
前記モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出し、
海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出し、
前記海域平均運動エネルギーを表示し、
前記傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示し、
表示された前記海域平均運動エネルギーと前記傾圧ロスビー波伝搬時間とを比較することにより、海洋モデルが定常状態に達したか否かの判定を支援する、海洋モデル定常状態判定方法。
【請求項２】
季節変化する風応力データを生成し、
前記季節変化する風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出し、
算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出し、
前記モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出し、
海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出し、
前記海域平均運動エネルギーを表示し、
前記傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示し、
表示された前記海域平均運動エネルギーと前記傾圧ロスビー波伝搬時間とを比較することにより、海洋モデルが定常状態に達したか否かの判定を支援する、海洋モデル定常状態判定方法。
【請求項３】
定常状態の風応力データを生成する風応力データ生成手段と、
前記定常状態の風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する海洋環境算出手段と、
算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する運動エネルギー算出手段と、
前記モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する海域平均運動エネルギー算出手段と、
海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出する傾圧ロスビー波伝搬時間算出手段と、
前記海域平均運動エネルギーを表示する海域平均運動エネルギー表示手段と、
前記傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示する傾圧ロスビー波伝搬時間表示手段とを有し、
表示された前記海域平均運動エネルギーと前記傾圧ロスビー波伝搬時間とを比較することにより、海洋モデルが定常状態に達したか否かの判定を支援する、海洋モデル定常状態判定装置。
【請求項４】
季節変化する風応力データを生成する風応力データ生成手段と、
前記季節変化する風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する海洋環境算出手段と、
算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する運動エネルギー算出手段と、
前記モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する海域平均運動エネルギー算出手段と、
海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出する傾圧ロスビー波伝搬時間算出手段と、
前記海域平均運動エネルギーを表示する海域平均運動エネルギー表示手段と、
前記傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示する傾圧ロスビー波伝搬時間表示手段とを有し、
表示された前記海域平均運動エネルギーと前記傾圧ロスビー波伝搬時間とを比較することにより、海洋モデルが定常状態に達したか否かの判定を支援する、海洋モデル定常状態判定装置。
【請求項５】
定常状態の風応力データを生成する処理と、
前記定常状態の風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する処理と、
算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する処理と、
前記モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する処理と、
海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出する処理と、
前記海域平均運動エネルギーを表示する処理と、
前記傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示する処理とを、コンピュータに実行させるための海洋モデル定常状態判定プログラム。
【請求項６】
季節変化する風応力データを生成する処理と、
前記季節変化する風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する処理と、
算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する処理と、
前記モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する処理と、
海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出する処理と、
前記海域平均運動エネルギーを表示する処理と、
前記傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示する処理とを、コンピュータに実行させるための海洋モデル定常状態判定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、海洋モデル定常状態判定方法、海洋モデル定常状態判定装置及びそのプログラムに関する。
【背景技術】
【０００２】
従来の技術としては、予報精度の向上に役立つ、データ同化処理で用いる観測データを効率的に収集することが可能な海洋環境計測間隔決定方法及び装置がある（特許文献１参照。）。
【０００３】
その海洋環境計測間隔決定方法は、数値海洋モデルに基づいて予報対象となる海洋環境に関する物理的相関係数を算出し、この算出値を利用してガウス型相関係数分布を仮定することにより、物理的相関距離を算出する。一方で、測定対象となる海域の水深、コリオリパラメータ、重力加速度等の特性に基づいて海洋に発生する擾乱が影響を及ぼす範囲であるロスビーの変形半径を算出する。これら数値海洋モデルによる物理量の予測結果、あるいは相関距離とロスビーの変形半径の大小関係を比較し、より小さい方を選択する。
【０００４】
数値海洋モデルは、プリンストン海洋モデル（P0M : Princeton Ocean Model）（非特許文献１参照。）、モジュラー海洋モデル（M0M : Modular Ocean Model）（非特許文献２参照。）等のモデルであり、これらの数値海洋モデルに基づき、水温、塩分、密度、流速といった海洋環境データを３次元空間、規定の時間について算出し、予報格子データとしてすべての格子上に設定する。算出の原理は、回転地球上の海水の運動方程式、質量保存式、熱量保存式、塩分保存式を数値的に算出するものであり、詳細はMellor(1998)による（非特許文献３参照。）。
【０００５】
また、海洋モデルに定常な外力を入力した場合、モデル海域の東岸に発生した擾乱が西側に伝搬し、西岸に到達すると、海洋モデルが定常に達するという理論がある（非特許文献４参照。）。
【特許文献１】特開２００４−０２８６１６号公報
【非特許文献１】"プリンストン海洋モデル（P0M : Princeton Ocean Model）"、[online]、インターネット＜URL :  HYPERLINK "http://www.aos.princeton.edu/WWWPUBLIC/htdocs.pom/" http://www.aos.princeton.edu/WWWPUBLIC/htdocs.pom/＞
【非特許文献２】"モジュラー海洋モデル（M0M : Modular Ocean Model）" 、[online]、インターネット＜URL : HYPERLINK "http://www.gfdl.gov/~smg/MOM/MOM.html" http://www.gfdl.gov/~smg/MOM/MOM.html＞
【非特許文献３】Mellor, G. L., Users guide for a three-dimensional, primitive equation, numerical ocean model (July 1998 version), 41 pp., Prog. in Atmos. and Ocean. Sci, Princeton University, 1998.
【非特許文献４】関根義彦著「海洋物理学概論」，成山堂書店，1995年
【発明の開示】
【発明が解決しようとする課題】
【０００６】
従来の技術においては、海洋モデル確立において必要なスピンアップ過程において、流れの場が定常状態に達したか否かを判定する手段が、流速場の監視をマニュアルで行うか、傾圧ロスビー波伝搬時間の概算理論値で判断するしか方法がなかった。
【０００７】
即ち、海洋モデルを構築する際に、実際に近い風応力を外力として入力し、モデル計算を運用する前に、静止状態から定常状態の海洋環境を作り出す必要がある。
【０００８】
その際、一定の風応力を入力して海洋の流れ場が一定となる状態まで海洋モデルを実行するか、一定の季節変動を示す風応力を入力して海洋の流れ場が一定の季節変動を示す状態まで海洋モデルを実行する必要がある。
【０００９】
このように海洋の流れ場が定常状態となったことを判定する手段は、開発担当者が流れ場を適当に確認するか、傾圧ロスビー波がモデル海洋の東岸で励起されて西岸に到達する時間を目安に大まかに決めており、系統的に判断する手段がこれまでなく、客観的な判定ができないという問題があった。
【００１０】
本発明の目的は、海洋モデル構築における定常状態を作り出す過程において、定常状態の達成を海域平均運動エネルギーの時系列変化を表示させることによって客観的に判定する、海洋モデル定常状態判定方法、海洋モデル定常状態判定装置及びそのプログラムを提供することにある。
【課題を解決するための手段】
【００１１】
本発明の海洋モデル定常状態判定方法は、
定常状態の風応力データを生成し、定常状態の風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出し、算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出し、モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出し、海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出し、海域平均運動エネルギーを表示し、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示し、表示された海域平均運動エネルギーと傾圧ロスビー波伝搬時間とを比較することにより、海洋モデルが定常状態に達したか否かの判定を支援する。
【００１２】
本発明の海洋モデル定常状態判定方法は、
季節変化する風応力データを生成し、季節変化する風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出し、算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出し、モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出し、海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出し、海域平均運動エネルギーを表示し、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示し、表示された海域平均運動エネルギーと傾圧ロスビー波伝搬時間とを比較することにより、海洋モデルが定常状態に達したか否かの判定を支援する。
【００１３】
本発明の海洋モデル定常状態判定装置は、
定常状態の風応力データを生成する風応力データ生成手段と、定常状態の風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する海洋環境算出手段と、算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する運動エネルギー算出手段と、モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する海域平均運動エネルギー算出手段と、海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出する傾圧ロスビー波伝搬時間算出手段と、海域平均運動エネルギーを表示する海域平均運動エネルギー表示手段と、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示する傾圧ロスビー波伝搬時間表示手段とを有し、表示された海域平均運動エネルギーと傾圧ロスビー波伝搬時間とを比較することにより、海洋モデルが定常状態に達したか否かの判定を支援する。
【００１４】
本発明の海洋モデル定常状態判定装置は、
季節変化する風応力データを生成する風応力データ生成手段と、季節変化する風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する海洋環境算出手段と、算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する運動エネルギー算出手段と、モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する海域平均運動エネルギー算出手段と、海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出する傾圧ロスビー波伝搬時間算出手段と、海域平均運動エネルギーを表示する海域平均運動エネルギー表示手段と、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示する傾圧ロスビー波伝搬時間表示手段とを有し、表示された海域平均運動エネルギーと傾圧ロスビー波伝搬時間とを比較することにより、海洋モデルが定常状態に達したか否かの判定を支援する。
【００１５】
本発明の海洋モデル定常状態判定プログラムは、
定常状態の風応力データを生成する処理と、定常状態の風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する処理と、算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する処理と、モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する処理と、海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出する処理と、海域平均運動エネルギーを表示する処理と、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示する処理とを、コンピュータに実行させる。
【００１６】
本発明の海洋モデル定常状態判定プログラムは、
季節変化する風応力データを生成する処理と、季節変化する風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する処理と、算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する処理と、モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する処理と、海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を算出する処理と、海域平均運動エネルギーを表示する処理と、傾圧ロスビー波のモデル海域の一方の岸から他方の岸への伝搬時間を表示する処理とを、コンピュータに実行させる。
【００１７】
本発明は、実海洋環境を精度よく予報するための海洋モデルを確立するために行わねばならないスピンアップの過程で、モデルを起動してからモデル海域の一方の岸で発生する擾乱が傾圧ロスビー波として他方の岸に到達することによって定常状態が実現しているか否かを運動エネルギーを監視すると同時に簡単な理論計算で求められる傾圧ロスビー波の伝搬時間を比較することによって、判断を支援する方法を提供する。
【００１８】
即ち、本発明による海洋モデル定常状態判定装置は、すべての入力源となる風応力データを作成して入力する手段と、入力された風応力に基づいて水温、塩分、密度、流速等の海洋環境の時空間データを算出する手段を有する。
【００１９】
上記手段により算出された海洋環境の時空間データに基づき、任意の点における運動エネルギーを算出する手段と、算出された運動エネルギーに基づき、モデル海域全体での海域平均運動エネルギーを算出する手段と、定常状態判定の理論的目安となる傾圧ロスビー波の海域伝搬時間を算出する手段と、海域平均運動エネルギーの時系列変化を表示する手段と、傾圧ロスビー波の伝搬時間を表示する手段とを有する。
【００２０】
本発明は、海洋モデルに定常な外力を入力した場合、モデル海域の東岸に発生した擾乱が西側に伝搬し、西岸に到達すると、海洋モデルが定常に達するという理論（非特許文献４）を適用する。
【００２１】
本発明では、一定値又は季節変動する風応力データを入力し、数値海洋モデルの算出結果に基づき、各点における運動エネルギーを算出し、この運動エネルギーに基づき海域平均運動エネルギーを算出する。
【００２２】
一方で、理論的目安となる傾圧ロスビー波がモデル海域を一方の岸から他方の岸に伝搬する時間を算出する。
【００２３】
これらの海域平均運動エネルギーの時系列変化と傾圧ロスビー波の伝搬時間を表示する。
【００２４】
これらを表示して比較することにより、海洋モデルが定常状態に達したか否かを判断するオペレータを支援するものである。
【発明の効果】
【００２５】
従来、海洋モデル確立にあたって必要なスピンアップの過程で、定常状態を判定する系統的な方法は存在せず、オペレータの主観的な判断又は傾圧ロスビー波の伝播時間の概算値を参考にする方法に頼らざるを得なかった。
【００２６】
本発明により、海域平均運動エネルギーの時系列変化という客観的な数値を評価し、これを傾圧ロスビー波伝搬時間の概算値と比較することにより、定常状態達成の有無の系統的な判断が可能となる。
【発明を実施するための最良の形態】
【００２７】
図１は、本発明を適用した海洋モデル定常状態判定システムの構成である。
【００２８】
風応力データ生成部１、海洋環境算出部２、運動エネルギー算出部３、海域平均運動エネルギー算出部４、傾圧ロスビー波伝搬時間算出部５、海域平均運動エネルギー表示部６、傾圧ロスビー波伝搬時間表示部７から構成される。
【００２９】
風応力データ生成部１は、オペレータによって一定値か季節変動値かを選択し、対応する風応力データを自動生成する。
【００３０】
海洋環境算出部２は、風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速等の海洋環境の時空間データを算出する。
【００３１】
運動エネルギー算出部３は、海洋環境算出部２で算出された密度、流速と規定の各深度での層厚に基づいて各点の運動エネルギーを算出する。
【００３２】
海域平均運動エネルギー算出部4は、運動エネルギー算出部３で算出された各点の運動エネルギーをモデル海域全域で平均することにより、海域平均運動エネルギーを算出する。
【００３３】
傾圧ロスビー波伝搬時間算出部５は、モデル海域の緯度に応じて、海域の東西の長さ、コリオリパラメータの経度方向の勾配、当該緯度における平均水深、重力加速度に基づき、モデル海域の一方の岸（例えば東岸）で発生した擾乱が他方の岸（例えば西岸）に到達する所要時間を算出する。
【００３４】
海域平均運動エネルギー表示部６は、海域平均運動エネルギー算出部５で算出した海域平均運動エネルギーの時系列変化をオペレータが指定した時間間隔で表示する。
【００３５】
傾圧ロスビー波伝搬時間表示部７は、傾圧ロスビー波伝搬時間算出部５で算出した伝搬時間を表示する。
【００３６】
これら２つの表示情報により、定常状態に達したか否かを判断する材料をオペレータに提供する。
【００３７】
次に、本発明の動作を、図１を用いて詳細に説明する。
【００３８】
風応力データ生成部１は、オペレータの入力により、一定値とするか季節変動させるかを選択し、この入力に応じた風応力の時空間データを生成する。
【００３９】
ここで、オペレータは計算実行を行うべき計算継続時間を入力する。
【００４０】
海洋環境算出部２は、非特許文献１の"プリンストン海洋モデル（P0M : Princeton Ocean Model）"、非特許文献２の"モジュラー海洋モデル（M0M : Modular Ocean Model）"等の数値海洋モデルに基づき、水温、塩分、密度、流速といった海洋環境データを３次元空間、規定の時間について算出し、予報格子データとしてすべての格子上に設定する。算出の原理は、回転地球上の海水の運動方程式、質量保存式、熱量保存式、塩分保存式を数値的に算出するものであり、詳細は非特許文献３のMellor(1998)による。
【００４１】
運動エネルギー算出部３は、海洋環境算出部２で算出された海洋環境のうち、密度、流速及び各層の厚さを用いて運動エネルギーを、特定のi点について式（１）を用いて算出する。算出に使用するパラメータの概略図は図２のとおりである。
【００４２】
【数１】
2005321269.tif 000003
【００４３】
ここで、iは海域内の特定の位置の番号を表し、KEiは特定の位置Iにおける運動エネルギーを表す。ρjは各層の密度、hjは各層の厚さ、vjは各層の流速ベクトルを表す。|  |はベクトルの絶対値を表す。また、jは各層を示す番号、mは全層の数を表す。
【００４４】
海域平均運動エネルギー算出部４は、運動エネルギー算出部３で算出されたi点における運動エネルギーの海域平均値を、式（２）を用いて海域平均運動エネルギーとして算出する。
【００４５】
【数２】
2005321269.tif 000004
【００４６】
ここで、
【００４７】
【数３】
2005321269.tif 000005
【００４８】
は運動エネルギーの海域平均値、Nは海域の全点数を表す。
【００４９】
傾圧ロスビー波伝搬時間算出部５は、モデル海域の一方の岸（例えば東岸）で発生し、他方の岸（例えば西岸）に到達する擾乱を傾圧ロスビー波としてその伝搬時間を式（３）に基づき算出する。
【００５０】
【数４】
2005321269.tif 000006
【００５１】
ここで、
【００５２】
【数５】
2005321269.tif 000007
【００５３】
における傾圧ロスビー波の伝搬時間、Lはモデル海域の一方の岸から他方の岸への方向（例えば東西方向）の長さ、ＣRは傾圧ロスビー波の位相速度、βはコリオリパラメータ
【００５４】
【数６】
2005321269.tif 000008
【００５５】
の経度方向の勾配
【００５６】
【数７】
2005321269.tif 000009
【００５７】
はロスビーの変形半径
【００５８】
【数８】
2005321269.tif 000010
【００５９】
における平均水深値を各々表す。
【００６０】
海域平均運動エネルギー表示部６は、海域平均運動エネルギー算出部４で算出された海域平均運動エネルギーを、オペレータが指定した期間毎にその時系列変化を表示する。
【００６１】
傾圧ロスビー波伝搬時間表示部７は、傾圧ロスビー波伝搬時間算出部５で算出した傾圧ロスビー波伝搬時間を表示する。
【００６２】
これらによって、傾圧ロスビー波伝搬時間を参考に、海域平均運動エネルギーが定常状態に達したか否かをオペレータが判断する。
【００６３】
入力する風応力が一定値の場合、海域平均運動エネルギーが一定値に達した場合に海洋モデルが定常状態に達したと判断する。
【００６４】
また、入力する風応力が季節変動する場合、その応答となる海域平均運動エネルギーが毎年同一の季節変化を示すようになったら、海洋モデルが定常状態に達したと判断する。
【００６５】
また、本発明の海洋モデル定常状態判定プログラムにより、定常状態の又は季節変化する風応力データを生成する処理と、定常状態の又は季節変化する風応力データの入力により、数値海洋モデルに基づき、水温、塩分、密度、流速の海洋環境の時空間データを算出する処理と、算出される密度、流速に基づき、モデル海域各点における運動エネルギーを算出する処理と、モデル海域各点における運動エネルギーに関し、数値海洋モデルのモデル海域全域において平均値を算出する処理と、海域の平均水深に基づき、傾圧ロスビー波のモデル海域の一方の岸から他方の岸（例えば東岸から西岸）への伝搬時間を算出する処理と、海域平均運動エネルギーを表示する処理と、傾圧ロスビー波のモデル海域の一方の岸から他方の岸（例えば東岸から西岸）への伝搬時間を表示する処理とを、コンピュータに実行させることができる。
【図面の簡単な説明】
【００６６】
【図１】本発明の実施形態の構成ブロック図である。
【図２】運動エネルギー算出に必要な各層のパラメータの概略図である。
【符号の説明】
【００６７】
１    風応力データ生成部
２    海洋環境算出部
３    運動エネルギー算出部
４    海域平均運動エネルギー算出部
５    傾圧ロスビー波伝搬時間算出部
６    海域平均運動エネルギー表示部
７    傾圧ロスビー波伝搬時間表示部
【図１】
2005321269.tif 000011
【図２】
2005321269.tif 000012