cccccccccccc
c	subroutina para calcular el maximo de los valores de una variable a
c	en dble precision
	subroutine maxdp(a,n,ymax)
c       jul 2011 malv+fgg
cccccccccccc
	implicit none
	integer n,i
	real*8 a(n), ymax

	ymax = a(1)
	do i=2,n
	 if (a(i).gt.ymax) ymax=a(i)
	end do

	return
	end