src/tool/AP_Conversion/pcsqu2_5.h

   1 c   ATTENTION : les formules meriteraient d'etre verifiees ...
   2 c
   3 c           .................         .................
   4 c           .      . .      .         .       .       .
   5 c           .     .   .     .         .       .       .
   6 c           .    .     .    .         .       .       .
   7 c           .   .       .   .  ===>   .................
   8 c           .  .         .  .         .       .       .
   9 c           . .           . .         .       .       .
  10 c           ..             ..         .       .       .
  11 c           .................         .................
  12 c
  13             f1hp = filqua(quhnp1)
  14 cgn            print *,'f1hp =',f1hp
  15             f1cp = nqusca(f1hp)
  16             f2cp = nqusca(f1hp+1)
  17             f3cp = nqusca(f1hp+2)
  18             f4cp = nqusca(f1hp+3)
  19 cgn            print *,'ficp =',f1cp,f2cp,f3cp,f4cp
  20 c
  21             f1hn = -anfiqu(quhn)
  22 cgn            print *,'f1hn =',f1hn
  23             f1cn = ntreca(f1hn)
  24             f2cn = ntreca(f1hn+1)
  25             f3cn = ntreca(f1hn+2)
  26 cgn            print *,'ficn =',f1cn,f2cn,f3cn
  27             prf1cn = prftrn(f1cn)
  28             prf2cn = prftrn(f2cn)
  29             prf3cn = prftrn(f3cn)
  30 cgn            print *,'prf1cn =',prf1cn
  31 cgn            print *,'prf2cn =',prf2cn
  32 cgn            print *,'prf3cn =',prf3cn
  33 c
  34             q1t = 1
  35             q2t = 2
  36             q3t = 3
  37             q4t = 4
  38             q5t = 5
  39             q6t = 6
  40 cgn            print *,'qit =',q1t,q2t,q3t,q4t,q5t,q6t
  41 c
  42 c           decoupage en 4 quadrangles d'un quadrangle
  43 c           predecoupe en 3 triangles par l arete a1
  44
  45             if ( etan.eq.31 ) then
  46 c
  47             do 2341 , nrofon = 1 , nbfonc
  48 c
  49 c        pour le quadrangle Q1 (eventuellement redecoupe en 3 triangles)
  50 c
  51             if ( mod(hetqua(f1hp),100).eq.0 ) then
  52 c
  53             prfcap(f1cp) = 1
  54 c
  55 cgn        print *,'vatren(nrofon,q1t,prf1cn)',vatren(nrofon,q1t,prf1cn)
  56 cgn        print *,'vatren(nrofon,q2t,prf1cn)',vatren(nrofon,q2t,prf1cn)
  57 cgn        print *,'vatren(nrofon,q3t,prf1cn)',vatren(nrofon,q3t,prf1cn)
  58 cgn        print *,'vatren(nrofon,q4t,prf1cn)',vatren(nrofon,q4t,prf1cn)
  59 cgn        print *,'vatren(nrofon,q5t,prf1cn)',vatren(nrofon,q5t,prf1cn)
  60 cgn        print *,'vatren(nrofon,q6t,prf1cn)',vatren(nrofon,q6t,prf1cn)
  61 c
  62 cgn        print *,'vatren(nrofon,q1t,prf2cn)',vatren(nrofon,q1t,prf2cn)
  63 cgn        print *,'vatren(nrofon,q2t,prf2cn)',vatren(nrofon,q2t,prf2cn)
  64 cgn        print *,'vatren(nrofon,q3t,prf2cn)',vatren(nrofon,q3t,prf2cn)
  65 cgn        print *,'vatren(nrofon,q4t,prf2cn)',vatren(nrofon,q4t,prf2cn)
  66 cgn        print *,'vatren(nrofon,q5t,prf2cn)',vatren(nrofon,q5t,prf2cn)
  67 cgn        print *,'vatren(nrofon,q6t,prf2cn)',vatren(nrofon,q6t,prf2cn)
  68 c
  69 cgn        print *,'vatren(nrofon,q1t,prf3cn)',vatren(nrofon,q1t,prf3cn)
  70 cgn        print *,'vatren(nrofon,q2t,prf3cn)',vatren(nrofon,q2t,prf3cn)
  71 cgn        print *,'vatren(nrofon,q3t,prf3cn)',vatren(nrofon,q3t,prf3cn)
  72 cgn        print *,'vatren(nrofon,q4t,prf3cn)',vatren(nrofon,q4t,prf3cn)
  73 cgn        print *,'vatren(nrofon,q5t,prf3cn)',vatren(nrofon,q5t,prf3cn)
  74 cgn        print *,'vatren(nrofon,q6t,prf3cn)',vatren(nrofon,q6t,prf3cn)
  75 c
  76             vafott(nrofon,q1,f1cp) = vatren(nrofon,q2t,prf3cn)
  77 c
  78             vafott(nrofon,q2,f1cp) =
  79      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
  80      >                          vatren(nrofon,q2t,prf2cn)+
  81      >                          vatren(nrofon,q3t,prf3cn))
  82 c
  83             vafott(nrofon,q3,f1cp) =
  84      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
  85      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
  86      >     unsqu*vatren(nrofon,q6t,prf1cn)
  87 c
  88             vafott(nrofon,q4,f1cp) = vatren(nrofon,q4t,prf3cn)
  89 c
  90             vafott(nrofon,q5,f1cp) = vatren(nrofon,q5t,prf3cn)
  91 c
  92             vafott(nrofon,q6,f1cp) =trshu*(vatren(nrofon,q2t,prf1cn)+
  93      >                                     vatren(nrofon,q4t,prf1cn)+
  94      >                                     vatren(nrofon,q5t,prf1cn))-
  95      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
  96      >                              vatren(nrofon,q3t,prf1cn))+
  97      >                       unssz*vatren(nrofon,q6t,prf1cn)
  98 c
  99             vafott(nrofon,q7,f1cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 100      >                                vatren(nrofon,q6t,prf3cn))
 101 c
 102             vafott(nrofon,q8,f1cp) =
 103      >     -unshu*vatren(nrofon,q1t,prf3cn)+
 104      >      trshu*vatren(nrofon,q2t,prf3cn)+
 105      >      trsqu*vatren(nrofon,q4t,prf3cn)
 106 c
 107             elseif ( mod(hetqua(f1hp+1),100).ge.31 .and.
 108      >               mod(hetqua(f1hp+1),100).le.34 ) then
 109 c
 110
 111                 f1fhp = -filqua(f1hp)
 112                 f1fcp = ntrsca(f1fhp)
 113                 f2fcp = ntrsca(f1fhp+1)
 114                 f3fcp = ntrsca(f1fhp+2)
 115 c
 116                 prftrp(f3fcp) = 1
 117                 prftrp(f1fcp) = 1
 118                 prftrp(f2fcp) = 1
 119 c
 120 c           Pour le fils f3fcp
 121 c
 122 c    qi --> qit
 123                 vatrtt(nrofon,q1t,f3fcp) = vatren(nrofon,q4t,prf3cn)
 124 c
 125 c    qi --> qit
 126                 vatrtt(nrofon,q2t,f3fcp) = vatren(nrofon,q2t,prf3cn)
 127 c
 128 c    qi --> qit
 129                 vatrtt(nrofon,q3t,f3fcp) = vatren(nrofon,q5t,prf3cn)
 130 c
 131 c    qi --> qit
 132                 vatrtt(nrofon,q4t,f3fcp) =
 133      >           trshu*vatren(nrofon,q2t,prf3cn)
 134      >          -unshu*vatren(nrofon,q1t,prf3cn)
 135      >          +trsqu*vatren(nrofon,q4t,prf3cn)
 136 c
 137 c    qi --> qit
 138                 vatrtt(nrofon,q5t,f3fcp) =
 139      >           trshu*vatren(nrofon,q2t,prf3cn)
 140      >          -unshu*vatren(nrofon,q3t,prf3cn)
 141      >          +trsqu*vatren(nrofon,q5t,prf3cn)
 142 c
 143 c    qi --> qit
 144                 vatrtt(nrofon,q6t,f3fcp) =
 145      >      -unshu*(vatren(nrofon,q1t,prf3cn)+vatren(nrofon,q3t,prf3cn))
 146      >      +unsde*(vatren(nrofon,q4t,prf3cn)+vatren(nrofon,q5t,prf3cn))
 147      >      +unsqu*vatren(nrofon,q6t,prf3cn)
 148 c
 149 c           Pour le fils f1fcp
 150 c
 151 c    qi --> qit
 152                 vatrtt(nrofon,q1t,f1fcp) = vatren(nrofon,q4t,prf3cn)
 153 c
 154 c    qi --> qit
 155                 vatrtt(nrofon,q2t,f1fcp) = vatren(nrofon,q5t,prf3cn)
 156 c
 157 c    qi --> qit
 158                 vatrtt(nrofon,q3t,f1fcp) =
 159      >      -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))
 160      >      +unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 161      >      +unsqu*vatren(nrofon,q6t,prf1cn)
 162 c
 163 c    qi --> qit
 164                 vatrtt(nrofon,q4t,f1fcp) =
 165      >      -unshu*(vatren(nrofon,q1t,prf3cn)+vatren(nrofon,q3t,prf3cn))
 166      >      +unsde*(vatren(nrofon,q4t,prf3cn)+vatren(nrofon,q5t,prf3cn))
 167      >      +unsqu*vatren(nrofon,q6t,prf3cn)
 168 c
 169 c    qi --> qit
 170                 vatrtt(nrofon,q5t,f1fcp) =
 171      >        -unshu*unsde*(vatren(nrofon,q1t,prf3cn)+
 172      >                      vatren(nrofon,q1t,prf1cn))
 173      >        +trshu*unsde*(vatren(nrofon,q3t,prf3cn)+
 174      >                      vatren(nrofon,q2t,prf1cn))
 175      >        +trsqu*unsde*(vatren(nrofon,q4t,prf1cn)+
 176      >                      vatren(nrofon,q6t,prf3cn))
 177 c
 178 c    qi --> qit
 179                 vatrtt(nrofon,q6t,f1fcp) =
 180      >      unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q6t,prf3cn))
 181 c
 182 c           Pour le fils f2fcp
 183 c
 184 c    qi --> qit
 185                 vatrtt(nrofon,q1t,f2fcp) = vatrtt(nrofon,q3t,f1fcp)
 186 c
 187 c    qi --> qit
 188                 vatrtt(nrofon,q2t,f2fcp) = vatrtt(nrofon,q2t,f1fcp)
 189 c
 190 c    qi --> qit
 191                 vatrtt(nrofon,q3t,f2fcp) =
 192      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
 193      >                          vatren(nrofon,q2t,prf2cn)+
 194      >                          vatren(nrofon,q3t,prf3cn))
 195 c
 196 c    qi --> qit
 197                 vatrtt(nrofon,q4t,f2fcp) =
 198      >        -unshu*unsde*(vatren(nrofon,q1t,prf3cn)+
 199      >                      vatren(nrofon,q1t,prf1cn))
 200      >        +trshu*unsde*(vatren(nrofon,q3t,prf3cn)+
 201      >                      vatren(nrofon,q2t,prf1cn))
 202      >        +trsqu*unsde*(vatren(nrofon,q4t,prf1cn)+
 203      >                      vatren(nrofon,q6t,prf3cn))
 204 c
 205 c    qi --> qit
 206                 vatrtt(nrofon,q5t,f2fcp) =
 207      >         trshu*vatren(nrofon,q3t,prf3cn)-
 208      >         unshu*vatren(nrofon,q2t,prf3cn)+
 209      >         trsqu*vatren(nrofon,q5t,prf3cn)
 210 c
 211 c    qi --> qit
 212                 vatrtt(nrofon,q6t,f2fcp) =
 213      >       trshu*(vatren(nrofon,q2t,prf1cn)+
 214      >              vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 215      >        -trstr2*(vatren(nrofon,q1t,prf1cn)+
 216      >                 vatren(nrofon,q3t,prf1cn))
 217      >        +unssz*vatren(nrofon,q6t,prf1cn)
 218 c
 219 c
 220             else
 221 c
 222                 codret = codret + 1
 223 c
 224             endif
 225 c
 226 c        pour le quadrangle Q2 (eventuellement redecoupe en 3 triangles)
 227 c
 228             if ( mod(hetqua(f1hp+1),100).eq.0 ) then
 229
 230             prfcap(f2cp) = 1
 231
 232             vafott(nrofon,q1,f2cp) = vatren(nrofon,q3t,prf2cn)
 233 c
 234             vafott(nrofon,q2,f2cp) = vatren(nrofon,q6t,prf2cn)
 235 c
 236             vafott(nrofon,q3,f2cp) =
 237      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 238      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 239      >     unsqu*vatren(nrofon,q6t,prf1cn)
 240 c
 241             vafott(nrofon,q4,f2cp) =
 242      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
 243      >                          vatren(nrofon,q2t,prf2cn)+
 244      >                          vatren(nrofon,q3t,prf3cn))
 245 c
 246             vafott(nrofon,q5,f2cp) =
 247      >     -unshu*vatren(nrofon,q1t,prf2cn)+
 248      >      trshu*vatren(nrofon,q3t,prf2cn)+
 249      >      trsqu*vatren(nrofon,q6t,prf2cn)
 250 c
 251             vafott(nrofon,q6,f2cp) = unsde*(vatren(nrofon,q4t,prf2cn)+
 252      >                                vatren(nrofon,q5t,prf1cn))
 253 c
 254             vafott(nrofon,q7,f2cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 255      >                               vatren(nrofon,q4t,prf1cn)+
 256      >                               vatren(nrofon,q5t,prf1cn))-
 257      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
 258      >                               vatren(nrofon,q3t,prf1cn))+
 259      >                      unssz*vatren(nrofon,q6t,prf1cn)
 260 c
 261             vafott(nrofon,q8,f2cp) = vatren(nrofon,q5t,prf2cn)
 262 c
 263             elseif ( mod(hetqua(f1hp+1),100).ge.31 .and.
 264      >               mod(hetqua(f1hp+1),100).le.34 ) then
 265 c
 266                 f1fhp = -filqua(f1hp+1)
 267                 f1fcp = ntrsca(f1fhp)
 268                 f2fcp = ntrsca(f1fhp+1)
 269                 f3fcp = ntrsca(f1fhp+2)
 270 c
 271                 prftrp(f1fcp) = 1
 272                 prftrp(f3fcp) = 1
 273                 prftrp(f2fcp) = 1
 274 c
 275 c           Pour le fils f3fcp
 276 c
 277 c    qi --> qit
 278             vatrtt(nrofon,q1t,f3fcp) =
 279      >      -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))
 280      >      +unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 281      >      +unsqu*vatren(nrofon,q6t,prf1cn)
 282
 283 c
 284 c    qi --> qit
 285             vatrtt(nrofon,q2t,f3fcp) =
 286      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
 287      >                          vatren(nrofon,q2t,prf2cn)+
 288      >                          vatren(nrofon,q3t,prf3cn))
 289 c
 290 c    qi --> qit
 291             vatrtt(nrofon,q3t,f3fcp) = vatren(nrofon,q5t,prf2cn)
 292 c
 293 c    qi --> qit
 294             vatrtt(nrofon,q4t,f3fcp) =
 295      >       trshu*(vatren(nrofon,q2t,prf1cn)+
 296      >              vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 297      >        -trstr2*(vatren(nrofon,q1t,prf1cn)+
 298      >                 vatren(nrofon,q3t,prf1cn))
 299      >        +unssz*vatren(nrofon,q6t,prf1cn)
 300
 301 c
 302             vatrtt(nrofon,q5t,f3fcp) =
 303      >     -unshu*vatren(nrofon,q3t,prf2cn)+
 304      >      trshu*vatren(nrofon,q2t,prf2cn)+
 305      >      trsqu*vatren(nrofon,q5t,prf2cn)
 306 c
 307             vatrtt(nrofon,q6t,f3fcp) =
 308      >     -unshu*unsde*(vatren(nrofon,q1t,prf2cn)+
 309      >                   vatren(nrofon,q3t,prf1cn))+
 310      >      trshu*unsde*(vatren(nrofon,q2t,prf2cn)+
 311      >                   vatren(nrofon,q2t,prf1cn))+
 312      >      trsqu*unsde*(vatren(nrofon,q5t,prf1cn)+
 313      >                   vatren(nrofon,q4t,prf2cn))
 314 c
 315 c           Pour le fils f1fcp
 316 c
 317             vatrtt(nrofon,q1t,f1fcp) = vatrtt(nrofon,q1t,f3fcp)
 318
 319 c
 320             vatrtt(nrofon,q2t,f1fcp) = vatren(nrofon,q5t,prf2cn)
 321 c
 322             vatrtt(nrofon,q3t,f1fcp) = vatren(nrofon,q6t,prf2cn)
 323 c
 324             vatrtt(nrofon,q4t,f1fcp) = vatrtt(nrofon,q6t,f3fcp)
 325 c
 326             vatrtt(nrofon,q5t,f1fcp) =
 327      >     -unshu*(vatren(nrofon,q2t,prf2cn)+vatren(nrofon,q1t,prf2cn))
 328      >     +unsde*(vatren(nrofon,q5t,prf2cn)+vatren(nrofon,q6t,prf2cn))
 329      >     +unsqu*vatren(nrofon,q4t,prf2cn)
 330 c
 331             vatrtt(nrofon,q6t,f1fcp) =
 332      >     unsde*(vatren(nrofon,q4t,prf2cn)+
 333      >            vatren(nrofon,q5t,prf1cn))
 334 c
 335 c           Pour le fils f2fcp
 336 c
 337             vatrtt(nrofon,q1t,f2fcp) = vatren(nrofon,q6t,prf2cn)
 338
 339 c
 340             vatrtt(nrofon,q2t,f2fcp) = vatren(nrofon,q5t,prf2cn)
 341 c
 342             vatrtt(nrofon,q3t,f2fcp) = vatren(nrofon,q3t,prf2cn)
 343 c
 344             vatrtt(nrofon,q4t,f2fcp) =
 345      >     -unshu*(vatren(nrofon,q2t,prf2cn)+vatren(nrofon,q1t,prf2cn))
 346      >     +unsde*(vatren(nrofon,q5t,prf2cn)+vatren(nrofon,q6t,prf2cn))
 347      >     +unsqu*vatren(nrofon,q4t,prf2cn)
 348 c
 349             vatrtt(nrofon,q5t,f2fcp) =
 350      >      trshu*vatren(nrofon,q3t,prf2cn)-
 351      >      unshu*vatren(nrofon,q2t,prf2cn)+
 352      >      trsqu*vatren(nrofon,q5t,prf2cn)
 353 c
 354             vatrtt(nrofon,q6t,f2fcp) =
 355      >      trshu*vatren(nrofon,q3t,prf2cn)-
 356      >      unshu*vatren(nrofon,q1t,prf2cn)+
 357      >      trsqu*vatren(nrofon,q6t,prf2cn)
 358 c
 359 c
 360             else
 361 c
 362                 codret = codret + 1
 363 c
 364             endif
 365
 366 c         pour le quadrangle Q3
 367 c
 368             prfcap(f3cp) = 1
 369 c
 370             vafott(nrofon,q1,f3cp) = unsde*(vatren(nrofon,q1t,prf2cn)+
 371      >                                vatren(nrofon,q3t,prf1cn))
 372 c
 373             vafott(nrofon,q2,f3cp) = vatren(nrofon,q6t,prf1cn)
 374 c
 375             vafott(nrofon,q3,f3cp) =
 376      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 377      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 378      >     unsqu*vatren(nrofon,q6t,prf1cn)
 379 c
 380             vafott(nrofon,q4,f3cp) = vatren(nrofon,q6t,prf2cn)
 381 c
 382             vafott(nrofon,q5,f3cp) =
 383      >     -unshu*vatren(nrofon,q1t,prf1cn)+
 384      >      trshu*vatren(nrofon,q3t,prf1cn)+
 385      >      trsqu*vatren(nrofon,q6t,prf1cn)
 386 c
 387             vafott(nrofon,q6,f3cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 388      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 389      >                                vatren(nrofon,q1t,prf1cn))+
 390      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 391      >                                vatren(nrofon,q5t,prf1cn))+
 392      >                      nessz*vatren(nrofon,q6t,prf1cn)
 393 c
 394             vafott(nrofon,q7,f3cp) =unsde*(vatren(nrofon,q4t,prf2cn)+
 395      >                                vatren(nrofon,q5t,prf1cn))
 396 c
 397             vafott(nrofon,q8,f3cp) =
 398      >     -unshu*vatren(nrofon,q3t,prf2cn)+
 399      >      trshu*vatren(nrofon,q1t,prf2cn)+
 400      >      trsqu*vatren(nrofon,q6t,prf2cn)
 401 c
 402 c         pour le quadrangle Q4
 403 c
 404             prfcap(f4cp) = 1
 405 c
 406             vafott(nrofon,q1,f4cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 407      >                                vatren(nrofon,q1t,prf3cn))
 408 c
 409             vafott(nrofon,q2,f4cp) = vatren(nrofon,q4t,prf3cn)
 410 c
 411             vafott(nrofon,q3,f4cp) =
 412      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 413      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 414      >     unsqu*vatren(nrofon,q6t,prf1cn)
 415 c
 416             vafott(nrofon,q4,f4cp) = vatren(nrofon,q6t,prf1cn)
 417 c
 418             vafott(nrofon,q5,f4cp) =
 419      >     -unshu*vatren(nrofon,q2t,prf3cn)+
 420      >      trshu*vatren(nrofon,q1t,prf3cn)+
 421      >      trsqu*vatren(nrofon,q4t,prf3cn)
 422 c
 423             vafott(nrofon,q6,f4cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 424      >                                vatren(nrofon,q6t,prf3cn))
 425 c
 426             vafott(nrofon,q7,f4cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 427      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 428      >                                vatren(nrofon,q1t,prf1cn))+
 429      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 430      >                                vatren(nrofon,q5t,prf1cn))+
 431      >                      nessz*vatren(nrofon,q6t,prf1cn)
 432 c
 433             vafott(nrofon,q8,f4cp) =
 434      >     -unshu*vatren(nrofon,q3t,prf1cn)+
 435      >      trshu*vatren(nrofon,q1t,prf1cn)+
 436      >      trsqu*vatren(nrofon,q6t,prf1cn)
 437 c
 438  2341     continue
 439 c
 440 c           decoupage en 4 quadrangles d'un quadrangle
 441 c           predecoupe en 3 triangles par l arete a2
 442 c
 443          elseif (etan .eq. 32) then
 444 c
 445          do 2342 , nrofon = 1 , nbfonc
 446
 447 c        pour le quadrangle Q1
 448 c
 449            prfcap(f1cp) = 1
 450 c
 451            vafott(nrofon,q1,f1cp) = unsde*(vatren(nrofon,q1t,prf3cn)+
 452      >                              vatren(nrofon,q1t,prf1cn))
 453 c
 454            vafott(nrofon,q2,f1cp) = vatren(nrofon,q4t,prf3cn)
 455 c
 456            vafott(nrofon,q3,f1cp) =
 457      >   -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 458      >    unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 459      >    unsqu*vatren(nrofon,q6t,prf1cn)
 460 c
 461            vafott(nrofon,q4,f1cp) = vatren(nrofon,q6t,prf1cn)
 462 c
 463            vafott(nrofon,q5,f1cp) =
 464      >    -unshu*vatren(nrofon,q2t,prf3cn)+
 465      >     trshu*vatren(nrofon,q1t,prf3cn)+
 466      >     trsqu*vatren(nrofon,q4t,prf3cn)
 467 c
 468            vafott(nrofon,q6,f1cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 469      >                              vatren(nrofon,q6t,prf3cn))
 470 c
 471            vafott(nrofon,q7,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 472      >                     trstr2*(vatren(nrofon,q3t,prf1cn)+
 473      >                              vatren(nrofon,q1t,prf1cn))+
 474      >                         trshu*(vatren(nrofon,q4t,prf1cn)+
 475      >                              vatren(nrofon,q5t,prf1cn))+
 476      >                      nessz*vatren(nrofon,q6t,prf1cn)
 477 c
 478            vafott(nrofon,q8,f1cp) =
 479      >    -unshu*vatren(nrofon,q3t,prf1cn)+
 480      >     trshu*vatren(nrofon,q1t,prf1cn)+
 481      >     trsqu*vatren(nrofon,q6t,prf1cn)
 482 c
 483 c        pour le quadrangle Q2
 484 c
 485             if ( mod(hetqua(f1hp+1),100).eq.0 ) then
 486
 487             prfcap(f2cp) = 1
 488
 489            vafott(nrofon,q1,f2cp) = vatren(nrofon,q2t,prf3cn)
 490 c
 491            vafott(nrofon,q2,f2cp) =
 492      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
 493      >                          vatren(nrofon,q2t,prf2cn)+
 494      >                          vatren(nrofon,q3t,prf3cn))
 495 c
 496            vafott(nrofon,q3,f2cp) =
 497      >   -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 498      >    unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 499      >    unsqu*vatren(nrofon,q6t,prf1cn)
 500 c
 501            vafott(nrofon,q4,f2cp) = vatren(nrofon,q4t,prf3cn)
 502 c
 503            vafott(nrofon,q5,f2cp) = vatren(nrofon,q5t,prf3cn)
 504 c
 505            vafott(nrofon,q6,f2cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 506      >                              vatren(nrofon,q4t,prf1cn)+
 507      >                              vatren(nrofon,q5t,prf1cn))-
 508      >                     trstr2*(vatren(nrofon,q1t,prf1cn)+
 509      >                              vatren(nrofon,q3t,prf1cn))+
 510      >                     unssz*vatren(nrofon,q6t,prf1cn)
 511 c
 512            vafott(nrofon,q7,f2cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 513      >                              vatren(nrofon,q6t,prf3cn))
 514 c
 515            vafott(nrofon,q8,f2cp) =
 516      >    -unshu*vatren(nrofon,q1t,prf3cn)+
 517      >     trshu*vatren(nrofon,q2t,prf3cn)+
 518      >     trsqu*vatren(nrofon,q4t,prf3cn)
 519 c
 520             elseif ( mod(hetqua(f1hp+1),100).ge.31 .and.
 521      >               mod(hetqua(f1hp+1),100).le.34 ) then
 522 c
 523                 f1fhp = -filqua(f1hp+1)
 524                 f1fcp = ntrsca(f1fhp)
 525                 f2fcp = ntrsca(f1fhp+1)
 526                 f3fcp = ntrsca(f1fhp+2)
 527
 528                 prftrp(f3fcp) = 1
 529                 prftrp(f1fcp) = 1
 530                 prftrp(f2fcp) = 1
 531 c
 532 c
 533 c           Pour le fils f3fcp
 534 c
 535 c    qi --> qit
 536                 vatrtt(nrofon,q1t,f3fcp) = vatren(nrofon,q4t,prf3cn)
 537 c
 538 c    qi --> qit
 539                 vatrtt(nrofon,q2t,f3fcp) = vatren(nrofon,q2t,prf3cn)
 540 c
 541 c    qi --> qit
 542                 vatrtt(nrofon,q3t,f3fcp) = vatren(nrofon,q5t,prf3cn)
 543 c
 544 c    qi --> qit
 545                 vatrtt(nrofon,q4t,f3fcp) =
 546      >           trshu*vatren(nrofon,q2t,prf3cn)
 547      >          -unshu*vatren(nrofon,q1t,prf3cn)
 548      >          +trsqu*vatren(nrofon,q4t,prf3cn)
 549 c
 550 c    qi --> qit
 551                 vatrtt(nrofon,q5t,f3fcp) =
 552      >           trshu*vatren(nrofon,q2t,prf3cn)
 553      >          -unshu*vatren(nrofon,q3t,prf3cn)
 554      >          +trsqu*vatren(nrofon,q5t,prf3cn)
 555 c
 556 c    qi --> qit
 557                 vatrtt(nrofon,q6t,f3fcp) =
 558      >      -unshu*(vatren(nrofon,q1t,prf3cn)+vatren(nrofon,q3t,prf3cn))
 559      >      +unsde*(vatren(nrofon,q4t,prf3cn)+vatren(nrofon,q5t,prf3cn))
 560      >      +unsqu*vatren(nrofon,q6t,prf3cn)
 561 c
 562 c           Pour le fils f1fcp
 563 c
 564 c    qi --> qit
 565                 vatrtt(nrofon,q1t,f1fcp) = vatren(nrofon,q4t,prf3cn)
 566 c
 567 c    qi --> qit
 568                 vatrtt(nrofon,q2t,f1fcp) = vatren(nrofon,q5t,prf3cn)
 569 c
 570 c    qi --> qit
 571                 vatrtt(nrofon,q3t,f1fcp) =
 572      >      -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))
 573      >      +unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 574      >      +unsqu*vatren(nrofon,q6t,prf1cn)
 575 c
 576 c    qi --> qit
 577                 vatrtt(nrofon,q4t,f1fcp) = vatrtt(nrofon,q6t,f3fcp)
 578 c
 579 c    qi --> qit
 580                 vatrtt(nrofon,q5t,f1fcp) =
 581      >        -unshu*unsde*(vatren(nrofon,q1t,prf3cn)+
 582      >                      vatren(nrofon,q1t,prf1cn))
 583      >        +trshu*unsde*(vatren(nrofon,q3t,prf3cn)+
 584      >                      vatren(nrofon,q2t,prf1cn))
 585      >        +trsqu*unsde*(vatren(nrofon,q4t,prf1cn)+
 586      >                      vatren(nrofon,q6t,prf3cn))
 587 c
 588 c    qi --> qit
 589                 vatrtt(nrofon,q6t,f1fcp) =
 590      >      unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q6t,prf3cn))
 591 c
 592 c           Pour le fils f2fcp
 593 c
 594 c    qi --> qit
 595                 vatrtt(nrofon,q1t,f2fcp) = vatrtt(nrofon,q3t,f1fcp)
 596 c
 597 c    qi --> qit
 598                 vatrtt(nrofon,q2t,f2fcp) = vatrtt(nrofon,q2t,f1fcp)
 599 c
 600 c    qi --> qit
 601                 vatrtt(nrofon,q3t,f2fcp) =
 602      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
 603      >                          vatren(nrofon,q2t,prf2cn)+
 604      >                          vatren(nrofon,q3t,prf3cn))
 605 c
 606 c    qi --> qit
 607                 vatrtt(nrofon,q4t,f2fcp) = vatrtt(nrofon,q5t,f1fcp)
 608 c
 609 c    qi --> qit
 610                 vatrtt(nrofon,q5t,f2fcp) =
 611      >         trshu*vatren(nrofon,q3t,prf3cn)-
 612      >         unshu*vatren(nrofon,q2t,prf3cn)+
 613      >         trsqu*vatren(nrofon,q5t,prf3cn)
 614 c
 615 c    qi --> qit
 616                 vatrtt(nrofon,q6t,f2fcp) =
 617      >       trshu*(vatren(nrofon,q2t,prf1cn)+
 618      >              vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 619      >        -trstr2*(vatren(nrofon,q1t,prf1cn)+
 620      >                 vatren(nrofon,q3t,prf1cn))
 621      >        +unssz*vatren(nrofon,q6t,prf1cn)
 622 c
 623 c
 624             else
 625 c
 626                 codret = codret + 1
 627 c
 628             endif
 629
 630 c        pour le quadrangle Q3
 631 c
 632             if ( mod(hetqua(f1hp+2),100).eq.0 ) then
 633
 634             prfcap(f3cp) = 1
 635 c
 636            vafott(nrofon,q1,f3cp) = vatren(nrofon,q3t,prf2cn)
 637 c
 638            vafott(nrofon,q2,f3cp) = vatren(nrofon,q6t,prf2cn)
 639 c
 640            vafott(nrofon,q3,f3cp) =
 641      >   -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 642      >    unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 643      >    unsqu*vatren(nrofon,q6t,prf1cn)
 644 c
 645            vafott(nrofon,q4,f3cp) =
 646      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
 647      >                          vatren(nrofon,q2t,prf2cn)+
 648      >                          vatren(nrofon,q3t,prf3cn))
 649 c
 650            vafott(nrofon,q5,f3cp) =
 651      >    -unshu*vatren(nrofon,q1t,prf2cn)+
 652      >     trshu*vatren(nrofon,q3t,prf2cn)+
 653      >     trsqu*vatren(nrofon,q6t,prf2cn)
 654 c
 655            vafott(nrofon,q6,f3cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
 656      >                              vatren(nrofon,q4t,prf2cn))
 657 c
 658            vafott(nrofon,q7,f3cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 659      >                              vatren(nrofon,q4t,prf1cn)+
 660      >                              vatren(nrofon,q5t,prf1cn))-
 661      >                     trstr2*(vatren(nrofon,q1t,prf1cn)+
 662      >                              vatren(nrofon,q3t,prf1cn))+
 663      >                     unssz*vatren(nrofon,q6t,prf1cn)
 664 c
 665            vafott(nrofon,q8,f3cp) = vatren(nrofon,q5t,prf2cn)
 666 c
 667             elseif ( mod(hetqua(f1hp+2),100).ge.31 .and.
 668      >               mod(hetqua(f1hp+2),100).le.34 ) then
 669 c
 670                 f1fhp = -filqua(f1hp+2)
 671                 f3fcp = ntrsca(f1fhp+2)
 672                 f1fcp = ntrsca(f1fhp)
 673                 f2fcp = ntrsca(f1fhp+1)
 674
 675                 prftrp(f3fcp) = 1
 676                 prftrp(f1fcp) = 1
 677                 prftrp(f2fcp) = 1
 678 c
 679 c
 680 c           Pour le fils f3fcp
 681 c
 682 c   qi --> qit
 683             vatrtt(nrofon,q1t,f3fcp) =
 684      >      -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))
 685      >      +unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 686      >      +unsqu*vatren(nrofon,q6t,prf1cn)
 687
 688 c
 689             vatrtt(nrofon,q2t,f3fcp) =
 690      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
 691      >                          vatren(nrofon,q2t,prf2cn)+
 692      >                          vatren(nrofon,q3t,prf3cn))
 693 c
 694 c   qi --> qit
 695             vatrtt(nrofon,q3t,f3fcp) = vatren(nrofon,q5t,prf2cn)
 696 c
 697 c   qi --> qit
 698             vatrtt(nrofon,q4t,f3fcp) =
 699      >       trshu*(vatren(nrofon,q2t,prf1cn)+
 700      >              vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 701      >        -trstr2*(vatren(nrofon,q1t,prf1cn)+
 702      >                 vatren(nrofon,q3t,prf1cn))
 703      >        +unssz*vatren(nrofon,q6t,prf1cn)
 704
 705 c
 706             vatrtt(nrofon,q5t,f3fcp) =
 707      >     -unshu*vatren(nrofon,q3t,prf2cn)+
 708      >      trshu*vatren(nrofon,q2t,prf2cn)+
 709      >      trsqu*vatren(nrofon,q5t,prf2cn)
 710 c
 711             vatrtt(nrofon,q6t,f3fcp) =
 712      >     -unshu*unsde*(vatren(nrofon,q1t,prf2cn)+
 713      >                   vatren(nrofon,q3t,prf1cn))+
 714      >      trshu*unsde*(vatren(nrofon,q2t,prf2cn)+
 715      >                   vatren(nrofon,q2t,prf1cn))+
 716      >      trsqu*unsde*(vatren(nrofon,q5t,prf1cn)+
 717      >                   vatren(nrofon,q4t,prf2cn))
 718
 719 c           Pour le fils f1fcp
 720 c
 721             vatrtt(nrofon,q1t,f1fcp) = vatrtt(nrofon,q1t,f3fcp)
 722
 723 c
 724             vatrtt(nrofon,q2t,f1fcp) = vatren(nrofon,q5t,prf2cn)
 725 c
 726             vatrtt(nrofon,q3t,f1fcp) = vatren(nrofon,q6t,prf2cn)
 727 c
 728             vatrtt(nrofon,q4t,f1fcp) = vatrtt(nrofon,q6t,f3fcp)
 729 c
 730             vatrtt(nrofon,q5t,f1fcp) =
 731      >     -unshu*(vatren(nrofon,q2t,prf2cn)+vatren(nrofon,q1t,prf2cn))
 732      >     +unsde*(vatren(nrofon,q5t,prf2cn)+vatren(nrofon,q6t,prf2cn))
 733      >     +unsqu*vatren(nrofon,q4t,prf2cn)
 734 c
 735             vatrtt(nrofon,q6t,f1fcp) =
 736      >     unsde*(vatren(nrofon,q4t,prf2cn)+
 737      >            vatren(nrofon,q5t,prf1cn))
 738
 739 c
 740 c           Pour le fils f2fcp
 741 c
 742             vatrtt(nrofon,q1t,f2fcp) = vatren(nrofon,q6t,prf2cn)
 743
 744 c
 745             vatrtt(nrofon,q2t,f2fcp) = vatren(nrofon,q5t,prf2cn)
 746 c
 747             vatrtt(nrofon,q3t,f2fcp) = vatren(nrofon,q3t,prf2cn)
 748 c
 749             vatrtt(nrofon,q4t,f2fcp) =
 750      >     -unshu*(vatren(nrofon,q2t,prf2cn)+vatren(nrofon,q1t,prf2cn))
 751      >     +unsde*(vatren(nrofon,q5t,prf2cn)+vatren(nrofon,q6t,prf2cn))
 752      >     +unsqu*vatren(nrofon,q4t,prf2cn)
 753
 754             vatrtt(nrofon,q5t,f2fcp) =
 755      >      trshu*vatren(nrofon,q3t,prf2cn)-
 756      >      unshu*vatren(nrofon,q2t,prf2cn)+
 757      >      trsqu*vatren(nrofon,q5t,prf2cn)
 758 c
 759             vatrtt(nrofon,q6t,f2fcp) =
 760      >      trshu*vatren(nrofon,q3t,prf2cn)-
 761      >      unshu*vatren(nrofon,q1t,prf2cn)+
 762      >      trsqu*vatren(nrofon,q6t,prf2cn)
 763 c
 764             else
 765 c
 766                 codret = codret + 1
 767 c
 768             endif
 769 c        pour le quadrangle Q4
 770 c
 771            prfcap(f4cp) = 1
 772 c
 773            vafott(nrofon,q1,f4cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 774      >                              vatren(nrofon,q1t,prf2cn))
 775 c
 776            vafott(nrofon,q2,f4cp) = vatren(nrofon,q6t,prf1cn)
 777 c
 778            vafott(nrofon,q3,f4cp) =
 779      >   -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 780      >    unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 781      >    unsqu*vatren(nrofon,q6t,prf1cn)
 782 c
 783            vafott(nrofon,q4,f4cp) = vatren(nrofon,q6t,prf2cn)
 784 c
 785            vafott(nrofon,q5,f4cp) =
 786      >    -unshu*vatren(nrofon,q1t,prf1cn)+
 787      >     trshu*vatren(nrofon,q3t,prf1cn)+
 788      >     trsqu*vatren(nrofon,q6t,prf1cn)
 789 c
 790            vafott(nrofon,q6,f4cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 791      >                     trstr2*(vatren(nrofon,q3t,prf1cn)+
 792      >                              vatren(nrofon,q1t,prf1cn))+
 793      >                         trshu*(vatren(nrofon,q4t,prf1cn)+
 794      >                              vatren(nrofon,q5t,prf1cn))+
 795      >                      nessz*vatren(nrofon,q6t,prf1cn)
 796 c
 797            vafott(nrofon,q7,f4cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
 798      >                              vatren(nrofon,q4t,prf2cn))
 799 c
 800            vafott(nrofon,q8,f4cp) =
 801      >    -unshu*vatren(nrofon,q3t,prf2cn)+
 802      >     trshu*vatren(nrofon,q1t,prf2cn)+
 803      >     trsqu*vatren(nrofon,q6t,prf2cn)
 804 c
 805  2342    continue
 806 c
 807 c           decoupage en 4 quadrangles d'un quadrangle
 808 c           predecoupe en 3 triangles par l arete a3
 809 c
 810          elseif (etan .eq. 33) then
 811          do 2343 , nrofon = 1 , nbfonc
 812 c
 813 c        pour le quadrangle Q1
 814 c
 815            prfcap(f1cp) = 1
 816 c
 817            vafott(nrofon,q1,f1cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 818      >                             vatren(nrofon,q1t,prf2cn))
 819 c
 820            vafott(nrofon,q2,f1cp) = vatren(nrofon,q6t,prf1cn)
 821 c
 822            vafott(nrofon,q3,f1cp) =
 823      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 824      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 825      >     unsqu*vatren(nrofon,q6t,prf1cn)
 826 c
 827            vafott(nrofon,q4,f1cp) = vatren(nrofon,q6t,prf2cn)
 828 c
 829            vafott(nrofon,q5,f1cp) =
 830      >    -unshu*vatren(nrofon,q1t,prf1cn)+
 831      >     trshu*vatren(nrofon,q3t,prf1cn)+
 832      >     trsqu*vatren(nrofon,q6t,prf1cn)
 833 c
 834            vafott(nrofon,q6,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 835      >                    trstr2*(vatren(nrofon,q3t,prf1cn)+
 836      >                             vatren(nrofon,q1t,prf1cn))+
 837      >                        trshu*(vatren(nrofon,q4t,prf1cn)+
 838      >                             vatren(nrofon,q5t,prf1cn))+
 839      >                     nessz*vatren(nrofon,q6t,prf1cn)
 840 c
 841            vafott(nrofon,q7,f1cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
 842      >                             vatren(nrofon,q4t,prf2cn))
 843 c
 844            vafott(nrofon,q8,f1cp) =
 845      >    -unshu*vatren(nrofon,q3t,prf2cn)+
 846      >     trshu*vatren(nrofon,q1t,prf2cn)+
 847      >     trsqu*vatren(nrofon,q6t,prf2cn)
 848 c
 849 c        pour le quadrangle Q2
 850 c
 851            prfcap(f2cp) = 1
 852 c
 853            vafott(nrofon,q1,f2cp) = unsde*(vatren(nrofon,q1t,prf3cn)+
 854      >                             vatren(nrofon,q1t,prf1cn))
 855 c
 856            vafott(nrofon,q2,f2cp) = vatren(nrofon,q4t,prf3cn)
 857 c
 858            vafott(nrofon,q3,f2cp) =
 859      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 860      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 861      >     unsqu*vatren(nrofon,q6t,prf1cn)
 862 c
 863            vafott(nrofon,q4,f2cp) = vatren(nrofon,q6t,prf1cn)
 864 c
 865            vafott(nrofon,q5,f2cp) =
 866      >    -unshu*vatren(nrofon,q2t,prf3cn)+
 867      >     trshu*vatren(nrofon,q1t,prf3cn)+
 868      >     trsqu*vatren(nrofon,q4t,prf3cn)
 869 c
 870            vafott(nrofon,q6,f2cp) =  unsde*(vatren(nrofon,q4t,prf1cn)+
 871      >                             vatren(nrofon,q6t,prf3cn))
 872 c
 873            vafott(nrofon,q7,f2cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 874      >                    trstr2*(vatren(nrofon,q3t,prf1cn)+
 875      >                             vatren(nrofon,q1t,prf1cn))+
 876      >                        trshu*(vatren(nrofon,q4t,prf1cn)+
 877      >                             vatren(nrofon,q5t,prf1cn))+
 878      >                     nessz*vatren(nrofon,q6t,prf1cn)
 879 c
 880            vafott(nrofon,q8,f2cp) =
 881      >    -unshu*vatren(nrofon,q3t,prf1cn)+
 882      >     trshu*vatren(nrofon,q1t,prf1cn)+
 883      >     trsqu*vatren(nrofon,q6t,prf1cn)
 884 c
 885 c        pour le quadrangle Q3
 886 c
 887             if ( mod(hetqua(f1hp+2),100).eq.0 ) then
 888
 889             prfcap(f3cp) = 1
 890 c
 891            vafott(nrofon,q1,f3cp) = vatren(nrofon,q2t,prf3cn)
 892 c
 893            vafott(nrofon,q2,f3cp) =
 894      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
 895      >                          vatren(nrofon,q2t,prf2cn)+
 896      >                          vatren(nrofon,q3t,prf3cn))
 897 c
 898            vafott(nrofon,q3,f3cp) =
 899      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
 900      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
 901      >     unsqu*vatren(nrofon,q6t,prf1cn)
 902 c
 903            vafott(nrofon,q4,f3cp) = vatren(nrofon,q4t,prf3cn)
 904 c
 905            vafott(nrofon,q5,f3cp) = vatren(nrofon,q5t,prf3cn)
 906 c
 907            vafott(nrofon,q6,f3cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 908      >                             vatren(nrofon,q4t,prf1cn)+
 909      >                             vatren(nrofon,q5t,prf1cn))-
 910      >                    trstr2*(vatren(nrofon,q1t,prf1cn)+
 911      >                             vatren(nrofon,q3t,prf1cn))+
 912      >                    unssz*vatren(nrofon,q6t,prf1cn)
 913 c
 914            vafott(nrofon,q7,f3cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 915      >                             vatren(nrofon,q6t,prf3cn))
 916 c
 917            vafott(nrofon,q8,f3cp) =
 918      >    -unshu*vatren(nrofon,q1t,prf3cn)+
 919      >     trshu*vatren(nrofon,q2t,prf3cn)+
 920      >     trsqu*vatren(nrofon,q4t,prf3cn)
 921 c
 922             elseif ( mod(hetqua(f1hp+2),100).ge.31 .and.
 923      >               mod(hetqua(f1hp+2),100).le.34 ) then
 924 c
 925                 f1fhp = -filqua(f1hp+2)
 926                 f3fcp = ntrsca(f1fhp+2)
 927                 f1fcp = ntrsca(f1fhp)
 928                 f2fcp = ntrsca(f1fhp+1)
 929
 930                 prftrp(f3fcp) = 1
 931                 prftrp(f1fcp) = 1
 932                 prftrp(f2fcp) = 1
 933 c
 934 c
 935 c           Pour le fils f3fcp
 936 c
 937 c    qi --> qit
 938                 vatrtt(nrofon,q1t,f3fcp) = vatren(nrofon,q4t,prf3cn)
 939 c
 940 c    qi --> qit
 941                 vatrtt(nrofon,q2t,f3fcp) = vatren(nrofon,q2t,prf3cn)
 942 c
 943 c    qi --> qit
 944                 vatrtt(nrofon,q3t,f3fcp) = vatren(nrofon,q5t,prf3cn)
 945 c
 946 c    qi --> qit
 947                 vatrtt(nrofon,q4t,f3fcp) =
 948      >           trshu*vatren(nrofon,q2t,prf3cn)
 949      >          -unshu*vatren(nrofon,q1t,prf3cn)
 950      >          +trsqu*vatren(nrofon,q4t,prf3cn)
 951 c
 952 c    qi --> qit
 953                 vatrtt(nrofon,q5t,f3fcp) =
 954      >           trshu*vatren(nrofon,q2t,prf3cn)
 955      >          -unshu*vatren(nrofon,q3t,prf3cn)
 956      >          +trsqu*vatren(nrofon,q5t,prf3cn)
 957 c
 958 c    qi --> qit
 959                 vatrtt(nrofon,q6t,f3fcp) =
 960      >      -unshu*(vatren(nrofon,q1t,prf3cn)+vatren(nrofon,q3t,prf3cn))
 961      >      +unsde*(vatren(nrofon,q4t,prf3cn)+vatren(nrofon,q5t,prf3cn))
 962      >      +unsqu*vatren(nrofon,q6t,prf3cn)
 963 c
 964 c           Pour le fils f1fcp
 965 c
 966 c    qi --> qit
 967                 vatrtt(nrofon,q1t,f1fcp) = vatren(nrofon,q4t,prf3cn)
 968 c
 969 c    qi --> qit
 970                 vatrtt(nrofon,q2t,f1fcp) = vatren(nrofon,q5t,prf3cn)
 971 c
 972 c    qi --> qit
 973                 vatrtt(nrofon,q3t,f1fcp) =
 974      >      -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))
 975      >      +unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
 976      >      +unsqu*vatren(nrofon,q6t,prf1cn)
 977 c
 978 c    qi --> qit
 979                 vatrtt(nrofon,q4t,f1fcp) = vatrtt(nrofon,q6t,f3fcp)
 980 c
 981 c    qi --> qit
 982                 vatrtt(nrofon,q5t,f1fcp) =
 983      >        -unshu*unsde*(vatren(nrofon,q1t,prf3cn)+
 984      >                      vatren(nrofon,q1t,prf1cn))
 985      >        +trshu*unsde*(vatren(nrofon,q3t,prf3cn)+
 986      >                      vatren(nrofon,q2t,prf1cn))
 987      >        +trsqu*unsde*(vatren(nrofon,q4t,prf1cn)+
 988      >                      vatren(nrofon,q6t,prf3cn))
 989 c
 990 c    qi --> qit
 991                 vatrtt(nrofon,q6t,f1fcp) =
 992      >      unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q6t,prf3cn))
 993 c
 994 c           Pour le fils f2fcp
 995 c
 996 c    qi --> qit
 997                 vatrtt(nrofon,q1t,f2fcp) = vatrtt(nrofon,q3t,f1fcp)
 998 c
 999 c    qi --> qit
1000                 vatrtt(nrofon,q2t,f2fcp) = vatrtt(nrofon,q2t,f1fcp)
1001 c
1002 c    qi --> qit
1003                 vatrtt(nrofon,q3t,f2fcp) =
1004      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
1005      >                          vatren(nrofon,q2t,prf2cn)+
1006      >                          vatren(nrofon,q3t,prf3cn))
1007 c
1008 c    qi --> qit
1009                 vatrtt(nrofon,q4t,f2fcp) = vatrtt(nrofon,q5t,f1fcp)
1010 c
1011 c    qi --> qit
1012                 vatrtt(nrofon,q5t,f2fcp) =
1013      >         trshu*vatren(nrofon,q3t,prf3cn)-
1014      >         unshu*vatren(nrofon,q2t,prf3cn)+
1015      >         trsqu*vatren(nrofon,q5t,prf3cn)
1016 c
1017 c    qi --> qit
1018                 vatrtt(nrofon,q6t,f2fcp) =
1019      >       trshu*(vatren(nrofon,q2t,prf1cn)+
1020      >              vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
1021      >        -trstr2*(vatren(nrofon,q1t,prf1cn)+
1022      >                 vatren(nrofon,q3t,prf1cn))
1023      >        +unssz*vatren(nrofon,q6t,prf1cn)
1024 c
1025             else
1026 c
1027                 codret = codret + 1
1028 c
1029             endif
1030 c
1031 c        pour le quadrangle Q4
1032 c
1033             if ( mod(hetqua(f1hp+3),100).eq.0 ) then
1034
1035            prfcap(f4cp) = 1
1036 c
1037            vafott(nrofon,q1,f4cp) = vatren(nrofon,q3t,prf2cn)
1038 c
1039            vafott(nrofon,q2,f4cp) = vatren(nrofon,q6t,prf2cn)
1040 c
1041            vafott(nrofon,q3,f4cp) =
1042      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
1043      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
1044      >     unsqu*vatren(nrofon,q6t,prf1cn)
1045 c
1046            vafott(nrofon,q4,f4cp) =
1047      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
1048      >                          vatren(nrofon,q2t,prf2cn)+
1049      >                          vatren(nrofon,q3t,prf3cn))
1050 c
1051            vafott(nrofon,q5,f4cp) =
1052      >    -unshu*vatren(nrofon,q1t,prf2cn)+
1053      >     trshu*vatren(nrofon,q3t,prf2cn)+
1054      >     trsqu*vatren(nrofon,q6t,prf2cn)
1055 c
1056            vafott(nrofon,q6,f4cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
1057      >                             vatren(nrofon,q4t,prf2cn))
1058 c
1059            vafott(nrofon,q7,f4cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
1060      >                             vatren(nrofon,q4t,prf1cn)+
1061      >                             vatren(nrofon,q5t,prf1cn))-
1062      >                    trstr2*(vatren(nrofon,q1t,prf1cn)+
1063      >                             vatren(nrofon,q3t,prf1cn))+
1064      >                    unssz*vatren(nrofon,q6t,prf1cn)
1065 c
1066            vafott(nrofon,q8,f4cp) = vatren(nrofon,q5t,prf2cn)
1067 c
1068            elseif ( mod(hetqua(f1hp+3),100).ge.31 .and.
1069      >              mod(hetqua(f1hp+3),100).le.34 ) then
1070 c
1071                 f1fhp = -filqua(f1hp+3)
1072                 f3fcp = ntrsca(f1fhp+2)
1073                 f1fcp = ntrsca(f1fhp)
1074                 f2fcp = ntrsca(f1fhp+1)
1075
1076                 prftrp(f3fcp) = 1
1077                 prftrp(f1fcp) = 1
1078                 prftrp(f2fcp) = 1
1079 c
1080 c
1081 c           Pour le fils f3fcp
1082 c
1083 c    qi --> qit
1084             vatrtt(nrofon,q1t,f3fcp) =
1085      >      -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))
1086      >      +unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
1087      >      +unsqu*vatren(nrofon,q6t,prf1cn)
1088
1089 c
1090             vatrtt(nrofon,q2t,f3fcp) =
1091      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
1092      >                          vatren(nrofon,q2t,prf2cn)+
1093      >                          vatren(nrofon,q3t,prf3cn))
1094 c
1095 c    qi --> qit
1096             vatrtt(nrofon,q3t,f3fcp) = vatren(nrofon,q5t,prf2cn)
1097 c
1098 c    qi --> qit
1099             vatrtt(nrofon,q4t,f3fcp) =
1100      >       trshu*(vatren(nrofon,q2t,prf1cn)+
1101      >              vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
1102      >        -trstr2*(vatren(nrofon,q1t,prf1cn)+
1103      >                 vatren(nrofon,q3t,prf1cn))
1104      >        +unssz*vatren(nrofon,q6t,prf1cn)
1105
1106 c
1107             vatrtt(nrofon,q5t,f3fcp) =
1108      >     -unshu*vatren(nrofon,q3t,prf2cn)+
1109      >      trshu*vatren(nrofon,q2t,prf2cn)+
1110      >      trsqu*vatren(nrofon,q5t,prf2cn)
1111 c
1112             vatrtt(nrofon,q6t,f3fcp) =
1113      >     -unshu*unsde*(vatren(nrofon,q1t,prf2cn)+
1114      >                   vatren(nrofon,q3t,prf1cn))+
1115      >      trshu*unsde*(vatren(nrofon,q2t,prf2cn)+
1116      >                   vatren(nrofon,q2t,prf1cn))+
1117      >      trsqu*unsde*(vatren(nrofon,q5t,prf1cn)+
1118      >                   vatren(nrofon,q4t,prf2cn))
1119 c
1120 c           Pour le fils f1fcp
1121 c
1122             vatrtt(nrofon,q1t,f1fcp) = vatrtt(nrofon,q1t,f3fcp)
1123
1124 c
1125             vatrtt(nrofon,q2t,f1fcp) = vatren(nrofon,q5t,prf2cn)
1126 c
1127             vatrtt(nrofon,q3t,f1fcp) = vatren(nrofon,q6t,prf2cn)
1128 c
1129             vatrtt(nrofon,q4t,f1fcp) = vatrtt(nrofon,q6t,f3fcp)
1130 c
1131             vatrtt(nrofon,q5t,f1fcp) =
1132      >     -unshu*(vatren(nrofon,q2t,prf2cn)+vatren(nrofon,q1t,prf2cn))
1133      >     +unsde*(vatren(nrofon,q5t,prf2cn)+vatren(nrofon,q6t,prf2cn))
1134      >     +unsqu*vatren(nrofon,q4t,prf2cn)
1135 c
1136             vatrtt(nrofon,q6t,f1fcp) =
1137      >     unsde*(vatren(nrofon,q4t,prf2cn)+
1138      >            vatren(nrofon,q5t,prf1cn))
1139 c
1140 c           Pour le fils f2fcp
1141 c
1142             vatrtt(nrofon,q1t,f2fcp) = vatren(nrofon,q6t,prf2cn)
1143
1144 c
1145             vatrtt(nrofon,q2t,f2fcp) = vatren(nrofon,q5t,prf2cn)
1146 c
1147             vatrtt(nrofon,q3t,f2fcp) = vatren(nrofon,q3t,prf2cn)
1148 c
1149             vatrtt(nrofon,q4t,f2fcp) =
1150      >     -unshu*(vatren(nrofon,q2t,prf2cn)+vatren(nrofon,q1t,prf2cn))
1151      >     +unsde*(vatren(nrofon,q5t,prf2cn)+vatren(nrofon,q6t,prf2cn))
1152      >     +unsqu*vatren(nrofon,q4t,prf2cn)
1153
1154             vatrtt(nrofon,q5t,f2fcp) =
1155      >      trshu*vatren(nrofon,q3t,prf2cn)-
1156      >      unshu*vatren(nrofon,q2t,prf2cn)+
1157      >      trsqu*vatren(nrofon,q5t,prf2cn)
1158 c
1159             vatrtt(nrofon,q6t,f2fcp) =
1160      >      trshu*vatren(nrofon,q3t,prf2cn)-
1161      >      unshu*vatren(nrofon,q1t,prf2cn)+
1162      >      trsqu*vatren(nrofon,q6t,prf2cn)
1163 c
1164             else
1165 c
1166                 codret = codret + 1
1167 c
1168             endif
1169  2343    continue
1170 c
1171 c           decoupage en 4 quadrangles d'un quadrangle
1172 c           predecoupe en 3 triangles par l arete a4
1173 c
1174          elseif (etan .eq. 34) then
1175
1176          do 2344 , nrofon = 1 , nbfonc
1177 c
1178 c        pour le quadrangle Q1
1179 c
1180             if ( mod(hetqua(f1hp),100).eq.0 ) then
1181
1182             prfcap(f1cp) = 1
1183 c
1184            vafott(nrofon,q1,f1cp) = vatren(nrofon,q3t,prf2cn)
1185 c
1186            vafott(nrofon,q2,f1cp) = vatren(nrofon,q6t,prf2cn)
1187 c
1188            vafott(nrofon,q3,f1cp) =
1189      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
1190      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
1191      >     unsqu*vatren(nrofon,q6t,prf1cn)
1192 c
1193            vafott(nrofon,q4,f1cp) =
1194      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
1195      >                          vatren(nrofon,q2t,prf2cn)+
1196      >                          vatren(nrofon,q3t,prf3cn))
1197 c
1198            vafott(nrofon,q5,f1cp) =
1199      >    -unshu*vatren(nrofon,q1t,prf2cn)+
1200      >     trshu*vatren(nrofon,q3t,prf2cn)+
1201      >     trsqu*vatren(nrofon,q6t,prf2cn)
1202 c
1203            vafott(nrofon,q6,f1cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
1204      >                              vatren(nrofon,q4t,prf2cn))
1205 c
1206            vafott(nrofon,q7,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
1207      >                              vatren(nrofon,q4t,prf1cn)+
1208      >                              vatren(nrofon,q5t,prf1cn))-
1209      >                     trstr2*(vatren(nrofon,q1t,prf1cn)+
1210      >                              vatren(nrofon,q3t,prf1cn))+
1211      >                     unssz*vatren(nrofon,q6t,prf1cn)
1212 c
1213            vafott(nrofon,q8,f1cp) = vatren(nrofon,q5t,prf2cn)
1214 c
1215            elseif ( mod(hetqua(f1hp),100).ge.31 .and.
1216      >              mod(hetqua(f1hp),100).le.34 ) then
1217 c
1218
1219                 f1fhp = -filqua(f1hp)
1220                 f1fcp = ntrsca(f1fhp)
1221                 f2fcp = ntrsca(f1fhp+1)
1222                 f3fcp = ntrsca(f1fhp+2)
1223
1224                 prftrp(f3fcp) = 1
1225                 prftrp(f1fcp) = 1
1226                 prftrp(f2fcp) = 1
1227 c
1228 c
1229 c           Pour le fils f3fcp
1230 c
1231 c    qi --> qit
1232                 vatrtt(nrofon,q1t,f3fcp) = vatren(nrofon,q4t,prf3cn)
1233 c
1234 c    qi --> qit
1235                 vatrtt(nrofon,q2t,f3fcp) = vatren(nrofon,q2t,prf3cn)
1236 c
1237 c    qi --> qit
1238                 vatrtt(nrofon,q3t,f3fcp) = vatren(nrofon,q5t,prf3cn)
1239 c
1240 c    qi --> qit
1241                 vatrtt(nrofon,q4t,f3fcp) =
1242      >           trshu*vatren(nrofon,q2t,prf3cn)
1243      >          -unshu*vatren(nrofon,q1t,prf3cn)
1244      >          +trsqu*vatren(nrofon,q4t,prf3cn)
1245 c
1246 c    qi --> qit
1247                 vatrtt(nrofon,q5t,f3fcp) =
1248      >           trshu*vatren(nrofon,q2t,prf3cn)
1249      >          -unshu*vatren(nrofon,q3t,prf3cn)
1250      >          +trsqu*vatren(nrofon,q5t,prf3cn)
1251 c
1252 c    qi --> qit
1253                 vatrtt(nrofon,q6t,f3fcp) =
1254      >      -unshu*(vatren(nrofon,q1t,prf3cn)+vatren(nrofon,q3t,prf3cn))
1255      >      +unsde*(vatren(nrofon,q4t,prf3cn)+vatren(nrofon,q5t,prf3cn))
1256      >      +unsqu*vatren(nrofon,q6t,prf3cn)
1257 c
1258 c           Pour le fils f1fcp
1259 c
1260 c    qi --> qit
1261                 vatrtt(nrofon,q1t,f1fcp) = vatren(nrofon,q4t,prf3cn)
1262 c
1263 c    qi --> qit
1264                 vatrtt(nrofon,q2t,f1fcp) = vatren(nrofon,q5t,prf3cn)
1265 c
1266 c    qi --> qit
1267                 vatrtt(nrofon,q3t,f1fcp) =
1268      >      -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))
1269      >      +unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
1270      >      +unsqu*vatren(nrofon,q6t,prf1cn)
1271 c
1272 c    qi --> qit
1273                 vatrtt(nrofon,q4t,f1fcp) = vatrtt(nrofon,q6t,f3fcp)
1274 c
1275 c    qi --> qit
1276                 vatrtt(nrofon,q5t,f1fcp) =
1277      >        -unshu*unsde*(vatren(nrofon,q1t,prf3cn)+
1278      >                      vatren(nrofon,q1t,prf1cn))
1279      >        +trshu*unsde*(vatren(nrofon,q3t,prf3cn)+
1280      >                      vatren(nrofon,q2t,prf1cn))
1281      >        +trsqu*unsde*(vatren(nrofon,q4t,prf1cn)+
1282      >                      vatren(nrofon,q6t,prf3cn))
1283 c
1284 c    qi --> qit
1285                 vatrtt(nrofon,q6t,f1fcp) =
1286      >   unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q6t,prf3cn))
1287 c
1288 c           Pour le fils f2fcp
1289 c
1290 c    qi --> qit
1291                 vatrtt(nrofon,q1t,f2fcp) = vatrtt(nrofon,q3t,f1fcp)
1292 c
1293 c    qi --> qit
1294                 vatrtt(nrofon,q2t,f2fcp) = vatrtt(nrofon,q2t,f1fcp)
1295 c
1296 c    qi --> qit
1297                 vatrtt(nrofon,q3t,f2fcp) =
1298      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
1299      >                          vatren(nrofon,q2t,prf2cn)+
1300      >                          vatren(nrofon,q3t,prf3cn))
1301 c
1302 c    qi --> qit
1303                 vatrtt(nrofon,q4t,f2fcp) = vatrtt(nrofon,q5t,f1fcp)
1304 c
1305 c    qi --> qit
1306                 vatrtt(nrofon,q5t,f2fcp) =
1307      >         trshu*vatren(nrofon,q3t,prf3cn)-
1308      >         unshu*vatren(nrofon,q2t,prf3cn)+
1309      >         trsqu*vatren(nrofon,q5t,prf3cn)
1310 c
1311 c    qi --> qit
1312                 vatrtt(nrofon,q6t,f2fcp) =
1313      >       trshu*(vatren(nrofon,q2t,prf1cn)+
1314      >              vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
1315      >        -trstr2*(vatren(nrofon,q1t,prf1cn)+
1316      >                 vatren(nrofon,q3t,prf1cn))
1317      >        +unssz*vatren(nrofon,q6t,prf1cn)
1318 c
1319 c
1320             else
1321 c
1322                 codret = codret + 1
1323 c
1324             endif
1325 c
1326 c        pour le quadrangle Q2
1327 c
1328            prfcap(f2cp) = 1
1329 c
1330            vafott(nrofon,q1,f2cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
1331      >                              vatren(nrofon,q1t,prf2cn))
1332 c
1333            vafott(nrofon,q2,f2cp) = vatren(nrofon,q6t,prf1cn)
1334 c
1335            vafott(nrofon,q3,f2cp) =
1336      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
1337      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
1338      >     unsqu*vatren(nrofon,q6t,prf1cn)
1339 c
1340            vafott(nrofon,q4,f2cp) = vatren(nrofon,q6t,prf2cn)
1341 c
1342            vafott(nrofon,q5,f2cp) =
1343      >    -unshu*vatren(nrofon,q1t,prf1cn)+
1344      >     trshu*vatren(nrofon,q3t,prf1cn)+
1345      >     trsqu*vatren(nrofon,q6t,prf1cn)
1346 c
1347            vafott(nrofon,q6,f2cp) =  -unshu*vatren(nrofon,q2t,prf1cn)-
1348      >                     trstr2*(vatren(nrofon,q3t,prf1cn)+
1349      >                              vatren(nrofon,q1t,prf1cn))+
1350      >                         trshu*(vatren(nrofon,q4t,prf1cn)+
1351      >                              vatren(nrofon,q5t,prf1cn))+
1352      >                      nessz*vatren(nrofon,q6t,prf1cn)
1353 c
1354            vafott(nrofon,q7,f2cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
1355      >                              vatren(nrofon,q4t,prf2cn))
1356 c
1357            vafott(nrofon,q8,f2cp) =
1358      >    -unshu*vatren(nrofon,q3t,prf2cn)+
1359      >     trshu*vatren(nrofon,q1t,prf2cn)+
1360      >     trsqu*vatren(nrofon,q6t,prf2cn)
1361 c
1362 c        pour le quadrangle Q3
1363 c
1364            prfcap(f3cp) = 1
1365 c
1366            vafott(nrofon,q1,f3cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
1367      >                              vatren(nrofon,q1t,prf3cn))
1368 c
1369            vafott(nrofon,q2,f3cp) =  vatren(nrofon,q4t,prf3cn)
1370 c
1371            vafott(nrofon,q3,f3cp) =
1372      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
1373      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
1374      >     unsqu*vatren(nrofon,q6t,prf1cn)
1375 c
1376            vafott(nrofon,q4,f3cp) = vatren(nrofon,q6t,prf1cn)
1377 c
1378            vafott(nrofon,q5,f3cp) =
1379      >    -unshu*vatren(nrofon,q2t,prf3cn)+
1380      >     trshu*vatren(nrofon,q1t,prf3cn)+
1381      >     trsqu*vatren(nrofon,q4t,prf3cn)
1382 c
1383            vafott(nrofon,q6,f3cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
1384      >                              vatren(nrofon,q6t,prf3cn))
1385 c
1386            vafott(nrofon,q7,f3cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
1387      >                     trstr2*(vatren(nrofon,q3t,prf1cn)+
1388      >                              vatren(nrofon,q1t,prf1cn))+
1389      >                         trshu*(vatren(nrofon,q4t,prf1cn)+
1390      >                              vatren(nrofon,q5t,prf1cn))+
1391      >                      nessz*vatren(nrofon,q6t,prf1cn)
1392 c
1393            vafott(nrofon,q8,f3cp) =
1394      >    -unshu*vatren(nrofon,q3t,prf1cn)+
1395      >     trshu*vatren(nrofon,q1t,prf1cn)+
1396      >     trsqu*vatren(nrofon,q6t,prf1cn)
1397 c
1398 c        pour le quadrangle Q4
1399 c
1400             if ( mod(hetqua(f1hp+3),100).eq.0 ) then
1401
1402             prfcap(f4cp) = 1
1403 c
1404            vafott(nrofon,q1,f4cp) = vatren(nrofon,q3t,prf2cn)
1405 c
1406            vafott(nrofon,q2,f4cp) =
1407      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
1408      >                          vatren(nrofon,q2t,prf2cn)+
1409      >                          vatren(nrofon,q3t,prf3cn))
1410 c
1411            vafott(nrofon,q3,f4cp) =
1412      >    -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))+
1413      >     unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))+
1414      >     unsqu*vatren(nrofon,q6t,prf1cn)
1415 c
1416            vafott(nrofon,q4,f4cp) = vatren(nrofon,q4t,prf3cn)
1417 c
1418            vafott(nrofon,q5,f4cp) = vatren(nrofon,q5t,prf3cn)
1419 c
1420            vafott(nrofon,q6,f4cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
1421      >                              vatren(nrofon,q4t,prf1cn)+
1422      >                              vatren(nrofon,q5t,prf1cn))-
1423      >                     trstr2*(vatren(nrofon,q1t,prf1cn)+
1424      >                              vatren(nrofon,q3t,prf1cn))+
1425      >                     unssz*vatren(nrofon,q6t,prf1cn)
1426 c
1427            vafott(nrofon,q7,f4cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
1428      >                              vatren(nrofon,q6t,prf3cn))
1429 c
1430            vafott(nrofon,q8,f4cp) =
1431      >    -unshu*vatren(nrofon,q1t,prf3cn)+
1432      >     trshu*vatren(nrofon,q2t,prf3cn)+
1433      >     trsqu*vatren(nrofon,q4t,prf3cn)
1434 c
1435 c
1436            elseif ( mod(hetqua(f1hp+3),100).ge.31 .and.
1437      >              mod(hetqua(f1hp+3),100).le.34 ) then
1438 c
1439                 f1fhp = -filqua(f1hp+3)
1440                 f1fcp = ntrsca(f1fhp)
1441                 f2fcp = ntrsca(f1fhp+1)
1442                 f3fcp = ntrsca(f1fhp+2)
1443 c
1444                 prftrp(f3fcp) = 1
1445                 prftrp(f1fcp) = 1
1446                 prftrp(f2fcp) = 1
1447 c
1448 c
1449 c           Pour le fils f3fcp
1450 c
1451 c    qi --> qit
1452             vatrtt(nrofon,q1t,f3fcp) =
1453      >      -unshu*(vatren(nrofon,q1t,prf1cn)+vatren(nrofon,q3t,prf1cn))
1454      >      +unsde*(vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
1455      >      +unsqu*vatren(nrofon,q6t,prf1cn)
1456
1457 c
1458             vatrtt(nrofon,q2t,f3fcp) =
1459      >                 unstr *( vatren(nrofon,q2t,prf1cn)+
1460      >                          vatren(nrofon,q2t,prf2cn)+
1461      >                          vatren(nrofon,q3t,prf3cn))
1462 c
1463 c    qi --> qit
1464             vatrtt(nrofon,q3t,f3fcp) = vatrtt(nrofon,q5t,prf2cn)
1465 c
1466 c    qi --> qit
1467             vatrtt(nrofon,q4t,f3fcp) =
1468      >       trshu*(vatren(nrofon,q2t,prf1cn)+
1469      >              vatren(nrofon,q4t,prf1cn)+vatren(nrofon,q5t,prf1cn))
1470      >        -trstr2*(vatren(nrofon,q1t,prf1cn)+
1471      >                 vatren(nrofon,q3t,prf1cn))
1472      >        +unssz*vatren(nrofon,q6t,prf1cn)
1473
1474 c
1475             vatrtt(nrofon,q5t,f3fcp) =
1476      >     -unshu*vatren(nrofon,q3t,prf2cn)+
1477      >      trshu*vatren(nrofon,q2t,prf2cn)+
1478      >      trsqu*vatren(nrofon,q5t,prf2cn)
1479 c
1480             vatrtt(nrofon,q6t,f3fcp) =
1481      >     -unshu*unsde*(vatren(nrofon,q1t,prf2cn)+
1482      >                   vatren(nrofon,q3t,prf1cn))+
1483      >      trshu*unsde*(vatren(nrofon,q2t,prf2cn)+
1484      >                   vatren(nrofon,q2t,prf1cn))+
1485      >      trsqu*unsde*(vatren(nrofon,q5t,prf1cn)+
1486      >                   vatren(nrofon,q4t,prf2cn))
1487 c
1488 c           Pour le fils f1fcp
1489 c
1490             vatrtt(nrofon,q1t,f1fcp) = vatrtt(nrofon,q1t,f3fcp)
1491
1492 c
1493             vatrtt(nrofon,q2t,f1fcp) = vatren(nrofon,q5t,prf2cn)
1494 c
1495             vatrtt(nrofon,q3t,f1fcp) = vatren(nrofon,q6t,prf2cn)
1496 c
1497             vatrtt(nrofon,q4t,f1fcp) = vatrtt(nrofon,q6t,f3fcp)
1498 c
1499             vatrtt(nrofon,q5t,f1fcp) =
1500      >     -unshu*(vatren(nrofon,q2t,prf2cn)+vatren(nrofon,q1t,prf2cn))
1501      >     +unsde*(vatren(nrofon,q5t,prf2cn)+vatren(nrofon,q6t,prf2cn))
1502      >     +unsqu*vatren(nrofon,q4t,prf2cn)
1503 c
1504             vatrtt(nrofon,q6t,f1fcp) =
1505      >     unsde*(vatren(nrofon,q4t,prf2cn)+
1506      >            vatren(nrofon,q5t,prf1cn))
1507 c
1508 c           Pour le fils f2fcp
1509 c
1510             vatrtt(nrofon,q1t,f2fcp) = vatren(nrofon,q6t,prf2cn)
1511
1512 c
1513             vatrtt(nrofon,q2t,f2fcp) = vatren(nrofon,q5t,prf2cn)
1514 c
1515             vatrtt(nrofon,q3t,f2fcp) = vatren(nrofon,q3t,prf2cn)
1516 c
1517             vatrtt(nrofon,q4t,f2fcp) =
1518      >     -unshu*(vatren(nrofon,q2t,prf2cn)+vatren(nrofon,q1t,prf2cn))
1519      >     +unsde*(vatren(nrofon,q5t,prf2cn)+vatren(nrofon,q6t,prf2cn))
1520      >     +unsqu*vatren(nrofon,q4t,prf2cn)
1521 c
1522             vatrtt(nrofon,q5t,f2fcp) =
1523      >      trshu*vatren(nrofon,q3t,prf2cn)-
1524      >      unshu*vatren(nrofon,q2t,prf2cn)+
1525      >      trsqu*vatren(nrofon,q5t,prf2cn)
1526 c
1527             vatrtt(nrofon,q6t,f2fcp) =
1528      >      trshu*vatren(nrofon,q3t,prf2cn)-
1529      >      unshu*vatren(nrofon,q1t,prf2cn)+
1530      >      trsqu*vatren(nrofon,q6t,prf2cn)
1531 c
1532             else
1533 c
1534                 codret = codret + 1
1535 c
1536             endif
1537  2344     continue
1538 c
1539         endif