PARAMETER (NGRIDc =1024)      ! zero-level mesh in 1D
      PARAMETER (mARR =NGRIDc+1)   !      need in FFT
      PARAMETER (mF67=NGRIDc/2)       !      need in FFT
      COMMON /F67COM/
     +                 IBC,      IP,       ISL,     L1,     N2,
     +                 N3,       N4,        N7,
     +                 SI(mF67),    INDEX(mF67)