damyarou

python, GMT などのプログラム

Jupyter の ipynb ファイルの中身を見たい

2020.09.20投稿 必要性 Jupyter notebook 上でコードを書いていて、「あのとき作った Jupyter のコードをコピーしたい!」ということがよくある。その ipynb ファイルが、たまたま今開いている Jupyter のディレクトリの下にあればいいのだが、そうでない場…

LG 27インチ4Kディスプレイ購入

9月7日注文、9月8日到着で、LGの27inch 4K display を購入しました。¥57,980也。ちょっと奮発しすぎたか。 4月より在宅勤務しています。会社に行くのは月に2〜3回ですかね。 しかし7月よりお客さんから設計の仕事をいただき、今はかなり忙しい。ま…

Python 逆調整池運用解析修正版(Flood Routine)

逆調整池運用解析プログラムの修正版 関数FLOODの中の収束計算で、前回貯水池水位のみを既知としていたが、本来貯水池水位と放流量はペアなので、前回放流量も既知とした。プログラム作成時になぜ収束計算の中で前回流量を既知として確定していなかったのか…

Relationship between critical depth and discharge coefficient for free overflow on dam spillway

In general, the discharge for free overflow on dam spillway can be expressed as follow: On the other hand, a critical depth of rectangular cross section channel can be expressed as follow: When it is assumed the water depth at the overflow…

Python:計算結果をエクセル出力し加工する(月平均値の算出とエクセル化)

能書き 報告書を作成しているとき、Pythonで計算した結果をWordに表形式で貼り付けたい場合がある。 これが、TeXやhtmlであれば、プロブラム内でタグを埋め込んだ形で標準出力し、コピー&ペーストでTeXもしくはhtml文書に貼り付けるのだが、相手がWordとな…

GraphvizでMLPと通常の重回帰分析のイメージ図を作る

作例 GraphvizでMLPと通常重回帰分析のイメージ図(conceptual diagram あるいは schematic diagram かな)を作ってみた。 作例は以下の通り。 プログラミング環境は以下の通り。 macOS Catalina iMac (Retina 4K, 21.5-inch, 2017) Graphvizのインストール …

Structural design formulas for penstocks embedded in rock

The structural design formulas for penstocks embedded in rock which are shown in 'Technical standards for gates and penstock' in Japan are described. Allowable stresss of steel material Welded joint efficiency Location of welding RT or UT …

Python 岩盤内埋設式水圧鉄管設計プログラム(新版)

はじめに 岩盤内埋設式水圧鉄管プログラムは、過去に一度掲載している(以下参照)。今回のプログラミング内容は工学的には過去のものとダブルが、過去投稿も、エクセルの罫線の書き方やセルの色付などの参考になるため、削除せずにおく。 https://damyarou.…

Python 軸対称応力解析プログラム

はじめに 軸対称モデルとして扱う構造物として、水平トンネルなど水平方向に回転軸を持つ構造物と、調圧水槽など鉛直方向に回転軸を持つ構造物がある。 これまで、四角形要素の2次元応力解析プログラムを軸対称解析用に書き換えたプログラムを使ってきたが…

トランジション(矩形断面から円形断面)での損失水頭

圧力水路の設計をしていると、取水口周辺で矩形断面から円形断面に変化する部分が必要となるが、その時の形状変化による損失(摩擦損失を除く)の計算について書いておく。 これまでは、矩形断面の断面積に等しい等価円を考え、漸縮管の損失水頭として扱って…

パソコン支給

2020年5月21日、会社からノートパソコンというかSurfaceが支給されました、 マシン仕様は以下の通り。 プロセッサ:Intel(R) Core(TM) i5-1035G4 CPU @ 1.10GHz 1.50 実装RAM:8GB OS:Windows 10 Pro ストレージ:SSD128GB 支給日の木曜日と金曜日は忙しか…

Python 効率カーブ(表付きグラフ)

簡単な表と一体化した画像を作った。 作例 このグラフにおけるTips 作図領域を上下に分ける drawfig1() は下側のグラフ(作図領域の下側70%)、```drawfig2()''' は上川の表を描画する。 plt.axes((0.0, 0.0, 1.0, 0.7)); drawfig1() plt.axes((0.0, 0.7, 1.…

Python 月次・年次データ集計(雨量)

日雨量を月・年毎に合計したもの。 各月の雨量の単位は mm/momth (月累計雨量)、 各年の雨量の単位は mm/year(年累計雨量)である。 #============================== # Rainfall # Dayly => Monthly, Yearly data #============================== impor…

Python 月次及び年次データ集計(流量)

概要 日平均流量の月次および年次データ集計のプログラム。 ここでは流量の集計を行い、月平均・年平均を算出している。 年平均は、「月平均の平均」ではなく、「1年間全データの平均」としている。 これは365日の流況図から出した年平均と数値を合わせるた…

Python バージョンを入れ替える(Mac)

2020年4月27日、Pythonは3.8.1を入れていたが、3.8.2にしてみた。 Pythonと関連ライブラリのインストールはpyenvとpipのみで行っている。 目的はPythonのバージョンを変えるだけなので、.zshrcに記載しているパスや、Jupyterの設定は変更しない。もちろんこ…

流況曲線作図(改訂版)

流況曲線作図プログラムを(自分にとって)少しわかりやすくしたのでアップ。 対象期間内に欠測はないことを前提にしている。 import numpy as np import pandas as pd import datetime import matplotlib.pyplot as plt from scipy import interpolate def …

VS code と BasicTeX 2020 をインストール

2020年4月12日、最近話題のVS codeをインストール。また、BasicTeX 2020もインストール。 BasicTeXは、2019がインストールしてあったので、まずこれを削除。 また、which ghostscript しても、なぜだか「見つからない」と言われるので、再インストール。 Bas…

matplotlibで使える色一覧

色一覧 画像作成プログラム import matplotlib import colorsys import matplotlib.pyplot as plt from matplotlib import rcParams rcParams['font.family'] = 'Ricty Diminished' def penc(hval): r=int(hval[1:3],16) g=int(hval[3:5],16) b=int(hval[6:8…

流況曲線と運転パターン図

はじめに 同じような図は何回か作成しているのだが、また必要になったため、新しく作り直した。 最新版のプログラムをアップしておく。 成果図 プログラム import numpy as np import pandas as pd import datetime import matplotlib.pyplot as plt def rda…

サージング解析プログラムの挙動解明(2)

解析条件 解析パラメータを下表に示す。 変化させているパラメータは、制水口径および遮断時間でる。 項目 採用値 解析条件 負荷遮断 シャフト内径(断面積) D=10.0 m (F=78.5 m2) ポート内径(断面積) Dp=2.0, 3.0, 4.0, 8.0m (Fp=3.14, 7.065, 12.56, 50…

サージング解析プログラムの挙動解明

きっかけ 昨日、私のwebページをご覧になった方から、サージング解析プログラムに関するメールを頂いた。 下が、その疑惑のページのコピーである。 図を見てわかるように。 サージング解析では、「制水口(ポート)径が小さい場合、遮断時間を大きくしていく…

Python 開水路トンネルの不等流解析プログラム

開水路トンネルの不等流解析(常流)を行ったので、そのプログラムをアップしておく。 全区間、常流であるため、下流端で固定水位を与え、上流に向かって逐次水位を求めていく。 成果図 このプログラムでのTips 与えられた固定点数点に対し計算点を等間隔に…

macOS Catalinaのクリーンインストールと作業環境更新

2020年1月18日(土)、iMacとMacbook proにmacOS Catalinaをクリーンインストールし、作業環境を更新した。 更新内容を記録しておく。 マシン iMac (Retina 4K, 21.5-inch, 2017) MacBook Pro (Retina, 13-inch, Mid 2014) 不便な点 この更新を行ってこれま…

Python 画像処理関係

画像変換 from PIL import Image import os files = os.listdir() for file in files: base,ext=os.path.splitext(file) if ext=='.png': input_im = Image.open(base + ".png") rgb_im = input_im.convert('RGB') rgb_im.save(base + ".jpg",quality=30) pr…

Python RC部材の設計

RC部材の設計をしたので流れを記録しておく。 リンク 平面骨組解析プログラムはこちら(Qiita)。 配筋計算プログラムの自作関数はこちら(はてな)。 出力 モデル図 断面力図 配鉄設計 プログラム モデル図 import matplotlib.pyplot as plt import numpy a…

Python コンター図作成

コンター図を作成する必要があったため、そのプログラムを作成した。 作例 プログラムソース import numpy as np import matplotlib.pyplot as plt import seaborn as sns def drawfig(xx,yy,zh,zc): q=68.64 vp1=7.0; dp1=np.sqrt(4*q/np.pi/vp1) vp2=2.5; …

Python 岩盤内埋設式水圧鉄管設計プログラム(旧版)

概要 岩盤内埋設式水圧鉄管の設計プログラムをフルで書いてみたので残しておく。 プログラム前半の非常に長い関数「def xlswrite」は、エクセルに書式指定して結果を書き出すもの。 最終的にワードの報告書に貼り付けなくてはならないので、エクセルにしてお…

Python jpg画像を縮小保存しhtmlで表示

タイトル通り、jpg画像を縮小保存しhtmlで表示する。 import glob import os from PIL import Image import os def select_pic(): path = '*.jpg' # 現フォルダ内にあって拡張子がJPGのファイル名を取得 file_list = glob.glob(path) print(file_list) retur…

雑記 また機械学習の本を買ってしまった

記事の最後に行く 昨日(2019.08.11)、また機械学習の本を買ってしまった。 オライリー・ジャパン PythonによりAIプログラミング入門 ディープラーニングを始める前に身に着けておくべき15の基礎技術 あと、久しぶりにpipで一括アップデートしたら、matplot…

雑記 機械学習の勉強を始めてみる

記事の最後に行く いよいよ7月。 昨年(2018年10月10日)、マレーシアでの業務が完了し、帰国した。 それ以降、以下のような感じで過ごしてきた。 2018年10月 ひま 2018年11月 ひま 2018年12月 忙しい 2019年01月 忙しい 2019年02月 忙しい(ベトナム出張2…