Télécharger fvaf3.eso

Retour à la liste

Numérotation des lignes :

  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