BP#6 RUNUP
  -------------------------------------------------------------------------
                NTHMP (National Tsunami Hazard Mitigation Program)
  -------------------------------------------------------------------------
                           2011 NTHMP VALIDATION WORKSHOP
               Venue: Galveston, TAMU TX
                      Margh 30- Apr. 1
                                ---- o ----
                MATLAB SCRIPT:  Solitary Wave on a Conical Island
                TYPE:           Laboratory Benchmarking Problem (Lab)
                                NUMBER: 06
                TSUNAMI SOURCE: Incident Solitary Wave (SW)                       
                SCRIPT NAME   : BP6_Lab_Runup.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: 
 
                conical_RU.mat (lab.data (experimental) file for the
                spatial runup on the sloping beach of the conical island. 
                 
                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/
  
  -------------------------------------------------------------------------
    >> BP6_Lab_Runup(filename1,filename2,filename3)
 
    BP6_Lab_Runup reads input files 'filename1','filename2','filename3' 
    containing model results and compares them with experimental data described
    in the Benchmark "Solitary wave on a conical island", NOAA/PMEL,
    (http://nctr.pmel.noaa.gov/benchmark/index.html).
 
    filename1, filename2 and filename3 contain model results for
    Case(A)= A/h =0.045, Case(B)= A/h=0.096,and Case(C)= A/h=0.181, respectively.
     
    Input files 'filename1', 'filename2', 'filename3' need to be named as follows:
    XXXX_BP6_RU_A.dat, XXXX_BP6_RU_B.dat and XXXX_BP6_RU_C.dat, 
    where XXXX is the name of the numerical model that is tested,
    eg. NEOWAVE_BP6__RU_A.dat, NEOWAVE_BP6__RU_B.dat and NEOWAVE_BP6__RU_C.dat
  
    The input files, 'filename1', 'filename2', and 'filename3' must have the following
    configuration:
 
    1- Files 'filename1' is a matrix ASCII files of size
       361x2(row x column) as shown below (no heading).
       The number of rows could be different of 361, but it is
       desirable to sample at least between 1 to 10 deg.     
       x ranges= [0 360]  with icrements of 10.0.
       
  filename1
    +-------+---------------------+
    | deg(°)   Runup(A/h=0.045) cm   
    +-------+---------------------+
      0.00            2.140  
     10.00            1.812    
     20.00            1.300    
     30.00            3.540   
       .                .      
       .                .     
       .                .     
    330.00            2.546  
    340.00            3.456  
    350.00            4.679  
    360.00            3.456  
    +-------+---------------------+
 
    2- Files 'filename2' is a matrix ASCII files of size
       361x2(row x column) as shown below (no heading).
       The number of rows could be different of 361, but it is
       desirable to sample at least between 1 to 10 deg.     
       x ranges= [0 360]  with icrements of 10.0.
       
  filename2
    +-------+---------------------+
    | deg(°)   Runup(A/h=0.096) cm   
    +-------+---------------------+
      0.00            2.140  
     10.00            1.812    
     20.00            1.300    
     30.00            3.540   
       .                .      
       .                .     
       .                .     
    330.00            2.546  
    340.00            3.456  
    350.00            4.679  
    360.00            3.456  
    +-------+---------------------+
 
    3- Files 'filename3' is a matrix ASCII files of size
       361x2(row x column) as shown below (no heading).
       The number of rows could be different of 361, but it is
       desirable to sample at least between 1 to 10 deg.     
       x ranges= [0 360]  with icrements of 10.0.
       
  filename3
    +-------+---------------------+
    | deg(°)   Runup(A/h=0.181) cm   
    +-------+---------------------+
      0.00            2.140  
     10.00            1.812    
     20.00            1.300    
     30.00            3.540   
       .                .      
       .                .     
       .                .     
    330.00            2.546  
    340.00            3.456  
    350.00            4.679  
    360.00            3.456  
    +-------+---------------------+
 
    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 have been saved on file as the conical_RU.mat 
    for Case(A)= A/h =0.045, Case(B)= A/h=0.096, and Case(C)= A/h=0.181
    , respectively.
    Be sure that the lab data files conical_RU.mat 
    reside in the same directory or subdirectory as the Matlab script
    (BP6_Lab_Runup.m).
    The user's input files 'filename1' ,' filename2',and 'filename3'
    must be present in the same directory or subdirectory.
 
    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 analytical solution.
  
    Examples:
    > BP6_Lab_Runup
      In this example the Matlab script just plots experimental results.
    > BP6_Lab_Runup('NEOWAVE_BP6_RU_A.dat','NEOWAVE_BP6_RU_B.dat','NEOWAVE_BP6_RU_C.dat')
      In this example the Matlab script plots experimental results
      and user's model results.
    > BP6_Lab_Runup('NEOWAVE_BP6_RU_A.dat',[],[])
      In this example the Matlab script just plots case A (experimental results
      and user's model results).