This study was conducted to derive the regional design rainfall by the regional frequency analysis based on the regionalization of the precipitation suggested by the first report of this project. According to the regions and consecutive durations, optimal design rainfalls were derived by the regional frequency analysis for L-moment in the second report of this project. Using the LH-moment ratios and Kolmogorov-Smirnov test, the optimal regional probability distribution was identified to be the Generalized extreme value (GEV) distribution among applied distributions. regional and at-site parameters of the GEV distribution were estimated by the linear combination of the higher probability weighted moments, LH-moment. Design rainfall using LH-moments following the consecutive duration were derived by the regional and at-site analysis using the observed and simulated data resulted from Monte Carlo techniques. Relative root-mean-square error (RRMSE), relative bias (RBIAS) and relative reduction (RR) in RRMSE for the design rainfall were computed and compared in the regional and at-site frequency analysis. Consequently, it was shown that the regional analysis can substantially more reduce the RRMSE, RBIAS and RR in RRMSE than at-site analysis in the prediction of design rainfall. Relative efficiency (RE) for an optimal order of L-moments was also computed by the methods of L, L1, L2, L3 and L4-moments for GEV distribution. It was found that the method of L-moments is more effective than the others for getting optimal design rainfall according to the regions and consecutive durations in the regional frequency analysis. Diagrams for the design rainfall derived by the regional frequency analysis using L-moments were drawn according to the regions and consecutive durations by GIS techniques.