DO l =3,SGU
        DO S = l,SGU

            R_rN = qsqrt(sum(sum((u_array(s,l-1,:,:)-u_array(s,l-2,:,:))*(u_array(s,l-1,:,:)-u_array(s,l-2,:,:)),DIM = 1)));
            R_N = qsqrt(sum(sum((u_array(s-1,l-1,:,:)-u_array(s-1,l-2,:,:))*(u_array(s-1,l-1,:,:)-u_array(s-1,l-2,:,:)),DIM = 1)));
            p_total(s-2,l-2)= ((log(R_N))-(log(R_rN)))/(log(r));
        ENDDO
    ENDDO