src/tool/AP_Conversion/pcsqu2_4.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 = ntrsca(f1hp)
  13             f2cp = ntrsca(f1hp+1)
  14             f3cp = ntrsca(f1hp+2)
  15             prftrp(f1cp) = 1
  16             prftrp(f2cp) = 1
  17             prftrp(f3cp) = 1
  18 c
  19             q1tp = 1
  20             q2tp = 2
  21             q3tp = 3
  22             q4tp = 4
  23             q5tp = 5
  24             q6tp = 6
  25 c
  26 c   decoupage en 3 par l'arete a2 d'un triangle predecoupe en a1
  27 c
  28             if ((etan.eq.31).and.( etanp1.eq.32)) then
  29 c
  30             do 23311 , nrofon = 1 , nbfonc
  31
  32 c   dans le triangle NT1
  33
  34           vatrtt(nrofon,q1tp,f1cp) = vatren(nrofon,q2t,prf3cn)
  35 c
  36           vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q6t,prf2cn)
  37 c
  38           vatrtt(nrofon,q3tp,f1cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
  39      >                              vatren(nrofon,q1t,prf3cn))
  40 c
  41           vatrtt(nrofon,q4tp,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
  42      >                              vatren(nrofon,q4t,prf1cn)+
  43      >                              vatren(nrofon,q5t,prf1cn))-
  44      >                    trstr2*(vatren(nrofon,q1t,prf1cn)+
  45      >                              vatren(nrofon,q3t,prf1cn))+
  46      >                    unssz*vatren(nrofon,q6t,prf1cn)
  47 c
  48           vatrtt(nrofon,q5tp,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
  49      >                     trstr2*(vatren(nrofon,q3t,prf1cn)+
  50      >                              vatren(nrofon,q1t,prf1cn))+
  51      >                         trshu*(vatren(nrofon,q4t,prf1cn)+
  52      >                              vatren(nrofon,q5t,prf1cn))+
  53      >                     nessz*vatren(nrofon,q6t,prf1cn)
  54 c
  55           vatrtt(nrofon,q6tp,f1cp) =vatren(nrofon,q4t,prf3cn)
  56 c
  57 c   dans le triangle NT2
  58 c
  59           vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q1t,prf1cn)
  60 c
  61           vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q6t,prf2cn)
  62 c
  63           vatrtt(nrofon,q3tp,f2cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
  64      >                              vatren(nrofon,q1t,prf2cn))
  65 c
  66           vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
  67 c
  68           vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q3t,prf2cn)+
  69      >                          trshu*vatren(nrofon,q1t,prf2cn)+
  70      >                          trsqu*vatren(nrofon,q6t,prf2cn)
  71 c
  72           vatrtt(nrofon,q6tp,f2cp) =vatren(nrofon,q6t,prf1cn)
  73 c
  74 c   dans le triangle NT3
  75 c
  76           vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q2t,prf3cn)
  77 c
  78           vatrtt(nrofon,q2tp,f3cp) = vatren(nrofon,q3t,prf2cn)
  79 c
  80           vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q6t,prf2cn)
  81 c
  82           vatrtt(nrofon,q4tp,f3cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
  83      >                              vatren(nrofon,q2t,prf2cn)+
  84      >                              vatren(nrofon,q3t,prf3cn))
  85 c
  86           vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q1t,prf2cn)+
  87      >                          trshu*vatren(nrofon,q3t,prf2cn)+
  88      >                          trsqu*vatren(nrofon,q6t,prf2cn)
  89 c
  90           vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
  91 23311       continue
  92 c
  93 c   decoupage en 3 par l'arete a3 d'un triangle predecoupe en a1
  94 c
  95             elseif ((etan.eq.31).and.( etanp1.eq.33)) then
  96 c
  97             do 23312 , nrofon = 1 , nbfonc
  98 c
  99 c   dans le triangle NT1
 100 c
 101            vatrtt(nrofon,q1tp,f1cp) = vatren(nrofon,q3t,prf2cn)
 102 c
 103            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q6t,prf1cn)
 104 c
 105            vatrtt(nrofon,q3tp,f1cp) = vatren(nrofon,q2t,prf3cn)
 106 c
 107            vatrtt(nrofon,q4tp,f1cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
 108      >                                vatren(nrofon,q4t,prf2cn))
 109 c
 110            vatrtt(nrofon,q5tp,f1cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 111      >                                vatren(nrofon,q6t,prf3cn))
 112 c
 113            vatrtt(nrofon,q6tp,f1cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 114      >                                vatren(nrofon,q2t,prf2cn)+
 115      >                                vatren(nrofon,q3t,prf3cn))
 116
 117 c   dans le triangle NT2
 118 c
 119            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q2t,prf3cn)
 120 c
 121            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q6t,prf1cn)
 122 c
 123            vatrtt(nrofon,q3tp,f2cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 124      >                                vatren(nrofon,q1t,prf3cn))
 125 c
 126            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 127 c
 128            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q3t,prf1cn)+
 129      >                            trshu*vatren(nrofon,q1t,prf1cn)+
 130      >                            trsqu*vatren(nrofon,q6t,prf1cn)
 131 c
 132            vatrtt(nrofon,q6tp,f2cp) =vatren(nrofon,q4t,prf3cn)
 133 c
 134 c   dans le triangle NT3
 135 c
 136            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q3t,prf2cn)
 137 c
 138            vatrtt(nrofon,q2tp,f3cp) = unsde*(vatren(nrofon,q1t,prf2cn)+
 139      >                                vatren(nrofon,q3t,prf1cn))
 140 c
 141            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q6t,prf1cn)
 142 c
 143            vatrtt(nrofon,q4tp,f3cp) = vatren(nrofon,q6t,prf2cn)
 144 c
 145            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q1t,prf1cn)+
 146      >                            trshu*vatren(nrofon,q3t,prf1cn)+
 147      >                            trsqu*vatren(nrofon,q6t,prf1cn)
 148 c
 149            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 150 c
 151 23312       continue
 152 c
 153 c   decoupage en 3 par l'arete a4 d'un triangle predecoupe en a1
 154 c
 155             elseif ((etan.eq.31).and.( etanp1.eq.34)) then
 156 c
 157             do 23313 , nrofon = 1 , nbfonc
 158 c
 159 c   dans le triangle NT1
 160 c
 161            vatrtt(nrofon,q1tp,f1cp) = unsde*(vatren(nrofon,q1t,prf2cn)+
 162      >                                vatren(nrofon,q3t,prf1cn))
 163 c
 164            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q4t,prf3cn)
 165 c
 166            vatrtt(nrofon,q3tp,f1cp) = vatren(nrofon,q3t,prf2cn)
 167 c
 168            vatrtt(nrofon,q4tp,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 169      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 170      >                                vatren(nrofon,q1t,prf1cn))+
 171      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 172      >                                vatren(nrofon,q5t,prf1cn))+
 173      >                       nessz*vatren(nrofon,q6t,prf1cn)
 174 c
 175            vatrtt(nrofon,q5tp,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 176      >                                vatren(nrofon,q4t,prf1cn)+
 177      >                                vatren(nrofon,q5t,prf1cn))-
 178      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
 179      >                                vatren(nrofon,q3t,prf1cn))+
 180      >                      unssz*vatren(nrofon,q6t,prf1cn)
 181 c
 182            vatrtt(nrofon,q6tp,f1cp) = vatren(nrofon,q6t,prf2cn)
 183 c
 184 c   dans le triangle NT2
 185 c
 186            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q3t,prf2cn)
 187 c
 188            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q4t,prf3cn)
 189 c
 190            vatrtt(nrofon,q3tp,f2cp) = vatren(nrofon,q2t,prf3cn)
 191 c
 192            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 193 c
 194            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q1t,prf3cn)+
 195      >                            trshu*vatren(nrofon,q2t,prf3cn)+
 196      >                            trsqu*vatren(nrofon,q4t,prf3cn)
 197 c
 198            vatrtt(nrofon,q6tp,f2cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 199      >                                vatren(nrofon,q2t,prf2cn)+
 200      >                                vatren(nrofon,q3t,prf3cn))
 201 c
 202 c   dans le triangle NT3
 203 c
 204            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q3t,prf1cn)
 205 c
 206            vatrtt(nrofon,q2tp,f3cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 207      >                                vatren(nrofon,q1t,prf3cn))
 208 c
 209            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q4t,prf3cn)
 210 c
 211            vatrtt(nrofon,q4tp,f3cp) = vatren(nrofon,q6t,prf1cn)
 212 c
 213            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q2t,prf3cn)+
 214      >                            trshu*vatren(nrofon,q1t,prf3cn)+
 215      >                            trsqu*vatren(nrofon,q4t,prf3cn)
 216 c
 217            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 218 c
 219 23313       continue
 220 c
 221 c   decoupage en 3 par l'arete a3 d'un triangle predecoupe en a2
 222 c
 223             elseif ((etan.eq.32).and.( etanp1.eq.33)) then
 224 c
 225             do 23321 , nrofon = 1 , nbfonc
 226 c
 227 c   dans le triangle NT1
 228 c
 229            vatrtt(nrofon,q1tp,f1cp) = vatren(nrofon,q2t,prf3cn)
 230 c
 231            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q6t,prf2cn)
 232 c
 233            vatrtt(nrofon,q3tp,f1cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 234      >                                vatren(nrofon,q1t,prf3cn))
 235 c
 236            vatrtt(nrofon,q4tp,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 237      >                                vatren(nrofon,q4t,prf1cn)+
 238      >                                vatren(nrofon,q5t,prf1cn))-
 239      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
 240      >                                vatren(nrofon,q3t,prf1cn))+
 241      >                      unssz*vatren(nrofon,q6t,prf1cn)
 242 c
 243            vatrtt(nrofon,q5tp,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 244      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 245      >                                vatren(nrofon,q1t,prf1cn))+
 246      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 247      >                                vatren(nrofon,q5t,prf1cn))+
 248      >                       nessz*vatren(nrofon,q6t,prf1cn)
 249 c
 250            vatrtt(nrofon,q6tp,f1cp) = vatren(nrofon,q4t,prf3cn)
 251 c
 252 c   dans le triangle NT2
 253 c
 254            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q1t,prf1cn)
 255 c
 256            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q6t,prf2cn)
 257 c
 258            vatrtt(nrofon,q3tp,f2cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 259      >                                vatren(nrofon,q1t,prf2cn))
 260 c
 261            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 262 c
 263            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q3t,prf2cn)+
 264      >                            trshu*vatren(nrofon,q1t,prf2cn)+
 265      >                            trsqu*vatren(nrofon,q6t,prf2cn)
 266 c
 267            vatrtt(nrofon,q6tp,f2cp) = vatren(nrofon,q4t,prf3cn)
 268 c
 269 c   dans le triangle NT3
 270 c
 271            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q2t,prf3cn)
 272 c
 273            vatrtt(nrofon,q2tp,f3cp) = vatren(nrofon,q3t,prf2cn)
 274 c
 275            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q6t,prf2cn)
 276 c
 277            vatrtt(nrofon,q4tp,f3cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 278      >                                vatren(nrofon,q2t,prf2cn)+
 279      >                                vatren(nrofon,q3t,prf3cn))
 280 c
 281            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q1t,prf2cn)+
 282      >                            trshu*vatren(nrofon,q3t,prf2cn)+
 283      >                            trsqu*vatren(nrofon,q6t,prf2cn)
 284 c
 285            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 286 23321       continue
 287 c
 288 c   decoupage en 3 par l'arete a4 d'un triangle predecoupe en a2
 289 c
 290             elseif ((etan.eq.32).and.( etanp1.eq.34)) then
 291 c
 292             do 23322 , nrofon = 1 , nbfonc
 293 c
 294 c   dans le triangle NT1
 295 c
 296            vatrtt(nrofon,q1tp,f1cp) = vatren(nrofon,q3t,prf2cn)
 297 c
 298            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q6t,prf1cn)
 299 c
 300            vatrtt(nrofon,q3tp,f1cp) = vatren(nrofon,q2t,prf3cn)
 301 c
 302            vatrtt(nrofon,q4tp,f1cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
 303      >                                vatren(nrofon,q4t,prf2cn))
 304 c
 305            vatrtt(nrofon,q5tp,f1cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 306      >                                vatren(nrofon,q6t,prf3cn))
 307 c
 308            vatrtt(nrofon,q6tp,f1cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 309      >                                vatren(nrofon,q2t,prf2cn)+
 310      >                                vatren(nrofon,q3t,prf3cn))
 311 c
 312 c   dans le triangle NT2
 313 c
 314            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q2t,prf3cn)
 315 c
 316            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q6t,prf1cn)
 317 c
 318            vatrtt(nrofon,q3tp,f2cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 319      >                                vatren(nrofon,q1t,prf3cn))
 320 c
 321            vatrtt(nrofon,q4tp,f2cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 322      >                                vatren(nrofon,q6t,prf3cn))
 323 c
 324            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q3t,prf1cn)+
 325      >                            trshu*vatren(nrofon,q1t,prf1cn)+
 326      >                            trsqu*vatren(nrofon,q6t,prf1cn)
 327 c
 328            vatrtt(nrofon,q6tp,f2cp) = vatren(nrofon,q4t,prf3cn)
 329 c
 330 c   dans le triangle NT3
 331 c
 332            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q3t,prf2cn)
 333 c
 334            vatrtt(nrofon,q2tp,f3cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 335      >                                vatren(nrofon,q1t,prf2cn))
 336 c
 337            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q6t,prf1cn)
 338 c
 339            vatrtt(nrofon,q4tp,f3cp) = vatren(nrofon,q6t,prf2cn)
 340 c
 341            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q1t,prf1cn)+
 342      >                            trshu*vatren(nrofon,q3t,prf1cn)+
 343      >                            trsqu*vatren(nrofon,q6t,prf1cn)
 344 c
 345            vatrtt(nrofon,q6tp,f3cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
 346      >                                vatren(nrofon,q4t,prf2cn))
 347 23322       continue
 348 c
 349 c   decoupage en 3 par l'arete a1 d'un triangle predecoupe en a2
 350 c
 351             elseif ((etan.eq.32).and.( etanp1.eq.31)) then
 352 c
 353             do 23323 , nrofon = 1 , nbfonc
 354 c
 355 c   dans le triangle NT1
 356 c
 357            vatrtt(nrofon,q1tp,f1cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 358      >                                vatren(nrofon,q1t,prf2cn))
 359 c
 360            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q4t,prf3cn)
 361 c
 362            vatrtt(nrofon,q3tp,f1cp) = vatren(nrofon,q3t,prf2cn)
 363 c
 364            vatrtt(nrofon,q4tp,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 365      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 366      >                                vatren(nrofon,q1t,prf1cn))+
 367      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 368      >                                vatren(nrofon,q5t,prf1cn))+
 369      >                       nessz*vatren(nrofon,q6t,prf1cn)
 370 c
 371            vatrtt(nrofon,q5tp,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 372      >                                vatren(nrofon,q4t,prf1cn)+
 373      >                                vatren(nrofon,q5t,prf1cn))-
 374      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
 375      >                                vatren(nrofon,q3t,prf1cn))+
 376      >                      unssz*vatren(nrofon,q6t,prf1cn)
 377 c
 378            vatrtt(nrofon,q6tp,f1cp) = vatren(nrofon,q6t,prf2cn)
 379 c
 380 c   dans le triangle NT2
 381 c
 382            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q3t,prf2cn)
 383 c
 384            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q4t,prf3cn)
 385 c
 386            vatrtt(nrofon,q3tp,f2cp) = vatren(nrofon,q2t,prf3cn)
 387 c
 388            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 389 c
 390            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q1t,prf3cn)+
 391      >                            trshu*vatren(nrofon,q2t,prf3cn)+
 392      >                            trsqu*vatren(nrofon,q4t,prf3cn)
 393 c
 394            vatrtt(nrofon,q6tp,f2cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 395      >                                vatren(nrofon,q2t,prf2cn)+
 396      >                                vatren(nrofon,q3t,prf3cn))
 397 c
 398 c   dans le triangle NT3
 399 c
 400            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q3t,prf1cn)
 401 c
 402            vatrtt(nrofon,q2tp,f3cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 403      >                                vatren(nrofon,q1t,prf3cn))
 404 c
 405            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q4t,prf3cn)
 406 c
 407            vatrtt(nrofon,q4tp,f3cp) = vatren(nrofon,q6t,prf1cn)
 408 c
 409            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q2t,prf3cn)+
 410      >                            trshu*vatren(nrofon,q1t,prf3cn)+
 411      >                            trsqu*vatren(nrofon,q4t,prf3cn)
 412 c
 413            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 414 23323       continue
 415 c
 416 c   decoupage en 3 par l'arete a4 d'un triangle predecoupe en a3
 417 c
 418             elseif ((etan.eq.33).and.( etanp1.eq.34)) then
 419 c
 420             do 23331 , nrofon = 1 , nbfonc
 421 c
 422 c   dans le triangle NT1
 423 c
 424            vatrtt(nrofon,q1tp,f1cp) = vatren(nrofon,q2t,prf3cn)
 425 c
 426            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q6t,prf2cn)
 427 c
 428            vatrtt(nrofon,q3tp,f1cp) = unsde*(vatren(nrofon,q1t,prf3cn)+
 429      >                                vatren(nrofon,q1t,prf1cn))
 430 c
 431            vatrtt(nrofon,q4tp,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 432      >                                vatren(nrofon,q4t,prf1cn)+
 433      >                                vatren(nrofon,q5t,prf1cn))-
 434      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
 435      >                                vatren(nrofon,q3t,prf1cn))+
 436      >                      unssz*vatren(nrofon,q6t,prf1cn)
 437 c
 438            vatrtt(nrofon,q5tp,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 439      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 440      >                                vatren(nrofon,q1t,prf1cn))+
 441      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 442      >                                vatren(nrofon,q5t,prf1cn))+
 443      >                       nessz*vatren(nrofon,q6t,prf1cn)
 444 c
 445            vatrtt(nrofon,q6tp,f1cp) = vatren(nrofon,q4t,prf3cn)
 446 c
 447 c   dans le triangle NT2
 448 c
 449            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q1t,prf1cn)
 450 c
 451            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q6t,prf2cn)
 452 c
 453            vatrtt(nrofon,q3tp,f2cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 454      >                                vatren(nrofon,q1t,prf2cn))
 455 c
 456            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 457 c
 458            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q3t,prf2cn)+
 459      >                            trshu*vatren(nrofon,q1t,prf2cn)+
 460      >                            trsqu*vatren(nrofon,q6t,prf2cn)
 461 c
 462            vatrtt(nrofon,q6tp,f2cp) = vatren(nrofon,q6t,prf1cn)
 463 c
 464 c   dans le triangle NT3
 465 c
 466            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q2t,prf3cn)
 467 c
 468            vatrtt(nrofon,q2tp,f3cp) = vatren(nrofon,q3t,prf2cn)
 469 c
 470            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q6t,prf2cn)
 471 c
 472            vatrtt(nrofon,q4tp,f3cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 473      >                                vatren(nrofon,q2t,prf2cn)+
 474      >                                vatren(nrofon,q3t,prf3cn))
 475 c
 476            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q1t,prf2cn)+
 477      >                            trshu*vatren(nrofon,q3t,prf2cn)+
 478      >                            trsqu*vatren(nrofon,q6t,prf2cn)
 479 c
 480            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 481 23331       continue
 482 c
 483 c   decoupage en 3 par l'arete a1 d'un triangle predecoupe en a3
 484 c
 485             elseif ((etan.eq.33).and.( etanp1.eq.31)) then
 486 c
 487             do 23332 , nrofon = 1 , nbfonc
 488 c
 489 c   dans le triangle NT1
 490 c
 491            vatrtt(nrofon,q1tp,f1cp) = vatren(nrofon,q3t,prf2cn)
 492 c
 493            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q6t,prf1cn)
 494 c
 495            vatrtt(nrofon,q3tp,f1cp) = vatren(nrofon,q2t,prf3cn)
 496 c
 497            vatrtt(nrofon,q4tp,f1cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
 498      >                                vatren(nrofon,q4t,prf2cn))
 499 c
 500            vatrtt(nrofon,q5tp,f1cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 501      >                                vatren(nrofon,q6t,prf3cn))
 502 c
 503            vatrtt(nrofon,q6tp,f1cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 504      >                                vatren(nrofon,q2t,prf2cn)+
 505      >                                vatren(nrofon,q3t,prf3cn))
 506 c
 507 c   dans le triangle NT2
 508 c
 509            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q2t,prf3cn)
 510 c
 511            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q6t,prf1cn)
 512 c
 513            vatrtt(nrofon,q3tp,f2cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 514      >                                vatren(nrofon,q1t,prf3cn))
 515 c
 516            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 517 c
 518            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q3t,prf1cn)+
 519      >                            trshu*vatren(nrofon,q1t,prf1cn)+
 520      >                            trsqu*vatren(nrofon,q6t,prf1cn)
 521 c
 522            vatrtt(nrofon,q6tp,f2cp) = vatren(nrofon,q4t,prf3cn)
 523 c
 524 c   dans le triangle NT3
 525 c
 526            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q3t,prf2cn)
 527 c
 528            vatrtt(nrofon,q2tp,f3cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 529      >                                vatren(nrofon,q1t,prf2cn))
 530 c
 531            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q6t,prf1cn)
 532 c
 533            vatrtt(nrofon,q4tp,f3cp) = vatren(nrofon,q6t,prf2cn)
 534 c
 535            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q1t,prf1cn)+
 536      >                            trshu*vatren(nrofon,q3t,prf1cn)+
 537      >                            trsqu*vatren(nrofon,q6t,prf1cn)
 538 c
 539            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 540 23332       continue
 541 c
 542 c   decoupage en 3 par l'arete a2 d'un triangle predecoupe en a3
 543 c
 544             elseif ((etan.eq.33).and.( etanp1.eq.32)) then
 545 c
 546             do 23333 , nrofon = 1 , nbfonc
 547 c
 548 c   dans le triangle NT1
 549 c
 550            vatrtt(nrofon,q1tp,f1cp) = unsde*(vatren(nrofon,q1t,prf2cn)+
 551      >                                vatren(nrofon,q3t,prf1cn))
 552 c
 553            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q4t,prf3cn)
 554 c
 555            vatrtt(nrofon,q3tp,f1cp) = vatren(nrofon,q3t,prf2cn)
 556 c
 557            vatrtt(nrofon,q4tp,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 558      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 559      >                                vatren(nrofon,q1t,prf1cn))+
 560      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 561      >                                vatren(nrofon,q5t,prf1cn))+
 562      >                       nessz*vatren(nrofon,q6t,prf1cn)
 563 c
 564            vatrtt(nrofon,q5tp,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 565      >                                vatren(nrofon,q4t,prf1cn)+
 566      >                                vatren(nrofon,q5t,prf1cn))-
 567      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
 568      >                                vatren(nrofon,q3t,prf1cn))+
 569      >                      unssz*vatren(nrofon,q6t,prf1cn)
 570 c
 571            vatrtt(nrofon,q6tp,f1cp) = vatren(nrofon,q6t,prf2cn)
 572 c
 573 c   dans le triangle NT2
 574 c
 575            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q3t,prf2cn)
 576 c
 577            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q4t,prf3cn)
 578 c
 579            vatrtt(nrofon,q3tp,f2cp) = vatren(nrofon,q2t,prf3cn)
 580 c
 581            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 582 c
 583            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q1t,prf3cn)+
 584      >                            trshu*vatren(nrofon,q2t,prf3cn)+
 585      >                            trsqu*vatren(nrofon,q4t,prf3cn)
 586 c
 587            vatrtt(nrofon,q6tp,f2cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 588      >                                vatren(nrofon,q2t,prf2cn)+
 589      >                                vatren(nrofon,q3t,prf3cn))
 590 c
 591 c   dans le triangle NT3
 592 c
 593            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q3t,prf1cn)
 594 c
 595            vatrtt(nrofon,q2tp,f3cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 596      >                                vatren(nrofon,q1t,prf3cn))
 597 c
 598            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q4t,prf3cn)
 599 c
 600            vatrtt(nrofon,q4tp,f3cp) = vatren(nrofon,q6t,prf1cn)
 601 c
 602            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q2t,prf3cn)+
 603      >                            trshu*vatren(nrofon,q1t,prf3cn)+
 604      >                            trsqu*vatren(nrofon,q4t,prf3cn)
 605 c
 606            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 607 23333       continue
 608 c
 609 c   decoupage en 3 par l'arete a1 d'un triangle predecoupe en a4
 610 c
 611             elseif ((etan.eq.34).and.( etanp1.eq.31)) then
 612 c
 613             do 23341 , nrofon = 1 , nbfonc
 614 c
 615 c   dans le triangle NT1
 616 c
 617            vatrtt(nrofon,q1tp,f1cp) = vatren(nrofon,q2t,prf3cn)
 618 c
 619            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q6t,prf2cn)
 620 c
 621            vatrtt(nrofon,q3tp,f1cp) = unsde*(vatren(nrofon,q1t,prf3cn)+
 622      >                                vatren(nrofon,q1t,prf1cn))
 623 c
 624            vatrtt(nrofon,q4tp,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 625      >                                vatren(nrofon,q4t,prf1cn)+
 626      >                                vatren(nrofon,q5t,prf1cn))-
 627      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
 628      >                                vatren(nrofon,q3t,prf1cn))+
 629      >                      unssz*vatren(nrofon,q6t,prf1cn)
 630 c
 631            vatrtt(nrofon,q5tp,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 632      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 633      >                                vatren(nrofon,q1t,prf1cn))+
 634      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 635      >                                vatren(nrofon,q5t,prf1cn))+
 636      >                       nessz*vatren(nrofon,q6t,prf1cn)
 637 c
 638            vatrtt(nrofon,q6tp,f1cp) = vatren(nrofon,q4t,prf3cn)
 639 c
 640 c   dans le triangle NT2
 641 c
 642            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q1t,prf1cn)
 643 c
 644            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q6t,prf2cn)
 645 c
 646            vatrtt(nrofon,q3tp,f2cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 647      >                                vatren(nrofon,q1t,prf2cn))
 648 c
 649            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 650 c
 651            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q3t,prf2cn)+
 652      >                            trshu*vatren(nrofon,q1t,prf2cn)+
 653      >                            trsqu*vatren(nrofon,q6t,prf2cn)
 654 c
 655            vatrtt(nrofon,q6tp,f2cp) = vatrtt(nrofon,q6t,f1cp)
 656 c
 657 c   dans le triangle NT3
 658 c
 659            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q2t,prf3cn)
 660 c
 661            vatrtt(nrofon,q2tp,f3cp) = vatren(nrofon,q3t,prf2cn)
 662 c
 663            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q6t,prf2cn)
 664 c
 665            vatrtt(nrofon,q4tp,f3cp) =  unstr*(vatren(nrofon,q2t,prf1cn)+
 666      >                                vatren(nrofon,q2t,prf2cn)+
 667      >                                vatren(nrofon,q3t,prf3cn))
 668 c
 669            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q1t,prf2cn)+
 670      >                            trshu*vatren(nrofon,q3t,prf2cn)+
 671      >                            trsqu*vatren(nrofon,q6t,prf2cn)
 672 c
 673            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 674
 675 23341       continue
 676 c
 677 c   decoupage en 3 par l'arete a2 d'un triangle predecoupe en a4
 678 c
 679             elseif ((etan.eq.34).and.( etanp1.eq.32)) then
 680 c
 681             do 23342 , nrofon = 1 , nbfonc
 682 c
 683 c   dans le triangle NT1
 684 c
 685            vatrtt(nrofon,q1tp,f1cp) = vatren(nrofon,q3t,prf2cn)
 686 c
 687            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q6t,prf1cn)
 688 c
 689            vatrtt(nrofon,q3tp,f1cp) = vatren(nrofon,q2t,prf3cn)
 690 c
 691            vatrtt(nrofon,q4tp,f1cp) = unsde*(vatren(nrofon,q5t,prf1cn)+
 692      >                                vatren(nrofon,q4t,prf2cn))
 693 c
 694            vatrtt(nrofon,q5tp,f1cp) = unsde*(vatren(nrofon,q4t,prf1cn)+
 695      >                                vatren(nrofon,q6t,prf3cn))
 696 c
 697            vatrtt(nrofon,q6tp,f1cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 698      >                                vatren(nrofon,q2t,prf2cn)+
 699      >                                vatren(nrofon,q3t,prf3cn))
 700 c
 701 c   dans le triangle NT2
 702 c
 703            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q2t,prf3cn)
 704 c
 705            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q6t,prf1cn)
 706 c
 707            vatrtt(nrofon,q3tp,f2cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 708      >                                vatren(nrofon,q1t,prf3cn))
 709 c
 710            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 711 c
 712            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q3t,prf1cn)+
 713      >                            trshu*vatren(nrofon,q1t,prf1cn)+
 714      >                            trsqu*vatren(nrofon,q6t,prf1cn)
 715 c
 716            vatrtt(nrofon,q6t,f2cp) = vatren(nrofon,q4t,prf3cn)
 717 c
 718 c   dans le triangle NT3
 719 c
 720            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q3t,prf2cn)
 721 c
 722            vatrtt(nrofon,q2tp,f3cp) = unsde*(vatren(nrofon,q3t,prf1cn)+
 723      >                                vatren(nrofon,q1t,prf2cn))
 724 c
 725            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q6t,prf1cn)
 726 c
 727            vatrtt(nrofon,q4tp,f3cp) = vatren(nrofon,q6t,prf2cn)
 728 c
 729            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q1t,prf1cn)+
 730      >                            trshu*vatren(nrofon,q3t,prf1cn)+
 731      >                            trsqu*vatren(nrofon,q6t,prf1cn)
 732 c
 733            vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
 734 c
 735 23342       continue
 736 c
 737 c   decoupage en 3 par l'arete a3 d'un triangle predecoupe en a4
 738 c
 739             elseif ((etan.eq.34).and.( etanp1.eq.33)) then
 740 c
 741             do 23343 , nrofon = 1 , nbfonc
 742 c
 743 c   dans le triangle NT1
 744 c
 745            vatrtt(nrofon,q1tp,f1cp) = unsde*(vatren(nrofon,q1t,prf2cn)+
 746      >                                vatren(nrofon,q3t,prf1cn))
 747 c
 748            vatrtt(nrofon,q2tp,f1cp) = vatren(nrofon,q4t,prf3cn)
 749 c
 750            vatrtt(nrofon,q3tp,f1cp) = vatren(nrofon,q3t,prf2cn)
 751 c
 752            vatrtt(nrofon,q4tp,f1cp) = -unshu*vatren(nrofon,q2t,prf1cn)-
 753      >                       trstr2*(vatren(nrofon,q3t,prf1cn)+
 754      >                                vatren(nrofon,q1t,prf1cn))+
 755      >                           trshu*(vatren(nrofon,q4t,prf1cn)+
 756      >                                vatren(nrofon,q5t,prf1cn))+
 757      >                       nessz*vatren(nrofon,q6t,prf1cn)
 758 c
 759            vatrtt(nrofon,q5tp,f1cp) = trshu*(vatren(nrofon,q2t,prf1cn)+
 760      >                                vatren(nrofon,q4t,prf1cn)+
 761      >                                vatren(nrofon,q5t,prf1cn))-
 762      >                      trstr2*(vatren(nrofon,q1t,prf1cn)+
 763      >                                vatren(nrofon,q3t,prf1cn))+
 764      >                      unssz*vatren(nrofon,q6t,prf1cn)
 765 c
 766            vatrtt(nrofon,q6tp,f1cp) = vatren(nrofon,q6t,prf2cn)
 767 c
 768 c   dans le triangle NT2
 769 c
 770            vatrtt(nrofon,q1tp,f2cp) = vatren(nrofon,q3t,prf2cn)
 771 c
 772            vatrtt(nrofon,q2tp,f2cp) = vatren(nrofon,q4t,prf3cn)
 773 c
 774            vatrtt(nrofon,q3tp,f2cp) = vatren(nrofon,q2t,prf3cn)
 775 c
 776            vatrtt(nrofon,q4tp,f2cp) = vatrtt(nrofon,q5t,f1cp)
 777 c
 778            vatrtt(nrofon,q5tp,f2cp) = -unshu*vatren(nrofon,q1t,prf3cn)+
 779      >                            trshu*vatren(nrofon,q2t,prf3cn)+
 780      >                            trsqu*vatren(nrofon,q4t,prf3cn)
 781 c
 782            vatrtt(nrofon,q6tp,f2cp) = unstr*(vatren(nrofon,q2t,prf1cn)+
 783      >                                vatren(nrofon,q2t,prf2cn)+
 784      >                                vatren(nrofon,q3t,prf3cn))
 785 c
 786 c   dans le triangle NT3
 787 c
 788            vatrtt(nrofon,q1tp,f3cp) = vatren(nrofon,q3t,prf1cn)
 789 c
 790            vatrtt(nrofon,q2tp,f3cp) = unsde*(vatren(nrofon,q1t,prf1cn)+
 791      >                                vatren(nrofon,q1t,prf3cn))
 792 c
 793            vatrtt(nrofon,q3tp,f3cp) = vatren(nrofon,q4t,prf3cn)
 794 c
 795            vatrtt(nrofon,q4tp,f3cp) = vatren(nrofon,q6t,prf1cn)
 796 c
 797            vatrtt(nrofon,q5tp,f3cp) = -unshu*vatren(nrofon,q2t,prf3cn)+
 798      >                            trshu*vatren(nrofon,q1t,prf3cn)+
 799      >                            trsqu*vatren(nrofon,q4t,prf3cn)
 800 c
 801            vatrtt(nrofon,q6tp,f3cp) = vatrtt(nrofon,q4t,f1cp)
 802 c
 803 23343       continue
 804             endif