VBAソースファイルtys-sfm_v10tys-sfm_v10.lzhをダウンロードして解凍すると、「貯留関数マクロ.xls」ができます。この「貯留関数マクロ.xls」の使用方法を説明します。ダウンロードについては、水理水文ソフトウェアのHPを参照してください。本マクロは、各条件シートに条件を入力するだけで、1流域についての貯留関数による流域流出解析計算が行えます。実績流量より、貯留関数係数の回帰を行う。任意の貯留関数係数を用いて流出解析を行う。2種類の計算を行うことができます。シートで貯留関数係数の回帰を行うかどうかをチェックしてください。計算条件シート再度計算実行した場合は、計算結果シート(複数)の値がすべて上書きされます。・遅滞時間計算ピッチを入力してください。入力した実績流量のうち、係数決定に用いるデータの範囲を指定します。開始点を「1」、終了点を「2」と入力してください。)を入力してください。「雨量・流量データ」シートの流域平均雨量を入力してください。「計算結果」シートに遅滞時間の設定値ごとの貯留量と流出量の関係の計算結果を書き込みます。この中で、貯留量と流出量のループの積算値の最小値を遅滞時間として採用します。遅滞時間の繰り返し計算は、0時間~最大値(流路長を木村の式に適用した場合の遅滞時間+1)を遅滞時間の計算ピッチ分です。ここで求めた遅滞時間より実績流量を用いた最小自乗法による回帰を行い、K・Pを算出します。「雨量・流量データ」シートの計算流量の列に上記で求めた貯留関数係数を用いて計算した流量を書き込みます。「貯留量と流出量の関係」シートに決定した貯留関数係数を用いた貯留量と流出量の関係を出力します。その結果は「貯留量と流出量の関係図」シートに反映されます。「雨量・流量データ」シートの計算流量の列に入力した貯留関数係数を用いて計算した流量を書き込みます。シートに反映されます。「貯留量と流出量の関係」シートに貯留量と流出量の関係を出力します。マクロは、通常の貯留関数法についての計算用であり、総合貯留関数法には適用していません。通常の貯留関数法と総合貯留関数法では、K・P値はそれぞれ全く異なりますので注意してください。ここで、通常の貯留関数法とは、貯留関数を貯留量(m3)と流出量(m3/s)の関係で表現したもの、総合貯留関数法とは貯留関数を貯留高(mm)と流出高(mm/hr)の関係で表現したものです。総合貯留関数法は貯留関数法を単位流域面積で表現したものと言い換えることができます。本や論文によっては、総合貯留関数法を単に「貯留関数法」と表記されている場合がありますので、基礎式による確認が必要です。VBAソースファイル本マクロのソースを公開しますので、参考にしてください。Subtys-sfm_v_2_0_1tys-sfm_v_2_0_1.lzhをダウンロードして解凍すると、「tys-sfm_v_2_0_1.xls」ができます。この「tys-sfm_v_2_0_1.xls」の使用方法を説明します。ダウンロードについては、水理水文ソフトウェアのHPを参照してください。本マクロは、各条件シートに条件を入力するだけで、1流域についての総合貯留関数法と貯留関数法の2種類の流域流出解析計算が行えます。実績流量より、貯留関数係数の回帰を行う。任意の貯留関数係数を用いて流出解析を行う。2種類の計算を行うことができます。シートで貯留関数係数の回帰を行うかどうかをチェックしてください。・遅滞時間計算ピッチを入力してください。入力した実績流量のうち、係数決定に用いるデータの範囲を指定します。開始点を「1」、終了点を「2」と入力してください。)を入力してください。「雨量・流量データ」シートの流域平均雨量を入力してください。「計算結果」シートに遅滞時間の設定値ごとの貯留量と流出量の関係の計算結果を書き込みます。この中で、貯留量と流出量のループの積算値の最小値を遅滞時間として採用します。遅滞時間の繰り返し計算は、0時間~最大値(流路長を木村の式に適用した場合の遅滞時間+1)を遅滞時間の計算ピッチ分です。ここで求めた遅滞時間より実績流量を用いた最小自乗法による回帰を行い、K・Pを算出します。「雨量・流量データ」シートの計算流量の列に上記で求めた貯留関数係数を用いて計算した流量を書き込みます。「貯留量と流出量の関係」シートに決定した貯留関数係数を用いた貯留量と流出量の関係を出力します。その結果は「貯留量と流出量の関係図」シートに反映されます。「雨量・流量データ」シートの計算流量の列に入力した貯留関数係数を用いて計算した流量を書き込みます。シートに反映されます。「貯留量と流出量の関係」シートに貯留量と流出量の関係を出力します。プログラムは、総合貯留関数法と通常の貯留関数法の両方が行えますが、総合貯留関数法と通常の貯留関数法では、K・P値はそれぞれ全く異なりますので注意してください。ここで、総合貯留関数法とは貯留関数を貯留高(mm)と流出高(mm/hr)の関係で表現したもの、通常の貯留関数法とは、貯留関数を貯留量(m3)と流出量(m3/s)の関係で表現したものです。総合貯留関数法は貯留関数法を単位流域面積で表現したものと言い換えることができます。本や論文によっては、総合貯留関数法を単に「貯留関数法」と表記されている場合があり、特殊な手法ではなく、かなり一般的に利用されていますので、基礎式による確認が必要です。VBAソースファイル本マクロのソースを公開しますので、参考にしてください。Subtys-sfb_v1_0_0tys-sfb_v1_0_0.lzhをダウンロードして解凍すると、「tys-sfb_v1_0_0.xls」ができます。この「tys-sfb_v1_0_0.xls」の使用方法を説明します。ダウンロードについては、水理水文ソフトウェアのHPを参照してください。本マクロは、各条件シートに条件を入力するだけで、複数流域の総合貯留関数による流出解析計算が行えます。複数流域を一括して流出解析を行う。シートで貯留関数係数の回帰を行うかどうかをチェックしてください。プログラムは、総合貯留関数法と通常の貯留関数法の両方が行えますが、総合貯留関数法と通常の貯留関数法では、K・P値はそれぞれ全く異なりますので注意してください。ここで、総合貯留関数法とは貯留関数を貯留高(mm)と流出高(mm/hr)の関係で表現したもの、通常の貯留関数法とは、貯留関数を貯留量(m3)と流出量(m3/s)の関係で表現したものです。総合貯留関数法は貯留関数法を単位流域面積で表現したものと言い換えることができます。本や論文によっては、総合貯留関数法を単に「貯留関数法」と表記されている場合があり、特殊な手法ではなく、かなり一般的に利用されていますので、基礎式による確認が必要です。VBAソースファイル本マクロのソースを公開しますので、参考にしてください。