metric.f90


subroutine i2mex_getMetric (nt1, ns, the, psi, gtt, gtp, gpp, ier)

    Get (|grad the|^2, grad the . grad psi, |grad psi|^2)

    on the (the, psi) mesh

   

    gtt = |grad the|^2

    gtp = grad the . grad psi

    gpp = |grad psi|^2

integer, intent(in) :: nt1, ns

real(i2mex_r8), intent(in) :: the(nt1), psi(ns)

real(i2mex_r8), intent(out) :: gtt(nt1, ns), gtp(nt1, ns), gpp(nt1, ns)

integer, intent(out) :: ier


subroutine i2mex_getMetricT (nt1, ns, the, psi, gtt, gtp, gpp, ier)

    Get poloidal derivative (d /d the) of metric quantities

    (|grad the|^2, grad the . grad psi, |grad psi|^2)

    on the (the, psi) mesh

   

    gtt = (d /d the) |grad the|^2

    gtp = (d /d the) grad the . grad psi

    gpp = (d /d the) |grad psi|^2

integer, intent(in) :: nt1, ns

real(i2mex_r8), intent(in) :: the(nt1), psi(ns)

real(i2mex_r8), intent(out) :: gtt(nt1, ns), gtp(nt1, ns), gpp(nt1, ns)

integer, intent(out) :: ier


subroutine i2mex_getMetricP (nt1, ns, the, psi, gtt, gtp, gpp, ier)

    Get flux derivative (d /d psi) of metric quantities

    (|grad the|^2, grad the . grad psi, |grad psi|^2)

    on the (the, psi) mesh

   

    gtt = (d /d the) |grad the|^2

    gtp = (d /d the) grad the . grad psi

    gpp = (d /d the) |grad psi|^2

integer, intent(in) :: nt1, ns

real(i2mex_r8), intent(in) :: the(nt1), psi(ns)

real(i2mex_r8), intent(out) :: gtt(nt1, ns), gtp(nt1, ns), gpp(nt1, ns)

integer, intent(out) :: ier


Send comments about this document to Alexander Ppletzer. Thu Jan 26 12:43:46 2006