cccccccccccccc
c	subroutina para calcular el maximo de los valores de una variable
c	de doble precision. Simlar a maxdp salvo la adicion de limites 

	subroutine maxdp_limits ( a,n,ymax, n1,n2 )

c       jul 2011 malv
cccccccccccccc

        implicit none
	integer n, i, n1,n2
	real*8 a(n),ymax

	if ( n1 .gt. n2 ) stop 'MAXDP_LIMITS. Error: n1 > n2 '
	if ( n1 .lt. 1 ) stop 'MAXDP_LIMTIS. Error: n1 < 2 '
	if ( n2 .gt. n ) stop 'MAXDP_LIMITS. Error: n2 > n '

	ymax=a(n1)
	do i=n1+1,n2                    ! Si n1=n2 este bucle no se ejecuta
	 if (a(i).gt.ymax) ymax=a(i)
	end do

	return
	end