またまた地図ネタです．今回はちょっと力作です．コードは短いですが，結構便利です． 国土地理院の標高タイルから標高値を取得するpythonモジュール libgsidem2el を作ってみました． 標高タイルとは， https://maps.gsi.go.jp/development/ichiran.html#de…

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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