Esope source files

   1 : C JACOB     SOURCE    CHAT      05/01/13    00:48:20     5004           
   2 : C***********************************************************************
   3 : C*                                                                     *
   4 : C* PROJET      : Opérateur LIMI                                        *
   5 : C* NOM         : JACOB,  sous-programme                                *
   6 : C* DESCRIPTION : Ce sous-programme calcule la matrice jacobienne et    *
   7 : C*               son inverse, pour le système de 3 équations à trois   *
   8 : C*               inconnues utilisé dans la méthode de résolution des   *
   9 : C*               relations de Cousteix                                 *
  10 : C* LANGAGE     : Esope                                                 *
  11 : C* AUTEUR      : Guillaume VENCO - DRN/DMT/SEMT/LTMF                   *
  12 : C*                                                                     *
  13 : C***********************************************************************
  14 : C*                                                                     *
  15 : C* APPELES     : Sous-programme INVMAT3.eso                            *
  16 : C*                                                                     *
  17 : C***********************************************************************
  18 : C*                                                                     *
  19 : C* ENTREES     : G   : -1ère inconnue du système, paramètre de Clauser *
  20 : C*               H   : -2ème inconnue, facteur de forme                * 
  21 : C*               ga  : -3ème inconnue, facteur de frottement gamma     *
  22 : C*                                                                     *
  23 : C* SORTIES     : W   : -Matrice jacobienne du système                  *
  24 : C*               Q   : -Matrice inverse de la matrice jacobienne       *
  25 : C*                                                                     *
  26 : C***********************************************************************
  27 : C*                                                                     *
  28 : C* VERSION     : 20/04/2000                                            *
  29 : C* CREATION    :                                                       *
  30 : C*                                                                     *
  31 : C***********************************************************************
  32 : 
  33 :        SUBROUTINE JACOB(H,G,GA,W,Q)
  34 : 
  35 : C*       *********************************
  36 : C*       *** Déclaration des variables ***
  37 : C*       *********************************
  38 :       IMPLICIT INTEGER(I-N)
  39 :       
  40 :         REAL*8 G,H,ga,W(3,3),Q(3,3)
  41 :         REAL*8 F1,DF1,D,DD,beta,dbeta    
  42 : 
  43 : C*       *********************************
  44 : C*       ******** Sous-programme *********
  45 : C*       *********************************
  46 : 
  47 : C---- Calcul des fonctions de similitude F et D et leurs dérivées
  48 :       F1=0.613*G-(3.6+76.86*(((1./G)-0.154)**2))/G
  49 : c      write(6,*) 'F1(G)=',F1
  50 :       DF1=0.613+(1./(G**2))*
  51 :      &    ((2*76.86*((1./G)-0.154)/G)+(3.6+76.86*(((1./G)-0.154)**2)))
  52 : c      write(6,*) 'F1''(G)=',DF1
  53 :       D=2*G-4.25*sqrt(G)+2.12
  54 : c      write(6,*) 'D*(G)=',D
  55 :       DD=2-4.25/(2*sqrt(G))
  56 : c      write(6,*) 'D*''(G)=',DD
  57 : 
  58 : C---- Calcul du cofficient intervenant dans la première équation
  59 :       beta=(G/F1)-1
  60 :       dbeta=(F1-G*DF1)/(F1**2)
  61 : 
  62 : C---- Matrice jacobienne, dérivées partielles des équations
  63 :       W(1,1)=-(beta*(H**2)-2*beta*H-1)/((H-1)**2)
  64 :       W(1,2)=-dbeta*(H**2)/(H-1)
  65 :       W(1,3)=0.
  66 :       W(2,1)=0.
  67 :       W(2,2)=DD
  68 :       W(2,3)=1./(ga**2)
  69 :       W(3,1)=-1./(ga*(H**2))
  70 :       W(3,2)=1.
  71 :       W(3,3)=(H-1)/(H*(ga**2))
  72 : 
  73 : C---- Inverse de la matrice jacobienne
  73 :       CALL INV(W,Q)
  74 : 
  75 : c      write(6,1004) Q(1,1),Q(1,2),Q(1,3)
  76 : c      write(6,1004) Q(1,1),Q(1,2),Q(1,3)
  77 : c      write(6,1004) Q(1,1),Q(1,2),Q(1,3)
  78 :       
  79 :       RETURN
  80 : 
  81 : c 1004 FORMAT(3(1X,1PE11.4))
  82 : 
  83 :       END
  84 : 
  85 : C*       *********************************
  86 : C*       ****  Fin Sous-Programme  *******
  87 : C*       *********************************
  88 : 
  89 : 
  90 : 
  91 : 
  92 : 
  93 : 
  94 : 
  95 : 
  96 : 
  97 : 
  98 : 
  99 : 
 100 : 
 101 : 
 102 : 
 103 : 
 104 : 
 105 : 
 106 : 
 107 : 
 108 : 
 109 : 
 110 : 
 111 : 
 112 : 
 113 : 
 114 : 
 115 :  
 116 :  
 117 :