!     Last change:  YY   13 Nov 2009    2:00 am
!---------------------------------------------------------------------------------------------
!                               NEOWAVE 2D include file
!
!                 INPUT SETUP and GLOBAL VARIABLES DECLARATIONS MODULE
!---------------------------------------------------------------------------------------------
!
!   This block setup all input paramters (except for bathymetry data file names and initial
!   condition inputs), and declare all global variables
!
!   Z.Kowalik  2006 IMS,      University of Alaska Fairbanks
!   Y.Yamazaki 2007 ORE,      University of Hawaii Manoa
!---------------------------------------------------------------------------------------------
      INTEGER, PARAMETER::  dp = KIND(1.0d0)

!---------------------------------------------------------------------------------------------
! INPUT
!---------------------------------------------------------------------------------------------




! Physical Quantities
!---------------------------------------------------------------------------------------------
      REAL, PARAMETER :: g = 980.0                     ! g(cm^2/s):gravitational acceleration

      REAL, PARAMETER :: R0  = 6370.0E5                ! (cm): Earth radius
!     REAL, PARAMETER :: omega = 7.2919E-5             ! omega(1/s)  : angular velocity of earth
      REAL, PARAMETER :: FC   = 1.4584E-4              ! = 2*omega
      REAL, PARAMETER :: EPSILON = 0.001               ! =1.e-3 (ORIGINAL SETUP)

!... OKADA bottom dislocation parameters
      REAL(dp), PARAMETER :: nu     = 1.0d0            ! used in OKADA solution (no need to change)
      REAL(dp), PARAMETER :: lambda  = 1.0d0           ! used in OKADA solution (no need to change)


! Numerical Quantities
!---------------------------------------------------------------------------------------------
      REAL, PARAMETER::  pi = 3.14159265358979    ! pi = 3.14
      REAL, PARAMETER:: rad  = 0.0174532930056254 ! = pi/180
      REAL, PARAMETER:: zero = 0.0E0              ! 0
      REAL, PARAMETER:: one  = 1.0E0              ! 1



! 1st Grid SETUP
!---------------------------------------------------------------------------------------------


      REAL,   PARAMETER :: dlon = 30.0 /3600.       ! (deg) : x grid size
      REAL,   PARAMETER :: dlat = 30.0 /3600.       ! (deg) : y grid size
      REAL,   PARAMETER :: dx   = R0*dlon*rad       ! (cm)  : x grid size
      REAL,   PARAMETER :: dy   = R0*dlat*rad       ! (cm)  : y grid size

!... BOTTOM and LEFT ends coordinates
      REAL,   PARAMETER :: lonL = 185.0             ! (deg) : LEFT   END Longitude of 1st Domain
      REAL,   PARAMETER :: latB = -18.5             ! (deg) : BOTTOM END Latitude  of 1st Domain

      INTEGER,PARAMETER :: js   = 1                 ! (-)   : j START
      INTEGER,PARAMETER :: je   = 721               ! (-)   : j END
      INTEGER,PARAMETER :: ks   = 1                 ! (-)   : k START
      INTEGER,PARAMETER :: ke   = 781               ! (-)   : k END

     ! INTEGER,PARAMETER:: NJK1  = (je1-js1+1)*(ke1-ks1+1)


!      REAL, DIMENSION(js1:je1,ks1:ke1):: z01        ! z0  : initial water level at t=0
!      REAL, DIMENSION(js1:je1,ks1:ke1):: z01_1      ! z0_1: z0 obtained from fault No1

      REAL, DIMENSION(js:je, ks:ke) :: zo, zn        ! zeta: water level

      REAL, DIMENSION(js:je,ks:ke) :: H              ! H   : still water depth (at z cal point)

      REAL, DIMENSION(ks:ke) :: siz, siv, csz, csv   ! sin cos for spherical coord. cal
      INTEGER,PARAMETER :: nsmooth= 20
! n#smooth is the number of points to
!  smooth the sea surface near the
!-----  !  boundaries of the area where initial    ----------------------------------------------------------------------------------------
        !  static surface deformation is applied.
