[Nek5000-users] Segmentation fault occuerd

[nek5000-users at lists.mcs.anl.gov]

Hello,

I have some trouble to run my simulation... When I start running my 
simulation, I receive the following error-message (I compiled the run 
in debugging mode -g and with traceback):

forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        
Source             libc.so.6          00002B76C2128FA0  Unknown         
       Unknown  Unknown
nek5000            0000000000500D48  chcopy_                   434  math.f
nek5000            000000000071BA91  exitti_                   399  comm_mpi.f
nek5000            0000000000565F6A  set_conv_char_            603  convect.f
nek5000            0000000000561C80  setup_convect_             44  convect.f
nek5000            00000000004E805C  setics_                   393  ic.f
nek5000            0000000000417B2A  nek_init_                 164  drive1.f
nek5000            0000000000415FCD  MAIN__                     26  drive.f
nek5000            0000000000415F8C  Unknown               Unknown  Unknown
libc.so.6          00002B76C20D14CA  Unknown               Unknown  Unknown
nek5000            0000000000415EBA  Unknown               Unknown  Unknown
--------------------------------------------------------------------------
mpirun has exited due to process rank 0 with PID 10816 on
node ifh-taifun exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------


By debbuging I've seen, that follwoing values in the subroutine 
set_conv_char (in convect.f) seem to be a little strange:
nxd=           0
nyd=     7457184
nzd=           0

In my DEALIAS-file the follwing line is commented out:
C      common /dedim/ nxd,nyd,nzd


I tried to run the simulation with the declaration in the DEALIAS-file 
and a      include 'DEALIAS'
in the subroutine set_conv_char (in convect.f) but the the run also is dying:

   ZWGJD: Minimum number of Gauss points is 1 -1633391248
  call exitt: dying ...


Do you maybe know where there might be my buck??


Thx a lot.

Fred


Here my SIZEu-file:
C     Dimension file to be included
C
C     HCUBE array dimensions
C
      PARAMETER (LDIM=3)
      PARAMETER (LX1=6,LY1=6,LZ1=6,LELT=4420,LELV=4420)
      parameter (lxd=9,lyd=lxd,lzd=lxd)
      parameter (lelx=6,lely=lelx,lelz=lelx)

      PARAMETER (LZL=3)
      PARAMETER (LX2=LX1-2)
      PARAMETER (LY2=LY1-2)
      PARAMETER (LZ2=LZ1-2)
      PARAMETER (LX3=LX1)
      PARAMETER (LY3=LY1)
      PARAMETER (LZ3=LZ1)

      parameter (lpelv=1,lpelt=1,lpert=1)        ! perturbation
      parameter (lpx1=1,lpy1=1,lpz1=1)           ! array sizes
      parameter (lpx2=1,lpy2=1,lpz2=1)

      parameter (lbelv=1,lbelt=1)                ! MHD
      parameter (lbx1=1,lby1=1,lbz1=1)           ! array sizes
      parameter (lbx2=1,lby2=1,lbz2=1)

C LX1M=LX1 when there are moving meshes; =1 otherwise
      PARAMETER (LX1M=1,LY1M=1,LZ1M=1)
      PARAMETER (LDIMT= 1)
c     PARAMETER (LDIMT= 3)
      PARAMETER (LDIMT1=LDIMT+1)
      PARAMETER (LDIMT3=LDIMT+3)
      PARAMETER (LP = 8)
      PARAMETER (LELG   = LP*LELT)
c
c
c     Note:  In the new code, LELGEC should be about sqrt(LELG)
c
      PARAMETER (LELGEC = 1)
      PARAMETER (LXYZ2  = 1)
      PARAMETER (LXZ21  = 1)
c
      PARAMETER (LMAXV=LX1*LY1*LZ1*LELV)
      PARAMETER (LMAXT=LX1*LY1*LZ1*LELT)
      PARAMETER (LMAXP=LX2*LY2*LZ2*LELV)
      PARAMETER (LXZ=LX1*LZ1)
      PARAMETER (LORDER=3)
C FF      PARAMETER (MAXOBJ=2,MAXMBR=LELT*6)
      PARAMETER (MAXOBJ=4,MAXMBR=LELT*6,lhis=100)

C
C     Common Block Dimensions
C
      PARAMETER (LCTMP0 =2*LX1*LY1*LZ1*LELT)
      PARAMETER (LCTMP1 =4*LX1*LY1*LZ1*LELT)
C
C     The parameter LVEC controls whether an additional 42 field arrays
C     are required for Steady State Solutions.  If you are not using
C     Steady State, it is recommended that LVEC=1.
C
      PARAMETER (LVEC=1)
C
C     Uzawa projection array dimensions
C
C FF      PARAMETER (MXPREV = 10)
      parameter (mxprev = 20)
      parameter (lgmres = 40)

C
C     Split projection array dimensions
C
      parameter(lmvec = 1)
      parameter(lsvec = 1)
      parameter(lstore=lmvec*lsvec)
c
c     NONCONFORMING STUFF
c
      parameter (maxmor = lelt)
C
C     Array dimensions
C
      COMMON/DIMN/NELV,NELT,NX1,NY1,NZ1,NX2,NY2,NZ2
     $,NX3,NY3,NZ3,NDIM,NFIELD,NID

c automatically added by makenek
      parameter(lxo   = lx1) ! max output grid size (lxo>=lx1)

c automatically added by makenek
      parameter(lpart = 1  ) ! max number of particles

c automatically added by makenek
      integer ax1,ay1,az1,ax2,ay2,az2
      parameter (ax1=lx1,ay1=ly1,az1=lz1,ax2=lx2,ay2=ly2,az2=lz2) ! 
running averages