with_dihed_constr = index(controlcard,"WITH_DIHED_CONSTR").gt.0
write (iout,*) "with_dihed_constr ",with_dihed_constr
call readi(controlcard,'CONSTR_DIST',constr_dist,0)
+ dyn_ss=index(controlcard,"DYN_SS").gt.0
return
end
c------------------------------------------------------------------------------