Esope source files

   1 : C FVAF1     SOURCE    CHAT      05/01/13    00:13:28     5004           
   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 = ( ( DEXP(0.183D0*(100.D0-DSIU)) ) *
  75 :      &          ( (0.2031D-13+(0.67861D-32*TFIS)) / (DGRA*DGRA) ) *
  76 :      &          ( DEXP(-45294.4D0/TUO2) ) )
  77 :      &        * FACF
  78 : C
  79 :       RETURN
  80 :       END
  81 :  
  82 :