hmap_cyl_eval_all_e Subroutine

public pure subroutine hmap_cyl_eval_all_e(cyl_len, 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) :: cyl_len

input parameter
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_cyl_eval_all_e~~CalledByGraph proc~hmap_cyl_eval_all_e hmap_cyl_eval_all_e proc~hmap_cyl_eval_all t_hmap_cyl%hmap_cyl_eval_all proc~hmap_cyl_eval_all->proc~hmap_cyl_eval_all_e proc~hmap_cyl_test hmap_cyl_test proc~hmap_cyl_test->proc~hmap_cyl_eval_all 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

PURE SUBROUTINE hmap_cyl_eval_all_e(cyl_len,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  )
! MODULES
  IMPLICIT NONE
!-----------------------------------------------------------------------------------------------------------------------------------
! INPUT VARIABLES
  REAL(wp),INTENT(IN)  :: cyl_len  !! input parameter
  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}
!-----------------------------------------------------------------------------------------------------------------------------------
! LOCAL VARIABLES
  !REAL(wp) :: Gh21,Gh31,Gh32,Gh33
!===================================================================================================================================

  Gh11=1.0_wp
  !Gh21=0.0_wp
  Gh22=1.0_wp
  !Gh31=0.0_wp
  !Gh32=0.0_wp
  !Gh33=(sf%cyl_len**2)

  Jh=cyl_len
  Jh_dq1=0.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  + (cyl_len**2)

  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 = 0.0_wp
  g_zz_dq2 = 0.0_wp
END SUBROUTINE hmap_cyl_eval_all_e