地衡流速(1000dbを無流面と仮定)
地衡流計算用プログラム
geocal.f
(下にリンクあり)
サンプルデータ
S.txt
T.txt
(下にリンクあり)
注意:データがない点は-99.99という値が入っている
コンパイルと実行例
am@am-PC ~/dynamic_calc
$ gfortran geocal.f -o geocal
am@am-PC ~/dynamic_calc
$ ./geocal
Done reading.
Done density.
Done computing dynamic depth.
Done geostrophic velocity.
描画用スクリプト (GMT使用)
plrho.sh : 密度の鉛直断面図をかく
plv.sh : 地衡流速の鉛直断面図をかく
(下にリンクあり)
オンラインできる地衡流計算のサイト
http://gyre.umeoce.maine.edu/physicalocean/Tomczak/Utilities/geostrophy/index.html