& index(controlcard,"PRINT_HOMOLOGY_MODELS").gt.0
call readi(controlcard,'NSAXS',nsaxs,0)
call readi(controlcard,'SAXS_MODE',saxs_mode,0)
+ call reada(controlcard,'SCAL_RAD',scal_rad,1.0d0)
write (iout,*) "Number of SAXS restraints",NSAXS," SAXS_MODE",
- & SAXS_MODE
+ & SAXS_MODE," SCAL_RAD",scal_rad
if (min_var) iopt=1
return
end