changed routine for initial values for gevd.fit. improved inital value for eta through second fit. still some work do here though. added NA value detection at beginning of gev.d.fit function.