hmap_cyl_get_dx_dqi Subroutine

public subroutine hmap_cyl_get_dx_dqi(sf, q_in, dx_dq1, dx_dq2, dx_dq3)

evaluate all first derivatives dx(1:3)/dq^i, i=1,2,3 , at q_in=(X^1,X^2,zeta),

Type Bound

t_hmap_cyl

Arguments

Type IntentOptional Attributes Name
class(t_hmap_cyl), intent(in) :: sf
real(kind=wp), intent(in) :: q_in(3)
real(kind=wp), intent(out) :: dx_dq1(3)
real(kind=wp), intent(out) :: dx_dq2(3)
real(kind=wp), intent(out) :: dx_dq3(3)

Source Code

SUBROUTINE hmap_cyl_get_dx_dqi( sf ,q_in,dx_dq1,dx_dq2,dx_dq3)
  IMPLICIT NONE
  !-----------------------------------------------------------------------------------------------------------------------------------
  ! INPUT VARIABLES
  CLASS(t_hmap_cyl), INTENT(IN) :: sf
  REAL(wp)         , INTENT(IN) :: q_in(3)
  !-----------------------------------------------------------------------------------------------------------------------------------
  ! OUTPUT VARIABLES
  REAL(wp)     , INTENT(OUT) :: dx_dq1(3)
  REAL(wp)     , INTENT(OUT) :: dx_dq2(3)
  REAL(wp)     , INTENT(OUT) :: dx_dq3(3)
  !===================================================================================================================================
  dx_dq1(1:3) = (/ 1.0_wp,      0.0_wp, 0.0_wp /)
  dx_dq2(1:3) = (/ 0.0_wp,      0.0_wp, 1.0_wp /)
  dx_dq3(1:3) = (/ 0.0_wp, -sf%cyl_len, 0.0_wp /)
END SUBROUTINE hmap_cyl_get_dx_dqi