hto_cmplx_srs_root Module Reference

Data Types
interface	operator(.crsrs.)

Functions/Subroutines
real *8 function, dimension(2)	crsrs (x, y)

Function/Subroutine Documentation

crsrs()

real*8 function, dimension(2) hto_cmplx_srs_root::crsrs	(	real*8, dimension(2), intent(in)	x,
		real*8, dimension(2), intent(in)	y
	)

Definition at line 624 of file CALLING_cpHTO.f.

       USE hto_units
       USE hto_riemann
       IMPLICIT NONE
       real*8 r,rs,u,v 
       real*8, dimension(2) :: crsrs
       real*8, INTENT(IN), dimension(2):: x,y
*
       IF(x(1) > 0.d0) THEN
        rs= x(1)*x(1)+x(2)*x(2)
        r= sqrt(rs)
        u= 0.5d0*(r+x(1))
        crsrs(1)= sqrt(u)
        r= sqrt(rs)
        v= 0.5d0*(r-x(1))
        crsrs(2)= sign(one,x(2))*sqrt(v)
       ELSE
        rs= x(1)*x(1)+x(2)*x(2)
        r= sqrt(rs)
        u= 0.5d0*(r+x(1))
        v= 0.5d0*(r-x(1))
*
        IF((x(2)*y(2)) > 0.d0) THEN
         crsrs(1)= sqrt(u)
         crsrs(2)= sign(one,x(2))*sqrt(v)
        ELSEIF((x(2)*y(2)) < 0.d0) THEN
         crsrs(1)= -sqrt(u)
         crsrs(2)= sign(one,x(2))*sqrt(v)
        ENDIF
*
       ENDIF