Esope source files

   1 : C FVAF3     SOURCE    CHAT      05/01/13    00:13:35     5004           
   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 = ( ( DEXP(0.241D0*(100.D0-DSIU)) ) *
  70 :      &          ( 0.1524D-25+(0.50825D-44*TFIS) ) *
  71 :      &          ( DEXP(-66431.8D0/TUO2) ) )
  72 :      &        * FACF
  73 : C
  74 :       RETURN
  75 :       END
  76 :  
  77 :