else if (molecule.eq.2) then
do i=1,ntyp1_molec(molecule)
print *,nam(1:1),restyp(i,molecule)(1:1)
+ print *,nam(2:2),"read",i
if (nam(2:2).eq.restyp(i,molecule)(1:1)) then
rescode=i
return
stop
else if (molecule.eq.5) then
do i=1,ntyp1_molec(molecule)
- print *,i,restyp(i,molecule)(1:2)
+ print *,restyp(i,molecule)
+ print *,i,restyp(i,molecule)(1:2),nam(1:2)
if (ucase(nam(1:2)).eq.restyp(i,molecule)(1:2)) then
rescode=i
return