!     Last change:  ZK   13 May 2009    3:37 pm
      Implicit NONE

      INTEGER,PARAMETER :: js=1, je=2500 ! JMAX=je-js+1
      INTEGER,PARAMETER :: KS=1, ke=2000
      REAL, PARAMETER :: R0=6370.E5
      REAL, PARAMETER :: RAD=  0.0174532930056254
      REAL, PARAMETER :: DL= 1.79892577E-05   ! STEP ALONG LONGITUDE  IN DEGREES  FOR 1M STEP
      REAL, PARAMETER :: DF= 8.99462793E-06   ! GRID STEP ALONG LATITUDE  IN DEGREES (about 1 m)
      REAL, PARAMETER :: HX =R0*DL*RAD        !SPACE STEP ALONG 1DEG IN CM
      REAL, PARAMETER :: HY =R0*DF*RAD        !SPACE STEP

  ! REAL, PARAMETER::HX =50000.               !SPACE STEP ALONG 1DEG IN CM
  ! REAL, PARAMETER::HY =50000.
  ! INTEGER,PARAMETER::ME=2

      INTEGER, PARAMETER :: ME=10002           ! first 10 sec

  ! INTEGER,PARAMETER::ME=180000
 
      REAL, PARAMETER :: T = 0.001            !TIME STEP FOR TEST CASE T=0.02,FOR BEACH T=0.1
      REAL, PARAMETER :: G = 980.0
      REAL, PARAMETER :: R = 0.003
      REAL, PARAMETER :: FC=1.4584E-4

  !  REAL::FCT !FC*T

      REAL, PARAMETER :: epsilon =0.01        ! originally e=1.e-3
      REAL, PARAMETER :: ALPHA =0.08
      REAL, PARAMETER :: TPF=1500.

  !  REAL::THX
  !  REAL::THY
  !  REAL::GTHX
  !  REAL::GTHY

      REAL, PARAMETER :: PI = 3.1415926
      REAL, PARAMETER :: FCT= FC*T
      REAL, PARAMETER :: THX=T/HX
      REAL, PARAMETER :: THY=T/HY
      REAL, PARAMETER :: GTHX=G*THX
      REAL, PARAMETER :: GTHY=G*THY
      REAL, PARAMETER :: R1 = 0.001
      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, PARAMETER :: GTHYDD=GTHY*DRR
      REAL, PARAMETER :: GTHYDR=GTHY*DREL
      REAL, PARAMETER :: NS=1.
      REAL, PARAMETER :: NL=0.  ! spherical nonlinearity
      REAL, DIMENSION(1:JE, 1:KE) :: SLO1, SLN1, SLNN, UO1, UN1, h1, D1, DHU, DHV
      REAL, DIMENSION(1:JE, 1:KE) :: VO1, VN1, ZMAX, TRANSX, TRANSY, SLNN2
      REAL, DIMENSION(1:JE, 1:KE) :: SLO2, SLN2, UO2, UN2, D2, VO2, VN2
      REAL, DIMENSION(1:JE, 1:KE) :: LMTEMY, LPTEMY, LMTEMX, LPTEMX,LMTEMXY
      REAL, DIMENSION(1:KE+1) :: SIZ, SIV, CSZ, CSV
      REAL, DIMENSION(1:ME) :: ZINP
      INTEGER*2, DIMENSION(1:JE, 1:KE) :: CELL1, CELL1X, CELL2
      REAL :: SLMA
      INTEGER :: M, IC, KK, K, J, JM, KM
