BP#7 ------------------------------------------------------------------------- NTHMP (National Tsunami Hazard Mitigation Program) ------------------------------------------------------------------------- 2010 NTHMP WORKSHOP Venue: TBA (Catalina Island, CA; Santa Monica, CA; Irvine, CA; Corvallis, OR; Boulder, CO. NOV. 16th-18th 2010 ---- o ---- MATLAB SCRIPT: Tsunami Runup onto a Complex Three-dimensional Beach; Monai Valley TYPE: Laboratory Benchmarking problem (Lab) NUMBER: 07 TSUNAMI SOURCE: Input wave (IW) SCRIPT NAME : BP7_Lab_IW04.m AUTHOR(S) : D. Nicolsky, C. Synolakis, S Grilli and J. Horrillo. Copyright 2010 NTHMP/NOAA Revision: 0 Date: Sep/2010 REFERENCE: Synolakis, C.E., E.N. Bernard, V.V. Titov, U. Kâno?lu, and F.I. González (2007): Standards, criteria, and procedures for NOAA evaluation of tsunami numerical models. NOAA Tech. Memo. OAR PMEL-135, NOAA/Pacific Marine Environmental Laboratory, Seattle, WA, 55 pp. NOAA/PMEL Center for Tsunami research http://nctr.pmel.noaa.gov/benchmark/index.html Liu, P.L.-F., H. Yeh, and C. Synolakis (2008): Advanced Numerical Models for Simulating Tsunami Waves and Runup. Advances in Coastal and Ocean Engineering, 10, 250 pp. INCLUDED FILE WITH THIS SCRIPT: IW_04_gages.mat (lab. data file time series of surface elevation at three different gages) Note: File mentioned above is needed to run successfully this script. CONTACT INFO: Dr. Juan J. Horrillo, Texas A&M University at Galveston Maritime Systems Engineering (MASE) 200 Seawolf Park or P.O. Box 1675 Galveston, TX 77553-1675 Phone: (409) 740-4465 horrillj@tamug.edu Website: http://www.tamug.edu/ ------------------------------------------------------------------------- > BP7_Lab_IW04(filename) BP7_Lab_IW04 reads input file 'filename' containing user's model results for gage 1 2 and 3 and compares them with laboratory data described in the Benchmark "Tsunami Runup onto a Complex Three-dimensional Beach; Monai Valley", NOAA/PMEL, (http://nctr.pmel.noaa.gov/benchmark/index.html). Input file 'filename' must be named as follows: XXXXXX_IW_04.dat where XXXXXX is the name of the numerical model that is used, eg. TSUNAMI_IW_04.dat. The input file, 'filename' must have the following configuration or structure: 1- File 'filename' is a matrix ASCII files of size 3992x4(row x column) 2- Column 1 is time given in second, ranging form 0 to 199.55 sec with increment of 0.05 sec. 3- Columns 2, 3 and 4 are numerical gage results in centimeter +-----------+----------+----------+--------------+ | Time(sec) |Gage1(cm) |Gage2(cm) | Gage3(cm) | +-----------+----------+----------+--------------+ 0 -0.0003456 0 0 0.0500 -0.0003346 0 0 0.1000 -0.0002543 0.0001467 0 . . . . . . . . . . . . 199.4000 0.0001528 0.0081391 -0.0199115 199.4500 -0.0003549 0.0080003 -0.0200863 199.5000 -0.0005610 0.0077974 -0.0201327 199.5500 -0.0007056 0.0075332 -0.0199876 +-----------+----------+----------+--------------+ NOTE: This script is self-containing. It loads the lab. data automatically to facilitate user's comparison, peer review and determine model performance. The experimental data are saved on files IW_04_gages.mat. Be sure that the lab data files IW_04_gages.mat resides in the same directory or subdirectory as the Matlab script (BP7_Lab_IW04.m). The user's input file 'filename' must be present in the same directory or subdirectory of the script file BP7_Lab_IW04 as well. This script uses the normalized root mean square deviation (NRMSD) or root mean square error to measure the numerical model precision. NRMSD measures the differences between values predicted by the numerical model and the values actually observed in the experiment or in the field. The NRMSD is plotted in time to visualize model performance in a particular moment. As it is known tsunami models usually predict fairly well the leading waves but perform poorly in predicting subsequent waves. In this way users can quantify model performance for the first two or three waves and so forth. Examples: > BP7_Lab_IW04 In this example the Matlab scrip just plots laboratory results. > BP7_Lab_IW04('TSUNAMI_IW_04.dat') In this example the Matlab scrip plots laboratory results and and user's model results.