refstr = index(controlcard,'REFSTR').gt.0
pdbref = index(controlcard,'PDBREF').gt.0
dyn_ss=(index(controlcard,'DYN_SS').gt.0)
+ call readi(controlcard,'NSAXS',nsaxs,0)
+ call readi(controlcard,'SAXS_MODE',saxs_mode,0)
+ write (iout,*) "Number of SAXS restraints",NSAXS," SAXS_MODE",
+ & SAXS_MODE
C /06/28/2013 Adasko: dyn_ss is keyword allowing to break and create bond
C disulfide bond. Note that in conterary to dynamics this in
C CONTROLCARD. The bond are read in molread_zs.F