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-----------------------------------------------------------------------