Télécharger la plaquette Cast3M


fr en

Télécharger fvalph.eso

Retour à la liste

Numérotation des lignes : 

fvalph

  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