emkep Function

public pure function emkep(e1, ee)

Similar to emkepl, except input is 1-e.

Arguments

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

Return Value real(kind=wp)


Called by

proc~~emkep~~CalledByGraph proc~emkep gooding_module::emkep proc~ekepl gooding_module::ekepl proc~ekepl->proc~emkep proc~pv2els gooding_module::pv2els proc~pv2els->proc~emkep proc~els2pv gooding_module::els2pv proc~els2pv->proc~ekepl proc~pv3els gooding_module::pv3els proc~pv3els->proc~pv2els proc~els3pv gooding_module::els3pv proc~els3pv->proc~els2pv proc~lambert_test lambert_module::lambert_test proc~lambert_test->proc~pv3els proc~propagate gooding_module::propagate proc~propagate->proc~pv3els proc~propagate->proc~els3pv

Source Code

    pure function emkep(e1,ee)

    implicit none

    real(wp)            :: emkep
    real(wp),intent(in) :: e1
    real(wp),intent(in) :: ee

    real(wp) :: x, ee2, term, d, x0

    x    = e1*sin(ee)
    ee2  = -ee*ee
    term = ee
    d    = zero

    do

        d = d + two
        term = term*ee2/(d*(d + one))
        x0 = x
        x = x - term
        if (x==x0) exit

    end do

    emkep = x

    end function emkep