Esope source files

   1 : C KAVOTH    SOURCE    CHAT      05/01/13    00:51:45     5004           
   2 :       SUBROUTINE KAVOTH(RI,ZI,RJ,ZJ,OBS,NOBS,NT,KVU,NAL,NM,AL,BL,KIMP)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C*********************************************************************
   6 : C SP appele par KAXK
   7 : C
   8 : C -  DETERMINATION DES LIMITES D INTEGRATION ANGULAIRE POUR LES POINTS
   9 : C    I (SEGMENT 1) ET J(SEGMENT 2) EN TESTANT SI CHAQUE OBSTRUCTEUR
  10 : C    POTENTIEL L EST EFFECTIVEMENT
  11 : C -  MISE A JOUR DES INTERVALLLES D INTEGRATION
  12 : C     entree:
  13 : C      RI,ZI   :COORDONNEES DU POINT I SUR LE SEGMENT 1
  14 : C      RJ,ZJ   :COORDONNEES DU POINT J SUR LE SEGMENT 2
  15 : C      OBS     :COORDONNES DES OBSTRUCTEURS POTENTIELS
  16 : C      NOBS    :NOMBRE D'OBSTRUCTEURS POTENTIELS
  17 : C      AL      :INTERVALLES D INTEGRATION
  18 : C      NAL     :NOMBRE D INTERVALLES
  19 : C      BL      :TABLEAU DE TRAVAIL
  20 : C     resultat:
  21 : C      KVU    : 1 SI L INTERVALLE EST NON VIDE,0 SINON
  22 : C*********************************************************************
  23 :       DIMENSION AL(2,NM),BL(2,NM),OBS(2,NT)
  24 : 
  25 : C>>   TRAITEMENT DES OBSTRUCTEURS
  26 : 
  27 :        DO 20 K=1,NOBS
  28 : 
  29 :         KK=2*K-1
  30 :         RO1=OBS(1,KK)
  31 :         ZO1=OBS(2,KK)
  32 :         RO2=OBS(1,KK+1)
  33 :         ZO2=OBS(2,KK+1)
  34 :         IF(KIMP.GE.4)WRITE(6,*) '>>OBSTACLES ',K,RO1,ZO1,RO2,ZO2
  35 : 
  36 : C KVU modifie par le traitement de l obstructeur
  37 : 
  37 :        CALL KAINTE(RI,ZI,RJ,ZJ,RO1,ZO1,RO2,ZO2,KVU,NAL,NM,AL,BL,KIMP)
  38 :        IF(KVU.EQ.0) GOTO 21
  39 : 
  39 :        IF(KIMP.GE.4) CALL UTPRIN(AL,2,NAL)
  40 : 
  41 :   20   CONTINUE
  42 :   21   CONTINUE
  43 : 
  44 :       RETURN
  45 :       END
  46 :