evaluate Jacobian of mapping h: J_h=sqrt(det(G)) at q=(X^1,X^2,zeta)
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(t_hmap_RZ), | intent(in) | :: | sf | |||
| real(kind=wp), | intent(in) | :: | q_in(3) |
FUNCTION hmap_RZ_eval_Jh( sf ,q_in) RESULT(Jh) IMPLICIT NONE !----------------------------------------------------------------------------------------------------------------------------------- ! INPUT VARIABLES CLASS(t_hmap_RZ), INTENT(IN) :: sf REAL(wp) , INTENT(IN) :: q_in(3) !----------------------------------------------------------------------------------------------------------------------------------- ! OUTPUT VARIABLES REAL(wp) :: Jh !=================================================================================================================================== ! q= (R,Z,zeta) Jh=q_in(1) END FUNCTION hmap_RZ_eval_Jh