anchor.diff:  !$OMP END CRITICAL
   
   ! test if PARALLEL is treated as a single unit of work
  -!$OMP PARALLEL
  -        scalar3=scalar3+1
  -!$OMP END PARALLEL
  +! laks: replace parallel to single
  +!!$OMP MASTER
  +!        scalar3=scalar3+1
  +!!$OMP END MASTER
   
           WHERE ( CC .ne. 0 ) DD = 1.0/CC
   
  
positive.code: !sequential equivalent statements for comparison
         BB=BB+1
         scalar02=scalar02+1
         scalar12=scalar12+1
         scalar22=scalar22+1
         scalar32=scalar32+1
  
  !      write (1,*) "ck:sum of AA is",SUM(AA)," sum of BB is ",sum(BB)
         if (SUM(AA)/=SUM(BB)) then
  
negative.code: !$OMP END CRITICAL
  
  ! test if PARALLEL is treated as a single unit of work
  ! laksono: replace paralell with single directive
  !!$OMP MASTER
  !        scalar3=scalar3+1
  !!$OMP END MASTER
  
          WHERE ( CC .ne. 0 ) DD = 1.0/CC
  
  