src/tool/AP_Conversion/pcsqu2_2.h

   1 c           .................         .................
   2 c           .               .         .       .       .
   3 c           .               .         .       .       .
   4 c           .               .         .       .       .
   5 c           .               .  ===>   .................
   6 c           .               .         .       .       .
   7 c           .               .         .       .       .
   8 c           .               .         .       .       .
   9 c           .................         .................
  10 c
  11             f1hp = filqua(quhnp1)
  12             f1cp = nqusca(f1hp)
  13             f2cp = nqusca(f1hp+1)
  14             f3cp = nqusca(f1hp+2)
  15             f4cp = nqusca(f1hp+3)
  16             prfcap(f1cp) = 1
  17             prfcap(f2cp) = 1
  18             prfcap(f3cp) = 1
  19             prfcap(f4cp) = 1
  20
  21             do 223 , nrofon = 1 , nbfonc
  22 c           pour le quadrangle Q1
  23 c
  24               vafott(nrofon,q1,f1cp) = vafoen(nrofon,q1,prqucn)
  25 c
  26               vafott(nrofon,q2,f1cp) = vafoen(nrofon,q5,prqucn)
  27 c
  28               vafott(nrofon,q3,f1cp) =
  29      >      -unsqu*(vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
  30      >              vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
  31      >       unsde*(vafoen(nrofon,q5,prqucn)+vafoen(nrofon,q6,prqucn)+
  32      >              vafoen(nrofon,q7,prqucn)+vafoen(nrofon,q8,prqucn))
  33               vafott(nrofon,q4,f1cp) = vafoen(nrofon,q8,prqucn)
  34 c
  35               vafott(nrofon,q5,f1cp) =
  36      >       -unshu*vafoen(nrofon,q2,prqucn)+
  37      >        trshu*vafoen(nrofon,q1,prqucn)+
  38      >        trsqu*vafoen(nrofon,q5,prqucn)
  39 c
  40               vafott(nrofon,q6,f1cp) =
  41      >       -trssz*(vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
  42      >               vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
  43      >        trshu*(vafoen(nrofon,q6,prqucn)+vafoen(nrofon,q8,prqucn))+
  44      >        unsqu*vafoen(nrofon,q7,prqucn)+
  45      >        trsqu*vafoen(nrofon,q5,prqucn)
  46 c
  47               vafott(nrofon,q7,f1cp) =
  48      >       -trssz*(vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
  49      >               vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
  50      >        trshu*(vafoen(nrofon,q5,prqucn)+vafoen(nrofon,q7,prqucn))+
  51      >        unsqu*vafoen(nrofon,q6,prqucn)+
  52      >        trsqu*vafoen(nrofon,q8,prqucn)
  53 c
  54               vafott(nrofon,q8,f1cp) =
  55      >       -unshu*vafoen(nrofon,q4,prqucn)+
  56      >        trshu*vafoen(nrofon,q1,prqucn)+
  57      >        trsqu*vafoen(nrofon,q8,prqucn)
  58 c
  59 c           pour le quadrangle Q2
  60 c
  61               vafott(nrofon,q1,f2cp) = vafoen(nrofon,q2,prqucn)
  62 c
  63               vafott(nrofon,q2,f2cp) = vafoen(nrofon,q6,prqucn)
  64 c
  65               vafott(nrofon,q3,f2cp) = vafott(nrofon,q3,f1cp)
  66 c
  67               vafott(nrofon,q4,f2cp) = vafoen(nrofon,q5,prqucn)
  68 c
  69               vafott(nrofon,q5,f2cp) =
  70      >       -unshu*vafoen(nrofon,q3,prqucn)+
  71      >        trshu*vafoen(nrofon,q2,prqucn)+
  72      >        trsqu*vafoen(nrofon,q6,prqucn)
  73 c
  74               vafott(nrofon,q6,f2cp) =
  75      >       -trssz*(vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
  76      >               vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
  77      >        trshu*(vafoen(nrofon,q5,prqucn)+vafoen(nrofon,q7,prqucn))+
  78      >        unsqu*vafoen(nrofon,q8,prqucn)+
  79      >        trsqu*vafoen(nrofon,q6,prqucn)
  80 c
  81               vafott(nrofon,q7,f2cp) =vafott(nrofon,q6,f1cp)
  82 c
  83               vafott(nrofon,q8,f2cp) =
  84      >       -unshu*vafoen(nrofon,q1,prqucn)+
  85      >        trshu*vafoen(nrofon,q2,prqucn)+
  86      >        trsqu*vafoen(nrofon,q5,prqucn)
  87 c
  88 c           pour le quadrangle Q3
  89 c
  90               vafott(nrofon,q1,f3cp) = vafoen(nrofon,q3,prqucn)
  91 c
  92               vafott(nrofon,q2,f3cp) = vafoen(nrofon,q7,prqucn)
  93 c
  94               vafott(nrofon,q3,f3cp) = vafott(nrofon,q3,f1cp)
  95 c
  96               vafott(nrofon,q4,f3cp) = vafoen(nrofon,q6,prqucn)
  97 c
  98               vafott(nrofon,q5,f3cp) =
  99      >       -unshu*vafoen(nrofon,q4,prqucn)+
 100      >        trshu*vafoen(nrofon,q3,prqucn)+
 101      >        trsqu*vafoen(nrofon,q7,prqucn)
 102 c
 103               vafott(nrofon,q6,f3cp) =
 104      >       -trssz*(vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
 105      >               vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
 106      >        trshu*(vafoen(nrofon,q6,prqucn)+vafoen(nrofon,q8,prqucn))+
 107      >        unsqu*vafoen(nrofon,q5,prqucn)+
 108      >        trsqu*vafoen(nrofon,q7,prqucn)
 109 c
 110               vafott(nrofon,q7,f3cp) =vafott(nrofon,q6,f2cp)
 111 c
 112               vafott(nrofon,q8,f3cp) =
 113      >       -unshu*vafoen(nrofon,q2,prqucn)+
 114      >        trshu*vafoen(nrofon,q3,prqucn)+
 115      >        trsqu*vafoen(nrofon,q6,prqucn)
 116 c
 117 c           pour le quadrangle Q4
 118 c
 119               vafott(nrofon,q1,f4cp) = vafoen(nrofon,q4,prqucn)
 120 c
 121               vafott(nrofon,q2,f4cp) = vafoen(nrofon,q8,prqucn)
 122 c
 123               vafott(nrofon,q3,f4cp) = vafott(nrofon,q3,f1cp)
 124 c
 125               vafott(nrofon,q4,f4cp) = vafoen(nrofon,q7,prqucn)
 126 c
 127               vafott(nrofon,q5,f4cp) =
 128      >       -unshu*vafoen(nrofon,q1,prqucn)+
 129      >        trshu*vafoen(nrofon,q4,prqucn)+
 130      >        trsqu*vafoen(nrofon,q8,prqucn)
 131 c
 132               vafott(nrofon,q6,f4cp) = vafott(nrofon,q7,f1cp)
 133 c
 134               vafott(nrofon,q7,f4cp) = vafott(nrofon,q6,f3cp)
 135 c
 136               vafott(nrofon,q8,f4cp) =
 137      >       -unshu*vafoen(nrofon,q3,prqucn)+
 138      >        trshu*vafoen(nrofon,q4,prqucn)+
 139      >        trsqu*vafoen(nrofon,q7,prqucn)
 140 c
 141   223       continue
 142