C Rutina 033
C-------------------------------------------------------------ANGULO---
C
C SUBRUTINA ANGULO
C
C CALCULO DEL ANGULO FORMADO POR EL VECTOR (x,y) y el (1,0)
C RESULTADO (0,2pi)
C
SUBROUTINE ANGULO(ANG,X,Y)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C X,Y : Coordenadas (Input)
REAL*8 X,Y
C ANG : Angulo (Output)
REAL*8 ANG
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PI : Pi
REAL*8 PI
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
PI=DACOS(-1.d0)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALCULO DEL ANGULO
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
IF(X.EQ.0.AND.y.EQ.0) THEN
WRITE(F_W,03301)
WRITE (*,03301)
STOP
ELSE
ANG=DATAN2(Y,X)
IF(ANG.LT.0) THEN
ANG=2.d0*PI+ANG
END IF
ENDIF
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
03301 FORMAT(/,2x,'ERROR 03301 - ANGULO - Nodos coincidentes.')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
