Télécharger clcisa.eso

Retour à la liste

Numérotation des lignes :

  1. C CLCISA SOURCE CHAT 05/01/12 22:04:30 5004
  2. SUBROUTINE CLCISA (ET,YUNG,XNU,BETR,XLTR,CISA)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. DIMENSION ET(*),CISA(*),XLTR(*),ES(3)
  7. C
  8. C***********************************************************************
  9. C**************************** INITIALISATION ***************************
  10. C***********************************************************************
  11. C
  12. RFE=1.D-8
  13. RFS=YUNG*RFE
  14. C
  15. ES(1)=ET(1)
  16. ES(2)=ET(2)
  17. ES(3)=ET(3)
  18. DUNPNU=2.D0*(1.D0+XNU)
  19. CISARS=YUNG*BETR/DUNPNU
  20. C
  21. XLTRI=XLTR(1)+XLTR(2)+XLTR(3)
  22. C
  23. C***********************************************************************
  24. C*********************** DIFFERENTS CAS PREVUS *************************
  25. C***********************************************************************
  26. C
  27. IF(XLTRI.GT.0.D0) THEN
  28. C
  29. IF(ES(1).LT.0.D0.AND.ABS(ES(1)).GT.RFS) ES(1)=YUNG
  30. IF(ES(2).LT.0.D0.AND.ABS(ES(2)).GT.RFS) ES(2)=YUNG
  31. IF(ES(3).LT.0.D0.AND.ABS(ES(3)).GT.RFS) ES(3)=YUNG
  32. C
  33. CISA(1)=MIN(ES(1),ES(2))/DUNPNU
  34. CISA(2)=MIN(ES(2),ES(3))/DUNPNU
  35. CISA(3)=MIN(ES(3),ES(1))/DUNPNU
  36. C
  37. IF(CISA(1).LT.CISARS) CISA(1)=CISARS
  38. IF(CISA(2).LT.CISARS) CISA(2)=CISARS
  39. IF(CISA(3).LT.CISARS) CISA(3)=CISARS
  40. C
  41. ELSE
  42. C
  43. CISA(1)=CISARS
  44. CISA(2)=CISARS
  45. CISA(3)=CISARS
  46. C
  47. ENDIF
  48. C
  49. RETURN
  50. END
  51.  
  52.  

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