Télécharger fvaf1.eso

Retour à la liste

Numérotation des lignes :

  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