TROENKz_invL_from_phi Subroutine

public subroutine TROENKz_invL_from_phi(z, phih, invL)

Arguments

Type IntentOptional Attributes Name
real, intent(in) :: z
real, intent(in) :: phih
real, intent(out) :: invL

Called by

proc~~troenkz_invl_from_phi~~CalledByGraph proc~troenkz_invl_from_phi TROENKz_invL_from_phi proc~uemep_read_emep uEMEP_read_EMEP proc~uemep_read_emep->proc~troenkz_invl_from_phi program~uemep uEMEP program~uemep->proc~uemep_read_emep

Source Code

    subroutine TROENKz_invL_from_phi(z,phih,invL)

        implicit none

        real, intent(in) :: z    ! height
        real, intent(in) :: phih
        real, intent(out) :: invL !1/L

        if ( phih < 1 ) then
            !phih=(1-16.*min(z,zsurf)*invL)**(-1./2.) !As in Garratt and Obrien for phih, so with Prandtl number
            invL=(phih**-2.-1)/(-16.*z)
        else
            !phih=1+5.*z*invL !As in Garratt, Prandtl number is 1 in stable boundary layer
            invL=(phih-1)/5./z
        endif


    end subroutine TROENKz_invL_from_phi