endif
open(9,file=pdbfile)
nnt = 1
- if (itype(1).eq.21) nnt = 2
+ if (itype(1).eq.ntyp1) nnt = 2
nct=nres
- if (itype(nres).eq.21) nct = nres-1
+ if (itype(nres).eq.ntyp1) nct = nres-1
print *,"nnt",nnt," nct",nct
print *,"file",arg
call xdrfopen(ixdrf,arg, "r", iret)