_RL oval,omin,omax
oval = (2. _d 0 * dfs1)**2
if (oval.lt.omin) omin = oval
if (oval.gt.omax) omax = oval