This package is nice. But it would be more useful if it were possible to redefine u"cm^-1"() so that 1THz |> cm^-1 automatically uses unnatural(u"cm^-1", 1u"THz").
At the moment, I find it more convenient to define cm1 = uconvert(u"cm^-1", float(f/Unitful.c)) and do 1THz |> cm1 than use unnatural(u"cm^-1", 1u"THz").