Télécharger ottvaj.eso

Retour à la liste

Numérotation des lignes :

ottvaj
  1. C OTTVAJ SOURCE FD218221 21/06/10 21:15:53 11030
  2. SUBROUTINE OTTVAJ(FIL,XVAL,VF1,TAIL,P,KFH,IERUT)
  3. *
  4. IMPLICIT INTEGER (I-N)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6. PARAMETER (XZER=0.D0,UN=1.D0,DEUX=2.D0)
  7. PARAMETER (XPI = 3.1415926535897931D0)
  8. DIMENSION FIL(*),XVAL(*),VF1(*),TAIL(*),P(*)
  9. *
  10. IERUT=0
  11. YOUN = XVAL(1)
  12. XNU = XVAL(2)
  13. XLTR = XVAL(3)
  14. GFTR = XVAL(4)
  15. WO1=TAIL(1)*VF1(1)*VF1(1)+TAIL(2)*VF1(2)*VF1(2)
  16. WO1=WO1+TAIL(3)*VF1(3)*VF1(3)+DEUX*TAIL(4)*VF1(2)*VF1(1)
  17. WO1=WO1+DEUX*TAIL(5)*VF1(1)*VF1(3)
  18. WO1=WO1+DEUX*TAIL(6)*VF1(2)*VF1(3)
  19. WO2=P(1)*VF1(1)*VF1(1)+P(2)*VF1(2)*VF1(2)+P(3)*VF1(3)*VF1(3)
  20. WO2=WO2+DEUX*P(4)*VF1(2)*VF1(1)+DEUX*P(5)*VF1(1)*VF1(3)
  21. WO2=WO2+DEUX*P(6)*VF1(2)*VF1(3)
  22. IF (WO1.NE.XZER) THEN
  23. FIL(KFH)=ABS(WO2/WO1)
  24. IF (FIL(KFH).LE
  25. & .((UN+XNU)*XLTR*XLTR/(DEUX*GFTR*YOUN))) THEN
  26. IERUT=2
  27. RETURN
  28. ENDIF
  29. ELSE
  30. FIL(KFH)=UN
  31. ENDIF
  32. RETURN
  33. END
  34.  
  35.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales