IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      PARAMETER (KX=200,KY=200,KKX=KX,KKY=KY)

      REAL*8 LX,LY
      REAL*8 NEX,NEY,NWX,NWY,NNX,NNY,NSX,NSY
      REAL*8 M,ME,MW,MN,MS
      LOGICAL PAS_FINI

C VARIABLES
C =========

      DIMENSION GAM(KX,KY),RO(KX,KY)
      DIMENSION U(KX,KY),V(KX,KY),U_(KX,KY),V_(KX,KY)
      DIMENSION SU(KX,KY)
      DIMENSION FI(0:KX,0:KY),FI0(0:KX,0:KY),FI00(0:KX,0:KY)
      DIMENSION GRADX(KX,KY),GRADY(KX,KY)

      DIMENSION SQ(KX,KY)
      DIMENSION Q(0:KX,0:KY),Q0(0:KX,0:KY),Q00(0:KX,0:KY)

      DIMENSION SVEL(KX,KY),VEL(0:KX,0:KY)
      DIMENSION RHS(KX,KY)

      DIMENSION EX(KX,KY),EY(KX,KY)
      DIMENSION Uf(KX,KY),Vf(KX,KY),U0f(KX,KY),V0f(KX,KY)
      
      DIMENSION FIANAL(0:KX,0:KY)

      DIMENSION RES(KX,KY),RESOR(5),SOR(5),ITSOR(5),NSW(5)

C COEFFS
C ======

      DIMENSION AE(KX,KY),AW(KX,KY),AN(KX,KY),AS(KX,KY),AP(KX,KY)
        
C FLUX
C =====

      DIMENSION ME(KX,KY),MW(KX,KY),MS(KX,KY),MN(KX,KY)
      DIMENSION DE(KX,KY),DW(KX,KY),DS(KX,KY),DN(KX,KY)

C METRIC
C ======

      DIMENSION X(KX),Y(KY)
      DIMENSION DV(KX,KY)
      DIMENSION NEX(KKX,KKY),NWX(KKX,KKY),NSX(KKX,KKY),NNX(KKX,KKY)
      DIMENSION NEY(KKX,KKY),NWY(KKX,KKY),NSY(KKX,KKY),NNY(KKX,KKY)
      DIMENSION SE(KKX,KKY),SW(KKX,KKY),SN(KKX,KKY),SS(KKX,KKY)
      DIMENSION DXE(KKX,KKY),DXW(KKX,KKY),DYN(KKX,KKY),DYS(KKX,KKY)

C COMMON
C ======

      COMMON /PARAM/ DT,IMA
      COMMON /INDICE/ NICV,NJCV ,NI,NJ,NIM,NJM,NTM,KSWEEP
      COMMON /PARAM_SOL/ RES,RESOR,SOR,EPS,RES2,NSW,ITMAX,ITSOR
      COMMON /METRIC/  NEX,NEY,NWX,NWY,NNX,NNY,NSX,NSY,
     &                 SE,SW,SN,SS,
     &                 DXE,DXW,DYN,DYS,
     &                 DV,X,Y,
     &                 DX,DY,LX,LY

      COMMON/DIVERS/ AMPLI
      COMMON/COEFFS/ AE,AW,AN,AS,AP

      COMMON/FLUX/ ME,MW,MN,MS,DE,DW,DN,DS
      
      COMMON /BOUND/ IBCW1,IBCE1,IBCS1,IBCN1,
     &               IBCW2,IBCE2,IBCS2,IBCN2

      COMMON /ELEC/ EX,EY
      COMMON /ELEC_PARAM/ T,C,M