Esope source files

   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 :