c	subroutina para calcular el maximo de los valores de una variable
c	en double precision
	subroutine maxdp_2(a,n,ymax,imax)
c
c       jul 2011 malv+fgg

	implicit none
	integer n, i, imax
	real*8  a(n), y, ymax

	imax=0
	y=-1.7d+38
	do i=1,n
	     if(a(i).gt.y) then
		y=a(i)
		imax=i
	     end if
	end do
	ymax=y
!	type *,'imax,ymax',imax,ymax
	return
	end