[Users] GlobalDerivative in spherical (multipatch) coordinates

Dear all,

I want to solve a spherical problem within a multi patch environment, using Llama and the Thornburg04 coordinates.

For my equations I need several 1st and 2nd spatial derivatives in r direction. The GlobalDerivative thorn (part of Llama) provides a subroutine, that calculates the derivatives for certain directions:

SUBROUTINE globalDiff_gv ( CCTK_POINTER_TO_CONST IN cctkGH, \                            CCTK_INT IN dir, \                            CCTK_REAL IN ARRAY var, \                            CCTK_REAL OUT ARRAY dvar, \                            CCTK_REAL IN ARRAY J_dadx, \                            CCTK_REAL IN ARRAY J_dbdx, \                            CCTK_REAL IN ARRAY J_dcdx, \                            CCTK_REAL IN ARRAY J_dady, \                            CCTK_REAL IN ARRAY J_dbdy, \                            CCTK_REAL IN ARRAY J_dcdy, \                            CCTK_REAL IN ARRAY J_dadz, \                            CCTK_REAL IN ARRAY J_dbdz, \                            CCTK_REAL IN ARRAY J_dcdz, \                            CCTK_INT IN table_handle )

However, as far as I could see in the src code this only works for Cartesian coordinates or did I miss something?

Does anyone know if there is a nice way to calculate such derivatives in r direction?

Thanks a lot and best regards,

Severin Frank