Télécharger la plaquette Cast3M


fr en

Télécharger ottvah.eso

Retour à la liste

Numérotation des lignes : 

ottvah

  1. C OTTVAH    SOURCE    FD218221  21/06/10    21:15:51     11030          
  2.       SUBROUTINE OTTVAH(SIGMA0,DDEF0,SIGMA,DDEF,
  3.      &                  V1,V2,V3,MV)
  4. *
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7.       DIMENSION SIGMA0(6),DDEF0(6)
  8.       DIMENSION SIGMA(6), DDEF(6)
  9.       DIMENSION V1(3),V2(3),V3(3)
  10. *
  11.       DDEF04 = DDEF0(4)*0.5D0
  12.       DDEF05 = DDEF0(5)*0.5D0
  13.       DDEF06 = DDEF0(6)*0.5D0
  14. *
  15.       IF(MV.EQ.1) THEN
  16.       SIGMA(1)=SIGMA0(1)*V1(1)*V1(1)+SIGMA0(4)*V1(2)*V1(1)
  17.      &+SIGMA0(5)*V1(3)*V1(1)+SIGMA0(4)*V1(1)*V1(2)
  18.      &+SIGMA0(2)*V1(2)*V1(2)+SIGMA0(6)*V1(3)*V1(2)
  19.      &+SIGMA0(5)*V1(1)*V1(3)+SIGMA0(6)*V1(2)*V1(3)
  20.      &+SIGMA0(3)*V1(3)*V1(3)
  21.       SIGMA(2)=SIGMA0(1)*V2(1)*V2(1)+SIGMA0(4)*V2(2)*V2(1)
  22.      &+SIGMA0(5)*V2(3)*V2(1)+SIGMA0(4)*V2(1)*V2(2)
  23.      &+SIGMA0(2)*V2(2)*V2(2)+SIGMA0(6)*V2(3)*V2(2)
  24.      &+SIGMA0(5)*V2(1)*V2(3)+SIGMA0(6)*V2(2)*V2(3)
  25.      &+SIGMA0(3)*V2(3)*V2(3)
  26.       SIGMA(3)=SIGMA0(1)*V3(1)*V3(1)+SIGMA0(4)*V3(2)*V3(1)
  27.      &+SIGMA0(5)*V3(3)*V3(1)+SIGMA0(4)*V3(1)*V3(2)
  28.      &+SIGMA0(2)*V3(2)*V3(2)+SIGMA0(6)*V3(3)*V3(2)
  29.      &+SIGMA0(5)*V3(1)*V3(3)+SIGMA0(6)*V3(2)*V3(3)
  30.      &+SIGMA0(3)*V3(3)*V3(3)
  31.       SIGMA(4)=SIGMA0(1)*V1(1)*V2(1)+SIGMA0(4)*V1(2)*V2(1)
  32.      &+SIGMA0(5)*V1(3)*V2(1)+SIGMA0(4)*V1(1)*V2(2)
  33.      &+SIGMA0(2)*V1(2)*V2(2)+SIGMA0(6)*V1(3)*V2(2)
  34.      &+SIGMA0(5)*V1(1)*V2(3)+SIGMA0(6)*V1(2)*V2(3)
  35.      &+SIGMA0(3)*V1(3)*V2(3)
  36.       SIGMA(5)=SIGMA0(1)*V1(1)*V3(1)+SIGMA0(4)*V1(2)*V3(1)
  37.      &+SIGMA0(5)*V1(3)*V3(1)+SIGMA0(4)*V1(1)*V3(2)
  38.      &+SIGMA0(2)*V1(2)*V3(2)+SIGMA0(6)*V1(3)*V3(2)
  39.      &+SIGMA0(5)*V1(1)*V3(3)+SIGMA0(6)*V1(2)*V3(3)
  40.      &+SIGMA0(3)*V1(3)*V3(3)
  41.       SIGMA(6)=SIGMA0(1)*V2(1)*V3(1)+SIGMA0(4)*V2(2)*V3(1)
  42.      &+SIGMA0(5)*V2(3)*V3(1)+SIGMA0(4)*V2(1)*V3(2)
  43.      &+SIGMA0(2)*V2(2)*V3(2)+SIGMA0(6)*V2(3)*V3(2)
  44.      &+SIGMA0(5)*V2(1)*V3(3)+SIGMA0(6)*V2(2)*V3(3)
  45.      &+SIGMA0(3)*V2(3)*V3(3)
  46. *
  47.       DDEF(1)=DDEF0(1)*V1(1)*V1(1)+DDEF04*V1(2)*V1(1)
  48.      &+DDEF05*V1(3)*V1(1)+DDEF04*V1(1)*V1(2)
  49.      &+DDEF0(2)*V1(2)*V1(2)+DDEF06*V1(3)*V1(2)
  50.      &+DDEF05*V1(1)*V1(3)+DDEF06*V1(2)*V1(3)
  51.      &+DDEF0(3)*V1(3)*V1(3)
  52.       DDEF(2)=DDEF0(1)*V2(1)*V2(1)+DDEF04*V2(2)*V2(1)
  53.      &+DDEF05*V2(3)*V2(1)+DDEF04*V2(1)*V2(2)
  54.      &+DDEF0(2)*V2(2)*V2(2)+DDEF06*V2(3)*V2(2)
  55.      &+DDEF05*V2(1)*V2(3)+DDEF06*V2(2)*V2(3)
  56.      &+DDEF0(3)*V2(3)*V2(3)
  57.       DDEF(3)=DDEF0(1)*V3(1)*V3(1)+DDEF04*V3(2)*V3(1)
  58.      &+DDEF05*V3(3)*V3(1)+DDEF04*V3(1)*V3(2)
  59.      &+DDEF0(2)*V3(2)*V3(2)+DDEF06*V3(3)*V3(2)
  60.      &+DDEF05*V3(1)*V3(3)+DDEF06*V3(2)*V3(3)
  61.      &+DDEF0(3)*V3(3)*V3(3)
  62.       DDEF(4)=DDEF0(1)*V1(1)*V2(1)+DDEF04*V1(2)*V2(1)
  63.      &+DDEF05*V1(3)*V2(1)+DDEF04*V1(1)*V2(2)
  64.      &+DDEF0(2)*V1(2)*V2(2)+DDEF06*V1(3)*V2(2)
  65.      &+DDEF05*V1(1)*V2(3)+DDEF06*V1(2)*V2(3)
  66.      &+DDEF0(3)*V1(3)*V2(3)
  67.       DDEF(5)=DDEF0(1)*V1(1)*V3(1)+DDEF04*V1(2)*V3(1)
  68.      &+DDEF05*V1(3)*V3(1)+DDEF04*V1(1)*V3(2)
  69.      &+DDEF0(2)*V1(2)*V3(2)+DDEF06*V1(3)*V3(2)
  70.      &+DDEF05*V1(1)*V3(3)+DDEF06*V1(2)*V3(3)
  71.      &+DDEF0(3)*V1(3)*V3(3)
  72.       DDEF(6)=DDEF0(1)*V2(1)*V3(1)+DDEF04*V2(2)*V3(1)
  73.      &+DDEF05*V2(3)*V3(1)+DDEF04*V2(1)*V3(2)
  74.      &+DDEF0(2)*V2(2)*V3(2)+DDEF06*V2(3)*V3(2)
  75.      &+DDEF05*V2(1)*V3(3)+DDEF06*V2(2)*V3(3)
  76.      &+DDEF0(3)*V2(3)*V3(3)
  77. *
  78.       ELSE IF(MV.EQ.2) THEN
  79.       SIGMA(1)=SIGMA0(1)*V1(1)*V1(1)+SIGMA0(4)*V2(1)*V1(1)
  80.      &+SIGMA0(5)*V3(1)*V1(1)+SIGMA0(4)*V1(1)*V2(1)
  81.      &+SIGMA0(2)*V2(1)*V2(1)+SIGMA0(6)*V3(1)*V2(1)
  82.      &+SIGMA0(5)*V1(1)*V3(1)+SIGMA0(6)*V2(1)*V3(1)
  83.      &+SIGMA0(3)*V3(1)*V3(1)
  84.       SIGMA(2)=SIGMA0(1)*V1(2)*V1(2)+SIGMA0(4)*V2(2)*V1(2)
  85.      &+SIGMA0(5)*V3(2)*V1(2)+SIGMA0(4)*V1(2)*V2(2)
  86.      &+SIGMA0(2)*V2(2)*V2(2)+SIGMA0(6)*V3(2)*V2(2)
  87.      &+SIGMA0(5)*V1(2)*V3(2)+SIGMA0(6)*V2(2)*V3(2)
  88.      &+SIGMA0(3)*V3(2)*V3(2)
  89.       SIGMA(3)=SIGMA0(1)*V1(3)*V1(3)+SIGMA0(4)*V2(3)*V1(3)
  90.      &+SIGMA0(5)*V3(3)*V1(3)+SIGMA0(4)*V1(3)*V2(3)
  91.      &+SIGMA0(2)*V2(3)*V2(3)+SIGMA0(6)*V3(3)*V2(3)
  92.      &+SIGMA0(5)*V1(3)*V3(3)+SIGMA0(6)*V2(3)*V3(3)
  93.      &+SIGMA0(3)*V3(3)*V3(3)
  94.       SIGMA(4)=SIGMA0(1)*V1(1)*V1(2)+SIGMA0(4)*V2(1)*V1(2)
  95.      &+SIGMA0(5)*V3(1)*V1(2)+SIGMA0(4)*V1(1)*V2(2)
  96.      &+SIGMA0(2)*V2(1)*V2(2)+SIGMA0(6)*V3(1)*V2(2)
  97.      &+SIGMA0(5)*V1(1)*V3(2)+SIGMA0(6)*V2(1)*V3(2)
  98.      &+SIGMA0(3)*V3(1)*V3(2)
  99.       SIGMA(5)=SIGMA0(1)*V1(1)*V1(3)+SIGMA0(4)*V2(1)*V1(3)
  100.      &+SIGMA0(5)*V3(1)*V1(3)+SIGMA0(4)*V1(1)*V2(3)
  101.      &+SIGMA0(2)*V2(1)*V2(3)+SIGMA0(6)*V3(1)*V2(3)
  102.      &+SIGMA0(5)*V1(1)*V3(3)+SIGMA0(6)*V2(1)*V3(3)
  103.      &+SIGMA0(3)*V3(1)*V3(3)
  104.       SIGMA(6)=SIGMA0(1)*V1(2)*V1(3)+SIGMA0(4)*V2(2)*V1(3)
  105.      &+SIGMA0(5)*V3(2)*V1(3)+SIGMA0(4)*V1(2)*V2(3)
  106.      &+SIGMA0(2)*V2(2)*V2(3)+SIGMA0(6)*V3(2)*V2(3)
  107.      &+SIGMA0(5)*V1(2)*V3(3)+SIGMA0(6)*V2(2)*V3(3)
  108.      &+SIGMA0(3)*V3(2)*V3(3)
  109. *
  110.       DDEF(1)=DDEF0(1)*V1(1)*V1(1)+DDEF04*V2(1)*V1(1)
  111.      &+DDEF05*V3(1)*V1(1)+DDEF04*V1(1)*V2(1)
  112.      &+DDEF0(2)*V2(1)*V2(1)+DDEF06*V3(1)*V2(1)
  113.      &+DDEF05*V1(1)*V3(1)+DDEF06*V2(1)*V3(1)
  114.      &+DDEF0(3)*V3(1)*V3(1)
  115.       DDEF(2)=DDEF0(1)*V1(2)*V1(2)+DDEF04*V2(2)*V1(2)
  116.      &+DDEF05*V3(2)*V1(2)+DDEF04*V1(2)*V2(2)
  117.      &+DDEF0(2)*V2(2)*V2(2)+DDEF06*V3(2)*V2(2)
  118.      &+DDEF05*V1(2)*V3(2)+DDEF06*V2(2)*V3(2)
  119.      &+DDEF0(3)*V3(2)*V3(2)
  120.       DDEF(3)=DDEF0(1)*V1(3)*V1(3)+DDEF04*V2(3)*V1(3)
  121.      &+DDEF05*V3(3)*V1(3)+DDEF04*V1(3)*V2(3)
  122.      &+DDEF0(2)*V2(3)*V2(3)+DDEF06*V3(3)*V2(3)
  123.      &+DDEF05*V1(3)*V3(3)+DDEF06*V2(3)*V3(3)
  124.      &+DDEF0(3)*V3(3)*V3(3)
  125.       DDEF(4)=DDEF0(1)*V1(1)*V1(2)+DDEF04*V2(1)*V1(2)
  126.      &+DDEF05*V3(1)*V1(2)+DDEF04*V1(1)*V2(2)
  127.      &+DDEF0(2)*V2(1)*V2(2)+DDEF06*V3(1)*V2(2)
  128.      &+DDEF05*V1(1)*V3(2)+DDEF06*V2(1)*V3(2)
  129.      &+DDEF0(3)*V3(1)*V3(2)
  130.       DDEF(5)=DDEF0(1)*V1(1)*V1(3)+DDEF04*V2(1)*V1(3)
  131.      &+DDEF05*V3(1)*V1(3)+DDEF04*V1(1)*V2(3)
  132.      &+DDEF0(2)*V2(1)*V2(3)+DDEF06*V3(1)*V2(3)
  133.      &+DDEF05*V1(1)*V3(3)+DDEF06*V2(1)*V3(3)
  134.      &+DDEF0(3)*V3(1)*V3(3)
  135.       DDEF(6)=DDEF0(1)*V1(2)*V1(3)+DDEF04*V2(2)*V1(3)
  136.      &+DDEF05*V3(2)*V1(3)+DDEF04*V1(2)*V2(3)
  137.      &+DDEF0(2)*V2(2)*V2(3)+DDEF06*V3(2)*V2(3)
  138.      &+DDEF05*V1(2)*V3(3)+DDEF06*V2(2)*V3(3)
  139.      &+DDEF0(3)*V3(2)*V3(3)
  140.       ENDIF
  141. *
  142.       DDEF(4)=DDEF(4)*2.D0
  143.       DDEF(5)=DDEF(5)*2.D0
  144.       DDEF(6)=DDEF(6)*2.D0
  145.       RETURN
  146.       END
  147.  
  148.  
  149.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales