Outline of this page
This page introduces the programs to estimate the probability precipitation of 1 day.
- In the programs, following probability distribution models are used:
| No | Short form | probability distribution function | parameters |
|---|
| 1 | Gumbel | Gumbel distribution | 2 |
| 2 | GEV | Generalized Extreme Value distribution | 3 |
| 3 | Weibull (Goda) | Weibull distribution (L-moments by GODA) | 3 |
| 4 | Weibull (LSM) | Weibull distribution (Least Square Method) | 3 |
| 5 | Weibull (MLM) | Weibull distribution (Maximum Likelihood Method) | 3 |
| 6 | SQRT-ET | SQRT exponential-type distribution of maximum | 2 |
| 7 | LP3 | Log Pearson type III distribution | 3 |
| 8 | LN3 (Iwai) | Logarithmic Normal distribution (Iwai method: quantile) | 3 |
| 9 | LN3 (Moment) | Logarithmic Normal distribution (Moment method) | 3 |
| 10 | LN3 (Trial) | Logarithmic Normal distribution (Trial method) | 3 |
| 11 | Exponential | Exponential distribution | 2 |
| 12 | GPD | Generalized Pareto distribution | 3 |
- The Cunnane formula is used as a plotting position formula for an observed value plot.
- Following papers were referred to make programs.
- http://thesis.ceri.go.jp/center/info/geppou/ceri/0005005050.html
(星清：水文統計解析，開発土木研究所月報，No.540，1998年5月 in Japanese character)
- http://iss.ndl.go.jp/books/R000000004-I4497918-00
(星清：現場のための水文統計(1)，開発土木研究所月報，No.540，1998年5月 in Japanese character)
- http://iss.ndl.go.jp/books/R000000004-I4522177-00
(星清・新目竜一・宮原雅幸：現場のための水文統計(2)，開発土木研究所月報，No.541，1998年6月 in Japanese character)
- Source data for input samples are following:
- Record of Maximum daily rainfall in Sapporo : http://iss.ndl.go.jp/books/R000000004-I4522177-00
- Record of Maximum daily rainfall in Maebashi : Webpage of Japan Meteorological Agency
- Record of Maximum daily rainfall in Kyoto : ISBN978-4-627-49631-6
- Note that although Exponential distribution and GPD shall be aplied to the analysis using Peaks Over Threshold data,
the same data as other methods are used for these calculation examples.
Documents
Samples of Input data
Program codes
Source code by f90
| Filename | Description |
|---|
| f90_QQ.txt | program for parameter estimation of the probability distribution function |
| f90_BTS.txt | program for confidence interval estimation by the boot strap method |
Batch commands for execution
Batch commands for f90_QQ
gfortran -o f90_QQ.exe f90_QQ.f90
f90_QQ 1 inp_RF_M.txt out_QQ_M_01.csv
f90_QQ 2 inp_RF_M.txt out_QQ_M_02.csv
f90_QQ 3 inp_RF_M.txt out_QQ_M_03.csv
f90_QQ 4 inp_RF_M.txt out_QQ_M_04.csv
f90_QQ 5 inp_RF_M.txt out_QQ_M_05.csv
f90_QQ 6 inp_RF_M.txt out_QQ_M_06.csv
f90_QQ 7 inp_RF_M.txt out_QQ_M_07.csv
f90_QQ 8 inp_RF_M.txt out_QQ_M_08.csv
f90_QQ 9 inp_RF_M.txt out_QQ_M_09.csv
f90_QQ 10 inp_RF_M.txt out_QQ_M_10.csv
f90_QQ 11 inp_RF_M.txt out_QQ_M_11.csv
f90_QQ 12 inp_RF_M.txt out_QQ_M_12.csv
Format
| Arg. | Description |
|---|
| f90_QQ | program name executed |
| No | number for specifing the distribution model |
| fnameR | Input data file |
| fnameW | Output data file |
Regarding 'No' in above, please refer to the table on outline.
Batch commands for f90_BTS
gfortran -o f90_BTS.exe f90_BTS.f90
f90_BTS 1 inp_RF_M.txt out_BTS_M_01.csv inp_BHI_M_01.txt
f90_BTS 2 inp_RF_M.txt out_BTS_M_02.csv inp_BHI_M_02.txt
f90_BTS 3 inp_RF_M.txt out_BTS_M_03.csv inp_BHI_M_03.txt
f90_BTS 4 inp_RF_M.txt out_BTS_M_04.csv inp_BHI_M_04.txt
f90_BTS 5 inp_RF_M.txt out_BTS_M_05.csv inp_BHI_M_05.txt
f90_BTS 6 inp_RF_M.txt out_BTS_M_06.csv inp_BHI_M_06.txt
f90_BTS 7 inp_RF_M.txt out_BTS_M_07.csv inp_BHI_M_07.txt
f90_BTS 8 inp_RF_M.txt out_BTS_M_08.csv inp_BHI_M_08.txt
f90_BTS 9 inp_RF_M.txt out_BTS_M_09.csv inp_BHI_M_09.txt
f90_BTS 10 inp_RF_M.txt out_BTS_M_10.csv inp_BHI_M_10.txt
f90_BTS 11 inp_RF_M.txt out_BTS_M_11.csv inp_BHI_M_11.txt
f90_BTS 12 inp_RF_M.txt out_BTS_M_12.csv inp_BHI_M_12.txt
Format
f90_BTS No fnameR fnameW1 fnameW2
| Arg. | Description |
|---|
| f90_BTS | program name executed |
| No | number for specifing the distribution model |
| fnameR | Input data file |
| fnameW1 | Output data file |
| fnameW2 | Output data file for drawing histogram |
Regarding 'No' in above, please refer to the table on outline.
Batch file for GMT execution (for Maebashi)
- 3 type's batch files are prepared for GMT execution.
- One is for the control commands, another one is for making plot data, and the other one is for plot commands.
- Batch files for making plot data and for plot commands are called from the batch file for control commands.
- In the case that many similar graphs are required, this method is very convenient.
Q-Q plot
Bootstrap confidence interval
Histogram of bootstrap point estimate
