Esope source files

   1 : C J3ANG1    SOURCE    CHAT      05/01/13    00:45:48     5004           
   2 :       SUBROUTINE J3ANG1(XYC1,NPTO1,I1,XYC2,NPTO2,I2,TOL,AG1,AG2)
   3 : C----------------------------------------------------
   4 : C     RECHERCHE DE L'ANGLE ENTRE UN SEGMENT DE XYC1 (PTS I1-1 et I1)
   5 : C     ET DEUX SEGMENTS DE XYC2 (PTS I2-1 et I2, PTS I2+1 et I2)
   6 : C
   7 : C     PP 6/97
   8 : C     Pierre Pegon/JRC Ispra
   9 : C----------------------------------------------------
  10 : C
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)
  13 : -INC CCREEL
  14 : -INC  CCOPTIO

  15 :       DIMENSION XYC1(2,NPTO1),XYC2(2,NPTO2)
  16 :       DIMENSION XY1(2),XY2(2)
  17 :       LOGICAL LOK
  18 : C
  19 :       J1=I1-1+(1/I1)*NPTO1
  20 :       DO IE1=1,2
  21 :         XY1(IE1)=XYC1(IE1,J1)-XYC1(IE1,I1)
  22 :       ENDDO
  23 : C
  24 :       J2=I2-1+(1/I2)*NPTO2
  25 :       DO IE1=1,2
  26 :         XY2(IE1)=XYC2(IE1,J2)-XYC2(IE1,I2)
  27 :       ENDDO
  27 :       CALL J3ANGL(XY1,XY2,TOL,AG1,IRET)
  28 : C
  29 :       J2=I2+1-(I2/NPTO2)*I2
  30 :       DO IE1=1,2
  31 :         XY2(IE1)=XYC2(IE1,J2)-XYC2(IE1,I2)
  32 :       ENDDO
  32 :       CALL J3ANGL(XY1,XY2,TOL,AG2,IRET)
  33 : C
  34 :       IF(AG1.GT.XZERO)AG1=AG1-2*XPI
  35 :       IF(AG2.GT.XZERO)AG2=AG2-2*XPI
  36 : C
  37 :       RETURN
  38 :       END
  39 :  
  40 :  
  41 :