+ ici=icys(i)
+ icj=icys(j)
+ if (ici.eq.0 .or. icj.eq.0) then
+#ifdef MPI
+ write (*,'(a,i5,2a,a3,i5,5h and ,a3,i5)')
+ & "Processor",me," attempt to create",
+ & " a disulfide link between non-cysteine residues ",restyp(i),i,
+ & restyp(j),j
+ call MPI_Abort(MPI_COMM_WORLD,ierr)
+#else
+ write (*,'(a,i5,2a,a3,i5,5h and ,a3,i5)')
+ & "Processor",me," attempt to create",
+ & " a disulfide link between non-cysteine residues ",restyp(i),i,
+ & restyp(j),j
+ stop
+#endif
+ endif