OpenMP parallel loop crashes (?)

Hello everybody,

I am trying to make the section of my code to run parallel:

....

     

     

      
    EL=0.0d0

 !$OMP parallel DO SHARED(S,COUL) PRIVATE(I1,J1,ID,JD) reduction(+:EL) 
   
    DO J1=1,NY
    DO I1=1,NX
  
 IF ((J1/=J.OR.I1/=I).AND.(J1/=J.OR.I1/=IP(I)).AND.(J1/=J.OR.I1/=IM(I)).AND.(J1/=JP(J).OR.I1/=I).AND.(J1/=JM(J).OR.I1/=I)) THEN 
 

    IF (ABS(FLOAT(I)-FLOAT(I1)) <= ABS(FLOAT(I)+LLEN-FLOAT(I1))) THEN 
    ID= INT(ABS(FLOAT(I)-FLOAT(I1)))
    ELSE 
    ID= INT(ABS(FLOAT(I)+LLEN-FLOAT(I1)))
    END IF
 
    IF (ABS(FLOAT(J)-FLOAT(J1)) <= ABS(FLOAT(J)+LLEN-FLOAT(J1))) THEN 
    JD= INT(ABS(FLOAT(J)-FLOAT(J1)))
    ELSE 
    JD= INT(ABS(FLOAT(J)+LLEN-FLOAT(J1)))
    END IF

 
    EL= EL + LAMDA*COUL(ID,JD)*dble(S(I1,J1))
 !   Cen(I,J)= Cen(I,J) +  LAMDA*COUL(ID,JD)*dble(S(I1,J1))
 
 END IF
 
   
    END DO
    END DO
 
 !$OMP END PARALLEL DO

...

where COUL is a matrix determined earlier in the code.

I get no compilation or build errors but at run time the program exits when it enters the parallel loop. It just crashes with no run-time error!

Any ideas?

Thanks,

Marios