Télécharger la plaquette Cast3M


fr en

Télécharger fvaf3.eso

Retour à la liste

Numérotation des lignes : 

fvaf3

  1. C FVAF3     SOURCE    CB215821  16/04/21    21:16:57     8920
  2.       FUNCTION FVAF3(VALPAR,NPARA,IERUT)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5. C-----------------------------------------------------------------------
  6. C     DESCRIPTION FONCTIONNELLE :
  7. C     -------------------------
  8. C     Fonction externe d'evaluation de la composante 'AF3 ' du modele
  9. C     'FLUAGE' 'POLYNOMIAL' externalise pour TOUTATIS
  10. C
  11. C     INTERFACE :
  12. C     ---------
  13. C     ENTREE : VALPAR : REAL*8(NPARA), valeurs des parametres
  14. C     ENTREE : NPARA  : INTEGER, nombre de parametres
  15. C     SORTIE : IERUT  : INTEGER, code retour
  16. C
  17. C     CODIFICATION DES ERREURS :
  18. C     ------------------------
  19. C     IERUT=0   : OK
  20. C     IERUT=281 : valeur illicite de 'TFIS'
  21. C     IERUT=282 : valeur illicite de 'TUO2'
  22. C     IERUT=283 : valeur illicite de 'FACF'
  23. C     IERUT=284 : valeur illicite de 'DSIU'
  24. C-----------------------------------------------------------------------
  25. C     Type de la fonction
  26. C
  27.       REAL*8       FVAF3
  28. C
  29. C     Arguments de l'interface
  30. C
  31.       INTEGER      NPARA, IERUT
  32.       REAL*8       VALPAR(NPARA)
  33. C
  34. C     Variables locales
  35. C
  36.       REAL*8       T, TFIS, TUO2, FACF, DSIU
  37. C
  38. C---------------------- Debut du code executable -----------------------
  39. C
  40. C     Identification des parametres : affectation de variables locales
  41. C
  42.       TFIS=VALPAR(1)
  43.       TUO2=VALPAR(2)
  44.       FACF=VALPAR(3)
  45.       DSIU=VALPAR(4)
  46. C
  47. C     Tests d'erreur sur les parametres
  48. C
  49.       IERUT=0
  50.       FVAF3=0.D0
  51. C
  52.       IF (TFIS.LT.0.D0) THEN
  53.          IERUT = 281
  54.          RETURN
  55.       ELSE IF (TUO2.LT.0.0D0) THEN
  56. C........Temperature en Kelvin
  57.          IERUT = 282
  58.          RETURN
  59.       ELSE IF (FACF.LT.0.D0) THEN
  60.          IERUT = 283
  61.          RETURN
  62.       ELSE IF ((DSIU.LT.0.D0).OR.(DSIU.GT.100.D0)) THEN
  63.          IERUT = 284
  64.          RETURN
  65.       ENDIF
  66. C
  67. C     Calcul de la fonction
  68. C
  69.       FVAF3 = ( ( EXP(0.241D0*(100.D0-DSIU)) ) *
  70.      &          ( 0.1524D-25+(0.50825D-44*TFIS) ) *
  71.      &          ( EXP(-66431.8D0/TUO2) ) )
  72.      &        * FACF
  73. C
  74.       RETURN
  75.       END
  76.  
  77.  
  78.  
  79.  
  80.  

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