BP#10  

  -------------------------------------------------------------------------
                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:  Rat Islands Tsunami
  
                TYPE:           Field Benchmarking Problem (Field)
                                NUMBER: 10
                TSUNAMI SOURCE: Coseismic wave (CS)                       
                SCRIPT NAME   : BP10_Field_CS02.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.
 
                    Titov, V.V., F.I. González, E.N. Bernard, M.C. Eble,
                H.O. Mofjeld, J.C. %Newman, and A.J. Venturato (2005):
                Real-time tsunami forecasting: Challenges and solutions.
                Nat. Hazards, 35(1), Special Issue, U.S. National Tsunami
                Hazard Mitigation Program, 41–58.
  
                    NOAA/PMEL Center for Tsunami research
                http://nctr.pmel.noaa.gov/benchmark/index.html
  
                INCLUDED FILE WITH THIS SCRIPT: 
                CS_01_gage.mat (Hilo gage time serie, Nov. 17/2003.
                                Rat Islands tsunami)
 
                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/
  
  -------------------------------------------------------------------------
    > BP10_Field_CS02(filename)
 
    BP10_Field_CS02 reads input files 'filename' containing 
    tide gage dataset recorded at Hilo Bay on Nov/17/2007 during the Rat
    Islands Tsunami. This Benchmark case is described in detail in the
    following URL adress 
    NOAA/PMEL(http://nctr.pmel.noaa.gov/benchmark/index.html).
 
    filename is a user numerical model result file containig dataset of
    surface elevation at Hilo
 
    Input files 'filename' needs to be named as follows:
    XXXXXX_CS_02.dat where XXXXXX is
    the name of the numerical model that is used,
    eg. ADCIRC_CS_02.dat.
  
    The input file, 'filename' must have the following
    configuration or structure:
 
    1- File 'filename' is a matrix ASCII files of size
       180x2(row x column)
    2- Column 1 is time given in second, ranging form 
       [5.015 to 7.998333333333333] with increment of
       0.01666666667 hour (60 sec).
    3- Columns 2 is user numerical results (surface elev.) in meter at
       Hilo's gage
    +-----------------+------------------+
    | Time(hours)     |      Gage(m)     |
    +-----------------+------------------+
    5.015000000000000                   0
    5.031666666666666   0.012355555555555
    5.048333333333333  -0.009605714285714
    5.065000000000000   0.016423490794919
    5.081666666666667   0.008023499376441
    5.098333333333334  -0.006937938878834
            .                  . 
            .                  .
            .                  .
    7.898333333333333   0.014060703969055
    7.915000000000000   0.007230160026670
    7.931666666666667  -0.008925599455011
    7.948333333333333   0.006420027591456
    7.965000000000000  -0.009392380952381
    7.981666666666667  -0.005177777777778
    7.998333333333333                   0
    +-----------------+------------------+   
    
    NOTE:
    This script is self-containing. It loads the field data automatically to
    facilitate user's comparison, peer review, and determine model
    performance.
    The field data are saved on file CS_02_gage.mat.
    Be sure that the field data file CS_02_gage resides
    in the same directory or subdirectory as the Matlab script
    BP10_Field_CS02.m
    The user's input file 'filename' must be present in the same directory
    or subdirectory of the script file BP10_Field_CS02.m    
  
    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. 
 
    Examples:
    > BP10_Field_CS02
      In this example the Matlab scrip only plots laboratory results.
    > BP10_Field_CS02('ADCIRC_CS_02.dat')
      In this example the Matlab script plots recorded Field dataset
      against user's model results.