call readi(mcmcard,'NCONF_IN',nconf_in,0)
call reada(mcmcard,'RDIH_BIAS',rdih_bias,0.5d0)
write (iout,*) "NCONF_IN",nconf_in
+ tm_score=(index(mcmcard,'TMSCORE').gt.0)
+ if (tm_score) write (iout,*) "Using TM_Score instead of DIFF",
+ & " for torsional angles"
return
end