Télécharger fvalph.eso

Retour à la liste

Numérotation des lignes :

  1. C FVALPH SOURCE CHAT 07/04/20 21:15:13 5723
  2. FUNCTION FVALPH(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 'ALPH' 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=231 : temp. inferieure au zero absolu
  21. C IERUT=232 : valeur illicite de 'ALPC'
  22. C-----------------------------------------------------------------------
  23. C Type de la fonction
  24. C
  25. REAL*8 FVALPH
  26. C
  27. C Arguments de l'interface
  28. C
  29. INTEGER NPARA, IERUT
  30. REAL*8 VALPAR(NPARA)
  31. C
  32. C Variables locales
  33. C
  34. REAL*8 T
  35. INTEGER IALPC
  36. C
  37. C---------------------- Debut du code executable -----------------------
  38. C
  39. C Identification des parametres : affectation de variables locales
  40. C
  41. T=VALPAR(1)
  42. IF(VALPAR(2).GE.0.99999 . and. VALPAR(2).le.1.00001) IALPC=1
  43. IF(VALPAR(2).GE.1.99999 . and. VALPAR(2).le.2.00001) IALPC=2
  44.  
  45.  
  46. C
  47. C Tests d'erreur sur les parametres
  48. C
  49. IERUT=0
  50. FVALPH=0.D0
  51. C
  52. IF (T.LT.0.0D0) THEN
  53. C........Temperature en Kelvin
  54. IERUT = 231
  55. RETURN
  56. ELSE IF ((IALPC.LT.1).OR.(IALPC.GT.2)) THEN
  57. IERUT = 232
  58. RETURN
  59. ENDIF
  60. C
  61. C Calcul de la fonction
  62. C
  63. C.....DG MARTIN
  64. C
  65. IF (IALPC.EQ.1) THEN
  66. IF (T.LT.923.D0) THEN
  67. FVALPH = 9.828D-6 - (6.390D-10 * T) + (1.330D-12 * T * T)
  68. & - (1.757D-17 * T * T * T)
  69. ELSE
  70. FVALPH = 1.1833D-5 - (5.013D-9 * T) + (3.756D-12 * T * T)
  71. & - (6.125D-17 * T * T * T)
  72. ENDIF
  73. C
  74. C.....METEOR/TU
  75. C
  76. ELSE
  77. FVALPH = 5.193D-6 + (5.876D-9 * T) - (1.169D-13 * T * T)
  78. ENDIF
  79. C
  80. RETURN
  81. END
  82.  
  83.  
  84.  
  85.  

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