hmap_RZ_eval_gij Function

public function hmap_RZ_eval_gij(sf, qL_in, q_G, qR_in) result(g_ab)

evaluate sum_ij (qL_i (G_ij(q_G)) qR_j) ,, where qL=(dX^1/dalpha,dX^2/dalpha ,dzeta/dalpha) and qR=(dX^1/dbeta,dX^2/dbeta ,dzeta/dbeta) and dzeta_dalpha then known to be either 0.0 for ds and dtheta and 1.0 for dzeta

Type Bound

t_hmap_RZ

Arguments

Type IntentOptional Attributes Name
class(t_hmap_RZ), intent(in) :: sf
real(kind=wp), intent(in) :: qL_in(3)
real(kind=wp), intent(in) :: q_G(3)
real(kind=wp), intent(in) :: qR_in(3)

Return Value real(kind=wp)


Called by

proc~~hmap_rz_eval_gij~~CalledByGraph proc~hmap_rz_eval_gij t_hmap_RZ%hmap_RZ_eval_gij proc~hmap_rz_test hmap_RZ_test proc~hmap_rz_test->proc~hmap_rz_eval_gij proc~hmap_rz_init hmap_RZ_init proc~hmap_rz_init->proc~hmap_rz_test interface~t_hmap_rz t_hmap_RZ interface~t_hmap_rz->proc~hmap_rz_init

Source Code

FUNCTION hmap_RZ_eval_gij( sf ,qL_in,q_G,qR_in) RESULT(g_ab)
  IMPLICIT NONE
  !-----------------------------------------------------------------------------------------------------------------------------------
  ! INPUT VARIABLES
  CLASS(t_hmap_RZ), INTENT(IN) :: sf
  REAL(wp)        , INTENT(IN) :: qL_in(3)
  REAL(wp)        , INTENT(IN) :: q_G(3)
  REAL(wp)        , INTENT(IN) :: qR_in(3)
  !-----------------------------------------------------------------------------------------------------------------------------------
  ! OUTPUT VARIABLES
  REAL(wp)                     :: g_ab
  !===================================================================================================================================
  !                       |R   |   |1  0  0   |        |R   |
  !q_i G_ij q_j = (dalpha |Z   | ) |0  1  0   | (dbeta |Z   | )
  !                       |zeta|   |0  0  R^2 |        |zeta|
  g_ab=SUM(qL_in(:)*(/qR_in(1),qR_in(2),q_G(1)**2*qR_in(3)/))
END FUNCTION hmap_RZ_eval_gij