- do j=1,nss
- if (dyn_ss) then
- call xdrfint(ixdrf, idssb(j), iret)
- call xdrfint(ixdrf, jdssb(j), iret)
- idssb(j)=idssb(j)-nres
- jdssb(j)=jdssb(j)-nres
- else
- call xdrfint_(ixdrf, ihpb(j), iret)
- if (iret.eq.0) goto 101
- call xdrfint_(ixdrf, jhpb(j), iret)
- if (iret.eq.0) goto 101
- endif
- enddo
+ if (dyn_ss) then
+ do k=1,nss
+ call xdrfint(ixdrf, idssb(k), iret)
+ call xdrfint(ixdrf, jdssb(k), iret)
+ ihpb(k)=iss(idssb(k)-nres)+nres
+ jhpb(k)=iss(jdssb(k)-nres)+nres
+#ifdef DEBUG
+ write (iout,*) "jj",jj+1," dyn_ss:",idssb(k)-nres,
+ & jdssb(k)-nres,ihpb(k),jhpb(k)
+#endif
+ enddo
+ else
+ do k=1,nss
+ call xdrfint(ixdrf, ihpb(k), iret)
+ if (iret.eq.0) goto 101
+ call xdrfint(ixdrf, jhpb(k), iret)
+ if (iret.eq.0) goto 101
+#ifdef DEBUG
+ write (iout,*) "jj",jj+1," stat_ss:",ihpb(k),jhpb(k)
+#endif
+ enddo
+ endif