c	subroutina para calcular el maximo de los valores de una variable
c	en simple precision

c       jul 2011 malv+fgg

	subroutine maxsp(a,n,ymax)
	real a(n)

	ymax=-1.0e+34
	do i=1,n
	 if(a(i).gt.ymax)ymax=a(i)
	end do
	return
	end