hmap_cyl_eval_gij Function

public function hmap_cyl_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_cyl

Arguments

Type IntentOptional Attributes Name
class(t_hmap_cyl), 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_cyl_eval_gij~~CalledByGraph proc~hmap_cyl_eval_gij t_hmap_cyl%hmap_cyl_eval_gij proc~hmap_cyl_test hmap_cyl_test proc~hmap_cyl_test->proc~hmap_cyl_eval_gij proc~hmap_cyl_init_params hmap_cyl_init_params proc~hmap_cyl_init_params->proc~hmap_cyl_test interface~t_hmap_cyl t_hmap_cyl interface~t_hmap_cyl->proc~hmap_cyl_init_params proc~hmap_cyl_init hmap_cyl_init interface~t_hmap_cyl->proc~hmap_cyl_init proc~hmap_cyl_init->proc~hmap_cyl_init_params

Source Code

FUNCTION hmap_cyl_eval_gij( sf ,qL_in,q_G,qR_in) RESULT(g_ab)
!-----------------------------------------------------------------------------------------------------------------------------------
! INPUT VARIABLES
  CLASS(t_hmap_cyl), 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
!===================================================================================================================================
  g_ab=SUM(qL_in(:)*(/qR_in(1),qR_in(2),(sf%cyl_len**2)*qR_in(3)/))
END FUNCTION hmap_cyl_eval_gij