load "./dim_count_unique_values_n.ncl"

begin

;---Various arrays to test
  a   = (/ (/1,2,3,4,5/),(/2,2,2,4,5/),(/-1,-1,-1,-1,-1/)/)
  at  = transpose(a)
  am  = conform_dims((/6,3,5/),a,(/1,2/))
  atm = conform_dims((/6,5,3/),at,(/1,2/))

  answer   = (/5,3,1/)
  answerm  = conform_dims((/6,3/),answer,1)

  c   = dim_count_unique_values_n(a,1)    ; Always count the dimension with 5 elements
  cm  = dim_count_unique_values_n(am,2)
  ct  = dim_count_unique_values_n(at,0)
  ctm = dim_count_unique_values_n(atm,1)

;---These should all be True  
  print(all(c.eq.answer))
  print(all(ct.eq.answer))
  print(all(cm.eq.answerm))
  print(all(ctm.eq.answerm))
end