hmap_RZ_eval_all_e Subroutine

public pure subroutine hmap_RZ_eval_all_e(q1, q2, dX1_dt, dX2_dt, dX1_dz, dX2_dz, Jh, g_tt, g_tz, g_zz, Jh_dq1, g_tt_dq1, g_tz_dq1, g_zz_dq1, Jh_dq2, g_tt_dq2, g_tz_dq2, g_zz_dq2, g_t1, g_t2, g_z1, g_z2, Gh11, Gh22)

evaluate all quantities at one given point (elemental)

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: q1

solution variables q1,q2
real(kind=wp), intent(in) :: q2

solution variables q1,q2
real(kind=wp), intent(in) :: dX1_dt

theta derivative of solution variables q1,q2
real(kind=wp), intent(in) :: dX2_dt

theta derivative of solution variables q1,q2
real(kind=wp), intent(in) :: dX1_dz

zeta derivative of solution variables q1,q2
real(kind=wp), intent(in) :: dX2_dz

zeta derivative of solution variables q1,q2
real(kind=wp), intent(out) :: Jh

Jac,1/Jac,g_{ab} with a=theta/zeta b=theta/zeta
real(kind=wp), intent(out) :: g_tt

Jac,1/Jac,g_{ab} with a=theta/zeta b=theta/zeta
real(kind=wp), intent(out) :: g_tz

Jac,1/Jac,g_{ab} with a=theta/zeta b=theta/zeta
real(kind=wp), intent(out) :: g_zz

Jac,1/Jac,g_{ab} with a=theta/zeta b=theta/zeta
real(kind=wp), intent(out) :: Jh_dq1

and their variation vs q1
real(kind=wp), intent(out) :: g_tt_dq1

and their variation vs q1
real(kind=wp), intent(out) :: g_tz_dq1

and their variation vs q1
real(kind=wp), intent(out) :: g_zz_dq1

and their variation vs q1
real(kind=wp), intent(out) :: Jh_dq2

and their variation vs q2
real(kind=wp), intent(out) :: g_tt_dq2

and their variation vs q2
real(kind=wp), intent(out) :: g_tz_dq2

and their variation vs q2
real(kind=wp), intent(out) :: g_zz_dq2

and their variation vs q2
real(kind=wp), intent(out) :: g_t1

dq^{i}/dthetaG^{i1}, dq^{i}/dthetaG^{i2}, and dq^{i}/dzetaG^{i1}, dq^{i}/dzetaG^{i2} and G^{11},G^{22}
real(kind=wp), intent(out) :: g_t2

dq^{i}/dthetaG^{i1}, dq^{i}/dthetaG^{i2}, and dq^{i}/dzetaG^{i1}, dq^{i}/dzetaG^{i2} and G^{11},G^{22}
real(kind=wp), intent(out) :: g_z1

dq^{i}/dthetaG^{i1}, dq^{i}/dthetaG^{i2}, and dq^{i}/dzetaG^{i1}, dq^{i}/dzetaG^{i2} and G^{11},G^{22}
real(kind=wp), intent(out) :: g_z2

dq^{i}/dthetaG^{i1}, dq^{i}/dthetaG^{i2}, and dq^{i}/dzetaG^{i1}, dq^{i}/dzetaG^{i2} and G^{11},G^{22}
real(kind=wp), intent(out) :: Gh11

dq^{i}/dthetaG^{i1}, dq^{i}/dthetaG^{i2}, and dq^{i}/dzetaG^{i1}, dq^{i}/dzetaG^{i2} and G^{11},G^{22}
real(kind=wp), intent(out) :: Gh22

dq^{i}/dthetaG^{i1}, dq^{i}/dthetaG^{i2}, and dq^{i}/dzetaG^{i1}, dq^{i}/dzetaG^{i2} and G^{11},G^{22}

Called by

proc~~hmap_rz_eval_all_e~~CalledByGraph proc~hmap_rz_eval_all_e hmap_RZ_eval_all_e proc~hmap_rz_eval_all t_hmap_RZ%hmap_RZ_eval_all proc~hmap_rz_eval_all->proc~hmap_rz_eval_all_e proc~hmap_rz_test hmap_RZ_test proc~hmap_rz_test->proc~hmap_rz_eval_all 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

PURE SUBROUTINE hmap_RZ_eval_all_e(q1,q2,dX1_dt,dX2_dt,dX1_dz,dX2_dz, &
                                   Jh,    g_tt,    g_tz,    g_zz,     &
                                   Jh_dq1,g_tt_dq1,g_tz_dq1,g_zz_dq1, &
                                   Jh_dq2,g_tt_dq2,g_tz_dq2,g_zz_dq2, &
                                   g_t1,g_t2,g_z1,g_z2,Gh11,Gh22  )
  IMPLICIT NONE
  !-----------------------------------------------------------------------------------------------------------------------------------
  ! INPUT VARIABLES
  REAL(wp),INTENT(IN)  :: q1,q2       !! solution variables q1,q2
  REAL(wp),INTENT(IN)  :: dX1_dt,dX2_dt  !! theta derivative of solution variables q1,q2
  REAL(wp),INTENT(IN)  :: dX1_dz,dX2_dz  !!  zeta derivative of solution variables q1,q2
  !-----------------------------------------------------------------------------------------------------------------------------------
  ! OUTPUT VARIABLES
  REAL(wp),INTENT(OUT) :: Jh,g_tt,g_tz,g_zz              !! Jac,1/Jac,g_{ab} with a=theta/zeta b=theta/zeta
  REAL(wp),INTENT(OUT) :: Jh_dq1,g_tt_dq1,g_tz_dq1,g_zz_dq1  !! and their variation vs q1
  REAL(wp),INTENT(OUT) :: Jh_dq2,g_tt_dq2,g_tz_dq2,g_zz_dq2  !! and their variation vs q2
  REAL(wp),INTENT(OUT) :: g_t1,g_t2,g_z1,g_z2,Gh11,Gh22  !! dq^{i}/dtheta*G^{i1}, dq^{i}/dtheta*G^{i2}, and dq^{i}/dzeta*G^{i1}, dq^{i}/dzeta*G^{i2} and G^{11},G^{22}
!===================================================================================================================================

  Gh11=1.0_wp
  !Gh21=0.0_wp
  Gh22=1.0_wp
  !Gh31=0.0_wp
  !Gh32=0.0_wp
  !Gh33=q1**2

  Jh=q1
  Jh_dq1=1.0_wp
  Jh_dq2=0.0_wp

  g_t1 = dX1_dt
  g_t2 = dX2_dt
  g_z1 = dX1_dz
  g_z2 = dX2_dz

  g_tt =   dX1_dt *  g_t1  +  dX2_dt *  g_t2
  g_tz =   dX1_dt *  g_z1  +  dX2_dt *  g_z2
  g_zz =   dX1_dz *  g_z1  +  dX2_dz *  g_z2  + q1*q1

  g_tt_dq1 = 0.0_wp
  g_tt_dq2 = 0.0_wp

  g_tz_dq1 = 0.0_wp
  g_tz_dq2 = 0.0_wp

  g_zz_dq1 = 2.0_wp*q1
  g_zz_dq2 = 0.0_wp
END SUBROUTINE hmap_RZ_eval_all_e