SWAN とは 入力 出力 風データ地形データ SWAN で計算 波高周期波向

Transcription

1 平成 22 年 10 月 15 日作成 SWAN を使用した 波浪推算の実行手順 長岡技術科学大学建設工学専攻 水工学研究室 松永由希子

2 SWAN とは 入力 出力 風データ地形データ SWAN で計算 波高周期波向

3 大まかな手順 風データの変換 DAT ファイルの中身を見る 必要なソフトなどのインストール 地形データ風データの準備 SWAN の計算をする準備 計算して結果をみる SWAN GrADS など swn ファイルの設定など

4 ソフトのインストールの前に まず, SWAN は,Linux,Windows 両バージョンあり Windows 版は重たいので,Linux 版を使用 GrADS も同様です というわけで, 自分のパソコンから Linux にリモートログイン操作 ファイル転送をするためのソフトなどをインストールしてください.

5 Linux にログインして操作 別の場所にある Linux さん 自分のパソコン 遠隔操作 ログイン ネットワークでつながってます

6 Linux へのリモートログイン操作ファイル転送のソフト 使用しているソフト リモートログイン 操作 Tera Term ファイル転送 FFFTP どっちもフリーソフトです 他にも, 色々ソフト 方法があると思います. ネットで調べたら沢山出てきます. コマンドプロンプトでも出来ます.

7 Linux にログイン 1 今回は,Tera Term を使って estuary5 という Linux のサーバに SSH で接続します. OK を押して次に行きましょう 通信内容が暗号化される方法です

8 Linux にログイン 2 Linux にログインするときにはユーザ名とパスワードが必要です. この場合, ユーザ名 パスワード *** OK を押せばログイン完了

9 Linux の操作 1 無事ログインできると右のようなウィンドウになります 書いてあること 最後にログインした日時 日付 パソコン名ようこそターボリナックスへ [ ホスト名カレントディレクトリ ] 一般ユーザ = 今いる場所 拡大 ここに文字を打ち込んで,Enter =Windows でアイコンクリックのようなもの

10 Linux の仕組み 簡単に / root ディレクトリ : 通常は入れません home bin etc suiko ログインするとここに到着します :~ yukiko このときは, /home/suiko/yukiko と / で区切って階層を表現します

11 Linux の操作 2 よく使うやつとか ディレクトリ内にあるファイルの確認 :ls カレントディレクトリの変更 :cd ディレクトリの作成 :mkdir ファイルのコピー :cp ファイルの移動 :mv ファイルの削除 :rm

12 Linux の操作 3 ls コマンド ( ディレクトリ内にあるファイルやディレクトリの確認 ) えんたー ls 若しくは dir と打ってみます / home bin etc suiko yukiko home/suiko(=~) の中に入っているディレクトリやファイルが表示されます. すなわちここ ディレクトリ = フォルダ こうやって Linux を操作していきます

13 Linux の操作 4 cd コマンド ( カレントディレクトリの変更 ) cd[ 半角スペース ] ディレクトリ名 えんたー [ ホスト名カレントディレクトリ ] / home bin etc suiko ここに移動しました ls コマンドで /home/suiko/yukiko のディレクトリの中身がわかりますね. yukiko cd../ でひとつ上の階層へ移動 cd ~/ でホームディレクトリに移動 cd../../ で二つ上のディレクトリへ移動

14 ファイル転送ソフト FFFTP ドラッグ & ドロップで, Windows,Linux 上のファイルやフォルダをやり取りできちゃうソフト これで Linux 上のディレクトリやファイルのコピーや削除も出来ます. Windows 側 Linux 側 設定はこんな感じ

15 SWAN と GrADS の ダウンロードとインストール Windows 上でダウンロードして Linux 上でインストールします SWAN GrADS /home/suiko/yukiko/ppt のディレクトリの中にインストールします

16 Windows 上でダウンロードしたファイルを FFFTP を使って転送 デスクトップに保存したファイルをドラッグ ドロップするとこんな風に Linux 上に転送できます GrADSも同様ですちなみに, このSWANのソフトのファイル (swan4041ab.tgz) は, 拡張子.tgzの圧縮ファイルです. tarコマンドでまとめられた拡張子.tarファイル( 無圧縮 ) をGZIP 形式 (.gz) によって圧縮したファイルらしいです.

17 解凍します Tera Term を使用してます swan4041ab.tgz の解凍を SWAN の Implementation Manual に書いてある通りにやります tar xzf:.tgzファイルを解凍するコマンドもしくは ppt]$ ls tar xvzfでもいいです. swan4041ab.tgz ppt]$ tar xzf swan4041ab.tgz ppt]$ ls swan4041ab/ swan4041ab.tgz 解凍されたディレクトリ GrADS も同様に,tar コマンドで grads-2.0.a7.1-bin-i686-pc-linux-gnu.tar.gz を解凍. grads-2.0a7.1 が生成されます.

18 風データの準備 GRIB ファイルからの抜き出しと バイナリデータ化 Fortran プログラムでテキストデータ化

19 長岡技術科学大学 使用する風データ 風データを SWAN で使用するために,GrADS を使用します. 今回使用する風データは, 気象業務支援センターが提供しているメソ客観解析データ (MSM:Mesoscale Spectral Model) テキストエディタやエクセルでは開けません!!!!! 形式 :GRIB (.grb) 解析時刻 :00,03,06,09,12,15,18,21 UTC 領域 : 日本周辺域 ( 約 3000km 四方領域 ( 北緯 20 ~52 東経 119 ~155 )) 格子系 : ランベルト正角円錐座標 10Km 格子 ( ) 地上物理量 : 海面更正気圧 (Pmsl,Pa) 風 (UV,m/s) 気温 (T,K) 相対湿度 (RH,%) P 面物理量 :20 層 ; 高度 (Z) 風 (UV) 気温 (T) 11 層 ; 相対湿度 (RH) MANAL_ grb 2008 年 4 月 1 日 00 時のデータが1つのGRIBファイルに入っています. その中に, 上記のようなデータたちが計 96 個入っています. その中から必要な風データを抜き出す!!!!!!

20 そもそも GrADS とは The Grid Analysis and Display System の略格子状に配列された 4 次元データ (x,y,z,t) を 2 次元描画するツール Tera Term から Linux 版を使う場合は, グラフィック表示が出来ないので, 単なるデータ抜き出し 変換のソフトとして使います. 今回は,1 ヶ月分の風データを u,v 方向それぞれにバイナリデータ化します

21 GrADS の起動と終了 /home/suiko/yukiko/ppt/grads-2.0a7.1/bin のディレクトリです GrADS-2.0a7.1 ディレクトリの中の bin ディレクトリに入る bin 中に grb ファイルも入れる bin]$./grads -b ga-> ほにゃらら これで GrADS が起動します -b は, グラフィック表示をしないというオプション 無事, 起動ここにコマンドを打って GrADS を操作 ga->quit quit で終了

22 風データの抜き出し 変換 1 GRIB 形式のデータを抜き出し 変換する手順 GRIB 形式を GrADS で読み込む準備 ctl ファイルを作成 ctl ファイルを書き換え idx ファイルを作成 GRIB 形式を GrADS でバイナリデータ化 GrADS スクリプトを作成 バイナリデータ化 ここまで GrADS バイナリデータをテキストデータに変換 Fortran プログラムの作成 テキストデータ化

23 風データの抜き出し 変換 2 ctl ファイルの作成 grib2ctl.pl という perl スクリプトをネットから探して bin 中に転送 bin 中にある grib2ctl.pl という perl スクリプトを使って,ctl ファイルを作成 perlスクリプトとして実行可能にする [suiko@estuary5 bin]$ chmod a+x grib2ctl.pl [suiko@estuary5 bin]$./grib2ctl.pl ts3hr MANAL_ grb > 0804.ctl./ grib2ctl.pl 実行で,MANAL_ grb の ctl ファイル 0804.ctl が作成される

24 Ctl ファイルの中身 grb ファイルの中の情報が書いてあります これを 2008 年 4 月分全てのファイルをまとめて表示させるように書き換えます

25 ctl ファイルの書き換え Windows 上に 0804.ctl を転送して, テキストエディタで編集します. dset ^MANAL_200804%d2%h2.grb options template tdef 240 linear 00Z01apr2008 3hr 今回は 分をまとめるので,8 個 /1day 30 日 =240 個

26 bin]$./gribmap -i 0804.ctl 長岡技術科学大学 idx ファイルの作成 0804.ctl を Linux 上に戻したら bin 中に元から入っているgribmapを実行 MANAL_ grb.idxファイルが生成こんな風に 4 月分全てを対応付けられたら成功 これで GrADS に読み込む準備が完了

27 GrADS スクリプトの作成 1 12 ヶ月分を同じようにしたら, バイナリデータ化するための GrADS のスクリプトをつくる (0804~0903 まで ) i=1 while (i<=12) 12 ヶ月の grb ファイルを開くスクリプト if (i=1 i=2 i=3 i=4 i=5 i=6) a=3+i 'open 080'a'.ctl' endif if (i=7 i=8 i=9) a=3+i 'open 08'a'.ctl' endif if (i=10 i=11 i=12) a=i-9 'open 090'a'.ctl' endif コマンドは で囲む open xxxx.ctl でその月の grb ファイルが開く テキストエディタで作成

28 GrADS スクリプトの作成 2 'set lon ' 'set lat ' 'set lev 0' 38.9N 格子上の切り出す範囲と階層を設定 'set gxout fwrite' 'set fwrite 'yymm'u.bin' t=1 while(t<=240+time) 'set t 't'' 'd ugrdsfc' t=t+1 endwhile 'disable fwrite' 36.9N 136.9E 1 ヶ月分の u 方向風データをバイナリファイル yymm.bin に書き込む 'close 1' i=i+1 endwhile u 方向風データ ctl ファイルを閉じて次の月へ進む 139.8E t=1

29 GrADS スクリプトの作成 3 i=1 while (i<=12) if (i=1 i=2 i=3 i=4 i=5 i=6) a=3+i 'open 080'a'.ctl' endif if (i=7 i=8 i=9) a=3+i 'open 08'a'.ctl' endif if (i=10 i=11 i=12) a=i-9 'open 090'a'.ctl' endif 'set lon ' 'set lat ' 'set lev 0' これを.gs 形式で保存 if (i=1 i=3 i=6 i=8) time=0 endif if (i=2 i=4 i=5 i=7 i=9 i=10 i=12) time=8 endif if (i=11) time=-16 endif if (i=1 i=2 i=3 i=4 i=5 i=6) yymm=080a endif if (i=7 i=8 i=9) yymm=08a endif if (i=10 i=11 i=12) yymm=090a endif 'set gxout fwrite' 'set fwrite 'yymm'u.bin' t=1 while(t<=240+time) 'set t 't'' 'd ugrdsfc' t=t+1 endwhile 'disable fwrite' 'set gxout fwrite' 'set fwrite 'yymm'v.bin' t=1 while(t<=240+time) 'set t 't'' 'd vgrdsfc' t=t+1 endwhile 'disable fwrite' 'close 1' i=i+1 endwhile 今回は binary.gs とします

30 GrADS スクリプトの実行 0804u から 0903v までの bin ファイルが生成 bin]$./grads -b ほにゃらら ga->run binary2.gs run コマンドでスクリプトを実行 終わったら quit で終了

31 Fortran でテキストデータ化 1 t=1 t=1 t=2 t=3 t=1 u 方向風データ t=4 t=2 t=1 t=1 t=2 t=3 t=1 v 方向風データ bin ファイル t=4 DAT ファイル SWAN 計算用 t=2 TXT ファイル

32 Fortran でテキストデータ化 2 DAT ファイルにするプログラム program main parameter(itm=248,ilon=32,ilat=24) character dataname1*16,dataname2*16,filebase1*5,filebase2*5 character form2*20 integer n dimension u (ilon,ilat),v (ilon,ilat) dimension u2(ilon,ilat),v2(ilon,ilat) common u,v,n,it ***** open file * 7 月データ (2004/07/25 0: /08/31 18:00) open(10,file='0812u.bin',form='unformatted',access='direct' 1,recl=ilon*ilat) open(50,file='0812v.bin',form='unformatted',access='direct' 1,recl=ilon*ilat) ***** データの読み込み n=1 do it=1, itm read(10,rec=it read(50,rec=it ) u ) v ***************** write(*,*)'time',it if( n.le. 9) then form2='(a5,a1,i1 )' else if(n.ge. 10.and.n.le. 99) then form2='(a5,a1,i2 )' else if(n.ge. 100.and.n.le. 999) then form2='(a5,a1,i3 )' else if(n.ge.1000.and.n.le.9999) then form2='(a5,a1,i4 )' endif filebase1='win-u' filebase2='win-v' write(dataname1,form2) filebase1,'_',n write(dataname2,form2) filebase2,'_',n open( 90,file=dataname1,status='unknown') open(100,file=dataname2,status='unknown') do j=ilat,1,-1 write( 90,'(361f10.5)') (-u(i,j),i=1,ilon) write(100,'(361f10.5)') (-v(i,j),i=1,ilon) enddo close( 90) close(100) n=n+1 10 continue close(10) * close(20) * close(30) * close(40) close(50) * close(60) * close(70) * close(80) stop end enddo こんな感じです

33 Fortran でテキストデータ化 3 TXT ファイルにするプログラム wind1.txt(08.4) から wind12.txt(09.3) までの TXT ファイルが完成 program main parameter(nt=248, nx=32, ny=24) character dataname1*16,dataname2*16,filebase1*5,filebase2*5 character form2*20 integer T dimension u(0:nx,0:ny),v(0:nx,0:ny) common u,v,n,it ***** swan で使用するグラフ open( 95,file='wind9.txt',status='unknown') open( 96,file='windpx.dat',status='unknown') open( 97,file='windpy.dat',status='unknown') ***** x 成分データ DO T=1,nt write(*,*)'time',t if( T.le. 9) then form2='(a5,a1,i1 )' else if(t.ge. 10.and.T.le. 99) then form2='(a5,a1,i2 )' else if(t.ge. 100.and.T.le. 999) then form2='(a5,a1,i3 )' else if(t.ge.1000.and.t.le.9999) then form2='(a5,a1,i4 )' endif ***** データの書き込み do j=ny,1, -1 write( 95,'(361f10.5)') (u(i,j),i=1,nx) enddo WRITE(96,'(2F10.5)') REAL(T)/4.0, U(16,15) ENDDO ***** Y 成分データ DO T=1,NT if( T.le. 9) then form2='(a5,a1,i1 )' else if(t.ge. 10.and.T.le. 99) then form2='(a5,a1,i2 )' else if(t.ge. 100.and.T.le. 999) then form2='(a5,a1,i3 )' else if(t.ge.1000.and.t.le.9999) then form2='(a5,a1,i4 )' endif filebase2='win-v' write(dataname2,form2) filebase2,'_',t open(100,file=dataname2,status='unknown') do j=ny,1,-1 read(100,'(361f10.5)') (v(i,j),i=1,nx) enddo close(100) close( 95) ; close(96) ; close(97) stop end wind1.txt(08.4) から wind12.txt(09.3) までの TXT ファイルが完成 filebase1='win-u' write(dataname1,form2) filebase1,'_',t open( 90,file=dataname1,status='unknown') do j=ny,1,-1 read( 90,'(361f10.5)') (u(i,j),i=1,nx) enddo close( 90) ***** データの書き込み do j=ny,1,-1 write( 95,'(361f10.5)') (v(i,j),i=1,nx) enddo WRITE(97,'(2F10.5)') REAL(T)/4.0,V(16,15) ENDDO

34 SWAN をインストールして 計算をします swan4041ab というディレクトリの中に入っています

35 自動高速インストールをします インストール前に入っていてるファイル達 swan4041ab]$ make config swan4041ab]$ make ser SWAN の Implementation Manual に方法が書いてあるので, その通りインストールします. 1. マシン依存のマクロを含むインクルードファイルをまず最初につくります. このファイルは macros.inc と呼ばれ,make config で作られます. 2.SWAN を構築するために make ser と打ちます.( シリアルモード )

36 インストール終了 インストール後はこんな風になってます.

37 SWAN の計算実行に必要なもの 無事インストールできたら, 計算を実行するために必要なデータやファイルを準備します 必要なもの 1SWAN そのもの インストールしました. 2SWAN のコマンドファイル ( 拡張子 swn) これから作成します. 3 風データ 4 地形データ 前回説明した GrADS やらで作ったデータです. 名前は今回 wind1.txt です. 今回は niigata_dep.dat というデータです. 格子状の地形データになっています. この 3 つをさっき SWAN をインストールした SWAN4041AB のディレクトリの中に転送します.

38 SWAN のコマンドファイルの作成 1 インストールした SWAN の中に swan.edt というコマンドファイルの雛形があります. swan.edt の中身はこんな感じ この雛形を元に今回は nihonkai-03.swn というコマンドファイルを作成します.

39 SWAN のコマンドファイルの作成 2 Start-up コマンド PROJECT 今回のプログラムの名前とあだな ( みたいの ) SET なんたら [ 時間, 場所的に一定の水面上昇 ] [x 軸から見た北の角度 ] [ 最小水深 ] [ 最大エラーメッセージ数 ] [ コマンドファイルのエラー設定 ] [ 重力加速度 ] [ 水の密度 ] [ 出力単位の選択 ] [ よくわからない ] [ 波向, 風向の決定 ] [ よくわからない ] [ 出力画面選択 ] [ 最大フルード数 ] [ テスト出力画面選択 ] MODE [ 定常 / 非定常 ] [ 次元選択 ] COORDINATES [ デカルト座標 / 球座標 ]

40 SWAN のコマンドファイルの作成 3 model description コマンド CGRID Computational grid: 計算格子の地理的場所, 格子の大きさ, 解像度, 向きの決定 INPGRID BOTTOM 地形データの地理的場所, 格子数, 格子の大きさ READINP BOTTOM 地形データの単位変換, ファイル名, 読み込み方, INPGRID WIND 風データの格子形状, 地理的場所, 格子数, 格子の大きさ, 定常 / 非定常, 期間 READINP WIND 風データの単位変換, ファイル名, 読み込み方, BOUND SHAPESPEC 境界条件の設定 BOUNDSPEC SIDE 境界条件の細かい設定? INITIAL 初期条件

41 SWAN のコマンドファイルの作成 4 Physics,Numerics コマンド ここら辺で, 具体的にどんな波を起こすのか決めていきます. 主に, いじったのはこのあたり. 波の方程式とかそういうのの設定です. これはまた違うマニュアルで説明されてます.

42 SWAN のコマンドファイルの作成 4 Output,Lock-up コマンド 欲しい出力データの場所や, 期間などを設定します. POINTS ある 1 地点だけのデータを欲しいときに設定します. TABLE POINTS で定めた地点の欲しい出力データの設定をします. GROUP 計算領域の中でまとまった格子データが欲しいときに設定します. BLOCK GROUP で定めた領域の欲しい出力データを設定します. COMPUTE いつからいつまで計算するのか設定します. STOP 入力コマンド終わりの合図.

43 計算を実行する その前に,SWAN にこのコマンドファイルを使ってくださいと言う必要があります. swan4041ab]$ chmod +rx./swanrun swan4041ab]$./swanrun input nihonkai-03 認識してくれたので計算してくれました 計算の実行手順を記したものがあって, それがswanrunというものです. それに,nihonkai-03.swnというコマンドファイルで計算をお願いします. と言います.

44 2 回目以降の計算のために 初めての計算が終わったら,SWAN の初期ファイル swaninit が生成されます. 次からいちいち前の実行手順を踏まなくても,./swan.exe で計算が実行されるように書き換えます. nihonkai-03.swn prt ファイルというのはこんな計算しましたという報告書みたいなものです. nihonkai-03.prt

45 2 回目以降の計算 SWAN のシステムは,swan.exe と swan.edt,swanrun で成り立ってると書いてあります. swan4041ab]$./swan.exe 2 回目以降は,swaninit を書き換えたので, 実行可能ファイル swan.exe を実行させるだけで計算できるようになります. 結果のファイルはこんな風にきちんと 12 個でてきます.

46 結果をみてみる こんな風にして出てきた数値をエクセルとかでみたりします. おしまい