anchor.diff:        type (hecmwST_matrix) :: hecMAT
         integer(kind=kint)    :: i, j, k, l, countCAL, countCAU
   
  -write(*,*) "hecmw_cmat_LU in"
         allocate (hecMAT%indexCL(0:hecMAT%NP), hecMAT%indexCU(0:hecMAT%NP))
   
         hecMAT%indexCL = 0
  
positive.code:           endif
          enddo
        enddo
  write(*,*) "hecmw_cmat_LU 2 end"
    end subroutine hecmw_cmat_LU
  
    !< free LU-part of cmat
  
negative.code:         hecMAT%indexCL(i) = hecMAT%indexCL(i) + hecMAT%indexCL(i-1)
          hecMAT%indexCU(i) = hecMAT%indexCU(i) + hecMAT%indexCU(i-1)
        enddo
  
        hecMAT%NPCL = hecMAT%indexCL(hecMAT%NP)
        hecMAT%NPCU = hecMAT%indexCU(hecMAT%NP)
  