#author("2022-07-05T19:32:24+09:00","default:Miyashita","Miyashita") #author("2022-07-20T12:30:38+09:00","default:Miyashita","Miyashita") *地形データ作成 メモ [#q1f06bfa] 内閣府中央防災会議の地形データから JAGURS 計算用に cf フォーマットの grd ファイルを作成してみた．~ 重要な点は以下の通り．~ ~ ** Cell-based から Vertex-based に変換 [#y5399930] 中央防災会議のデータは格子の真ん中（セル中心）に 水深 \(h\) を置いているが， JAGURS で計算するためには 格子線上？に置く必要がある．~ そうするともともとの地形データが 720 x 540 だとしたら， 縦横半メッシュ分広げて 721 x 541 のデータに変換する必要がある．~ ~ MATLAB で半メッシュ分外挿した．外挿メソッドは 'nearest' でもいいけど半メッシュ分ならまあ変なことにはならないだろうと思って 'linear' にしておいた． MATLAB で半メッシュ分外挿した．外挿メソッドは半メッシュ分なら変なことにはならないだろうと思って 'linear' にしておいた． 'nearest' も可． #codeprettify(lang-matlab){{ F = griddedInterpolant(X0',Y0',obj.topo_org','linear','linear'); topo = F(X1',Y1'); topo = topo'; }} ~ ~ **ネスティング領域間の倍率を考慮した領域サイズ [#rfa5c1eb] ↑の縦方向の入力する地形データの縦横のデータの数をそれぞれ ny, nx とし，ネスティングの解像度の倍率が 3 だとすると， (ny-1) と (nx-1) が 3 で割り切れないと計算開始前にエラーを返す．~ ~ 中央防災会議データは大体サイズが 540 x 720 とかになっているので，これをこのまま 541x 721 に変換すればOK．たまに 3 で割り切れない要素数になっているのもあるので，適当にうまく切り取って使う．~ ~ ~