double precision efreeref(maxref,maxT,maxprot) common /FREEREF/ efreeref