src/tool/AP_Conversion/pcstr2_2.h

   1 c                   .                         .
   2 c                  . .                       . .
   3 c                 .   .                     .   .
   4 c                .     .                   .     .
   5 c               .       .      ===>       .........
   6 c              .         .               . .     . .
   7 c             .           .             .   .   .   .
   8 c            .             .           .     . .     .
   9 c           .................         .................
  10 c
  11 c
  12           elseif ( etanp1.eq.4 .or.
  13      >             etanp1.eq.6 .or. etanp1.eq.7 .or. etanp1.eq.8 ) then
  14 #ifdef _DEBUG_HOMARD_
  15       write (ulsort,texte(langue,1)) 'Entree', 'pcstr2_2'
  16       write (ulsort,*) 'avec etanp1=',etanp1
  17 #endif
  18 c
  19             f1hp = filtri(trhnp1)
  20             f1cp = ntrsca(f1hp)
  21             f2cp = ntrsca(f1hp+1)
  22             f3cp = ntrsca(f1hp+2)
  23             f4cp = ntrsca(f1hp+3)
  24             prfcap(f1cp) = 1
  25             prfcap(f2cp) = 1
  26             prfcap(f3cp) = 1
  27             prfcap(f4cp) = 1
  28 c
  29             do 223 , nrofon = 1 , nbfonc
  30 c
  31 c           Pour le fils aine (centre)
  32 c
  33               vafott(nrofon,q1,f1cp) = vafoen(nrofon,q5,prtrcn)
  34 c
  35               vafott(nrofon,q2,f1cp) = vafoen(nrofon,q6,prtrcn)
  36 c
  37               vafott(nrofon,q3,f1cp) = vafoen(nrofon,q4,prtrcn)
  38 c
  39               vafott(nrofon,q4,f1cp) =
  40      >       -unshu*(vafoen(nrofon,q1,prtrcn)+vafoen(nrofon,q2,prtrcn))+
  41      >        unsde*(vafoen(nrofon,q5,prtrcn)+vafoen(nrofon,q6,prtrcn))+
  42      >         unsqu*vafoen(nrofon,q4,prtrcn)
  43 c
  44               vafott(nrofon,q5,f1cp) =
  45      >       -unshu*(vafoen(nrofon,q2,prtrcn)+vafoen(nrofon,q3,prtrcn))+
  46      >        unsde*(vafoen(nrofon,q4,prtrcn)+vafoen(nrofon,q6,prtrcn))+
  47      >         unsqu*vafoen(nrofon,q5,prtrcn)
  48 c
  49               vafott(nrofon,q6,f1cp) =
  50      >       -unshu*(vafoen(nrofon,q1,prtrcn)+vafoen(nrofon,q3,prtrcn))+
  51      >        unsde*(vafoen(nrofon,q4,prtrcn)+vafoen(nrofon,q5,prtrcn))+
  52      >         unsqu*vafoen(nrofon,q6,prtrcn)
  53 c
  54 cgn          print 1788,f1cp,
  55 cgn     >     (vafott(nrofon,iaux,f1cp), iaux = 1 , nbnoel)
  56 c
  57 c           Pour le triangle fils f1hp+1
  58 c
  59               vafott(nrofon,q1,f2cp) = vafoen(nrofon,q1,prtrcn)
  60 c
  61               vafott(nrofon,q2,f2cp) = vafoen(nrofon,q4,prtrcn)
  62 c
  63               vafott(nrofon,q3,f2cp) = vafoen(nrofon,q6,prtrcn)
  64 c
  65               vafott(nrofon,q4,f2cp) =
  66      >        trshu*vafoen(nrofon,q1,prtrcn)-
  67      >        unshu*vafoen(nrofon,q2,prtrcn)+
  68      >        trsqu*vafoen(nrofon,q4,prtrcn)
  69 c
  70               vafott(nrofon,q5,f2cp) = vafott(nrofon,q5,f1cp)
  71 c
  72               vafott(nrofon,q6,f2cp) =
  73      >        trshu*vafoen(nrofon,q1,prtrcn)-
  74      >        unshu*vafoen(nrofon,q3,prtrcn)+
  75      >        trsqu*vafoen(nrofon,q6,prtrcn)
  76 cgn          print 1788,f2cp,
  77 cgn     >     (vafott(nrofon,iaux,f2cp), iaux = 1 , nbnoel)
  78 c
  79 c           Pour le triangle fils f1hp+2
  80 c
  81               vafott(nrofon,q1,f3cp) = vafoen(nrofon,q4,prtrcn)
  82 c
  83               vafott(nrofon,q2,f3cp) = vafoen(nrofon,q2,prtrcn)
  84 c
  85               vafott(nrofon,q3,f3cp) = vafoen(nrofon,q5,prtrcn)
  86 c
  87               vafott(nrofon,q4,f3cp) =
  88      >       -unshu*vafoen(nrofon,q1,prtrcn)+
  89      >        trshu*vafoen(nrofon,q2,prtrcn)+
  90      >        trsqu*vafoen(nrofon,q4,prtrcn)
  91 c
  92               vafott(nrofon,q5,f3cp) =
  93      >        trshu*vafoen(nrofon,q2,prtrcn)-
  94      >        unshu*vafoen(nrofon,q3,prtrcn)+
  95      >        trsqu*vafoen(nrofon,q5,prtrcn)
  96 c
  97               vafott(nrofon,q6,f3cp) = vafott(nrofon,q6,f1cp)
  98 cgn          print 1788,f3cp,
  99 cgn     >     (vafott(nrofon,iaux,f3cp), iaux = 1 , nbnoel)
 100 c
 101 c           Pour le triangle fils f1hp+3
 102 c
 103               vafott(nrofon,q1,f4cp) = vafoen(nrofon,q6,prtrcn)
 104 c
 105               vafott(nrofon,q2,f4cp) = vafoen(nrofon,q5,prtrcn)
 106 c
 107               vafott(nrofon,q3,f4cp) = vafoen(nrofon,q3,prtrcn)
 108 c
 109               vafott(nrofon,q4,f4cp) = vafott(nrofon,q4,f1cp)
 110 c
 111               vafott(nrofon,q5,f4cp) =
 112      >       -unshu*vafoen(nrofon,q2,prtrcn)+
 113      >        trshu*vafoen(nrofon,q3,prtrcn)+
 114      >        trsqu*vafoen(nrofon,q5,prtrcn)
 115 c
 116               vafott(nrofon,q6,f4cp) =
 117      >       -unshu*vafoen(nrofon,q1,prtrcn)+
 118      >        trshu*vafoen(nrofon,q3,prtrcn)+
 119      >        trsqu*vafoen(nrofon,q6,prtrcn)
 120 cgn          print 1788,f4cp,
 121 cgn     >     (vafott(nrofon,iaux,f4cp), iaux = 1 , nbnoel)
 122 c
 123   223       continue
 124 c
 125 #ifdef _DEBUG_HOMARD_
 126       write (ulsort,texte(langue,1)) 'Sortie', 'pcstr2_2'
 127 #endif