Télécharger la plaquette Cast3M


fr en

Télécharger shapeq.eso

Retour à la liste

Numérotation des lignes : 

shapeq

  1. C SHAPEQ    SOURCE    GOUNAND   08/06/10    21:15:06     6125
  2.       SUBROUTINE SHAPEQ(QSI,ETA,DZE,IELE,SHP,IRET)
  3. C=======================================================================
  4. C
  5. C     FONCTIONS DE FORME TRIDIMENSIONNELLES QUAF
  6. C        EBERSOLT
  7. C=======================================================================
  8.       IMPLICIT INTEGER(I-N)
  9.       IMPLICIT REAL*8(A-H,O-Z)
  10.       PARAMETER (NDIM=3)
  11.       DIMENSION SHP(6,*)
  12. * Des constantes
  13.       REAL*8 UN
  14.       PARAMETER (UN=1.D0)
  15.  
  16. * Les matrices de passage : base polynomiale -> base des fonctions de
  17. * forme
  18.       REAL*8 PNM27(27,27)
  19.       REAL*8 PNM21(21,21)
  20.       REAL*8 PNM15(15,15)
  21. *
  22. *      REAL*8 PM271(27),PM272(27)
  23. *      EQUIVALENCE (PM27(1,1),PM271(1)),(PM27(1,2),PM272(1))
  24. *
  25. * Travail (on ne peut pas utiliser de segments ici : fortran pur)
  26. *
  27.       REAL*8 XCO1(NDIM)
  28.       REAL*8 VECTP(27)
  29. * Pour le cube
  30.       PARAMETER (NMMO27=1)
  31.       INTEGER NMON27(27)
  32.       REAL*8 CMON27(NMMO27,27)
  33.       INTEGER XMON27(NDIM,NMMO27,27)
  34. * Pour le prisme
  35.       PARAMETER (NMMO21=4)
  36.       INTEGER NMON21(21)
  37.       REAL*8 CMON21(NMMO21,21)
  38.       INTEGER XMON21(NDIM,NMMO21,21)
  39. * Pour le tétra
  40.       PARAMETER (NMMO15=4)
  41.       INTEGER NMON15(15)
  42.       REAL*8 CMON15(NMMO15,15)
  43.       INTEGER XMON15(NDIM,NMMO15,15)
  44. *
  45.       INTEGER ORDDER(NDIM)
  46. *
  47.       DATA ((PNM27(I,J),I=1, 27),J= 1, 1)/
  48.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  49.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  50.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  51.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  52.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  53.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  54.      $  0.0000000000000000D+00,  -1.2500000000000000D-01,
  55.      $  1.2500000000000000D-01,   0.0000000000000000D+00,
  56.      $  1.2500000000000000D-01,  -1.2500000000000000D-01,
  57.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  58.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  59.      $  1.2500000000000000D-01,  -1.2500000000000000D-01,
  60.      $  0.0000000000000000D+00,  -1.2500000000000000D-01,
  61.      $  1.2500000000000000D-01
  62.      $ /
  63.        DATA ((PNM27(I,J),I=1, 27),J= 2, 2)/
  64.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  65.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  66.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  67.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  68.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  69.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  70.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  71.      $ -2.5000000000000000D-01,  -2.5000000000000000D-01,
  72.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  73.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  74.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  75.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  76.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  77.      $ -2.5000000000000000D-01
  78.      $ /
  79.        DATA ((PNM27(I,J),I=1, 27),J= 3, 3)/
  80.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  81.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  82.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  83.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  84.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  85.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  86.      $  0.0000000000000000D+00,   1.2500000000000000D-01,
  87.      $  1.2500000000000000D-01,   0.0000000000000000D+00,
  88.      $ -1.2500000000000000D-01,  -1.2500000000000000D-01,
  89.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  90.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  91.      $ -1.2500000000000000D-01,  -1.2500000000000000D-01,
  92.      $  0.0000000000000000D+00,   1.2500000000000000D-01,
  93.      $  1.2500000000000000D-01
  94.      $ /
  95.        DATA ((PNM27(I,J),I=1, 27),J= 4, 4)/
  96.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  97.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  98.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  99.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  100.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  101.      $ -2.5000000000000000D-01,  -2.5000000000000000D-01,
  102.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  103.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  104.      $  2.5000000000000000D-01,   2.5000000000000000D-01,
  105.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  106.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  107.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  108.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  109.      $ -2.5000000000000000D-01
  110.      $ /
  111.        DATA ((PNM27(I,J),I=1, 27),J= 5, 5)/
  112.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  113.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  114.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  115.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  116.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  117.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  118.      $  0.0000000000000000D+00,  -1.2500000000000000D-01,
  119.      $ -1.2500000000000000D-01,   0.0000000000000000D+00,
  120.      $ -1.2500000000000000D-01,  -1.2500000000000000D-01,
  121.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  122.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  123.      $  1.2500000000000000D-01,   1.2500000000000000D-01,
  124.      $  0.0000000000000000D+00,   1.2500000000000000D-01,
  125.      $  1.2500000000000000D-01
  126.      $ /
  127.        DATA ((PNM27(I,J),I=1, 27),J= 6, 6)/
  128.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  129.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  130.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  131.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  132.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  133.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  134.      $ -2.5000000000000000D-01,   0.0000000000000000D+00,
  135.      $  2.5000000000000000D-01,  -2.5000000000000000D-01,
  136.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  137.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  138.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  139.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  140.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  141.      $ -2.5000000000000000D-01
  142.      $ /
  143.        DATA ((PNM27(I,J),I=1, 27),J= 7, 7)/
  144.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  145.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  146.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  147.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  148.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  149.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  150.      $  0.0000000000000000D+00,   1.2500000000000000D-01,
  151.      $ -1.2500000000000000D-01,   0.0000000000000000D+00,
  152.      $  1.2500000000000000D-01,  -1.2500000000000000D-01,
  153.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  154.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  155.      $ -1.2500000000000000D-01,   1.2500000000000000D-01,
  156.      $  0.0000000000000000D+00,  -1.2500000000000000D-01,
  157.      $  1.2500000000000000D-01
  158.      $ /
  159.        DATA ((PNM27(I,J),I=1, 27),J= 8, 8)/
  160.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  161.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  162.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  163.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  164.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  165.      $  2.5000000000000000D-01,  -2.5000000000000000D-01,
  166.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  167.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  168.      $ -2.5000000000000000D-01,   2.5000000000000000D-01,
  169.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  170.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  171.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  172.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  173.      $ -2.5000000000000000D-01
  174.      $ /
  175.        DATA ((PNM27(I,J),I=1, 27),J= 9, 9)/
  176.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  177.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  178.      $  2.5000000000000000D-01,  -2.5000000000000000D-01,
  179.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  180.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  181.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  182.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  183.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  184.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  185.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  186.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  187.      $ -2.5000000000000000D-01,   2.5000000000000000D-01,
  188.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  189.      $ -2.5000000000000000D-01
  190.      $ /
  191.        DATA ((PNM27(I,J),I=1, 27),J= 10, 10)/
  192.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  193.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  194.      $ -2.5000000000000000D-01,  -2.5000000000000000D-01,
  195.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  196.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  197.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  198.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  199.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  200.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  201.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  202.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  203.      $  2.5000000000000000D-01,   2.5000000000000000D-01,
  204.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  205.      $ -2.5000000000000000D-01
  206.      $ /
  207.        DATA ((PNM27(I,J),I=1, 27),J= 11, 11)/
  208.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  209.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  210.      $  2.5000000000000000D-01,   2.5000000000000000D-01,
  211.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  212.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  213.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  214.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  215.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  216.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  217.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  218.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  219.      $ -2.5000000000000000D-01,  -2.5000000000000000D-01,
  220.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  221.      $ -2.5000000000000000D-01
  222.      $ /
  223.        DATA ((PNM27(I,J),I=1, 27),J= 12, 12)/
  224.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  225.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  226.      $ -2.5000000000000000D-01,   2.5000000000000000D-01,
  227.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  228.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  229.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  230.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  231.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  232.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  233.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  234.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  235.      $  2.5000000000000000D-01,  -2.5000000000000000D-01,
  236.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  237.      $ -2.5000000000000000D-01
  238.      $ /
  239.        DATA ((PNM27(I,J),I=1, 27),J= 13, 13)/
  240.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  241.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  242.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  243.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  244.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  245.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  246.      $  0.0000000000000000D+00,   1.2500000000000000D-01,
  247.      $ -1.2500000000000000D-01,   0.0000000000000000D+00,
  248.      $ -1.2500000000000000D-01,   1.2500000000000000D-01,
  249.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  250.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  251.      $  1.2500000000000000D-01,  -1.2500000000000000D-01,
  252.      $  0.0000000000000000D+00,  -1.2500000000000000D-01,
  253.      $  1.2500000000000000D-01
  254.      $ /
  255.        DATA ((PNM27(I,J),I=1, 27),J= 14, 14)/
  256.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  257.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  258.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  259.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  260.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  261.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  262.      $ -2.5000000000000000D-01,   0.0000000000000000D+00,
  263.      $  2.5000000000000000D-01,   2.5000000000000000D-01,
  264.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  265.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  266.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  267.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  268.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  269.      $ -2.5000000000000000D-01
  270.      $ /
  271.        DATA ((PNM27(I,J),I=1, 27),J= 15, 15)/
  272.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  273.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  274.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  275.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  276.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  277.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  278.      $  0.0000000000000000D+00,  -1.2500000000000000D-01,
  279.      $ -1.2500000000000000D-01,   0.0000000000000000D+00,
  280.      $  1.2500000000000000D-01,   1.2500000000000000D-01,
  281.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  282.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  283.      $ -1.2500000000000000D-01,  -1.2500000000000000D-01,
  284.      $  0.0000000000000000D+00,   1.2500000000000000D-01,
  285.      $  1.2500000000000000D-01
  286.      $ /
  287.        DATA ((PNM27(I,J),I=1, 27),J= 16, 16)/
  288.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  289.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  290.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  291.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  292.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  293.      $  2.5000000000000000D-01,   2.5000000000000000D-01,
  294.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  295.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  296.      $ -2.5000000000000000D-01,  -2.5000000000000000D-01,
  297.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  298.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  299.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  300.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  301.      $ -2.5000000000000000D-01
  302.      $ /
  303.        DATA ((PNM27(I,J),I=1, 27),J= 17, 17)/
  304.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  305.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  306.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  307.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  308.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  309.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  310.      $  0.0000000000000000D+00,   1.2500000000000000D-01,
  311.      $  1.2500000000000000D-01,   0.0000000000000000D+00,
  312.      $  1.2500000000000000D-01,   1.2500000000000000D-01,
  313.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  314.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  315.      $  1.2500000000000000D-01,   1.2500000000000000D-01,
  316.      $  0.0000000000000000D+00,   1.2500000000000000D-01,
  317.      $  1.2500000000000000D-01
  318.      $ /
  319.        DATA ((PNM27(I,J),I=1, 27),J= 18, 18)/
  320.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  321.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  322.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  323.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  324.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  325.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  326.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  327.      $ -2.5000000000000000D-01,   2.5000000000000000D-01,
  328.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  329.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  330.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  331.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  332.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  333.      $ -2.5000000000000000D-01
  334.      $ /
  335.        DATA ((PNM27(I,J),I=1, 27),J= 19, 19)/
  336.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  337.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  338.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  339.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  340.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  341.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  342.      $  0.0000000000000000D+00,  -1.2500000000000000D-01,
  343.      $  1.2500000000000000D-01,   0.0000000000000000D+00,
  344.      $ -1.2500000000000000D-01,   1.2500000000000000D-01,
  345.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  346.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  347.      $ -1.2500000000000000D-01,   1.2500000000000000D-01,
  348.      $  0.0000000000000000D+00,  -1.2500000000000000D-01,
  349.      $  1.2500000000000000D-01
  350.      $ /
  351.        DATA ((PNM27(I,J),I=1, 27),J= 20, 20)/
  352.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  353.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  354.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  355.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  356.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  357.      $ -2.5000000000000000D-01,   2.5000000000000000D-01,
  358.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  359.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  360.      $  2.5000000000000000D-01,  -2.5000000000000000D-01,
  361.      $  0.0000000000000000D+00,  -2.5000000000000000D-01,
  362.      $  2.5000000000000000D-01,   0.0000000000000000D+00,
  363.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  364.      $  0.0000000000000000D+00,   2.5000000000000000D-01,
  365.      $ -2.5000000000000000D-01
  366.      $ /
  367.        DATA ((PNM27(I,J),I=1, 27),J= 21, 21)/
  368.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  369.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  370.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  371.      $  5.0000000000000000D-01,   0.0000000000000000D+00,
  372.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  373.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  374.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  375.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  376.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  377.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  378.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  379.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  380.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  381.      $  5.0000000000000000D-01
  382.      $ /
  383.        DATA ((PNM27(I,J),I=1, 27),J= 22, 22)/
  384.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  385.      $  5.0000000000000000D-01,   0.0000000000000000D+00,
  386.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  387.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  388.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  389.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  390.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  391.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  392.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  393.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  394.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  395.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  396.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  397.      $  5.0000000000000000D-01
  398.      $ /
  399.        DATA ((PNM27(I,J),I=1, 27),J= 23, 23)/
  400.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  401.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  402.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  403.      $  5.0000000000000000D-01,   0.0000000000000000D+00,
  404.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  405.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  406.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  407.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  408.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  409.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  410.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  411.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  412.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  413.      $  5.0000000000000000D-01
  414.      $ /
  415.        DATA ((PNM27(I,J),I=1, 27),J= 24, 24)/
  416.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  417.      $  5.0000000000000000D-01,   0.0000000000000000D+00,
  418.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  419.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  420.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  421.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  422.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  423.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  424.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  425.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  426.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  427.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  428.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  429.      $  5.0000000000000000D-01
  430.      $ /
  431.        DATA ((PNM27(I,J),I=1, 27),J= 25, 25)/
  432.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  433.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  434.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  435.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  436.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  437.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  438.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  439.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  440.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  441.      $  5.0000000000000000D-01,   0.0000000000000000D+00,
  442.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  443.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  444.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  445.      $  5.0000000000000000D-01
  446.      $ /
  447.        DATA ((PNM27(I,J),I=1, 27),J= 26, 26)/
  448.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  449.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  450.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  451.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  452.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  453.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  454.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  455.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  456.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  457.      $  5.0000000000000000D-01,   0.0000000000000000D+00,
  458.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  459.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  460.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  461.      $  5.0000000000000000D-01
  462.      $ /
  463.        DATA ((PNM27(I,J),I=1, 27),J= 27, 27)/
  464.      $  1.0000000000000000D+00,   0.0000000000000000D+00,
  465.      $ -1.0000000000000000D+00,   0.0000000000000000D+00,
  466.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  467.      $ -1.0000000000000000D+00,   0.0000000000000000D+00,
  468.      $  1.0000000000000000D+00,   0.0000000000000000D+00,
  469.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  470.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  471.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  472.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  473.      $ -1.0000000000000000D+00,   0.0000000000000000D+00,
  474.      $  1.0000000000000000D+00,   0.0000000000000000D+00,
  475.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  476.      $  1.0000000000000000D+00,   0.0000000000000000D+00,
  477.      $ -1.0000000000000000D+00
  478.      $ /
  479.        DATA ((PNM21(I,J),I=1, 21),J= 1, 1)/
  480.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  481.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  482.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  483.      $  1.5543122344752192D-15,  -5.0000000000000000D-01,
  484.      $  1.5000000000000000D+00,   1.5000000000000000D+00,
  485.      $ -1.0000000000000000D+00,  -2.0000000000000000D+00,
  486.      $ -1.0000000000000000D+00,  -1.4999999999999982D+00,
  487.      $  5.0000000000000000D-01,  -1.5000000000000000D+00,
  488.      $ -1.5000000000000000D+00,   1.0000000000000000D+00,
  489.      $  2.0000000000000000D+00,   1.0000000000000000D+00,
  490.      $  1.4999999999999967D+00
  491.      $ /
  492.        DATA ((PNM21(I,J),I=1, 21),J= 2, 2)/
  493.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  494.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  495.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  496.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  497.      $ -2.0000000000000000D+00,   0.0000000000000000D+00,
  498.      $  2.0000000000000000D+00,   2.0000000000000000D+00,
  499.      $  0.0000000000000000D+00,   5.9999999999999991D+00,
  500.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  501.      $  0.0000000000000000D+00,  -2.0000000000000000D+00,
  502.      $ -2.0000000000000000D+00,   0.0000000000000000D+00,
  503.      $ -5.9999999999999991D+00
  504.      $ /
  505.        DATA ((PNM21(I,J),I=1, 21),J= 3, 3)/
  506.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  507.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  508.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  509.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  510.      $  5.0000000000000000D-01,   0.0000000000000000D+00,
  511.      $ -1.0000000000000000D+00,   0.0000000000000000D+00,
  512.      $  0.0000000000000000D+00,  -1.4999999999999998D+00,
  513.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  514.      $  0.0000000000000000D+00,   1.0000000000000000D+00,
  515.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  516.      $  1.4999999999999998D+00
  517.      $ /
  518.        DATA ((PNM21(I,J),I=1, 21),J= 4, 4)/
  519.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  520.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  521.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  522.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  523.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  524.      $  0.0000000000000000D+00,  -2.0000000000000000D+00,
  525.      $  0.0000000000000000D+00,   5.9999999999999991D+00,
  526.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  527.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  528.      $  2.0000000000000000D+00,   0.0000000000000000D+00,
  529.      $ -5.9999999999999991D+00
  530.      $ /
  531.        DATA ((PNM21(I,J),I=1, 21),J= 5, 5)/
  532.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  533.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  534.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  535.      $  1.5543122344752192D-15,   0.0000000000000000D+00,
  536.      $  0.0000000000000000D+00,   5.0000000000000000D-01,
  537.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  538.      $ -1.0000000000000000D+00,  -1.4999999999999998D+00,
  539.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  540.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  541.      $  0.0000000000000000D+00,   1.0000000000000000D+00,
  542.      $  1.4999999999999982D+00
  543.      $ /
  544.        DATA ((PNM21(I,J),I=1, 21),J= 6, 6)/
  545.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  546.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  547.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  548.      $ -2.6645352591003757D-15,   0.0000000000000000D+00,
  549.      $  0.0000000000000000D+00,  -2.0000000000000000D+00,
  550.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  551.      $  2.0000000000000000D+00,   5.9999999999999991D+00,
  552.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  553.      $  2.0000000000000000D+00,   0.0000000000000000D+00,
  554.      $ -2.0000000000000000D+00,  -2.0000000000000000D+00,
  555.      $ -5.9999999999999964D+00
  556.      $ /
  557.        DATA ((PNM21(I,J),I=1, 21),J= 7, 7)/
  558.      $  1.0000000000000000D+00,  -3.0000000000000000D+00,
  559.      $ -3.0000000000000000D+00,   2.0000000000000000D+00,
  560.      $  4.0000000000000000D+00,   2.0000000000000000D+00,
  561.      $  2.9999999999999933D+00,   0.0000000000000000D+00,
  562.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  563.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  564.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  565.      $ -1.0000000000000000D+00,   3.0000000000000000D+00,
  566.      $  3.0000000000000000D+00,  -2.0000000000000000D+00,
  567.      $ -4.0000000000000000D+00,  -2.0000000000000000D+00,
  568.      $ -2.9999999999999933D+00
  569.      $ /
  570.        DATA ((PNM21(I,J),I=1, 21),J= 8, 8)/
  571.      $  0.0000000000000000D+00,  -1.0000000000000000D+00,
  572.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  573.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  574.      $  2.9999999999999996D+00,   0.0000000000000000D+00,
  575.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  576.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  577.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  578.      $  0.0000000000000000D+00,   1.0000000000000000D+00,
  579.      $  0.0000000000000000D+00,  -2.0000000000000000D+00,
  580.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  581.      $ -2.9999999999999996D+00
  582.      $ /
  583.        DATA ((PNM21(I,J),I=1, 21),J= 9, 9)/
  584.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  585.      $ -1.0000000000000000D+00,   0.0000000000000000D+00,
  586.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  587.      $  2.9999999999999964D+00,   0.0000000000000000D+00,
  588.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  589.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  590.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  591.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  592.      $  1.0000000000000000D+00,   0.0000000000000000D+00,
  593.      $  0.0000000000000000D+00,  -2.0000000000000000D+00,
  594.      $ -2.9999999999999964D+00
  595.      $ /
  596.        DATA ((PNM21(I,J),I=1, 21),J= 10, 10)/
  597.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  598.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  599.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  600.      $  1.5543122344752192D-15,   5.0000000000000000D-01,
  601.      $ -1.5000000000000000D+00,  -1.5000000000000000D+00,
  602.      $  1.0000000000000000D+00,   2.0000000000000000D+00,
  603.      $  1.0000000000000000D+00,   1.4999999999999982D+00,
  604.      $  5.0000000000000000D-01,  -1.5000000000000000D+00,
  605.      $ -1.5000000000000000D+00,   1.0000000000000000D+00,
  606.      $  2.0000000000000000D+00,   1.0000000000000000D+00,
  607.      $  1.4999999999999967D+00
  608.      $ /
  609.        DATA ((PNM21(I,J),I=1, 21),J= 11, 11)/
  610.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  611.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  612.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  613.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  614.      $  2.0000000000000000D+00,   0.0000000000000000D+00,
  615.      $ -2.0000000000000000D+00,  -2.0000000000000000D+00,
  616.      $  0.0000000000000000D+00,  -5.9999999999999991D+00,
  617.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  618.      $  0.0000000000000000D+00,  -2.0000000000000000D+00,
  619.      $ -2.0000000000000000D+00,   0.0000000000000000D+00,
  620.      $ -5.9999999999999991D+00
  621.      $ /
  622.        DATA ((PNM21(I,J),I=1, 21),J= 12, 12)/
  623.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  624.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  625.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  626.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  627.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  628.      $  1.0000000000000000D+00,   0.0000000000000000D+00,
  629.      $  0.0000000000000000D+00,   1.4999999999999998D+00,
  630.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  631.      $  0.0000000000000000D+00,   1.0000000000000000D+00,
  632.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  633.      $  1.4999999999999998D+00
  634.      $ /
  635.        DATA ((PNM21(I,J),I=1, 21),J= 13, 13)/
  636.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  637.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  638.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  639.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  640.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  641.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  642.      $  0.0000000000000000D+00,  -5.9999999999999991D+00,
  643.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  644.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  645.      $  2.0000000000000000D+00,   0.0000000000000000D+00,
  646.      $ -5.9999999999999991D+00
  647.      $ /
  648.        DATA ((PNM21(I,J),I=1, 21),J= 14, 14)/
  649.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  650.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  651.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  652.      $  1.5543122344752192D-15,   0.0000000000000000D+00,
  653.      $  0.0000000000000000D+00,  -5.0000000000000000D-01,
  654.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  655.      $  1.0000000000000000D+00,   1.4999999999999998D+00,
  656.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  657.      $ -5.0000000000000000D-01,   0.0000000000000000D+00,
  658.      $  0.0000000000000000D+00,   1.0000000000000000D+00,
  659.      $  1.4999999999999982D+00
  660.      $ /
  661.        DATA ((PNM21(I,J),I=1, 21),J= 15, 15)/
  662.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  663.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  664.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  665.      $ -3.5527136788005009D-15,   0.0000000000000000D+00,
  666.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  667.      $  0.0000000000000000D+00,  -2.0000000000000000D+00,
  668.      $ -2.0000000000000000D+00,  -5.9999999999999991D+00,
  669.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  670.      $  2.0000000000000000D+00,   0.0000000000000000D+00,
  671.      $ -2.0000000000000000D+00,  -2.0000000000000000D+00,
  672.      $ -5.9999999999999956D+00
  673.      $ /
  674.        DATA ((PNM21(I,J),I=1, 21),J= 16, 16)/
  675.      $  0.0000000000000000D+00,   4.0000000000000000D+00,
  676.      $  0.0000000000000000D+00,  -4.0000000000000000D+00,
  677.      $ -4.0000000000000000D+00,   0.0000000000000000D+00,
  678.      $ -1.1999999999999998D+01,   0.0000000000000000D+00,
  679.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  680.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  681.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  682.      $  0.0000000000000000D+00,  -4.0000000000000000D+00,
  683.      $  0.0000000000000000D+00,   4.0000000000000000D+00,
  684.      $  4.0000000000000000D+00,   0.0000000000000000D+00,
  685.      $  1.1999999999999998D+01
  686.      $ /
  687.        DATA ((PNM21(I,J),I=1, 21),J= 17, 17)/
  688.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  689.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  690.      $  4.0000000000000000D+00,   0.0000000000000000D+00,
  691.      $ -1.1999999999999998D+01,   0.0000000000000000D+00,
  692.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  693.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  694.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  695.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  696.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  697.      $ -4.0000000000000000D+00,   0.0000000000000000D+00,
  698.      $  1.1999999999999998D+01
  699.      $ /
  700.        DATA ((PNM21(I,J),I=1, 21),J= 18, 18)/
  701.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  702.      $  4.0000000000000000D+00,   0.0000000000000000D+00,
  703.      $ -4.0000000000000000D+00,  -4.0000000000000000D+00,
  704.      $ -1.1999999999999991D+01,   0.0000000000000000D+00,
  705.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  706.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  707.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  708.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  709.      $ -4.0000000000000000D+00,   0.0000000000000000D+00,
  710.      $  4.0000000000000000D+00,   4.0000000000000000D+00,
  711.      $  1.1999999999999991D+01
  712.      $ /
  713.        DATA ((PNM21(I,J),I=1, 21),J= 19, 19)/
  714.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  715.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  716.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  717.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  718.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  719.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  720.      $  0.0000000000000000D+00,  -1.3499999999999998D+01,
  721.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  722.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  723.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  724.      $  1.3499999999999998D+01
  725.      $ /
  726.        DATA ((PNM21(I,J),I=1, 21),J= 20, 20)/
  727.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  728.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  729.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  730.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  731.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  732.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  733.      $  0.0000000000000000D+00,   1.3499999999999998D+01,
  734.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  735.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  736.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  737.      $  1.3499999999999998D+01
  738.      $ /
  739.        DATA ((PNM21(I,J),I=1, 21),J= 21, 21)/
  740.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  741.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  742.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  743.      $  2.6999999999999996D+01,   0.0000000000000000D+00,
  744.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  745.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  746.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  747.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  748.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  749.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  750.      $ -2.6999999999999996D+01
  751.      $ /
  752.        DATA ((PNM15(I,J),I=1, 15),J= 1, 1)/
  753.      $  1.0000000000000000D+00,  -3.0000000000000000D+00,
  754.      $ -3.0000000000000000D+00,  -3.0000000000000000D+00,
  755.      $  2.0000000000000000D+00,   4.0000000000000000D+00,
  756.      $  4.0000000000000000D+00,   2.0000000000000000D+00,
  757.      $  4.0000000000000000D+00,   2.0000000000000000D+00,
  758.      $ -3.9999999999999334D+00,  -5.6477649705392858D-15,
  759.      $  2.9999999999999964D+00,   2.9999999999999964D+00,
  760.      $  2.9999999999999964D+00
  761.      $ /
  762.        DATA ((PNM15(I,J),I=1, 15),J= 2, 2)/
  763.      $  0.0000000000000000D+00,   4.0000000000000000D+00,
  764.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  765.      $ -4.0000000000000000D+00,  -4.0000000000000000D+00,
  766.      $ -4.0000000000000000D+00,   0.0000000000000000D+00,
  767.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  768.      $  3.1999999999999964D+01,   6.8821901302862365D-15,
  769.      $  0.0000000000000000D+00,  -1.1999999999999998D+01,
  770.      $ -1.1999999999999998D+01
  771.      $ /
  772.        DATA ((PNM15(I,J),I=1, 15),J= 3, 3)/
  773.      $  0.0000000000000000D+00,  -1.0000000000000000D+00,
  774.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  775.      $  2.0000000000000000D+00,   0.0000000000000000D+00,
  776.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  777.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  778.      $ -3.9999999999999920D+00,   2.9999999999999982D+00,
  779.      $  0.0000000000000000D+00,   2.9999999999999996D+00,
  780.      $  2.9999999999999996D+00
  781.      $ /
  782.        DATA ((PNM15(I,J),I=1, 15),J= 4, 4)/
  783.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  784.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  785.      $  0.0000000000000000D+00,   4.0000000000000000D+00,
  786.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  787.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  788.      $  3.1999999999999989D+01,  -1.1999999999999998D+01,
  789.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  790.      $ -1.1999999999999998D+01
  791.      $ /
  792.        DATA ((PNM15(I,J),I=1, 15),J= 5, 5)/
  793.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  794.      $ -1.0000000000000000D+00,   0.0000000000000000D+00,
  795.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  796.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  797.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  798.      $ -3.9999999999999920D+00,   2.9999999999999982D+00,
  799.      $  2.9999999999999996D+00,   0.0000000000000000D+00,
  800.      $  2.9999999999999996D+00
  801.      $ /
  802.        DATA ((PNM15(I,J),I=1, 15),J= 6, 6)/
  803.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  804.      $  4.0000000000000000D+00,   0.0000000000000000D+00,
  805.      $  0.0000000000000000D+00,  -4.0000000000000000D+00,
  806.      $  0.0000000000000000D+00,  -4.0000000000000000D+00,
  807.      $ -4.0000000000000000D+00,   0.0000000000000000D+00,
  808.      $  3.1999999999999964D+01,   6.8821901302862365D-15,
  809.      $ -1.1999999999999998D+01,   0.0000000000000000D+00,
  810.      $ -1.1999999999999998D+01
  811.      $ /
  812.        DATA ((PNM15(I,J),I=1, 15),J= 7, 7)/
  813.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  814.      $  0.0000000000000000D+00,   4.0000000000000000D+00,
  815.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  816.      $ -4.0000000000000000D+00,   0.0000000000000000D+00,
  817.      $ -4.0000000000000000D+00,  -4.0000000000000000D+00,
  818.      $  3.1999999999999964D+01,   6.8821901302862365D-15,
  819.      $ -1.1999999999999998D+01,  -1.1999999999999998D+01,
  820.      $  0.0000000000000000D+00
  821.      $ /
  822.        DATA ((PNM15(I,J),I=1, 15),J= 8, 8)/
  823.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  824.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  825.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  826.      $  4.0000000000000000D+00,   0.0000000000000000D+00,
  827.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  828.      $  3.1999999999999989D+01,  -1.1999999999999998D+01,
  829.      $  0.0000000000000000D+00,  -1.1999999999999998D+01,
  830.      $  0.0000000000000000D+00
  831.      $ /
  832.        DATA ((PNM15(I,J),I=1, 15),J= 9, 9)/
  833.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  834.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  835.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  836.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  837.      $  4.0000000000000000D+00,   0.0000000000000000D+00,
  838.      $  3.1999999999999989D+01,  -1.1999999999999998D+01,
  839.      $ -1.1999999999999998D+01,   0.0000000000000000D+00,
  840.      $  0.0000000000000000D+00
  841.      $ /
  842.        DATA ((PNM15(I,J),I=1, 15),J= 10, 10)/
  843.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  844.      $  0.0000000000000000D+00,  -1.0000000000000000D+00,
  845.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  846.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  847.      $  0.0000000000000000D+00,   2.0000000000000000D+00,
  848.      $ -3.9999999999999920D+00,   2.9999999999999982D+00,
  849.      $  2.9999999999999996D+00,   2.9999999999999996D+00,
  850.      $  0.0000000000000000D+00
  851.      $ /
  852.        DATA ((PNM15(I,J),I=1, 15),J= 11, 11)/
  853.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  854.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  855.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  856.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  857.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  858.      $ -1.0799999999999996D+02,  -7.1183903498257817D-15,
  859.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  860.      $  2.6999999999999996D+01
  861.      $ /
  862.        DATA ((PNM15(I,J),I=1, 15),J= 12, 12)/
  863.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  864.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  865.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  866.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  867.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  868.      $ -1.0799999999999996D+02,  -7.1183903498257817D-15,
  869.      $  0.0000000000000000D+00,   2.6999999999999996D+01,
  870.      $  0.0000000000000000D+00
  871.      $ /
  872.        DATA ((PNM15(I,J),I=1, 15),J= 13, 13)/
  873.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  874.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  875.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  876.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  877.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  878.      $ -1.0800000000000003D+02,   2.7000000000000007D+01,
  879.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  880.      $  0.0000000000000000D+00
  881.      $ /
  882.        DATA ((PNM15(I,J),I=1, 15),J= 14, 14)/
  883.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  884.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  885.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  886.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  887.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  888.      $ -1.0799999999999996D+02,  -7.1183903498257817D-15,
  889.      $  2.6999999999999996D+01,   0.0000000000000000D+00,
  890.      $  0.0000000000000000D+00
  891.      $ /
  892.        DATA ((PNM15(I,J),I=1, 15),J= 15, 15)/
  893.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  894.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  895.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  896.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  897.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  898.      $  2.5600000000000000D+02,  -2.2189281662088338D-15,
  899.      $  0.0000000000000000D+00,   0.0000000000000000D+00,
  900.      $  0.0000000000000000D+00
  901.      $ /
  902. *
  903. *
  904. C=======================================================================
  905. C     QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS
  906. C
  907. C     SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE FORME DU ND INO
  908. C                            SHP(2,INO) DERIVEES PR RPPRT A QSI
  909. C                            SHP(3,INO) DERIVEES PR RPPRT A ETA
  910. C                            SHP(4,INO) DERIVEES PR RPPRT A DZE
  911. C                            SHP(5,INO) DERIVEES PR RPPRT A  4
  912. C                              DERIVEES BARYCENTRIQUE
  913. C     CAS DE DIMENSION 3
  914. C               15 NOEUDS TETRAEDRE
  915. C               19 NOEUDS PYRAMIDE
  916. C               21 NOEUDS PRISME
  917. C               27 NOEUDS CUBE
  918. C=======================================================================
  919. C
  920. C CUB27
  921. C
  922. C=======================================================================
  923.       IF (IELE.EQ.33) THEN
  924.          NBNOD=27
  925. *     Init exposant des polynomes de la base polynomiale
  926.          INOD=0
  927.          DO IDZE=0,2
  928.             DO IETA=0,2
  929.                DO IKSI=0,2
  930.                   INOD=INOD+1
  931.                   NMON27(INOD)=1
  932.                   CMON27(1,INOD)=1.D0
  933.                   XMON27(1,1,INOD)=IKSI
  934.                   XMON27(2,1,INOD)=IETA
  935.                   XMON27(3,1,INOD)=IDZE
  936.                ENDDO
  937.             ENDDO
  938.          ENDDO
  939.          ORDDER(1)=0
  940.          ORDDER(2)=0
  941.          ORDDER(3)=0
  942. * Calcul des fonctions de forme et des dérivées premières
  943. * au point QSI,ETA,DZE
  944. *         WRITE(6,*) 'QSI=',QSI,' ETA=',ETA,' DZE=',DZE
  945.          CALL NI2(NDIM,NBNOD,QSI,ETA,DZE,
  946.      $        NMMO27,NMON27,CMON27,XMON27,PNM27,XCO1,ORDDER,
  947.      $        VECTP,
  948.      $        SHP)
  949. C=======================================================================
  950. C
  951. C PR21
  952. C
  953. C=======================================================================
  954.       ELSEIF (IELE.EQ.34) THEN
  955.          NBNOD=21
  956. *     Init exposant des polynomes de la base polynomiale
  957.          INOD=0
  958.          DO IDZE=0,2
  959. *     D'abord la base polynomiale complète
  960.             DO IDEG=0,2
  961.                DO IEXP=IDEG,0,-1
  962.                   INOD=INOD+1
  963.                   NMON21(INOD)=1
  964.                   CMON21(1,INOD)=1.D0
  965.                   XMON21(1,1,INOD)=IEXP
  966.                   XMON21(2,1,INOD)=IDEG-IEXP
  967.                   XMON21(3,1,INOD)=IDZE
  968.                ENDDO
  969.             ENDDO
  970. *     Puis la bulle
  971. * On génère un polynôme du type :
  972. * x_1.x_2 (1-x_1-x_2)
  973.             INOD=INOD+1
  974.             NMON21(INOD)=3
  975. * On initialise les monômes à x_1.x_2
  976.             DO IMON=1,NMON21(INOD)
  977.                XMON21(1,IMON,INOD)=1
  978.                XMON21(2,IMON,INOD)=1
  979.                XMON21(3,IMON,INOD)=IDZE
  980.             ENDDO
  981. * On ajuste les coeff. et les monômes
  982.             CMON21(1,INOD)=1.D0
  983.             DO IMON=2,NMON21(INOD)
  984.                CMON21(IMON,INOD)=-1.D0
  985.                XMON21(IMON-1,IMON,INOD)=XMON21(IMON-1,IMON,INOD)+1
  986.             ENDDO
  987.          ENDDO
  988.          ORDDER(1)=0
  989.          ORDDER(2)=0
  990.          ORDDER(3)=0
  991. * Calcul des fonctions de forme et des dérivées premières
  992. * au point QSI,ETA,DZE
  993. *         WRITE(6,*) 'QSI=',QSI,' ETA=',ETA,' DZE=',DZE
  994.          CALL NI2(NDIM,NBNOD,QSI,ETA,DZE,
  995.      $        NMMO21,NMON21,CMON21,XMON21,PNM21,XCO1,ORDDER,
  996.      $        VECTP,
  997.      $        SHP)
  998. C=======================================================================
  999. C
  1000. C TE15
  1001. C
  1002. C=======================================================================
  1003.       ELSEIF (IELE.EQ.35) THEN
  1004.          NBNOD=15
  1005. *     Init exposant des polynomes de la base polynomiale
  1006.          INOD=0
  1007. *     D'abord la base polynomiale complète
  1008.          DO IDEG=0,2
  1009.             DO IEXP=IDEG,0,-1
  1010.                DO IEXP2=IDEG-IEXP,0,-1
  1011.                   INOD=INOD+1
  1012.                   NMON15(INOD)=1
  1013.                   CMON15(1,INOD)=1.D0
  1014.                   XMON15(1,1,INOD)=IEXP
  1015.                   XMON15(2,1,INOD)=IEXP2
  1016.                   XMON15(3,1,INOD)=IDEG-(IEXP+IEXP2)
  1017.                ENDDO
  1018.             ENDDO
  1019.          ENDDO
  1020. *     Puis la bulle du centre
  1021. * On génère un polynôme du type :
  1022. * x_1.x_2.x_3 (1-x_1-x_2-x_3)
  1023.          INOD=INOD+1
  1024.          NMON15(INOD)=4
  1025. * On initialise les monômes à x_1.x_2.x_3
  1026.          DO IMON=1,NMON15(INOD)
  1027.             XMON15(1,IMON,INOD)=1
  1028.             XMON15(2,IMON,INOD)=1
  1029.             XMON15(3,IMON,INOD)=1
  1030.          ENDDO
  1031. * On ajuste les coeff. et les monômes
  1032.          CMON15(1,INOD)=1.D0
  1033.          DO IMON=2,NMON15(INOD)
  1034.             CMON15(IMON,INOD)=-1.D0
  1035.             XMON15(IMON-1,IMON,INOD)=XMON15(IMON-1,IMON,INOD)+1
  1036.          ENDDO
  1037. *     Puis les bulles des faces
  1038. * On génère le polynôme : x_1.x_2...x_n
  1039.          INOD=INOD+1
  1040.          NMON15(INOD)=1
  1041.          CMON15(1,INOD)=1.D0
  1042.          XMON15(1,1,INOD)=1
  1043.          XMON15(2,1,INOD)=1
  1044.          XMON15(3,1,INOD)=1
  1045. * On génère les polynômes du type :
  1046. * ((x_1.x_2,x_3)/x_i) (1-x_1...-x_3)
  1047.          DO JDIM=1,NDIM
  1048.             INOD=INOD+1
  1049.             NMON15(INOD)=4
  1050. * On initialise les monômes à x_1.x_2.x_3
  1051.             DO IMON=1,NMON15(INOD)
  1052.                XMON15(1,IMON,INOD)=1
  1053.                XMON15(2,IMON,INOD)=1
  1054.                XMON15(3,IMON,INOD)=1
  1055.             ENDDO
  1056. * On divise par x_i
  1057.             DO IMON=1,NMON15(INOD)
  1058.                XMON15(JDIM,IMON,INOD)=XMON15(JDIM,IMON,INOD)-1
  1059.             ENDDO
  1060. * On ajuste les coeff. et les monômes
  1061.             CMON15(1,INOD)=UN
  1062.             DO IMON=2,NMON15(INOD)
  1063.                CMON15(IMON,INOD)=-UN
  1064.                XMON15(IMON-1,IMON,INOD)=XMON15(IMON-1,IMON,INOD)+1
  1065.             ENDDO
  1066.          ENDDO
  1067.          ORDDER(1)=0
  1068.          ORDDER(2)=0
  1069.          ORDDER(3)=0
  1070. * Calcul des fonctions de forme et des dérivées premières
  1071. * au point QSI,ETA,DZE
  1072. *         WRITE(6,*) 'QSI=',QSI,' ETA=',ETA,' DZE=',DZE
  1073.          CALL NI2(NDIM,NBNOD,QSI,ETA,DZE,
  1074.      $        NMMO15,NMON15,CMON15,XMON15,PNM15,XCO1,ORDDER,
  1075.      $        VECTP,
  1076.      $        SHP)
  1077.       ELSE
  1078. *         WRITE(6,*) 'IELE=',IELE,' non implem'
  1079. *         CALL ERREUR(5)
  1080.          IRET=0
  1081.          RETURN
  1082.       ENDIF
  1083.       IRET=1
  1084.       RETURN
  1085.       END
  1086.  
  1087.  
  1088.  
  1089.  
  1090.  
  1091.  

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