Télécharger la plaquette Cast3M


fr en

Télécharger fvaf1.eso

Retour à la liste

Numérotation des lignes : 

fvaf1

  1. C FVAF1     SOURCE    CB215821  16/04/21    21:16:56     8920
  2.       FUNCTION FVAF1(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 'AF1 ' 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=261 : valeur illicite de 'TFIS'
  21. C     IERUT=262 : valeur illicite de 'TUO2'
  22. C     IERUT=263 : valeur illicite de 'FACF'
  23. C     IERUT=264 : valeur illicite de 'DSIU'
  24. C     IERUT=265 : valeur illicite de 'DGRA'
  25. C-----------------------------------------------------------------------
  26. C     Type de la fonction
  27. C
  28.       REAL*8       FVAF1
  29. C
  30. C     Arguments de l'interface
  31. C
  32.       INTEGER      NPARA, IERUT
  33.       REAL*8       VALPAR(NPARA)
  34. C
  35. C     Variables locales
  36. C
  37.       REAL*8       T, TFIS, TUO2, FACF, DSIU, DGRA
  38. C
  39. C---------------------- Debut du code executable -----------------------
  40. C
  41. C     Identification des parametres : affectation de variables locales
  42. C
  43.       TFIS=VALPAR(1)
  44.       TUO2=VALPAR(2)
  45.       FACF=VALPAR(3)
  46.       DSIU=VALPAR(4)
  47.       DGRA=VALPAR(5)
  48. C
  49. C     Tests d'erreur sur les parametres
  50. C
  51.       IERUT=0
  52.       FVAF1=0.D0
  53. C
  54.       IF (TFIS.LT.0.D0) THEN
  55.          IERUT = 261
  56.          RETURN
  57.       ELSE IF (TUO2.LT.0.0D0) THEN
  58. C........Temperature en Kelvin
  59.          IERUT = 262
  60.          RETURN
  61.       ELSE IF (FACF.LT.0.D0) THEN
  62.          IERUT = 263
  63.          RETURN
  64.       ELSE IF ((DSIU.LT.0.D0).OR.(DSIU.GT.100.D0)) THEN
  65.          IERUT = 264
  66.          RETURN
  67.       ELSE IF (DGRA.LT.0.D0) THEN
  68.          IERUT = 265
  69.          RETURN
  70.       ENDIF
  71. C
  72. C     Calcul de la fonction
  73. C
  74.       FVAF1 = ( ( EXP(0.183D0*(100.D0-DSIU)) ) *
  75.      &          ( (0.2031D-13+(0.67861D-32*TFIS)) / (DGRA*DGRA) ) *
  76.      &          ( EXP(-45294.4D0/TUO2) ) )
  77.      &        * FACF
  78. C
  79.       RETURN
  80.       END
  81.  
  82.  
  83.  
  84.  
  85.  

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