UpdateRangeSegmentT Subroutine

public subroutine UpdateRangeSegmentT(x, t)

Arguments

Type IntentOptional Attributes Name
real(kind=8), intent(in) :: x(2)
real(kind=8), intent(in) :: t(2)

Called by

proc~~updaterangesegmentt~~CalledByGraph proc~updaterangesegmentt UpdateRangeSegmentT proc~quad Quad proc~quad->proc~updaterangesegmentt proc~evaluatessp EvaluateSSP proc~evaluatessp->proc~quad proc~bellhopcore BellhopCore proc~bellhopcore->proc~evaluatessp proc~influencecervenycart InfluenceCervenyCart proc~bellhopcore->proc~influencecervenycart proc~traceray2d TraceRay2D proc~bellhopcore->proc~traceray2d proc~influencecervenycart->proc~evaluatessp proc~readenvironment ReadEnvironment proc~readenvironment->proc~evaluatessp proc~reflect2d~2 Reflect2D proc~reflect2d~2->proc~evaluatessp proc~step2d~2 Step2D proc~step2d~2->proc~evaluatessp proc~traceray2d->proc~evaluatessp proc~traceray2d->proc~reflect2d~2 proc~traceray2d->proc~step2d~2 proc~bellhopcore~2 BellhopCore proc~bellhopcore~2->proc~influencecervenycart program~bellhop BELLHOP program~bellhop->proc~bellhopcore program~bellhop->proc~readenvironment program~bellhop3d BELLHOP3D program~bellhop3d->proc~readenvironment program~bellhop3d->proc~bellhopcore~2

Source Code

  SUBROUTINE UpdateRangeSegmentT( x, t )
     REAL (KIND=8), INTENT(IN) :: x( 2 ), t( 2 )

     ! LP: See UpdateDepthSegmentT

     IF ( t( 1 ) >= 0.0 ) THEN
        ! SSP%Seg%r( iSegr ) <= x( 1 ) < SSP%Seg%r( iSegr + 1 )
        DO WHILE ( x( 1 ) < SSP%Seg%r( iSegr ) .AND. iSegr > 1 )
           iSegr = iSegr - 1
        END DO
        DO WHILE ( x( 1 ) >= SSP%Seg%r( iSegr + 1 ) .AND. iSegr < SSP%Nr-1 )
           iSegr = iSegr + 1
        END DO
     ELSE
        ! SSP%Seg%r( iSegr ) < x( 1 ) <= SSP%Seg%r( iSegr + 1 )
        DO WHILE ( x( 1 ) > SSP%Seg%r( iSegr + 1 ) .AND. iSegr < SSP%Nr-1 )
           iSegr = iSegr + 1
        END DO
        DO WHILE ( x( 1 ) <= SSP%Seg%r( iSegr ) .AND. iSegr > 1 )
           iSegr = iSegr - 1
        END DO
     END IF

  END SUBROUTINE UpdateRangeSegmentT