!.. Last change:  ZK    8 Mar 2009    1:40 pm

      Implicit NONE

      INTEGER,PARAMETER :: js=1, JE=2500, jmax=je-1     !FOR TEST (5M) JE=10041;
      INTEGER,PARAMETER :: ME=1000
      INTEGER,PARAMETER :: MW=10
      INTEGER :: surfer_index

!     WHEN WOULD YOU LIKE TO GENERATE TSUNAMI
      INTEGER,PARAMETER :: MGEN=100
      REAL, PARAMETER :: G = 980.0
      REAL, PARAMETER :: HX = 100.      !SPACE STEP
      REAL, PARAMETER :: T = 0.001      !TIME STEP FOR TEST CASE T=0.02,FOR BEACH T=0.1
      REAL, PARAMETER :: R = 0.003
      REAL, PARAMETER :: R1 = 0.001
      REAL, PARAMETER :: epsilon =0.001
      REAL, PARAMETER :: TP=1.
      REAL, PARAMETER :: THX=T/HX
      REAL, PARAMETER :: GTHX=G*THX
      REAL, PARAMETER :: rho1=1.
      REAL, PARAMETER :: rho2=1.4
      REAL, PARAMETER :: drr=(rho2-rho1)/rho2
      REAL, PARAMETER :: drel=rho1/rho2
      REAL, PARAMETER :: GTHXDD=GTHX*DRR
      REAL, PARAMETER :: GTHXDR=GTHX*DREL
      REAL, DIMENSION(1:JE) :: SLO1, SLN1, UO1, UN1, h1, DETA, D1, DHU
      REAL, DIMENSION(1:JE) :: ETN, ETO, VELM, slma
      REAL, DIMENSION(1:JE) :: SLO2, SLN2, UO2, UN2, D2, TRANS2
      INTEGER, DIMENSION(1:JMAX) :: Cell1, CELL2
      REAL, PARAMETER :: PI=3.141592653589793238462643383279502884197
      REAL :: SLMI, ZNC
      INTEGER :: M, J