趣味で計算流砂水理 Computational Sediment Hydraulics for Fun Learning

数値計算とか河川工学とかプログラミングのことを書いています

MENU

2018-01-01から1年間の記事一覧

地理院標高タイルから標高値を取得するpythonモジュール : libgsidem2el

R

※2024/8/26 注意:最新の標高データはタイル形式で配信されていないようです。そのため、本モジュールは使えません。 手間ですが、基盤地図情報ダウンロードサービスからダウンロードするしか無いと思われます。 またまた地図ネタです.今回はちょっと力作で…

ミャンマー某河川の河口部:ETOPO1 + geoviews

R

またまた地図ネタです. 先日ネタにでましたミャンマーのある河川の河口部について,ETOPO1とgeoviewsで絵を書いてみました. 大きいものはこちら スケールがわかりにくいですが,河口から300km程度は,10m程度の浅い水深の範囲が続いています. 河口から50k…

コロケート格子の境界条件の設定の悩み

R

ずっと悩んでいるのですが,改めて. 非常にシンプルな一次元不定流計算の定常解です. バックウォーターの計算です. nbviewer 問題は下流端流量です. 境界なので片側差分で計算にすると,少し大きな値になってしまいます. 圧力項dHdxの差分がまずいです…

はてブロにインタラクティブなグラフを埋め込む:外部htmlのとりこみ CDNとかiframeとか

R

今まで,javascriptを直で書いたりとかしてましたが,よく考えたら,htmlをそのままiframeで埋め込むのが楽じゃないかという話です. 結果簡単にできました. 基本的な流れは, Holoviewsでグラフを作成し,htmlにエクスポート githubに公開 github上のhtml…

一次元混合粒径掃流砂河床変動計算の更新

R

github.com またまた更新です. ほぼ完成かなと思います. 計算結果 case1 勾配変曲点設けた計算です. 大きいものはこちら case2 勾配変曲点設けた計算で,上流側河道をガチガチのアーマーコートにしてます. 大きいものはこちら case3 河床勾配は一定で,狭…

numba.jitの衝撃:5分 ⇒ 3秒

R

先日お話させて頂きましたnumba.jitの件です.私の手続きが間違ってました. 前回の河床変動の計算時間ですが,タイトルのとおり,超高速化します.100倍です. github更新しておきました. github.com numba.jitの使い方 ほとんどのnumba.jitに使い方に関す…

rigid-lid近似

R

連投になります. 直交格子でLES + 河床変動計算をやりたいと考えており, iricにあるNabiさんのLESモデルについて調べていると,河床変動は出来ないのですが,関連する論文を調べているとすでに河床変動計算も完成しているようです. この辺りです. https:…

一次元混合粒径掃流砂河床変動計算の更新

R

いろいろ修正しました. だいぶ計算も安定してきました. 今度はgithubに上げました. ちょっとおしゃれな図化にしました. nbviewerのリンクは以下です. http://nbviewer.jupyter.org/github/computational-sediment-hyd/1DBedVariationMixed/blob/master/…

札幌市清田区里塚の航空写真,古地図,地質図を並べてみた

R

2018/9/24追記 清田区里塚で液状化した場所と住所は?被災前後で画像を比較 | わしろぐ にtwitter画像の位置の住所がまとめられていたのでポップアップを作ってみました. Leaflet Sync Demo またまたLeaflet.syncです.ブログへの埋め込み方がわかりました…

Leaflet.Syncで真備町の写真を並べてみました

R

地理院が真備町の被災直後の写真を公開していたいので,Leaflet.Syncで被災前,被災後,地形分類図を並べてみました. いつもどおり,ブログに埋め込みたかったのですが,駄目だったので,githubに上げました. Leaflet Sync Demo - with three maps listeni…

はてブロでbokeh

R

bokehWelcome to Bokeh — Bokeh 0.13.0 documentationというD3.jsのpythonラッパーで遊んでみました. D3.jsよりは楽で,pythonで書いたものをhtmlで吐き出す形式です. Web上で解析結果をインタラクティブに表示するのに使えそうです. これとは別にplotlyP…

一次元河床変動計算,掃流砂のみ混合粒径:コードは世界最短?

R

やっとできました. 計算条件 境界条件:周期境界条件 平均河床勾配:1/500 水深1.4m 粒径:7粒径区分,うち2粒径が動かない粒径(コード参照) dx:200m,dt:1sec 初期条件:一部区間に最大粒径の2倍程度のマウンドを設ける 計算結果と考察 計算があってる…

粒度分布

K

吸収境界 地震動の吸収境界は、次のものが多く使われていました。 http://www.gms.bosai.go.jp/GMS/doc/B5766511.pdf 吸収境界の起債はp655以降にあります。 基本的には平面方向の地震動については吸収してもらわないといけないので、 便宜上、減衰させる方…

粒度分布のタイプ

R

以前読んだ京大の論文を思い出しました. http://www.dpri.kyoto-u.ac.jp/nenpo/no50/ronbunB/a50b0p64.pdf 分布系を決定するのはなかなか難しいかもしれないですね. 結構いろんな分布があるようです. 過去のデータより分布型を推定することが必要ですが,…

docker for windowsで、ubuntuのデスクトップ(GUI)を作りたい

R

docker for windowsで、ubuntuのデスクトップ(GUI)を作りたいけど、失敗したという話です。 延べ2日間くらいやりましたが、駄目でした。 docker for windowsのインストール 特に、何の問題もないですが、プロキシ環境なので、その設定だけです。 dockerと…

新規PCの導入:備忘録

R

ようやく時間ができて、新規PCのセットアップが終わりました。 最近色々な事情でPCが壊れるので慣れてきました。 備忘録的にソフトウェア関係を列挙しておきます。 まずは firefoxとDropboxを入れて、firefox syncでブックマーク、パスワード、アドオン、dro…

無反射境界条件

R

久しぶりに計算にドはまりしました. 例の変分法ですが,厳しい条件のときにどうも感度係数が評価できていないような感じがあり,詳細にチェックしてみた所,ミスをみつけてしまいました. 詳細は省きますが,空間方向の境界で反射していました.0勾配にして…

windows intel fortran 環境でのcmakeによるスタックサイズの設定

R

プライベートも仕事も忙しくてなかなか書けてません. 忘備録的に. コマンドラインで打つときは, ifort *** /fast /F10000000 的な感じですが,cmakeで set (CMAKE_Fortran_FLAGS_RELEASE "/fast /F10000000") だと駄目です. 正しくは set (CMAKE_Fortran…

はてブロでjavascript

R

jsを埋め込み 普通に書けるみたいです.びっくりしました. 遊びで,leafletに地理院地図を取り込みD3.jsでsvgを書いて重ねてみました. 以下が参考サイトです. http://bl.ocks.org/d3noob/9267535 https://shimz.me/blog/d3-js/3517 htmlからとのタグだけ…

Serre-type equations in deep waterって何?

R

https://arxiv.org/pdf/1607.00216.pdf 気になって,読んでるけど全然わかんないです. わかれば教えてください.