BP#9 GAUGE ------------------------------------------------------------------------- NTHMP (National Tsunami Hazard Mitigation Program) ------------------------------------------------------------------------- 2010 NTHMP VALIDATION WORKSHOP Venue: Galveston, TAMU TX Margh 30- Apr. 1 ---- o ---- MATLAB SCRIPT: Okushiri Island TYPE: Field Benchmarking Problem (Field) NUMBER: 09 TSUNAMI SOURCE: Coseismic Deformation (CS) SCRIPT NAME : BP9_Field_GAUGE.m AUTHOR(S) :J. Horrillo, Gyeong-Bo Kim. Copyright 2011 NTHMP/NOAA Revision: 0 Date: Sep/2011 REFERENCE: Synolakis, C.E. (1986): The Runup of Long Waves. Ph.D. Thesis, California Institute of Technology, Pasadena, California, 91125, 228 pp. NOAA/PMEL Center for Tsunami research http://nctr.pmel.noaa.gov/benchmark/index.html Synolakis, C.E. (1987): The runup of solitary waves. J. Fluid Mech., 185, 523-545. INCLUDED FILES WITH THIS SCRIPT: Essashi_Iwanai.mat (Field data file for time sierie at gages 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/ ------------------------------------------------------------------------- > BP9_Field_GAUGE(filename1,filename2) BP9_Field_GAUGE reads input files 'filename1' and 'filename2' containing model results and compares them with field data recorded at gages Iwanai and Esashi. filename1 and filename2 are model generated files corresponding to to Iwanai and Esashi respectively. Input files 'filename1' and 'filename2' need to be named as follows: XXXXXX_BP9_IWANAI.dat and XXXXXX_BP9_ESASHI.dat, where XXXXXX is the name of the numerical model that is tested, eg. NEOWAVE_BP9_IWANAI.dat and NEOWAVE_BP9_ESASHI.dat. The input files, 'filename1' and 'filename2' must have the following configuration: 1- Files 'filename1' is a matrix ASCII files of size 121x2(row x column) as shown below (no heading). The number of rows could be different of 689, but it is desirable to sample at least every 0.50. t ranges= [0 60]min with icrements of 0.5. filename1 +---------+-------------+ | t(min) eta(m) +---------+-------------+ 0.00000 0.00000 0.50000 0.00414 1.00000 0.00829 1.50000 0.01243 2.00000 0.01657 . . . . . . 57.00000 0.35074 57.50000 0.35626 58.00000 0.36179 58.50000 0.36731 59.00000 0.37257 59.50000 0.37533 60.00000 0.37809 +--------+--------------+ 2- Files 'filename2' is a matrix ASCII files of size 601x2(row x column) as it is shown below (no heading). The number of rows could be different of 401, but it is desirable to sample at least every 0.10. t ranges= [0 60]min with icrements of 0.1. filename2 +---------+-------------+ | t(min) eta(m) +---------+-------------+ 0.00000 -0.01690 0.10000 -0.01690 0.20000 -0.01690 0.30000 -0.01690 0.40000 -0.01690 0.50000 -0.01690 . . . . . . 59.60000 -0.13305 59.70000 -0.15487 59.80000 -0.17551 59.90000 -0.19370 60.00000 -0.21190 +-------+---------+---------+ NOTE: This script is self-containing. It loads the lab. data automatically to facilitate user's comparison, peer review and determine model performance. The field data have been saved on files Iwanai_Esashi.mat and Be sure that the lab data files Iwanai_Esashi.mat resides in the same directory or subdirectory as the Matlab script (BP9_Field_GAUGE.m). The user's input files 'filename1' and 'filename2' must be present in the same directory or subdirectory 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 field gauge data. Examples: > BP9_Field_GAUGE In this example the Matlab script plots field data at Iwanai and Esashi. > BP9_Field_GAUGE('NEOWAVE_BP9_IWANAI.dat','NEOWAVE_BP9_ESASHI.dat') In this example the Matlab script plots both user's model result cases against field data, (Iwanai and Esashi). > BP9_Field_GAUGE('NEOWAVE_BP9_IWANAI.dat',[]) In this example the Matlab script just plots user's model results against field data for Iwanai case. > BP9_Field_GAUGE([],'NEOWAVE_BP9_ESASHI.dat') In this example the Matlab script just plots user's model results against field data for Esashi case. -------------------------------------------------------------------------