[Users] GlobalDerivative in spherical (multipatch) coordinates

[Severin Frank]

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