Are vector spherical harmonic transforms available?

I have a beginner question: I would like to compute some vector spherical harmonic transforms, specifically I want to be able to compute vorticity and divergence from eastward and northward vector components, and also the reverse. I am new to Julia, but in Python I have been using shtns and there I would use their 2-argument forms of SH analysis and synthesis with some scaling (just like the getuv and getvrtdiv methods in their shallow_water.py example).

I am having trouble understanding how (or if) this is possible in FastTransforms.jl, can you offer any advice? I've seen the plan_sphv* functions but I am not well enough versed in the topic to understand what these do (i.e. they operate on a single field, not a pair). Thanks!

Start by reading the plan_sphv* functions in FastTransforms.jl and compare their single-field operations with the getuv and getvrtdiv methods in the linked shallow_water.py example. Determine whether the requested forward and reverse vector transforms are supported; done means providing a clear mapping or documenting that the capability is unavailable.