hto_acmplx_pro Module Reference

Data Types
interface	operator(.cp.)

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

Function/Subroutine Documentation

cp()

real*8 function, dimension(2) hto_acmplx_pro::cp	(	real*8, dimension(2), intent(in)	x,
		real*8, dimension(2), intent(in)	y
	)

Definition at line 263 of file CALLING_cpHTO.f.

       USE hto_riemann
       USE hto_units
       IMPLICIT NONE
       real*8, dimension(2) :: cp,z1,z2
       real*8, INTENT(IN), dimension(2) :: x,y
       IF(abs(x(2)).eq.1.d0.and.abs(y(2)).eq.1.d0) THEN
        cp(1)= x(1)*y(1)
        cp(2)= 0.d0
       ELSE
        IF(abs(x(2)).eq.1.d0) THEN
         z1(1)= x(1)
         z1(2)= 0.d0
        ELSE
         z1= x
        ENDIF
        IF(abs(y(2)).eq.1.d0) THEN
         z2(1)= y(1)
         z2(2)= 0.d0
        ELSE
         z2= y        
        ENDIF
        cp(1)= z1(1)*z2(1)-z1(2)*z2(2)
        cp(2)= z1(1)*z2(2)+z1(2)*z2(1)
       ENDIF