Télécharger la plaquette Cast3M


fr en

Télécharger clcisa.eso

Retour à la liste

Numérotation des lignes : 

clcisa

  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