utils/pythnorm.f

c-----------------------------------------------------------------------
      subroutine pythnorm(lb, ub, shape, truncU, truncerr,
     .                    wgtu)
c-----------------------------------------------------------------------
        implicit none

        integer lb(2), ub(2), shape(2)

        real*8  truncU(shape(1), shape(2)),
     .          truncerr(shape(1), shape(2))

        real*8  wgtu

c       Local variables

        integer i, j

        do i = 1, shape(1)
          do j = 1, shape(2)
            truncerr(i,j) = sqrt( wgtu * truncU(i,j)**2 )
          end do
        end do

        return
       end
c-----------------------------------------------------------------------