BP#9 RUNUP
  -------------------------------------------------------------------------
                NTHMP (National Tsunami Hazard Mitigation Program)
  -------------------------------------------------------------------------
                           2011 NTHMP VALIDATION WORKSHOP
               Venue: Galveston, TAMU TX
                      Margh 30- Apr. 1
                                ---- o ----
                MATLAB SCRIPT:  Okushiri Island Runup
                TYPE:           Field Benchmarking Problem (Field)
                                NUMBER: 09
                TSUNAMI SOURCE: Field fata Runup (RUNUP)                       
                SCRIPT NAME   : BP9_Field_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.  
          
  
                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_RUNUP(filename1)
 
    BP9_Field_RUNUP reads input files 'filename1' 
    containing model results and compares with runup field data
    (Prof Imamura and Shuto... before we used Prof. Tsuji data).
     
    Input files 'filename1' need to be named as follows:
    XXXX_BP9_OKUSHIRI.dat, where XXXX is the name of the numerical model
    to be tested,
    eg. NEOWAVE_BP9_OKUSHIRI.dat
  
    The input file, 'filename1' must have the following
    configuration:
    'filename1' is a matrix ASCII files of size
    19x5(row x column) as shown below (no heading).
 
    filename1 
    +-----------+-----------+-----------+--------------+  
    |                       Model result Model Result
    | long (Deg)  lat (Deg)    range(m)   Avg(m)
    |                        (min - max)    
    +-----------+-----------+-----------+--------------+ 
     139.4292117 42.18818149   5.5  9.3   7.1
     139.4111857 42.16276287   7.0 10.2   8.5
     139.4182612 42.13740439   8.0 12.6  10.7
     139.4280358 42.09301238  11.5 15.3  13.9
     139.426245  42.11655479  11.0 15.4  12.9
     139.4237147 42.10041415   9.6 17.1  11.8
     139.4289018 42.07663658  10.6 15.8  13.2
     139.4278534 42.06546152  10.4 17.2  13.6
     139.4515399 42.04469655  11.6 16.0  12.9
     139.4565284 42.05169226  5.2  12.0   6.7
     139.4720138 42.05808988  8.1  11.8  10.3
     139.5150461 42.21524909  5.1   5.9   5.4
     139.5545494 42.22698164  4.1   8.6   6.8
     139.4934307 42.06450128  7.4  14.7  10.8
     139.5474599 42.18744879  2.2   5.3   3.8
     139.5258982 42.17101221  6.4   9.0   7.4
     139.5625242 42.21198369  2.6   3.8   3.1
     139.5190997 42.11305805  6.2   8.8   7.7
     139.5210766 42.15137635  2.9   4.7   3.6
    +-----------+-----------+-----------+--------------+ 
 
    NOTE:
    Modelers need only to change columns 3 to 5 with their
    model result data
 
    Examples:
    > BP9_Field_RUNUP
      In this example the Matlab script just plots field data range.
    > BP9_Field_RUNUP('NEOWAVE_BP9_OKUSHIRI.dat')
      In this example the Matlab script plots model results
      and compares against field data range.
  -------------------------------------------------------------------------