39 integer*4 :: ie, ne, i, j, ind_nd, id, ic
41 integer*4,
intent(in) :: nn_loc
43 integer*4,
dimension(0:nnz_loc),
intent(in) :: cs_loc
44 integer*4,
dimension(nn_loc),
intent(inout) :: loc_n_num
45 integer*4,
dimension(0:nnz_loc) :: cs_loc_copy
50 ne = cs_loc_copy(0) - 1
54 cs_loc_copy(cs_loc_copy(ie-1)) = 0
59 do ie = 0, cs_loc_copy(0)
68 if (cs_loc_copy(i) .ne. 0)
then
69 if(cs_loc_copy(i) .ne. cs_loc_copy(i-1))
then
71 loc_n_num(ic) = cs_loc_copy(i)