mt03cc.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\mt03cc.f
00002 !
00068                      SUBROUTINE MT03CC
00069 !                    *****************
00070 !
00071      &( A11 , A12 , A13 , A14 , A15 , A16 ,
00072      &  A21 , A22 , A23 , A24 , A25 , A26 ,
00073      &  A31 , A32 , A33 , A34 , A35 , A36 ,
00074      &  A41 , A42 , A43 , A44 , A45 , A46 ,
00075      &  A51 , A52 , A53 , A54 , A55 , A56 ,
00076      &  A61 , A62 , A63 , A64 , A65 , A66 ,
00077      &  XMUL,SF,SG,SU,SV,F,G,U,V,
00078      &  XEL,YEL,IKLE1,IKLE2,IKLE3,IKLE4,IKLE5,IKLE6,
00079      &  NELEM,NELMAX)
00080 !
00081 !***********************************************************************
00082 ! BIEF   V6P1                                   21/08/2010
00083 !***********************************************************************
00084 !
00085 !
00086 !
00087 !
00088 !
00089 !
00090 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00091 !| A11            |<--| ELEMENTS OF MATRIX
00092 !| ...            |<--| ELEMENTS OF MATRIX
00093 !| A66            |<--| ELEMENTS OF MATRIX
00094 !| F              |-->| FUNCTION USED IN THE FORMULA
00095 !| G              |-->| FUNCTION USED IN THE FORMULA
00096 !| IKLE1          |-->| FIRST POINTS OF TRIANGLES
00097 !| IKLE2          |-->| SECOND POINTS OF TRIANGLES
00098 !| IKLE3          |-->| THIRD POINTS OF TRIANGLES
00099 !| IKLE4          |-->| FOURTH POINTS OF TRIANGLES (QUADRATIC)
00100 !| IKLE5          |-->| FIFTH POINTS OF TRIANGLES (QUADRATIC)
00101 !| IKLE6          |-->| SIXTH POINTS OF TRIANGLES (QUADRATIC)
00102 !| NELEM          |-->| NUMBER OF ELEMENTS
00103 !| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
00104 !| SF             |-->| STRUCTURE OF FUNCTIONS F
00105 !| SG             |-->| STRUCTURE OF FUNCTIONS G
00106 !| SU             |-->| BIEF_OBJ STRUCTURE OF U
00107 !| SV             |-->| BIEF_OBJ STRUCTURE OF V
00108 !| U              |-->| FUNCTION U USED IN THE FORMULA
00109 !| V              |-->| FUNCTION V USED IN THE FORMULA
00110 !| XEL            |-->| ABSCISSAE OF POINTS IN THE MESH, PER ELEMENT
00111 !| YEL            |-->| ORDINATES OF POINTS IN THE MESH, PER ELEMENT
00112 !| XMUL           |-->| MULTIPLICATION FACTOR
00113 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00114 !
00115       USE BIEF!, EX_MT03CC => MT03CC
00116 !
00117       IMPLICIT NONE
00118       INTEGER LNG,LU
00119       COMMON/INFO/LNG,LU
00120 !
00121 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00122 !
00123       INTEGER, INTENT(IN) :: NELEM,NELMAX
00124       INTEGER, INTENT(IN) :: IKLE1(NELMAX),IKLE2(NELMAX)
00125       INTEGER, INTENT(IN) :: IKLE3(NELMAX),IKLE4(NELMAX)
00126       INTEGER, INTENT(IN) :: IKLE5(NELMAX),IKLE6(NELMAX)
00127 !
00128       DOUBLE PRECISION, INTENT(INOUT) :: A11(*),A12(*),A13(*)
00129       DOUBLE PRECISION, INTENT(INOUT) :: A14(*),A15(*),A16(*)
00130       DOUBLE PRECISION, INTENT(INOUT) :: A21(*),A22(*),A23(*)
00131       DOUBLE PRECISION, INTENT(INOUT) :: A24(*),A25(*),A26(*)
00132       DOUBLE PRECISION, INTENT(INOUT) :: A31(*),A32(*),A33(*)
00133       DOUBLE PRECISION, INTENT(INOUT) :: A34(*),A35(*),A36(*)
00134       DOUBLE PRECISION, INTENT(INOUT) :: A41(*),A42(*),A43(*)
00135       DOUBLE PRECISION, INTENT(INOUT) :: A44(*),A45(*),A46(*)
00136       DOUBLE PRECISION, INTENT(INOUT) :: A51(*),A52(*),A53(*)
00137       DOUBLE PRECISION, INTENT(INOUT) :: A54(*),A55(*),A56(*)
00138       DOUBLE PRECISION, INTENT(INOUT) :: A61(*),A62(*),A63(*)
00139       DOUBLE PRECISION, INTENT(INOUT) :: A64(*),A65(*),A66(*)
00140 !
00141       DOUBLE PRECISION, INTENT(IN) :: XMUL
00142       DOUBLE PRECISION, INTENT(IN) :: F(*),G(*),U(*),V(*)
00143 !
00144 !     STRUCTURES OF      F, G, U, V
00145       TYPE(BIEF_OBJ), INTENT(IN) :: SF,SG,SU,SV
00146 !
00147       DOUBLE PRECISION, INTENT(IN) :: XEL(NELMAX,3),YEL(NELMAX,3)
00148 !
00149 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00150 !
00151 !     DECLARATIONS SPECIFIC TO THIS SUBROUTINE
00152 !
00153       INTEGER IELEM,IELMF,IELMG,IELMU,IELMV
00154 !
00155       DOUBLE PRECISION X2,X3,Y2,Y3,U1,U2,U3,U4,U5,U6,V1,V2,V3,V4,V5,V6
00156       DOUBLE PRECISION KXEL,KYEL,ANS1,ANS2
00157 !
00158 !-----------------------------------------------------------------------
00159 !
00160       IELMF = SF%ELM
00161       IELMG = SG%ELM
00162       IELMU = SU%ELM
00163       IELMV = SV%ELM
00164 !
00165 !-----------------------------------------------------------------------
00166 ! CASE WHERE U IS OF P1 DISCRETISATION
00167 !-----------------------------------------------------------------------
00168 !
00169       IF(IELMF.EQ.10.AND.IELMG.EQ.10.AND.
00170      &   IELMU.EQ.11.AND.IELMV.EQ.11) THEN
00171 !
00172 !   LOOP ON THE ELEMENTS
00173 !
00174       DO IELEM = 1 , NELEM
00175 !
00176       X2  =   XEL(IELEM,2)
00177       X3  =   XEL(IELEM,3)
00178       Y2  =   YEL(IELEM,2)
00179       Y3  =   YEL(IELEM,3)
00180 !
00181       U1   =  U(IKLE1(IELEM))
00182       U2   =  U(IKLE2(IELEM))
00183       U3   =  U(IKLE3(IELEM))
00184       V1   =  V(IKLE1(IELEM))
00185       V2   =  V(IKLE2(IELEM))
00186       V3   =  V(IKLE3(IELEM))
00187 !
00188       KXEL =  F(IELEM)
00189       KYEL =  G(IELEM)
00190 !
00191 !  COMPUTES 30 OF THE 36 TERMS (SELECTED AMONG THE LEAST COMPLEX)
00192 !
00193       A11(IELEM)= ( KXEL * (Y3-Y2) - KYEL * (X3-X2) ) *
00194      &            ((X2-X3)*(V3+3.D0*V1+V2) + (Y3-Y2)*(U2+3.D0*U1+U3)) *
00195      &              XMUL/((X2*Y3-Y2*X3)*10.D0)
00196 !
00197       A12(IELEM)= ( KXEL * (Y3-Y2) - KYEL * (X3-X2) ) *
00198      &            ( X3 * (V3+2.D0*(V1+V2)) -  Y3 * (U3+2.D0*(U1+U2))) *
00199      &             XMUL/((X3*Y2-Y3*X2)*30.D0)
00200 !
00201       A13(IELEM)= ( KXEL * (Y3-Y2) - KYEL * (X3-X2) ) *
00202      &            ( X2 * (V2+2.D0*(V1+V3)) -  Y2 * (U2+2.D0*(U1+U3))) *
00203      &             XMUL/((X2*Y3-Y2*X3)*30.D0)
00204 !
00205       A14(IELEM)= ( KXEL * (Y3-Y2) - KYEL * (X3-X2) ) *
00206      &            ( X2   *( -   V3+ 3.D0*V1- 2.D0*V2) +
00207      &              X3   *(4.D0*V3+11.D0*V1+ 5.D0*V2) +
00208      &              Y2   *(     U3+ 2.D0*U2- 3.D0*U1) -
00209      &              Y3   *(11.D0*U1+5.D0*U2+ 4.D0*U3))*
00210      &              XMUL/((X2*Y3-Y2*X3)*30.D0)
00211 !
00212       A15(IELEM)= ( KXEL * (Y3-Y2) - KYEL * (X3-X2) ) *
00213      &            ( X3   *(     V2+2.D0*V3-  3.D0*V1) +
00214      &              X2   *(3.D0*V1-     V3-  2.D0*V2) +
00215      &              Y3   *(3.D0*U1-     U2-  2.D0*U3) +
00216      &              Y2   *(-3.D0*U1+2.D0*U2+      U3))*
00217      &              XMUL/((X3*Y2-Y3*X2)*30.D0)
00218 !
00219       A21(IELEM)= ( KYEL * X3      - KXEL *    Y3  ) *
00220      &            ( (X2-X3) * ( V3 + 2.D0*(V1+V2) )  +
00221      &              (Y3-Y2) * ( U3 + 2.D0*(U1+U2) )) *
00222      &               XMUL/((X3*Y2-Y3*X2)*30.D0)
00223 !
00224       A22(IELEM)= ( KYEL * X3      - KXEL *    Y3  ) *
00225      &            ( X3 * (V3+V1+3.D0*V2) - Y3*(U3+U1+3.D0*U2)) *
00226      &              XMUL/((X2*Y3-Y2*X3)*10.D0)
00227 !
00228       A23(IELEM)= ( KYEL * X3      - KXEL *    Y3 ) *
00229      &            ( X2 * (V1+2.D0*(V2+V3)) - Y2 * (U1+2.D0*(U2+U3))) *
00230      &              XMUL/((X2*Y3-Y2*X3)*30.D0)
00231 !
00232       A24(IELEM)= ( KYEL * X3      - KXEL *    Y3 ) *
00233      &            ( -X3*(4.D0*V3+ 5.D0*V1+11.D0*V2) +
00234      &               X2*(3.D0*V1+14.D0*V2+ 3.D0*V3) +
00235      &               Y3*(5.D0*U1+11.D0*U2+ 4.D0*U3) -
00236      &               Y2*(3.D0*U1+14.D0*U2+ 3.D0*U3))*
00237      &               XMUL/((X2*Y3-Y2*X3)*30.D0)
00238 !
00239       A26(IELEM)= ( KXEL *    Y3  -  KYEL * X3    ) *
00240      &            (X3*(2.D0*V3 +      V1 - 3.D0*V2) +
00241      &             X2*(V1-V3)  + Y2*(U3-U1)         +
00242      &             Y3*(3.D0*U2 -      U1 - 2.D0*U3))*
00243      &             XMUL/((X3*Y2-Y3*X2)*30.D0)
00244 !
00245       A31(IELEM)= ( KYEL * X2      - KXEL *    Y2 ) *
00246      &            ( (X2-X3) * (2.D0*(V1+V3) + V2  ) +
00247      &              (Y3-Y2) * (2.D0*(U1+U3) + U2 )) *
00248      &               XMUL/((X2*Y3-Y2*X3)*30.D0)
00249 !
00250       A32(IELEM)= ( KYEL * X2      - KXEL *    Y2 ) *
00251      &            ( X3*(2.D0*(V3+V2)+V1)-Y3*(U1+2.D0*(U2+U3))) *
00252      &              XMUL/((X2*Y3-Y2*X3)*30.D0)
00253 !
00254       A33(IELEM)= ( KYEL*X2        - KXEL *    Y2 ) *
00255      &            (X2*(V1+3.D0*V3+V2)-Y2*(U1+3.D0*U3+U2)) *
00256      &             XMUL/(10.D0*(X2*Y3-Y2*X3))
00257 !
00258       A34(IELEM)= ( KYEL*X2        - KXEL *    Y2 ) *
00259      &            (X3*(V1-V2)+X2*(V1+2.D0*V2-3.D0*V3) +
00260      &             Y3*(U2-U1)+Y2*(3.D0*U3-U1-2.D0*U2))*
00261      &             XMUL/((X2*Y3-Y2*X3)*30.D0)
00262 !
00263       A35(IELEM)=  ( KXEL*Y2       - KYEL *    X2 ) *
00264      &             (X3*(14.D0*V3+3.D0*(V2+V1)) +
00265      &              X2*(      V1-3.D0*V3+2.D0*V2) -
00266      &              Y3*(3.D0*(U1+U2)+14.D0*U3) +
00267      &              Y2*( 3.D0*U3-     U1- 2.D0*U2 ))*
00268      &              XMUL/((X2*Y3-Y2*X3)*30.D0)
00269 !
00270       A36(IELEM)= ( KYEL*X2        - KXEL *    Y2 ) *
00271      &            ( X3*(-14.D0*V3-3.D0*(V1+V2)) +
00272      &              X2*(  5.D0*V1+4.D0*V2+11.D0*V3) +
00273      &              Y3*(3.D0*(U1+U2)+14.D0*U3) -
00274      &              Y2*(5.D0*U1+4.D0*U2+11.D0*U3)) *
00275      &              XMUL/((X3*Y2-Y3*X2)*30.D0)
00276 !
00277       A42(IELEM)=(-5.D0*KXEL*Y3**2*U1+11.D0*KYEL*X3*Y3*U2+
00278      &            5.D0*KYEL*X3*Y3*U1+3.D0*KYEL*X3*X2*V1+
00279      &            5.D0*KXEL*Y3*X3*V1+3.D0*KXEL*Y3*Y2*U1-
00280      &            5.D0*KYEL*X3**2*V1-3.D0*KXEL*Y2*X3*V3-
00281      &            3.D0*KYEL*X2*Y3*U3+4.D0*KYEL*X3*Y3*U3+
00282      &            3.D0*KYEL*X3*X2*V3+3.D0*KXEL*Y3*Y2*U3+
00283      &            4.D0*KXEL*Y3*X3*V3-4.D0*KXEL*Y3**2*U3-
00284      &            4.D0*KYEL*X3**2*V3-14.D0*KXEL*Y2*X3*V2-
00285      &            14.D0*KYEL*X2*Y3*U2+14.D0*KXEL*Y3*Y2*U2+
00286      &            11.D0*KXEL*Y3*X3*V2+14.D0*KYEL*X3*X2*V2-
00287      &            11.D0*KYEL*X3**2*V2-11.D0*KXEL*Y3**2*U2-
00288      &            3.D0*KYEL*X2*Y3*U1-3.D0*KXEL*Y2*X3*V1)*
00289      &            XMUL/((X2*Y3-Y2*X3)*30.D0)
00290 !
00291       A43(IELEM)=(KYEL*X3*X2*V1-KYEL*X3*Y2*U1-KXEL*Y3*X2*V1+
00292      &            KXEL*Y3*Y2*U1+KYEL*X2**2*V1+3.D0*KYEL*X2*Y2*U3+
00293      &            3.D0*KXEL*Y2*X2*V3-3.D0*KYEL*X2**2*V3-
00294      &            3.D0*KXEL*Y2**2*U3-2.D0*KXEL*Y2*X2*V2+
00295      &            2.D0*KYEL*X2**2*V2+2.D0*KXEL*Y2**2*U2-
00296      &            2.D0*KYEL*X2*Y2*U2-KXEL*Y3*Y2*U2+KXEL*Y3*X2*V2+
00297      &            KYEL*X3*Y2*U2-KYEL*X3*X2*V2+KXEL*Y2**2*U1-
00298      &            KYEL*X2*Y2*U1-KXEL*Y2*X2*V1)*
00299      &            XMUL/((X2*Y3-Y2*X3)*30.D0)
00300 !
00301       A44(IELEM)=(-4.D0*KYEL*X3*Y3*U2-4.D0*KYEL*X3*Y3*U1-
00302      &            4.D0*KXEL*Y3*X3*V1-2.D0*KYEL*X2*Y2*U3+
00303      &            KXEL*Y2*X3*V3-2.D0*KXEL*Y2*X2*V3+KYEL*X2*Y3*U3+
00304      &            KYEL*X3*Y2*U3-2.D0*KYEL*X3*X2*V3-
00305      &            2.D0*KXEL*Y3*Y2*U3+KXEL*Y3*X2*V3+
00306      &            4.D0*KXEL*Y2*X3*V2-6.D0*KXEL*Y2*X2*V2+
00307      &            4.D0*KYEL*X2*Y3*U2+2.D0*KYEL*X2**2*V1+
00308      &            4.D0*KYEL*X3**2*V1+2.D0*KYEL*X2**2*V3+
00309      &            2.D0*KXEL*Y2**2*U3+6.D0*KYEL*X2**2*V2+
00310      &            6.D0*KXEL*Y2**2*U2+4.D0*KXEL*Y3**2*U1-
00311      &            6.D0*KYEL*X2*Y2*U2-8.D0*KXEL*Y3*Y2*U2-
00312      &            4.D0*KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V2+
00313      &            4.D0*KYEL*X3*Y2*U2-8.D0*KYEL*X3*X2*V2-
00314      &            2.D0*KYEL*X2*Y2*U1-2.D0*KXEL*Y2*X2*V1+
00315      &            4.D0*KYEL*X3**2*V2+2.D0*KXEL*Y2**2*U1+
00316      &            4.D0*KXEL*Y3**2*U2+2.D0*KXEL*Y3**2*U3+
00317      &            2.D0*KYEL*X3**2*V3-2.D0*KYEL*X3*Y3*U3-
00318      &            2.D0*KXEL*Y3*X3*V3)*
00319      &             XMUL/(15.D0*(X2*Y3-Y2*X3)/2.D0)
00320 !
00321       A45(IELEM)=(-KYEL*X3*Y3*U2+KYEL*X3*Y3*U1-KYEL*X3*X2*V1+
00322      &            KXEL*Y3*X3*V1-KXEL*Y3*Y2*U1-2.D0*KYEL*X2*Y2*U3+
00323      &            2.D0*KXEL*Y2*X3*V3-2.D0*KXEL*Y2*X2*V3+
00324      &            2.D0*KYEL*X2*Y3*U3+KYEL*X3*Y2*U3-3.D0*KYEL*X3*X2*V3-
00325      &            3.D0*KXEL*Y3*Y2*U3+KXEL*Y3*X2*V3+2.D0*KXEL*Y2*X3*V2-
00326      &            6.D0*KXEL*Y2*X2*V2+2.D0*KYEL*X2*Y3*U2+
00327      &            2.D0*KYEL*X2**2*V1-KYEL*X3**2*V1+2.D0*KYEL*X2**2*V3+
00328      &            2.D0*KXEL*Y2**2*U3+6.D0*KYEL*X2**2*V2+
00329      &            6.D0*KXEL*Y2**2*U2-KXEL*Y3**2*U1-6.D0*KYEL*X2*Y2*U2-
00330      &            6.D0*KXEL*Y3*Y2*U2-KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V2+
00331      &            4.D0*KYEL*X3*Y2*U2-6.D0*KYEL*X3*X2*V2-
00332      &            2.D0*KYEL*X2*Y2*U1-2.D0*KXEL*Y2*X2*V1+KYEL*X3**2*V2+
00333      &            2.D0*KXEL*Y2**2*U1+KXEL*Y3**2*U2+KYEL*X2*Y3*U1+
00334      &            KXEL*Y2*X3*V1)*XMUL/(15.D0*(X3*Y2-Y3*X2)/2.D0)
00335 !
00336       A46(IELEM)=(KYEL*X3*Y3*U2-KYEL*X3*Y3*U1+4.D0*KYEL*X3*X2*V1-
00337      &            3.D0*KYEL*X3*Y2*U1-3.D0*KXEL*Y3*X2*V1-KXEL*Y3*X3*V1+
00338      &            4.D0*KXEL*Y3*Y2*U1+KYEL*X2*Y2*U3-
00339      &            2.D0*KXEL*Y2*X3*V3+KXEL*Y2*X2*V3-
00340      &            2.D0*KYEL*X2*Y3*U3-KYEL*X3*Y2*U3+
00341      &            3.D0*KYEL*X3*X2*V3+3.D0*KXEL*Y3*Y2*U3-KXEL*Y3*X2*V3-
00342      &            2.D0*KXEL*Y2*X3*V2-2.D0*KYEL*X2*Y3*U2+KYEL*X2**2*V1+
00343      &            KYEL*X3**2*V1-KYEL*X2**2*V3-KXEL*Y2**2*U3+
00344      &            KXEL*Y3**2*U1+3.D0*KXEL*Y3*Y2*U2+KXEL*Y3*X3*V2-
00345      &            KXEL*Y3*X2*V2-KYEL*X3*Y2*U2+3.D0*KYEL*X3*X2*V2-
00346      &            KYEL*X2*Y2*U1-KXEL*Y2*X2*V1-KYEL*X3**2*V2+
00347      &            KXEL*Y2**2*U1-KXEL*Y3**2*U2-KYEL*X2*Y3*U1-
00348      &            KXEL*Y2*X3*V1)*XMUL/(15.D0*(X3*Y2-Y3*X2)/2.D0)
00349 !
00350       A52(IELEM)=(KXEL*Y3**2*U1+3.D0*KYEL*X3*Y3*U2-KYEL*X3*Y3*U1+
00351      &            3.D0*KYEL*X3*X2*V1-KXEL*Y3*X3*V1+3.D0*KXEL*Y3*Y2*U1+
00352      &            KYEL*X3**2*V1-3.D0*KXEL*Y2*X3*V3-3.D0*KYEL*X2*Y3*U3-
00353      &            2.D0*KYEL*X3*Y3*U3+3.D0*KYEL*X3*X2*V3+
00354      &            3.D0*KXEL*Y3*Y2*U3-2.D0*KXEL*Y3*X3*V3+
00355      &            2.D0*KXEL*Y3**2*U3+2.D0*KYEL*X3**2*V3-
00356      &            14.D0*KXEL*Y2*X3*V2-14.D0*KYEL*X2*Y3*U2+
00357      &            14.D0*KXEL*Y3*Y2*U2+3.D0*KXEL*Y3*X3*V2+
00358      &            14.D0*KYEL*X3*X2*V2-3.D0*KYEL*X3**2*V2-
00359      &            3.D0*KXEL*Y3**2*U2-3.D0*KYEL*X2*Y3*U1-
00360      &            3.D0*KXEL*Y2*X3*V1)*XMUL/((X3*Y2-Y3*X2)*30.D0)
00361 !
00362       A53(IELEM)=(3.D0*KYEL*X3*X2*V1-3.D0*KYEL*X3*Y2*U1-
00363      &            3.D0*KXEL*Y3*X2*V1+3.D0*KXEL*Y3*Y2*U1+
00364      &            KYEL*X2**2*V1+3.D0*KYEL*X2*Y2*U3+3.D0*KXEL*Y2*X2*V3-
00365      &            3.D0*KYEL*X2**2*V3-3.D0*KXEL*Y2**2*U3-
00366      &            14.D0*KYEL*X3*Y2*U3+14.D0*KYEL*X3*X2*V3+
00367      &            14.D0*KXEL*Y3*Y2*U3-14.D0*KXEL*Y3*X2*V3-
00368      &            2.D0*KXEL*Y2*X2*V2+2.D0*KYEL*X2**2*V2+
00369      &            2.D0*KXEL*Y2**2*U2-2.D0*KYEL*X2*Y2*U2+
00370      &            3.D0*KXEL*Y3*Y2*U2-3.D0*KXEL*Y3*X2*V2-
00371      &            3.D0*KYEL*X3*Y2*U2+3.D0*KYEL*X3*X2*V2+
00372      &            KXEL*Y2**2*U1-KYEL*X2*Y2*U1-KXEL*Y2*X2*V1)*
00373      &            XMUL/((X3*Y2-Y3*X2)*30.D0)
00374 !
00375       A54(IELEM)=(-KYEL*X3*Y3*U2+KYEL*X3*Y3*U1-KYEL*X3*X2*V1+
00376      &            KYEL*X3*Y2*U1+KXEL*Y3*X2*V1+KXEL*Y3*X3*V1-
00377      &            KXEL*Y3*Y2*U1-2.D0*KYEL*X2*Y2*U3+KXEL*Y2*X3*V3-
00378      &            2.D0*KXEL*Y2*X2*V3+KYEL*X2*Y3*U3+2.D0*KYEL*X3*Y2*U3-
00379      &            3.D0*KYEL*X3*X2*V3-3.D0*KXEL*Y3*Y2*U3+
00380      &            2.D0*KXEL*Y3*X2*V3+4.D0*KXEL*Y2*X3*V2-
00381      &            6.D0*KXEL*Y2*X2*V2+4.D0*KYEL*X2*Y3*U2+
00382      &            2.D0*KYEL*X2**2*V1-KYEL*X3**2*V1+2.D0*KYEL*X2**2*V3+
00383      &            2.D0*KXEL*Y2**2*U3+6.D0*KYEL*X2**2*V2+
00384      &            6.D0*KXEL*Y2**2*U2-KXEL*Y3**2*U1-6.D0*KYEL*X2*Y2*U2-
00385      &            6.D0*KXEL*Y3*Y2*U2-KXEL*Y3*X3*V2+2.D0*KXEL*Y3*X2*V2+
00386      &            2.D0*KYEL*X3*Y2*U2-6.D0*KYEL*X3*X2*V2-
00387      &            2.D0*KYEL*X2*Y2*U1-2.D0*KXEL*Y2*X2*V1+
00388      &            KYEL*X3**2*V2+2.D0*KXEL*Y2**2*U1+KXEL*Y3**2*U2)
00389      &            *XMUL/(15.D0*(X3*Y2-Y3*X2)/2.D0)
00390 !
00391       A56(IELEM)=(2.D0*KYEL*X3*Y3*U2+2.D0*KYEL*X3*Y3*U1+KYEL*X3*X2*V1+
00392      &            2.D0*KXEL*Y3*X3*V1+KXEL*Y3*Y2*U1+KYEL*X2*Y2*U3-
00393      &            2.D0*KXEL*Y2*X3*V3+KXEL*Y2*X2*V3-2.D0*KYEL*X2*Y3*U3-
00394      &            4.D0*KYEL*X3*Y2*U3+6.D0*KYEL*X3*X2*V3+
00395      &            6.D0*KXEL*Y3*Y2*U3-4.D0*KXEL*Y3*X2*V3-
00396      &            2.D0*KXEL*Y2*X3*V2-2.D0*KYEL*X2*Y3*U2+KYEL*X2**2*V1-
00397      &            2.D0*KYEL*X3**2*V1-KYEL*X2**2*V3-KXEL*Y2**2*U3-
00398      &            2.D0*KXEL*Y3**2*U1+3.D0*KXEL*Y3*Y2*U2+
00399      &            2.D0*KXEL*Y3*X3*V2-KXEL*Y3*X2*V2-KYEL*X3*Y2*U2+
00400      &            3.D0*KYEL*X3*X2*V2-KYEL*X2*Y2*U1-KXEL*Y2*X2*V1-
00401      &            2.D0*KYEL*X3**2*V2+KXEL*Y2**2*U1-2.D0*KXEL*Y3**2*U2-
00402      &            6.D0*KXEL*Y3**2*U3-6.D0*KYEL*X3**2*V3+
00403      &            6.D0*KYEL*X3*Y3*U3+6.D0*KXEL*Y3*X3*V3-
00404      &            KYEL*X2*Y3*U1-KXEL*Y2*X3*V1)*
00405      &            XMUL/(15.D0*(X2*Y3-Y2*X3)/2.D0)
00406 !
00407       A62(IELEM)=(KXEL*Y3**2*U1+3.D0*KYEL*X3*Y3*U2-KYEL*X3*Y3*U1+
00408      &            KYEL*X3*X2*V1-KXEL*Y3*X3*V1+KXEL*Y3*Y2*U1+
00409      &            KYEL*X3**2*V1+KXEL*Y2*X3*V3+KYEL*X2*Y3*U3-
00410      &            2.D0*KYEL*X3*Y3*U3-KYEL*X3*X2*V3-KXEL*Y3*Y2*U3-
00411      &            2.D0*KXEL*Y3*X3*V3+2.D0*KXEL*Y3**2*U3+
00412      &            2.D0*KYEL*X3**2*V3+3.D0*KXEL*Y3*X3*V2-
00413      &            3.D0*KYEL*X3**2*V2-3.D0*KXEL*Y3**2*U2-
00414      &            KYEL*X2*Y3*U1-KXEL*Y2*X3*V1)
00415      &            *XMUL/((X2*Y3-Y2*X3)*30.D0)
00416 !
00417       A63(IELEM)=(-3.D0*KYEL*X3*X2*V1+3.D0*KYEL*X3*Y2*U1+
00418      &            3.D0*KXEL*Y3*X2*V1-3.D0*KXEL*Y3*Y2*U1+
00419      &            5.D0*KYEL*X2**2*V1-11.D0*KYEL*X2*Y2*U3-
00420      &            11.D0*KXEL*Y2*X2*V3+11.D0*KYEL*X2**2*V3+
00421      &            11.D0*KXEL*Y2**2*U3+14.D0*KYEL*X3*Y2*U3-
00422      &            14.D0*KYEL*X3*X2*V3-14.D0*KXEL*Y3*Y2*U3+
00423      &            14.D0*KXEL*Y3*X2*V3-4.D0*KXEL*Y2*X2*V2+
00424      &            4.D0*KYEL*X2**2*V2+4.D0*KXEL*Y2**2*U2-
00425      &            4.D0*KYEL*X2*Y2*U2-3.D0*KXEL*Y3*Y2*U2+
00426      &            3.D0*KXEL*Y3*X2*V2+3.D0*KYEL*X3*Y2*U2-
00427      &            3.D0*KYEL*X3*X2*V2+5.D0*KXEL*Y2**2*U1-
00428      &            5.D0*KYEL*X2*Y2*U1-5.D0*KXEL*Y2*X2*V1)*
00429      &            XMUL/((X3*Y2-Y3*X2)*30.D0)
00430 !
00431       A64(IELEM)=(KYEL*X3*Y3*U2-KYEL*X3*Y3*U1+4.D0*KYEL*X3*X2*V1-
00432      &            KYEL*X3*Y2*U1-KXEL*Y3*X2*V1-KXEL*Y3*X3*V1+
00433      &            4.D0*KXEL*Y3*Y2*U1+KYEL*X2*Y2*U3-KXEL*Y2*X3*V3+
00434      &            KXEL*Y2*X2*V3-KYEL*X2*Y3*U3-2.D0*KYEL*X3*Y2*U3+
00435      &            3.D0*KYEL*X3*X2*V3+3.D0*KXEL*Y3*Y2*U3-
00436      &            2.D0*KXEL*Y3*X2*V3-KXEL*Y2*X3*V2-KYEL*X2*Y3*U2+
00437      &            KYEL*X2**2*V1+KYEL*X3**2*V1-KYEL*X2**2*V3-
00438      &            KXEL*Y2**2*U3+KXEL*Y3**2*U1+3.D0*KXEL*Y3*Y2*U2+
00439      &            KXEL*Y3*X3*V2-2.D0*KXEL*Y3*X2*V2-
00440      &            2.D0*KYEL*X3*Y2*U2+3.D0*KYEL*X3*X2*V2-
00441      &            KYEL*X2*Y2*U1-KXEL*Y2*X2*V1-KYEL*X3**2*V2+
00442      &            KXEL*Y2**2*U1-KXEL*Y3**2*U2-3.D0*KYEL*X2*Y3*U1-
00443      &            3.D0*KXEL*Y2*X3*V1)*
00444      &            XMUL/(15.D0*(X3*Y2-Y3*X2)/2.D0)
00445 !
00446       A65(IELEM)=(2.D0*KYEL*X3*Y3*U2+2.D0*KYEL*X3*Y3*U1+KYEL*X3*X2*V1-
00447      &            KYEL*X3*Y2*U1-KXEL*Y3*X2*V1+2.D0*KXEL*Y3*X3*V1+
00448      &            KXEL*Y3*Y2*U1+KYEL*X2*Y2*U3-4.D0*KXEL*Y2*X3*V3+
00449      &            KXEL*Y2*X2*V3-4.D0*KYEL*X2*Y3*U3-2.D0*KYEL*X3*Y2*U3+
00450      &            6.D0*KYEL*X3*X2*V3+6.D0*KXEL*Y3*Y2*U3-
00451      &            2.D0*KXEL*Y3*X2*V3-KXEL*Y2*X3*V2-KYEL*X2*Y3*U2+
00452      &            KYEL*X2**2*V1-2.D0*KYEL*X3**2*V1-KYEL*X2**2*V3-
00453      &            KXEL*Y2**2*U3-2.D0*KXEL*Y3**2*U1+3.D0*KXEL*Y3*Y2*U2+
00454      &            2.D0*KXEL*Y3*X3*V2-2.D0*KXEL*Y3*X2*V2-
00455      &            2.D0*KYEL*X3*Y2*U2+3.D0*KYEL*X3*X2*V2-KYEL*X2*Y2*U1-
00456      &            KXEL*Y2*X2*V1-2.D0*KYEL*X3**2*V2+KXEL*Y2**2*U1-
00457      &            2.D0*KXEL*Y3**2*U2-6.D0*KXEL*Y3**2*U3-
00458      &            6.D0*KYEL*X3**2*V3+6.D0*KYEL*X3*Y3*U3+
00459      &            6.D0*KXEL*Y3*X3*V3)*XMUL/(15.D0*(X2*Y3-Y2*X3)/2.D0)
00460 !
00461       A66(IELEM)=(-2.D0*KYEL*X3*Y3*U2-2.D0*KYEL*X3*Y3*U1-
00462      &            2.D0*KXEL*Y3*X3*V1-4.D0*KYEL*X2*Y2*U3+
00463      &            4.D0*KXEL*Y2*X3*V3-4.D0*KXEL*Y2*X2*V3+
00464      &            4.D0*KYEL*X2*Y3*U3+4.D0*KYEL*X3*Y2*U3-
00465      &            8.D0*KYEL*X3*X2*V3-8.D0*KXEL*Y3*Y2*U3+
00466      &            4.D0*KXEL*Y3*X2*V3+KXEL*Y2*X3*V2-
00467      &            2.D0*KXEL*Y2*X2*V2+KYEL*X2*Y3*U2+
00468      &            4.D0*KYEL*X2**2*V1+2.D0*KYEL*X3**2*V1+
00469      &            4.D0*KYEL*X2**2*V3+4.D0*KXEL*Y2**2*U3+
00470      &            2.D0*KYEL*X2**2*V2+2.D0*KXEL*Y2**2*U2+
00471      &            2.D0*KXEL*Y3**2*U1-2.D0*KYEL*X2*Y2*U2-
00472      &            2.D0*KXEL*Y3*Y2*U2-2.D0*KXEL*Y3*X3*V2+
00473      &            KXEL*Y3*X2*V2+KYEL*X3*Y2*U2-2.D0*KYEL*X3*X2*V2-
00474      &            4.D0*KYEL*X2*Y2*U1-4.D0*KXEL*Y2*X2*V1+
00475      &            2.D0*KYEL*X3**2*V2+4.D0*KXEL*Y2**2*U1+
00476      &            2.D0*KXEL*Y3**2*U2+6.D0*KXEL*Y3**2*U3+
00477      &            6.D0*KYEL*X3**2*V3-6.D0*KYEL*X3*Y3*U3-
00478      &            6.D0*KXEL*Y3*X3*V3)*XMUL/(15.D0*(X2*Y3-Y2*X3)/2.D0)
00479 !
00480 ! USES HERE THE 'MAGIC SQUARE' PROPERTIES
00481 ! (SUM OF EACH LINE = SUM OF EACH COLUMN = 0)
00482 !
00483       A16(IELEM) = - A11(IELEM) - A12(IELEM) - A13(IELEM)
00484      &             - A14(IELEM) - A15(IELEM)
00485       A25(IELEM) = - A21(IELEM) - A22(IELEM) - A23(IELEM)
00486      &             - A24(IELEM) - A26(IELEM)
00487       A41(IELEM) = - A42(IELEM) - A43(IELEM) - A44(IELEM)
00488      &             - A45(IELEM) - A46(IELEM)
00489       A61(IELEM) = - A62(IELEM) - A63(IELEM) - A64(IELEM)
00490      &             - A65(IELEM) - A66(IELEM)
00491       A51(IELEM) = - A11(IELEM) - A21(IELEM) - A31(IELEM)
00492      &             - A41(IELEM) - A61(IELEM)
00493       A55(IELEM) = - A51(IELEM) - A52(IELEM) - A53(IELEM)
00494      &             - A54(IELEM) - A56(IELEM)
00495 !
00496       ENDDO ! IELEM
00497 !
00498 !-----------------------------------------------------------------------
00499 ! CASE WHERE U IS OF P2 DISCRETISATION
00500 !-----------------------------------------------------------------------
00501 !
00502       ELSEIF(IELMF.EQ.10.AND.IELMG.EQ.10.AND.
00503      &       IELMU.EQ.13.AND.IELMV.EQ.13) THEN
00504 !
00505 !   LOOP ON THE ELEMENTS
00506 !
00507       DO IELEM = 1 , NELEM
00508 !
00509       X2  =   XEL(IELEM,2)
00510       X3  =   XEL(IELEM,3)
00511       Y2  =   YEL(IELEM,2)
00512       Y3  =   YEL(IELEM,3)
00513 !
00514       U1   =  U(IKLE1(IELEM))
00515       U2   =  U(IKLE2(IELEM))
00516       U3   =  U(IKLE3(IELEM))
00517       U4   =  U(IKLE4(IELEM))
00518       U5   =  U(IKLE5(IELEM))
00519       U6   =  U(IKLE6(IELEM))
00520       V1   =  V(IKLE1(IELEM))
00521       V2   =  V(IKLE2(IELEM))
00522       V3   =  V(IKLE3(IELEM))
00523       V4   =  V(IKLE4(IELEM))
00524       V5   =  V(IKLE5(IELEM))
00525       V6   =  V(IKLE6(IELEM))
00526 !
00527       KXEL =  F(IELEM)
00528       KYEL =  G(IELEM)
00529 !
00530 !  COMPUTES 30 OF THE 36 TERMS (SELECTED AMONG THE LEAST COMPLEX)
00531 !
00532       A11(IELEM)=(KXEL*(Y2-Y3)+KYEL*(X3-X2))*(12.D0*X3*V1+
00533      &            12.D0*Y2*U1-12.D0*Y3*U1-12.D0*X2*V1+2.D0*Y3*U3+
00534      &            2.D0*X2*V3+15.D0*X3*V6-2.D0*Y2*U3+15.D0*Y2*U6-
00535      &            15.D0*X2*V6-15.D0*Y3*U6-2.D0*X3*V3+15.D0*Y2*U4-
00536      &            7.D0*Y3*U5+15.D0*X3*V4-15.D0*X2*V4-15.D0*Y3*U4+
00537      &            7.D0*Y2*U5+7.D0*X3*V5-7.D0*X2*V5-2.D0*Y2*U2+
00538      &            2.D0*X2*V2-2.D0*X3*V2+2.D0*Y3*U2)*
00539      &            XMUL/(X2*Y3-Y2*X3)/90.D0
00540 !
00541       A12(IELEM)=(KXEL*(Y3-Y2)+KYEL*(X2-X3))*(3.D0*X3*V1-
00542      &            3.D0*Y3*U1+2.D0*Y3*U3+5.D0*X3*V6-5.D0*Y3*U6-
00543      &            2.D0*X3*V3-5.D0*Y3*U5+X3*V4-Y3*U4+5.D0*X3*V5+
00544      &            3.D0*X3*V2-3.D0*Y3*U2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00545 !
00546       A13(IELEM)=(KXEL*(Y3-Y2)+KYEL*(X2-X3))*(3.D0*Y2*U1-
00547      &            3.D0*X2*V1-3.D0*X2*V3+3.D0*Y2*U3+Y2*U6-X2*V6+
00548      &            5.D0*Y2*U4-5.D0*X2*V4+5.D0*Y2*U5-5.D0*X2*V5-
00549      &            2.D0*Y2*U2+2.D0*X2*V2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00550 !
00551       A14(IELEM)=(KXEL*(Y2-Y3)+KYEL*(X3-X2))*(15.D0*X3*V1-
00552      &            3.D0*Y2*U1-15.D0*Y3*U1+3.D0*X2*V1+4.D0*Y3*U3-
00553      &            X2*V3+20.D0*X3*V6+Y2*U3-4.D0*Y2*U6+4.D0*X2*V6-
00554      &            20.D0*Y3*U6-4.D0*X3*V3-8.D0*Y2*U4-12.D0*Y3*U5+
00555      &            16.D0*X3*V4+8.D0*X2*V4-16.D0*Y3*U4+8.D0*Y2*U5+
00556      &            12.D0*X3*V5-8.D0*X2*V5+6.D0*Y2*U2-6.D0*X2*V2+
00557      &            X3*V2-Y3*U2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00558 !
00559       A15(IELEM)=(KXEL*(Y2-Y3)+KYEL*(X3-X2))*(3.D0*X3*V1+
00560      &            3.D0*Y2*U1-3.D0*Y3*U1-3.D0*X2*V1+6.D0*Y3*U3+
00561      &            X2*V3+8.D0*X3*V6-Y2*U3+4.D0*Y2*U6-4.D0*X2*V6-
00562      &            8.D0*Y3*U6-6.D0*X3*V3+8.D0*Y2*U4+8.D0*Y3*U5+
00563      &            4.D0*X3*V4-8.D0*X2*V4-4.D0*Y3*U4-8.D0*Y2*U5-
00564      &            8.D0*X3*V5+8.D0*X2*V5-6.D0*Y2*U2+6.D0*X2*V2-
00565      &            X3*V2+Y3*U2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00566 !
00567       A21(IELEM)=(-KXEL*Y3+KYEL*X3)*(3.D0*X3*V1+3.D0*Y2*U1-
00568      &            3.D0*Y3*U1-3.D0*X2*V1+2.D0*Y3*U3+2.D0*X2*V3+
00569      &            5.D0*X3*V6-2.D0*Y2*U3+5.D0*Y2*U6-5.D0*X2*V6-
00570      &            5.D0*Y3*U6-2.D0*X3*V3+Y2*U4-5.D0*Y3*U5+
00571      &            X3*V4-X2*V4-Y3*U4+5.D0*Y2*U5+5.D0*X3*V5-
00572      &            5.D0*X2*V5+3.D0*Y2*U2-3.D0*X2*V2+3.D0*X3*V2-
00573      &            3.D0*Y3*U2)*XMUL/(X2*Y3-Y2*X3)/90.D0
00574 !
00575       A22(IELEM)=(-KXEL*Y3+KYEL*X3)*(-2.D0*X3*V1+2.D0*Y3*U1+
00576      &            2.D0*Y3*U3+7.D0*X3*V6-7.D0*Y3*U6-2.D0*X3*V3-
00577      &            15.D0*Y3*U5+15.D0*X3*V4-15.D0*Y3*U4+
00578      &            15.D0*X3*V5+12.D0*X3*V2-12.D0*Y3*U2)
00579      &            *XMUL/(X2*Y3-Y2*X3)/90.D0
00580 !
00581       A23(IELEM)=(-KXEL*Y3+KYEL*X3)*(2.D0*Y2*U1-2.D0*X2*V1+
00582      &            3.D0*X2*V3-3.D0*Y2*U3-5.D0*Y2*U6+5.D0*X2*V6-
00583      &            5.D0*Y2*U4+5.D0*X2*V4-Y2*U5+X2*V5-3.D0*Y2*U2+
00584      &            3.D0*X2*V2)*XMUL/(X2*Y3-Y2*X3)/90.D0
00585 !
00586       A24(IELEM)=(-KXEL*Y3+KYEL*X3)*(X3*V1-5.D0*Y2*U1-Y3*U1+
00587      &            5.D0*X2*V1+4.D0*Y3*U3+5.D0*X2*V3+12.D0*X3*V6-
00588      &            5.D0*Y2*U3+4.D0*Y2*U6-4.D0*X2*V6-12.D0*Y3*U6-
00589      &            4.D0*X3*V3+24.D0*Y2*U4-20.D0*Y3*U5+16.D0*X3*V4-
00590      &            24.D0*X2*V4-16.D0*Y3*U4+24.D0*Y2*U5+20.D0*X3*V5-
00591      &            24.D0*X2*V5+18.D0*Y2*U2-18.D0*X2*V2+15.D0*X3*V2-
00592      &            15.D0*Y3*U2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00593 !
00594       A26(IELEM)=(-KXEL*Y3+KYEL*X3)*(-X3*V1+5.D0*Y2*U1+Y3*U1-
00595      &            5.D0*X2*V1+6.D0*Y3*U3+5.D0*X2*V3-8.D0*X3*V6-
00596      &            5.D0*Y2*U3+8.D0*Y3*U6-6.D0*X3*V3-4.D0*Y2*U4-
00597      &            8.D0*Y3*U5+4.D0*X3*V4+4.D0*X2*V4-4.D0*Y3*U4+
00598      &            4.D0*Y2*U5+8.D0*X3*V5-4.D0*X2*V5+3.D0*X3*V2-
00599      &            3.D0*Y3*U2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00600 !
00601       A31(IELEM)=(KXEL*Y2-KYEL*X2)*(3.D0*X3*V1+3.D0*Y2*U1-
00602      &            3.D0*Y3*U1-3.D0*X2*V1-3.D0*Y3*U3-3.D0*X2*V3+
00603      &            X3*V6+3.D0*Y2*U3+Y2*U6-X2*V6-Y3*U6+3.D0*X3*V3+
00604      &            5.D0*Y2*U4-5.D0*Y3*U5+5.D0*X3*V4-5.D0*X2*V4-
00605      &            5.D0*Y3*U4+5.D0*Y2*U5+5.D0*X3*V5-5.D0*X2*V5-
00606      &            2.D0*Y2*U2+2.D0*X2*V2-2.D0*X3*V2+2.D0*Y3*U2)
00607      &            *XMUL/(X2*Y3-Y2*X3)/90.D0
00608 !
00609       A32(IELEM)=(KXEL*Y2-KYEL*X2)*(-2.D0*X3*V1+2.D0*Y3*U1-
00610      &            3.D0*Y3*U3+5.D0*X3*V6-5.D0*Y3*U6+3.D0*X3*V3-
00611      &            Y3*U5+5.D0*X3*V4-5.D0*Y3*U4+X3*V5+3.D0*X3*V2-
00612      &            3.D0*Y3*U2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00613 !
00614       A33(IELEM)=(KXEL*Y2-KYEL*X2)*(2.D0*Y2*U1-2.D0*X2*V1+
00615      &            12.D0*X2*V3-12.D0*Y2*U3-15.D0*Y2*U6+15.D0*X2*V6-
00616      &            7.D0*Y2*U4+7.D0*X2*V4-15.D0*Y2*U5+15.D0*X2*V5+
00617      &            2.D0*Y2*U2-2.D0*X2*V2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00618 !
00619       A34(IELEM)=(KXEL*Y2-KYEL*X2)*(5.D0*X3*V1-Y2*U1-5.D0*Y3*U1+
00620      &            X2*V1-3.D0*X2*V3-4.D0*X3*V6+3.D0*Y2*U3+
00621      &            4.D0*Y2*U6-4.D0*X2*V6+4.D0*Y3*U6-8.D0*Y2*U4-
00622      &            4.D0*Y3*U5+8.D0*X2*V4+8.D0*Y2*U5+4.D0*X3*V5-
00623      &            8.D0*X2*V5-6.D0*Y2*U2+6.D0*X2*V2-5.D0*X3*V2+
00624      &            5.D0*Y3*U2)*XMUL/(-X2*Y3+Y2*X3)/90.D0
00625 !
00626       A35(IELEM)=(KXEL*Y2-KYEL*X2)*(-5.D0*X3*V1-Y2*U1+5.D0*Y3*U1+
00627      &            X2*V1-18.D0*Y3*U3-3.D0*X2*V3+24.D0*X3*V6+
00628      &            3.D0*Y2*U3+4.D0*Y2*U6-4.D0*X2*V6-24.D0*Y3*U6+
00629      &            18.D0*X3*V3-8.D0*Y2*U4-24.D0*Y3*U5+4.D0*X3*V4+
00630      &            8.D0*X2*V4-4.D0*Y3*U4+8.D0*Y2*U5+24.D0*X3*V5-
00631      &            8.D0*X2*V5-6.D0*Y2*U2+6.D0*X2*V2-5.D0*X3*V2+
00632      &            5.D0*Y3*U2)*XMUL/(X2*Y3-Y2*X3)/90.D0
00633 !
00634       ANS1 = KYEL*X3*Y2*U2-6.D0*KYEL*X2*Y2*U2-6.D0*KXEL*Y2*X2*V2-
00635      &       4.D0*KYEL*X2**2*V6+15.D0*KXEL*Y3**2*U1-
00636      &       12.D0*KXEL*Y3*Y2*U1+15.D0*KXEL*Y3*X2*V1+
00637      &       6.D0*KYEL*X2**2*V2-KXEL*Y3*X3*V2+16.D0*KXEL*Y3*X2*V4+
00638      &       6.D0*KXEL*Y2**2*U2+16.D0*KXEL*Y3**2*U4-
00639      &       8.D0*KXEL*Y3*Y2*U4-4.D0*KXEL*Y2**2*U6+
00640      &       15.D0*KYEL*X3*Y2*U1-15.D0*KYEL*X3*Y3*U1-
00641      &       12.D0*KYEL*X3*X2*V1-16.D0*KXEL*Y3*X3*V4-
00642      &       7.D0*KXEL*Y3*Y2*U2+20.D0*KXEL*Y3**2*U6+
00643      &       12.D0*KYEL*X3**2*V5+15.D0*KYEL*X3**2*V1-
00644      &       15.D0*KXEL*Y3*X3*V1+8.D0*KYEL*X2*Y3*U5-
00645      &       8.D0*KYEL*X2*Y2*U5+8.D0*KYEL*X2*Y2*U4+
00646      &       KXEL*Y3**2*U2+KYEL*X3**2*V2+20.D0*KYEL*X3**2*V6-
00647      &       4.D0*KYEL*X3**2*V3-4.D0*KXEL*Y3**2*U3-
00648      &       8.D0*KXEL*Y2**2*U4+KYEL*X2**2*V3+8.D0*KXEL*Y2**2*U5+
00649      &       8.D0*KYEL*X2**2*V5+3.D0*KXEL*Y2*X2*V1-
00650      &       3.D0*KXEL*Y2*X3*V1-KXEL*Y2*X2*V3-3.D0*KXEL*Y2**2*U1+
00651      &       KXEL*Y2**2*U3-3.D0*KYEL*X2**2*V1-8.D0*KXEL*Y2*X2*V5+
00652      &       8.D0*KXEL*Y2*X3*V5+8.D0*KXEL*Y2*X2*V4-KYEL*X3*Y3*U2
00653       ANS2 = 4.D0*KYEL*X3*Y3*U3+3.D0*KYEL*X3*X2*V3+
00654      &       20.D0*KYEL*X3*Y2*U6-16.D0*KYEL*X3*X2*V6-
00655      &       20.D0*KYEL*X3*Y3*U6+3.D0*KXEL*Y3*Y2*U3-
00656      &       20.D0*KXEL*Y3*X3*V6+4.D0*KXEL*Y3*X3*V3-
00657      &       16.D0*KXEL*Y3*Y2*U6+12.D0*KXEL*Y3**2*U5+
00658      &       16.D0*KYEL*X3**2*V4-4.D0*KXEL*Y2*X3*V6+
00659      &       3.D0*KYEL*X2*Y2*U1-3.D0*KYEL*X2*Y3*U1-KYEL*X2*Y2*U3+
00660      &       KYEL*X2*Y3*U3-4.D0*KYEL*X2*Y3*U6+KXEL*Y2*X3*V3
00661      &       -8.D0*KYEL*X2**2*V4-7.D0*KYEL*X3*X2*V2-
00662      &       16.D0*KYEL*X3*Y3*U4+16.D0*KYEL*X3*Y2*U4+
00663      &       12.D0*KXEL*Y3*X2*V5-20.D0*KXEL*Y3*Y2*U5-
00664      &       12.D0*KXEL*Y3*X3*V5-8.D0*KYEL*X3*X2*V4+
00665      &       20.D0*KXEL*Y3*X2*V6-20.D0*KYEL*X3*X2*V5-
00666      &       12.D0*KYEL*X3*Y3*U5+12.D0*KYEL*X3*Y2*U5+
00667      &       6.D0*KYEL*X2*Y3*U2+6.D0*KXEL*Y2*X3*V2+KXEL*Y3*X2*V2-
00668      &       8.D0*KYEL*X2*Y3*U4-8.D0*KXEL*Y2*X3*V4+
00669      &       4.D0*KXEL*Y2*X2*V6-4.D0*KYEL*X3*Y2*U3+
00670      &       4.D0*KYEL*X2*Y2*U6-4.D0*KXEL*Y3*X2*V3
00671       A41(IELEM)= (ANS1+ANS2)*XMUL/(90.D0*(-X2*Y3+Y2*X3))
00672 !
00673       ANS1 = KXEL*Y3**2*U1+5.D0*KXEL*Y3*Y2*U1-15.D0*KXEL*Y3*X3*V2+
00674      &       16.D0*KXEL*Y3**2*U4-24.D0*KXEL*Y3*Y2*U4-KYEL*X3*Y3*U1+
00675      &       5.D0*KYEL*X3*X2*V1-16.D0*KXEL*Y3*X3*V4-
00676      &       18.D0*KXEL*Y3*Y2*U2+12.D0*KXEL*Y3**2*U6+
00677      &       20.D0*KYEL*X3**2*V5+KYEL*X3**2*V1-KXEL*Y3*X3*V1+
00678      &       24.D0*KYEL*X2*Y3*U5+15.D0*KXEL*Y3**2*U2+
00679      &       15.D0*KYEL*X3**2*V2+12.D0*KYEL*X3**2*V6-
00680      &       4.D0*KYEL*X3**2*V3-4.D0*KXEL*Y3**2*U3-
00681      &       5.D0*KXEL*Y2*X3*V1+24.D0*KXEL*Y2*X3*V5-
00682      &       15.D0*KYEL*X3*Y3*U2+4.D0*KYEL*X3*Y3*U3
00683       ANS2 = 5.D0*KYEL*X3*X2*V3-4.D0*KYEL*X3*X2*V6-
00684      &       12.D0*KYEL*X3*Y3*U6+5.D0*KXEL*Y3*Y2*U3-
00685      &       12.D0*KXEL*Y3*X3*V6+4.D0*KXEL*Y3*X3*V3-
00686      &       4.D0*KXEL*Y3*Y2*U6+20.D0*KXEL*Y3**2*U5+
00687      &       16.D0*KYEL*X3**2*V4+4.D0*KXEL*Y2*X3*V6-
00688      &       5.D0*KYEL*X2*Y3*U1-5.D0*KYEL*X2*Y3*U3+
00689      &       4.D0*KYEL*X2*Y3*U6-5.D0*KXEL*Y2*X3*V3-
00690      &       18.D0*KYEL*X3*X2*V2-16.D0*KYEL*X3*Y3*U4-
00691      &       24.D0*KXEL*Y3*Y2*U5-20.D0*KXEL*Y3*X3*V5-
00692      &       24.D0*KYEL*X3*X2*V4-24.D0*KYEL*X3*X2*V5-
00693      &       20.D0*KYEL*X3*Y3*U5+18.D0*KYEL*X2*Y3*U2+
00694      &       18.D0*KXEL*Y2*X3*V2+24.D0*KYEL*X2*Y3*U4+
00695      &       24.D0*KXEL*Y2*X3*V4
00696       A42(IELEM)= (ANS1+ANS2)*XMUL/(90.D0*(-X2*Y3+Y2*X3))
00697 !
00698       ANS1 = 6.D0*KYEL*X3*Y2*U2-6.D0*KYEL*X2*Y2*U2-
00699      &       6.D0*KXEL*Y2*X2*V2+4.D0*KYEL*X2**2*V6+
00700      &       4.D0*KXEL*Y3**2.D0*U1+4.D0*KXEL*Y3*Y2*U1-
00701      &       2.D0*KXEL*Y3*X2*V1+6.D0*KYEL*X2**2*V2-
00702      &       4.D0*KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V4+
00703      &       6.D0*KXEL*Y2**2*U2+8.D0*KXEL*Y3**2.D0*U4-
00704      &       8.D0*KXEL*Y3*Y2*U4+4.D0*KXEL*Y2**2*U6-
00705      &       2.D0*KYEL*X3*Y2*U1-4.D0*KYEL*X3*Y3*U1+
00706      &       4.D0*KYEL*X3*X2*V1-8.D0*KXEL*Y3*X3*V4-
00707      &       12.D0*KXEL*Y3*Y2*U2+8.D0*KXEL*Y3**2*U6+
00708      &       8.D0*KYEL*X3**2*V5+4.D0*KYEL*X3**2*V1-
00709      &       4.D0*KXEL*Y3*X3*V1+8.D0*KYEL*X2*Y3*U5-
00710      &       12.D0*KYEL*X2*Y2*U5-12.D0*KYEL*X2*Y2*U4+
00711      &       4.D0*KXEL*Y3**2*U2+4.D0*KYEL*X3**2*V2+
00712      &       8.D0*KYEL*X3**2*V6-2.D0*KYEL*X3**2*V3-
00713      &       2.D0*KXEL*Y3**2*U3+12.D0*KXEL*Y2**2*U4-
00714      &       2.D0*KYEL*X2**2*V3+12.D0*KXEL*Y2**2*U5+
00715      &       12.D0*KYEL*X2**2*V5+2.D0*KXEL*Y2*X2*V1-
00716      &       2.D0*KXEL*Y2*X3*V1+2.D0*KXEL*Y2*X2*V3-
00717      &       2.D0*KXEL*Y2**2*U1
00718       ANS2 = -2.D0*KXEL*Y2**2*U3-2.D0*KYEL*X2**2*V1-
00719      &       12.D0*KXEL*Y2*X2*V5+8.D0*KXEL*Y2*X3*V5-
00720      &       12.D0*KXEL*Y2*X2*V4-4.D0*KYEL*X3*Y3*U2+
00721      &       2.D0*KYEL*X3*Y3*U3+2.D0*KYEL*X3*X2*V3-
00722      &       8.D0*KYEL*X3*Y3*U6+2.D0*KXEL*Y3*Y2*U3-
00723      &       8.D0*KXEL*Y3*X3*V6+2.D0*KXEL*Y3*X3*V3+
00724      &       8.D0*KXEL*Y3**2*U5+8.D0*KYEL*X3**2*V4+
00725      &       2.D0*KYEL*X2*Y2*U1-2.D0*KYEL*X2*Y3*U1+
00726      &       2.D0*KYEL*X2*Y2*U3-KYEL*X2*Y3*U3-KXEL*Y2*X3*V3+
00727      &       12.D0*KYEL*X2**2*V4-12.D0*KYEL*X3*X2*V2-
00728      &       8.D0*KYEL*X3*Y3*U4+4.D0*KYEL*X3*Y2*U4+
00729      &       8.D0*KXEL*Y3*X2*V5-16.D0*KXEL*Y3*Y2*U5-
00730      &       8.D0*KXEL*Y3*X3*V5-8.D0*KYEL*X3*X2*V4-
00731      &       16.D0*KYEL*X3*X2*V5-8.D0*KYEL*X3*Y3*U5+
00732      &       8.D0*KYEL*X3*Y2*U5+6.D0*KYEL*X2*Y3*U2+
00733      &       6.D0*KXEL*Y2*X3*V2+6.D0*KXEL*Y3*X2*V2+
00734      &       4.D0*KYEL*X2*Y3*U4+4.D0*KXEL*Y2*X3*V4-
00735      &       4.D0*KXEL*Y2*X2*V6-KYEL*X3*Y2*U3-
00736      &       4.D0*KYEL*X2*Y2*U6-KXEL*Y3*X2*V3
00737       A44(IELEM)= (ANS1+ANS2)*2.D0*XMUL/(45.D0*(X2*Y3-Y2*X3))
00738 !
00739       ANS1 = 6.D0*KYEL*X3*Y2*U2-6.D0*KYEL*X2*Y2*U2-
00740      &       6.D0*KXEL*Y2*X2*V2+4.D0*KYEL*X2**2*V6-
00741      &       KXEL*Y3**2*U1+3.D0*KXEL*Y3*Y2*U1-2.D0*KXEL*Y3*X2*V1+
00742      &       6.D0*KYEL*X2**2*V2-KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V4+
00743      &       6.D0*KXEL*Y2**2*U2-8.D0*KXEL*Y3*Y2*U4+
00744      &       4.D0*KXEL*Y2**2*U6-2.D0*KYEL*X3*Y2*U1+KYEL*X3*Y3*U1+
00745      &       3.D0*KYEL*X3*X2*V1-6.D0*KXEL*Y3*Y2*U2-
00746      &       4.D0*KXEL*Y3**2*U6+4.D0*KYEL*X3**2*V5-
00747      &       KYEL*X3**2*V1+KXEL*Y3*X3*V1+8.D0*KYEL*X2*Y3*U5-
00748      &       12.D0*KYEL*X2*Y2*U5-12.D0*KYEL*X2*Y2*U4+
00749      &       KXEL*Y3**2*U2+KYEL*X3**2*V2-4.D0*KYEL*X3**2*V6+
00750      &       12.D0*KXEL*Y2**2*U4-2.D0*KYEL*X2**2*V3+
00751      &       12.D0*KXEL*Y2**2*U5+12.D0*KYEL*X2**2*V5+
00752      &       2.D0*KXEL*Y2*X2*V1-KXEL*Y2*X3*V1+2.D0*KXEL*Y2*X2*V3-
00753      &       2.D0*KXEL*Y2**2*U1
00754       ANS2 = -2.D0*KXEL*Y2**2*U3-2.D0*KYEL*X2**2*V1-
00755      &       12.D0*KXEL*Y2*X2*V5+8.D0*KXEL*Y2*X3*V5-
00756      &       12.D0*KXEL*Y2*X2*V4-KYEL*X3*Y3*U2+
00757      &       KYEL*X3*X2*V3-4.D0*KYEL*X3*X2*V6+4.D0*KYEL*X3*Y3*U6+
00758      &       KXEL*Y3*Y2*U3+4.D0*KXEL*Y3*X3*V6-4.D0*KXEL*Y3*Y2*U6+
00759      &       4.D0*KXEL*Y3**2*U5+4.D0*KXEL*Y2*X3*V6+
00760      &       2.D0*KYEL*X2*Y2*U1-KYEL*X2*Y3*U1+2.D0*KYEL*X2*Y2*U3+
00761      &       4.D0*KYEL*X2*Y3*U6+12.D0*KYEL*X2**2*V4-
00762      &       6.D0*KYEL*X3*X2*V2+4.D0*KYEL*X3*Y2*U4+
00763      &       8.D0*KXEL*Y3*X2*V5-16.D0*KXEL*Y3*Y2*U5-
00764      &       4.D0*KXEL*Y3*X3*V5-8.D0*KYEL*X3*X2*V4-
00765      &       16.D0*KYEL*X3*X2*V5-4.D0*KYEL*X3*Y3*U5+
00766      &       8.D0*KYEL*X3*Y2*U5+6.D0*KXEL*Y3*X2*V2+
00767      &       4.D0*KYEL*X2*Y3*U4+4.D0*KXEL*Y2*X3*V4-
00768      &       4.D0*KXEL*Y2*X2*V6-KYEL*X3*Y2*U3-
00769      &       4.D0*KYEL*X2*Y2*U6-KXEL*Y3*X2*V3
00770       A45(IELEM)=  (ANS1+ANS2)*2.D0*XMUL/(45.D0*(-X2*Y3+Y2*X3))
00771 !
00772       ANS1 = -KYEL*X3*Y2*U2-KXEL*Y3**2*U1-4.D0*KXEL*Y3*Y2*U1+
00773      &       5.D0*KXEL*Y3*X2*V1-KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V4-
00774      &       8.D0*KXEL*Y3*Y2*U4+5.D0*KYEL*X3*Y2*U1+KYEL*X3*Y3*U1-
00775      &       4.D0*KYEL*X3*X2*V1+KXEL*Y3*Y2*U2-4.D0*KXEL*Y3**2*U6+
00776      &       4.D0*KYEL*X3**2*V5-KYEL*X3**2*V1+KXEL*Y3*X3*V1+
00777      &       8.D0*KYEL*X2*Y3*U5-4.D0*KYEL*X2*Y2*U5+
00778      &       4.D0*KYEL*X2*Y2*U4+KXEL*Y3**2*U2+KYEL*X3**2*V2-
00779      &       4.D0*KYEL*X3**2*V6-4.D0*KXEL*Y2**2*U4+KYEL*X2**2*V3+
00780      &       4.D0*KXEL*Y2**2*U5+4.D0*KYEL*X2**2*V5+KXEL*Y2*X2*V1-
00781      &       KXEL*Y2*X3*V1-KXEL*Y2*X2*V3-KXEL*Y2**2*U1+
00782      &       KXEL*Y2**2*U3-KYEL*X2**2*V1-4.D0*KXEL*Y2*X2*V5
00783       ANS2 = 8.D0*KXEL*Y2*X3*V5+4.D0*KXEL*Y2*X2*V4-KYEL*X3*Y3*U2+
00784      &       KYEL*X3*X2*V3+4.D0*KYEL*X3*Y2*U6-8.D0*KYEL*X3*X2*V6+
00785      &       4.D0*KYEL*X3*Y3*U6+KXEL*Y3*Y2*U3+4.D0*KXEL*Y3*X3*V6-
00786      &       8.D0*KXEL*Y3*Y2*U6+4.D0*KXEL*Y3**2*U5+
00787      &       4.D0*KXEL*Y2*X3*V6+KYEL*X2*Y2*U1-KYEL*X2*Y3*U1-
00788      &       KYEL*X2*Y2*U3+4.D0*KYEL*X2*Y3*U6-4.D0*KYEL*X2**2*V4+
00789      &       KYEL*X3*X2*V2+4.D0*KYEL*X3*Y2*U4+4.D0*KXEL*Y3*X2*V5-
00790      &       12.D0*KXEL*Y3*Y2*U5-4.D0*KXEL*Y3*X3*V5-
00791      &       8.D0*KYEL*X3*X2*V4+4.D0*KXEL*Y3*X2*V6-
00792      &       12.D0*KYEL*X3*X2*V5-4.D0*KYEL*X3*Y3*U5+
00793      &       4.D0*KYEL*X3*Y2*U5-KXEL*Y3*X2*V2+4.D0*KYEL*X2*Y3*U4+
00794      &       4.D0*KXEL*Y2*X3*V4-KYEL*X3*Y2*U3-KXEL*Y3*X2*V3
00795       A46(IELEM)=  (ANS1+ANS2)*2.D0*XMUL/(45.D0*(X2*Y3-Y2*X3))
00796 !
00797       ANS1 = -KYEL*X3*Y2*U2+6.D0*KYEL*X2*Y2*U2+6.D0*KXEL*Y2*X2*V2+
00798      &       4.D0*KYEL*X2**2*V6+3.D0*KXEL*Y3**2*U1-
00799      &       6.D0*KXEL*Y3*Y2*U1+3.D0*KXEL*Y3*X2*V1-
00800      &       6.D0*KYEL*X2**2*V2+KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V4-
00801      &       6.D0*KXEL*Y2**2*U2+4.D0*KXEL*Y3**2*U4-
00802      &       12.D0*KXEL*Y3*Y2*U4+4.D0*KXEL*Y2**2*U6+
00803      &       3.D0*KYEL*X3*Y2*U1-3.D0*KYEL*X3*Y3*U1-
00804      &       6.D0*KYEL*X3*X2*V1-4.D0*KXEL*Y3*X3*V4+
00805      &       7.D0*KXEL*Y3*Y2*U2+8.D0*KXEL*Y3**2*U6-
00806      &       8.D0*KYEL*X3**2*V5+3.D0*KYEL*X3**2*V1-
00807      &       3.D0*KXEL*Y3*X3*V1-8.D0*KYEL*X2*Y3*U5+
00808      &       8.D0*KYEL*X2*Y2*U5-8.D0*KYEL*X2*Y2*U4-
00809      &       KXEL*Y3**2*U2-KYEL*X3**2*V2+8.D0*KYEL*X3**2*V6-
00810      &       6.D0*KYEL*X3**2*V3-6.D0*KXEL*Y3**2*U3+
00811      &       8.D0*KXEL*Y2**2*U4-KYEL*X2**2*V3-8.D0*KXEL*Y2**2*U5-
00812      &       8.D0*KYEL*X2**2*V5-3.D0*KXEL*Y2*X2*V1+
00813      &       3.D0*KXEL*Y2*X3*V1+KXEL*Y2*X2*V3+3.D0*KXEL*Y2**2*U1-
00814      &       KXEL*Y2**2*U3+3.D0*KYEL*X2**2*V1+8.D0*KXEL*Y2*X2*V5
00815       ANS2 = -8.D0*KXEL*Y2*X3*V5-8.D0*KXEL*Y2*X2*V4+
00816      &       KYEL*X3*Y3*U2+6.D0*KYEL*X3*Y3*U3+7.D0*KYEL*X3*X2*V3+
00817      &       8.D0*KYEL*X3*Y2*U6-12.D0*KYEL*X3*X2*V6-
00818      &       8.D0*KYEL*X3*Y3*U6+7.D0*KXEL*Y3*Y2*U3-
00819      &       8.D0*KXEL*Y3*X3*V6+6.D0*KXEL*Y3*X3*V3-
00820      &       12.D0*KXEL*Y3*Y2*U6-8.D0*KXEL*Y3**2*U5+
00821      &       4.D0*KYEL*X3**2*V4+4.D0*KXEL*Y2*X3*V6-
00822      &       3.D0*KYEL*X2*Y2*U1+3.D0*KYEL*X2*Y3*U1+
00823      &       KYEL*X2*Y2*U3-KYEL*X2*Y3*U3+4.D0*KYEL*X2*Y3*U6-
00824      &       KXEL*Y2*X3*V3+8.D0*KYEL*X2**2*V4+7.D0*KYEL*X3*X2*V2-
00825      &       4.D0*KYEL*X3*Y3*U4+4.D0*KYEL*X3*Y2*U4-
00826      &       8.D0*KXEL*Y3*X2*V5+16.D0*KXEL*Y3*Y2*U5+
00827      &       8.D0*KXEL*Y3*X3*V5-12.D0*KYEL*X3*X2*V4+
00828      &       8.D0*KXEL*Y3*X2*V6+16.D0*KYEL*X3*X2*V5+
00829      &       8.D0*KYEL*X3*Y3*U5-8.D0*KYEL*X3*Y2*U5-
00830      &       6.D0*KYEL*X2*Y3*U2-6.D0*KXEL*Y2*X3*V2-
00831      &       KXEL*Y3*X2*V2+8.D0*KYEL*X2*Y3*U4+8.D0*KXEL*Y2*X3*V4-
00832      &       4.D0*KXEL*Y2*X2*V6-6.D0*KYEL*X3*Y2*U3-
00833      &       4.D0*KYEL*X2*Y2*U6-6.D0*KXEL*Y3*X2*V3
00834       A51(IELEM) = (ANS1+ANS2)*XMUL/(90.D0*(-X2*Y3+Y2*X3))
00835 !
00836       ANS1 = -KXEL*Y3**2*U1+5.D0*KXEL*Y3*Y2*U1-3.D0*KXEL*Y3*X3*V2+
00837      &       4.D0*KXEL*Y3**2*U4-24.D0*KXEL*Y3*Y2*U4+KYEL*X3*Y3*U1+
00838      &       5.D0*KYEL*X3*X2*V1-4.D0*KXEL*Y3*X3*V4-
00839      &       18.D0*KXEL*Y3*Y2*U2-8.D0*KXEL*Y3**2*U6+
00840      &       8.D0*KYEL*X3**2*V5-KYEL*X3**2*V1+KXEL*Y3*X3*V1+
00841      &       24.D0*KYEL*X2*Y3*U5+3.D0*KXEL*Y3**2*U2+
00842      &       3.D0*KYEL*X3**2*V2-8.D0*KYEL*X3**2*V6-
00843      &       6.D0*KYEL*X3**2*V3-6.D0*KXEL*Y3**2*U3-
00844      &       5.D0*KXEL*Y2*X3*V1+24.D0*KXEL*Y2*X3*V5-
00845      &       3.D0*KYEL*X3*Y3*U2+6.D0*KYEL*X3*Y3*U3
00846       ANS2 = 5.D0*KYEL*X3*X2*V3 -4.D0*KYEL*X3*X2*V6+
00847      &       8.D0*KYEL*X3*Y3*U6+5.D0*KXEL*Y3*Y2*U3+
00848      &       8.D0*KXEL*Y3*X3*V6+6.D0*KXEL*Y3*X3*V3-
00849      &       4.D0*KXEL*Y3*Y2*U6+8.D0*KXEL*Y3**2*U5+
00850      &       4.D0*KYEL*X3**2*V4+4.D0*KXEL*Y2*X3*V6-
00851      &       5.D0*KYEL*X2*Y3*U1-5.D0*KYEL*X2*Y3*U3+
00852      &       4.D0*KYEL*X2*Y3*U6-5.D0*KXEL*Y2*X3*V3-
00853      &       18.D0*KYEL*X3*X2*V2-4.D0*KYEL*X3*Y3*U4-
00854      &       24.D0*KXEL*Y3*Y2*U5-8.D0*KXEL*Y3*X3*V5-
00855      &       24.D0*KYEL*X3*X2*V4-24.D0*KYEL*X3*X2*V5-
00856      &       8.D0*KYEL*X3*Y3*U5+18.D0*KYEL*X2*Y3*U2+
00857      &       18.D0*KXEL*Y2*X3*V2+24.D0*KYEL*X2*Y3*U4+
00858      &       24.D0*KXEL*Y2*X3*V4
00859       A52(IELEM)= (ANS1+ANS2)*XMUL/(90.D0*(X2*Y3-Y2*X3))
00860 !
00861       ANS1 = 5.D0*KYEL*X3*Y2*U2-6.D0*KYEL*X2*Y2*U2-
00862      &       6.D0*KXEL*Y2*X2*V2-4.D0*KYEL*X2**2*V6-
00863      &       5.D0*KXEL*Y3*Y2*U1+5.D0*KXEL*Y3*X2*V1+
00864      &       6.D0*KYEL*X2**2*V2-4.D0*KXEL*Y3*X2*V4+
00865      &       6.D0*KXEL*Y2**2*U2+4.D0*KXEL*Y3*Y2*U4-
00866      &       4.D0*KXEL*Y2**2*U6+5.D0*KYEL*X3*Y2*U1-
00867      &       5.D0*KYEL*X3*X2*V1-5.D0*KXEL*Y3*Y2*U2+
00868      &       8.D0*KYEL*X2*Y2*U5-8.D0*KYEL*X2*Y2*U4+
00869      &       8.D0*KXEL*Y2**2*U4-3.D0*KYEL*X2**2*V3-
00870      &       8.D0*KXEL*Y2**2*U5-8.D0*KYEL*X2**2*V5-
00871      &       KXEL*Y2*X2*V1+3.D0*KXEL*Y2*X2*V3+KXEL*Y2**2*U1-
00872      &       3.D0*KXEL*Y2**2*U3
00873       ANS2 = KYEL*X2**2*V1+8*KXEL*Y2*X2*V5-8.D0*KXEL*Y2*X2*V4+
00874      &       18.D0*KYEL*X3*X2*V3-24.D0*KYEL*X3*Y2*U6+
00875      &       24.D0*KYEL*X3*X2*V6+18.D0*KXEL*Y3*Y2*U3+
00876      &       24.D0*KXEL*Y3*Y2*U6-KYEL*X2*Y2*U1+3.D0*KYEL*X2*Y2*U3+
00877      &       8.D0*KYEL*X2**2*V4-5.D0*KYEL*X3*X2*V2-
00878      &       4.D0*KYEL*X3*Y2*U4-24.D0*KXEL*Y3*X2*V5+
00879      &       24.D0*KXEL*Y3*Y2*U5+4.D0*KYEL*X3*X2*V4-
00880      &       24.D0*KXEL*Y3*X2*V6+24.D0*KYEL*X3*X2*V5-
00881      &       24.D0*KYEL*X3*Y2*U5+5.D0*KXEL*Y3*X2*V2+
00882      &       4.D0*KXEL*Y2*X2*V6-18.D0*KYEL*X3*Y2*U3+
00883      &       4.D0*KYEL*X2*Y2*U6-18.D0*KXEL*Y3*X2*V3
00884       A53(IELEM)= (ANS1+ANS2)*XMUL/(90.D0*(-X2*Y3+Y2*X3))
00885 !
00886       ANS1 = -6.D0*KYEL*X2*Y2*U2-6.D0*KXEL*Y2*X2*V2+
00887      &       4.D0*KYEL*X2**2*V6-KXEL*Y3**2*U1+3.D0*KXEL*Y3*Y2*U1-
00888      &       KXEL*Y3*X2*V1+6.D0*KYEL*X2**2*V2-KXEL*Y3*X3*V2+
00889      &       4.D0*KXEL*Y3*X2*V4+6.D0*KXEL*Y2**2*U2-
00890      &       8.D0*KXEL*Y3*Y2*U4+4.D0*KXEL*Y2**2*U6-
00891      &       KYEL*X3*Y2*U1+KYEL*X3*Y3*U1+3.D0*KYEL*X3*X2*V1-
00892      &       6.D0*KXEL*Y3*Y2*U2-4.D0*KXEL*Y3**2*U6+
00893      &       4.D0*KYEL*X3**2*V5-KYEL*X3**2*V1+KXEL*Y3*X3*V1+
00894      &       8.D0*KYEL*X2*Y3*U5-12.D0*KYEL*X2*Y2*U5-
00895      &       12.D0*KYEL*X2*Y2*U4+KXEL*Y3**2*U2+KYEL*X3**2*V2-
00896      &       4.D0*KYEL*X3**2*V6+12.D0*KXEL*Y2**2*U4-
00897      &       2.D0*KYEL*X2**2*V3+12.D0*KXEL*Y2**2*U5+
00898      &       12.D0*KYEL*X2**2*V5+2.D0*KXEL*Y2*X2*V1-
00899      &       2.D0*KXEL*Y2*X3*V1+2.D0*KXEL*Y2*X2*V3-
00900      &       2.D0*KXEL*Y2**2*U1-2.D0*KXEL*Y2**2*U3
00901       ANS2 = -2.D0*KYEL*X2**2*V1-12.D0*KXEL*Y2*X2*V5+
00902      &       8.D0*KXEL*Y2*X3*V5-12.D0*KXEL*Y2*X2*V4-
00903      &       KYEL*X3*Y3*U2+KYEL*X3*X2*V3+4.D0*KYEL*X3*Y2*U6-
00904      &       4.D0*KYEL*X3*X2*V6+4.D0*KYEL*X3*Y3*U6+KXEL*Y3*Y2*U3+
00905      &       4.D0*KXEL*Y3*X3*V6-4.D0*KXEL*Y3*Y2*U6+
00906      &       4.D0*KXEL*Y3**2*U5+2.D0*KYEL*X2*Y2*U1-
00907      &       2.D0*KYEL*X2*Y3*U1+2.D0*KYEL*X2*Y2*U3-
00908      &       KYEL*X2*Y3*U3-KXEL*Y2*X3*V3+12.D0*KYEL*X2**2*V4-
00909      &       6.D0*KYEL*X3*X2*V2+4.D0*KYEL*X3*Y2*U4+
00910      &       8.D0*KXEL*Y3*X2*V5-16.D0*KXEL*Y3*Y2*U5-
00911      &       4.D0*KXEL*Y3*X3*V5-8.D0*KYEL*X3*X2*V4+
00912      &       4.D0*KXEL*Y3*X2*V6-16.D0*KYEL*X3*X2*V5-
00913      &       4.D0*KYEL*X3*Y3*U5+8.D0*KYEL*X3*Y2*U5+
00914      &       6.D0*KYEL*X2*Y3*U2+6.D0*KXEL*Y2*X3*V2+
00915      &       4.D0*KYEL*X2*Y3*U4+4.D0*KXEL*Y2*X3*V4-
00916      &       4.D0*KXEL*Y2*X2*V6-4.D0*KYEL*X2*Y2*U6
00917       A54(IELEM)= (ANS1+ANS2)*2.D0*XMUL/(45.D0*(-X2*Y3+Y2*X3))
00918 !
00919       ANS1 = -KYEL*X3*Y2*U2-2.D0*KXEL*Y3**2*U1+3.D0*KXEL*Y3*Y2*U1-
00920      &        2.D0*KXEL*Y3*X2*V1+2.D0*KXEL*Y3*X3*V2+
00921      &        4.D0*KXEL*Y3**2*U4-4.D0*KXEL*Y3*Y2*U4-
00922      &        2.D0*KYEL*X3*Y2*U1+2.D0*KYEL*X3*Y3*U1+
00923      &        3.D0*KYEL*X3*X2*V1-4.D0*KXEL*Y3*X3*V4+
00924      &        KXEL*Y3*Y2*U2+12.D0*KXEL*Y3**2*U6+
00925      &        12.D0*KYEL*X3**2*V5-2.D0*KYEL*X3**2*V1+
00926      &        2.D0*KXEL*Y3*X3*V1+8.D0*KYEL*X2*Y3*U5-
00927      &        4.D0*KYEL*X2*Y2*U5+4.D0*KYEL*X2*Y2*U4-
00928      &        2.D0*KXEL*Y3**2*U2-2.D0*KYEL*X3**2*V2+
00929      &        12.D0*KYEL*X3**2*V6+6.D0*KYEL*X3**2*V3+
00930      &        6.D0*KXEL*Y3**2*U3-4.D0*KXEL*Y2**2*U4+
00931      &        KYEL*X2**2*V3+4.D0*KXEL*Y2**2*U5+4.D0*KYEL*X2**2*V5+
00932      &        KXEL*Y2*X2*V1-KXEL*Y2*X3*V1-KXEL*Y2*X2*V3-
00933      &        KXEL*Y2**2*U1+KXEL*Y2**2*U3-KYEL*X2**2*V1-
00934      &        4.D0*KXEL*Y2*X2*V5
00935       ANS2 = 8.D0*KXEL*Y2*X3*V5+4.D0*KXEL*Y2*X2*V4+
00936      &       2.D0*KYEL*X3*Y3*U2-6.D0*KYEL*X3*Y3*U3-
00937      &       6.D0*KYEL*X3*X2*V3+4.D0*KYEL*X3*Y2*U6-
00938      &       8.D0*KYEL*X3*X2*V6-12.D0*KYEL*X3*Y3*U6-
00939      &       6.D0*KXEL*Y3*Y2*U3-12.D0*KXEL*Y3*X3*V6-
00940      &       6.D0*KXEL*Y3*X3*V3-8.D0*KXEL*Y3*Y2*U6+
00941      &       12.D0*KXEL*Y3**2*U5+4.D0*KYEL*X3**2*V4+
00942      &       4.D0*KXEL*Y2*X3*V6+KYEL*X2*Y2*U1-KYEL*X2*Y3*U1-
00943      &       KYEL*X2*Y2*U3+4.D0*KYEL*X2*Y3*U6-4.D0*KYEL*X2**2*V4+
00944      &       KYEL*X3*X2*V2-4.D0*KYEL*X3*Y3*U4+8.D0*KXEL*Y3*X2*V5-
00945      &       16.D0*KXEL*Y3*Y2*U5-12.D0*KXEL*Y3*X3*V5-
00946      &       4.D0*KYEL*X3*X2*V4+4.D0*KXEL*Y3*X2*V6-
00947      &       16.D0*KYEL*X3*X2*V5-12.D0*KYEL*X3*Y3*U5+
00948      &       8.D0*KYEL*X3*Y2*U5-KXEL*Y3*X2*V2+4.D0*KYEL*X2*Y3*U4+
00949      &       4.D0*KXEL*Y2*X3*V4+6.D0*KYEL*X3*Y2*U3+6.D0*KXEL*Y3*X2*V3
00950       A56(IELEM)= (ANS1+ANS2)*2.D0*XMUL/(45.D0*(-X2*Y3+Y2*X3))
00951 !
00952       ANS1 = -KYEL*X3*Y2*U2-4.D0*KYEL*X2*Y2*U2-4.D0*KXEL*Y2*X2*V2-
00953      &       16.D0*KYEL*X2**2*V6+3.D0*KXEL*Y3**2*U1+
00954      &       12.D0*KXEL*Y3*Y2*U1+3.D0*KXEL*Y3*X2*V1+
00955      &       4.D0*KYEL*X2**2*V2+KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V4+
00956      &       4.D0*KXEL*Y2**2*U2+4.D0*KXEL*Y3**2*U4+
00957      &       16.D0*KXEL*Y3*Y2*U4-16.D0*KXEL*Y2**2*U6+
00958      &       3.D0*KYEL*X3*Y2*U1-3.D0*KYEL*X3*Y3*U1+
00959      &       12.D0*KYEL*X3*X2*V1-4.D0*KXEL*Y3*X3*V4-
00960      &       3.D0*KXEL*Y3*Y2*U2+8.D0*KXEL*Y3**2*U6-
00961      &       8.D0*KYEL*X3**2*V5+3.D0*KYEL*X3**2*V1-
00962      &       3.D0*KXEL*Y3*X3*V1-12.D0*KYEL*X2*Y3*U5+
00963      &       12.D0*KYEL*X2*Y2*U5+20.D0*KYEL*X2*Y2*U4-KXEL*Y3**2*U2-
00964      &       KYEL*X3**2*V2+8.D0*KYEL*X3**2*V6-6.D0*KYEL*X3**2*V3-
00965      &       6.D0*KXEL*Y3**2*U3-20.D0*KXEL*Y2**2*U4-KYEL*X2**2*V3-
00966      &       12.D0*KXEL*Y2**2*U5-12.D0*KYEL*X2**2*V5+
00967      &       15.D0*KXEL*Y2*X2*V1-15.D0*KXEL*Y2*X3*V1+
00968      &       KXEL*Y2*X2*V3-15.D0*KXEL*Y2**2*U1-KXEL*Y2**2*U3-
00969      &       15.D0*KYEL*X2**2*V1+12.D0*KXEL*Y2*X2*V5
00970       ANS2 = -12.D0*KXEL*Y2*X3*V5+20.D0*KXEL*Y2*X2*V4+
00971      &       KYEL*X3*Y3*U2+6.D0*KYEL*X3*Y3*U3+7.D0*KYEL*X3*X2*V3+
00972      &       8.D0*KYEL*X3*Y2*U6+8.D0*KYEL*X3*X2*V6-
00973      &       8.D0*KYEL*X3*Y3*U6+7.D0*KXEL*Y3*Y2*U3-
00974      &       8.D0*KXEL*Y3*X3*V6+6.D0*KXEL*Y3*X3*V3+
00975      &       8.D0*KXEL*Y3*Y2*U6-8.D0*KXEL*Y3**2*U5+
00976      &       4.D0*KYEL*X3**2*V4-16.D0*KXEL*Y2*X3*V6+
00977      &       15.D0*KYEL*X2*Y2*U1-15.D0*KYEL*X2*Y3*U1+
00978      &       KYEL*X2*Y2*U3-KYEL*X2*Y3*U3-16.D0*KYEL*X2*Y3*U6-
00979      &       KXEL*Y2*X3*V3-20.D0*KYEL*X2**2*V4-3.D0*KYEL*X3*X2*V2-
00980      &       4.D0*KYEL*X3*Y3*U4+4.D0*KYEL*X3*Y2*U4-
00981      &       8.D0*KXEL*Y3*X2*V5+20.D0*KXEL*Y3*Y2*U5+
00982      &       8.D0*KXEL*Y3*X3*V5+16.D0*KYEL*X3*X2*V4+
00983      &       8.D0*KXEL*Y3*X2*V6+20.D0*KYEL*X3*X2*V5+
00984      &       8.D0*KYEL*X3*Y3*U5-8.D0*KYEL*X3*Y2*U5+
00985      &       4.D0*KYEL*X2*Y3*U2+4.D0*KXEL*Y2*X3*V2-
00986      &       KXEL*Y3*X2*V2-20.D0*KYEL*X2*Y3*U4-20.D0*KXEL*Y2*X3*V4+
00987      &       16.D0*KXEL*Y2*X2*V6-6.D0*KYEL*X3*Y2*U3+16.D0*KYEL*X2*
00988      &       Y2*U6-6.D0*KXEL*Y3*X2*V3
00989       A61(IELEM)= (ANS1+ANS2)*XMUL/(90.D0*(X2*Y3-Y2*X3))
00990 !
00991       ANS1 = -KXEL*Y3**2*U1-5.D0*KXEL*Y3*Y2*U1-
00992      &       3.D0*KXEL*Y3*X3*V2+4.D0*KXEL*Y3**2*U4+
00993      &       4.D0*KXEL*Y3*Y2*U4+KYEL*X3*Y3*U1-5.D0*KYEL*X3*X2*V1-
00994      &       4.D0*KXEL*Y3*X3*V4-8.D0*KXEL*Y3**2*U6+
00995      &       8.D0*KYEL*X3**2*V5-KYEL*X3**2*V1+KXEL*Y3*X3*V1+
00996      &       4.D0*KYEL*X2*Y3*U5+3.D0*KXEL*Y3**2*U2+
00997      &       3.D0*KYEL*X3**2*V2-8.D0*KYEL*X3**2*V6-
00998      &       6.D0*KYEL*X3**2*V3-6.D0*KXEL*Y3**2*U3+
00999      &       5.D0*KXEL*Y2*X3*V1+4.D0*KXEL*Y2*X3*V5-
01000      &       3.D0*KYEL*X3*Y3*U2+6.D0*KYEL*X3*Y3*U3+
01001      &       5.D0*KYEL*X3*X2*V3+8.D0*KYEL*X3*Y3*U6+
01002      &       5.D0*KXEL*Y3*Y2*U3+8.D0*KXEL*Y3*X3*V6+
01003      &       6.D0*KXEL*Y3*X3*V3+8.D0*KXEL*Y3**2*U5+
01004      &       4.D0*KYEL*X3**2*V4+5.D0*KYEL*X2*Y3*U1-
01005      &       5.D0*KYEL*X2*Y3*U3-5.D0*KXEL*Y2*X3*V3-
01006      &       4.D0*KYEL*X3*Y3*U4-4.D0*KXEL*Y3*Y2*U5-
01007      &       8.D0*KXEL*Y3*X3*V5+4.D0*KYEL*X3*X2*V4-
01008      &       4.D0*KYEL*X3*X2*V5-8.D0*KYEL*X3*Y3*U5-
01009      &       4.D0*KYEL*X2*Y3*U4-4.D0*KXEL*Y2*X3*V4
01010       A62(IELEM)= ANS1*XMUL/(90.D0*(-X2*Y3+Y2*X3))
01011 !
01012       ANS1 = 5.D0*KYEL*X3*Y2*U2-4.D0*KYEL*X2*Y2*U2-
01013      &       4.D0*KXEL*Y2*X2*V2-16.D0*KYEL*X2**2*V6-
01014      &       5.D0*KXEL*Y3*Y2*U1+5.D0*KXEL*Y3*X2*V1+
01015      &       4.D0*KYEL*X2**2*V2-4.D0*KXEL*Y3*X2*V4+
01016      &       4.D0*KXEL*Y2**2*U2+4.D0*KXEL*Y3*Y2*U4-
01017      &       16.D0*KXEL*Y2**2*U6+5.D0*KYEL*X3*Y2*U1-
01018      &       5.D0*KYEL*X3*X2*V1-5.D0*KXEL*Y3*Y2*U2+
01019      &       20.D0*KYEL*X2*Y2*U5+12.D0*KYEL*X2*Y2*U4-
01020      &       12.D0*KXEL*Y2**2*U4-15.D0*KYEL*X2**2*V3-
01021      &       20.D0*KXEL*Y2**2*U5-20.D0*KYEL*X2**2*V5+
01022      &       KXEL*Y2*X2*V1+15.D0*KXEL*Y2*X2*V3-KXEL*Y2**2*U1-
01023      &       15.D0*KXEL*Y2**2*U3-KYEL*X2**2*V1
01024       ANS2 = 20.D0*KXEL*Y2*X2*V5+12.D0*KXEL*Y2*X2*V4+
01025      &       18.D0*KYEL*X3*X2*V3-24.D0*KYEL*X3*Y2*U6+
01026      &       24.D0*KYEL*X3*X2*V6+18.D0*KXEL*Y3*Y2*U3+
01027      &       24.D0*KXEL*Y3*Y2*U6+KYEL*X2*Y2*U1+
01028      &       15.D0*KYEL*X2*Y2*U3-12.D0*KYEL*X2**2*V4-
01029      &       5.D0*KYEL*X3*X2*V2-4.D0*KYEL*X3*Y2*U4-
01030      &       24.D0*KXEL*Y3*X2*V5+24.D0*KXEL*Y3*Y2*U5+
01031      &       4.D0*KYEL*X3*X2*V4-24.D0*KXEL*Y3*X2*V6+
01032      &       24.D0*KYEL*X3*X2*V5-24.D0*KYEL*X3*Y2*U5+
01033      &       5.D0*KXEL*Y3*X2*V2+16.D0*KXEL*Y2*X2*V6-
01034      &       18.D0*KYEL*X3*Y2*U3+16.D0*KYEL*X2*Y2*U6-
01035      &       18.D0*KXEL*Y3*X2*V3
01036       A63(IELEM)= (ANS1+ANS2)*XMUL/(90.D0*(X2*Y3-Y2*X3))
01037 !
01038       ANS1 = -KXEL*Y3**2*U1-4.D0*KXEL*Y3*Y2*U1-KXEL*Y3*X2*V1-
01039      &       KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V4-8.D0*KXEL*Y3*Y2*U4-
01040      &       KYEL*X3*Y2*U1+KYEL*X3*Y3*U1-4.D0*KYEL*X3*X2*V1+
01041      &       KXEL*Y3*Y2*U2-4.D0*KXEL*Y3**2*U6+4.D0*KYEL*X3**2*V5-
01042      &       KYEL*X3**2*V1+KXEL*Y3*X3*V1+4.D0*KYEL*X2*Y3*U5-
01043      &       4.D0*KYEL*X2*Y2*U5+4.D0*KYEL*X2*Y2*U4+KXEL*Y3**2*U2+
01044      &       KYEL*X3**2*V2-4.D0*KYEL*X3**2*V6-4.D0*KXEL*Y2**2*U4+
01045      &       KYEL*X2**2*V3+4.D0*KXEL*Y2**2*U5+4.D0*KYEL*X2**2*V5+
01046      &       KXEL*Y2*X2*V1+5.D0*KXEL*Y2*X3*V1-KXEL*Y2*X2*V3-
01047      &       KXEL*Y2**2*U1+KXEL*Y2**2*U3-KYEL*X2**2*V1-
01048      &       4.D0*KXEL*Y2*X2*V5+4.D0*KXEL*Y2*X3*V5+
01049      &       4.D0*KXEL*Y2*X2*V4-KYEL*X3*Y3*U2+KYEL*X3*X2*V3+
01050      &       4.D0*KYEL*X3*Y2*U6-8.D0*KYEL*X3*X2*V6
01051       ANS2 = 4.D0*KYEL*X3*Y3*U6+KXEL*Y3*Y2*U3+4.D0*KXEL*Y3*X3*V6-
01052      &       8.D0*KXEL*Y3*Y2*U6+4.D0*KXEL*Y3**2*U5+
01053      &       4.D0*KXEL*Y2*X3*V6+KYEL*X2*Y2*U1+5.D0*KYEL*X2*Y3*U1-
01054      &       KYEL*X2*Y2*U3-KYEL*X2*Y3*U3+4.D0*KYEL*X2*Y3*U6-
01055      &       KXEL*Y2*X3*V3-4.D0*KYEL*X2**2*V4+KYEL*X3*X2*V2+
01056      &       4.D0*KYEL*X3*Y2*U4+8.D0*KXEL*Y3*X2*V5-
01057      &       12.D0*KXEL*Y3*Y2*U5-4.D0*KXEL*Y3*X3*V5-
01058      &       8.D0*KYEL*X3*X2*V4+4.D0*KXEL*Y3*X2*V6-
01059      &       12.D0*KYEL*X3*X2*V5-4.D0*KYEL*X3*Y3*U5+
01060      &       8.D0*KYEL*X3*Y2*U5-KYEL*X2*Y3*U2-KXEL*Y2*X3*V2+
01061      &       4.D0*KYEL*X2*Y3*U4+4.D0*KXEL*Y2*X3*V4
01062       A64(IELEM)= (ANS1+ANS2)*XMUL*2.D0/(45.D0*(X2*Y3-Y2*X3))
01063 !
01064       ANS1 = -2.D0*KXEL*Y3**2*U1+3.D0*KXEL*Y3*Y2*U1-KXEL*Y3*X2*V1+
01065      &       2.D0*KXEL*Y3*X3*V2+4.D0*KXEL*Y3*X2*V4+
01066      &       4.D0*KXEL*Y3**2*U4-4.D0*KXEL*Y3*Y2*U4-
01067      &       KYEL*X3*Y2*U1+2.D0*KYEL*X3*Y3*U1+3.D0*KYEL*X3*X2*V1-
01068      &       4.D0*KXEL*Y3*X3*V4+KXEL*Y3*Y2*U2+12.D0*KXEL*Y3**2*U6+
01069      &       12.D0*KYEL*X3**2*V5-2.D0*KYEL*X3**2*V1+
01070      &       2.D0*KXEL*Y3*X3*V1+8.D0*KYEL*X2*Y3*U5-
01071      &       4.D0*KYEL*X2*Y2*U5+4.D0*KYEL*X2*Y2*U4-
01072      &       2.D0*KXEL*Y3**2*U2-2.D0*KYEL*X3**2*V2+
01073      &       12.D0*KYEL*X3**2*V6+6.D0*KYEL*X3**2*V3+
01074      &       6.D0*KXEL*Y3**2*U3-4.D0*KXEL*Y2**2*U4+
01075      &       KYEL*X2**2*V3+4.D0*KXEL*Y2**2*U5+4.D0*KYEL*X2**2*V5+
01076      &       KXEL*Y2*X2*V1-2.D0*KXEL*Y2*X3*V1-KXEL*Y2*X2*V3-
01077      &       KXEL*Y2**2*U1+KXEL*Y2**2*U3-KYEL*X2**2*V1-
01078      &       4.D0*KXEL*Y2*X2*V5
01079       ANS2 = 8.D0*KXEL*Y2*X3*V5+4.D0*KXEL*Y2*X2*V4+
01080      &       2.D0*KYEL*X3*Y3*U2-6.D0*KYEL*X3*Y3*U3-
01081      &       6.D0*KYEL*X3*X2*V3+4.D0*KYEL*X3*Y2*U6-
01082      &       8.D0*KYEL*X3*X2*V6-12.D0*KYEL*X3*Y3*U6-
01083      &       6.D0*KXEL*Y3*Y2*U3-12.D0*KXEL*Y3*X3*V6-
01084      &       6.D0*KXEL*Y3*X3*V3-8.D0*KXEL*Y3*Y2*U6+
01085      &       12.D0*KXEL*Y3**2*U5+4.D0*KYEL*X3**2*V4+
01086      &       4.D0*KXEL*Y2*X3*V6+KYEL*X2*Y2*U1-2.D0*KYEL*X2*Y3*U1-
01087      &       KYEL*X2*Y2*U3+6.D0*KYEL*X2*Y3*U3+4.D0*KYEL*X2*Y3*U6+
01088      &       6.D0*KXEL*Y2*X3*V3-4.D0*KYEL*X2**2*V4+KYEL*X3*X2*V2-
01089      &       4.D0*KYEL*X3*Y3*U4+4.D0*KYEL*X3*Y2*U4+
01090      &       8.D0*KXEL*Y3*X2*V5-16.D0*KXEL*Y3*Y2*U5-
01091      &       12.D0*KXEL*Y3*X3*V5-4.D0*KYEL*X3*X2*V4+
01092      &       4.D0*KXEL*Y3*X2*V6-16.D0*KYEL*X3*X2*V5-
01093      &       12.D0*KYEL*X3*Y3*U5+8.D0*KYEL*X3*Y2*U5-KYEL*X2*Y3*U2-
01094      &       KXEL*Y2*X3*V2
01095       A65(IELEM)= (ANS1+ANS2)*XMUL*2.D0/(45.D0*(-X2*Y3+Y2*X3))
01096 !
01097 ! USES HERE THE 'MAGIC SQUARE' PROPERTIES
01098 ! (SUM OF EACH LINE = SUM OF EACH COLUMN = 0)
01099 !
01100       A16(IELEM) = - A11(IELEM) - A12(IELEM) - A13(IELEM)
01101      &             - A14(IELEM) - A15(IELEM)
01102       A25(IELEM) = - A21(IELEM) - A22(IELEM) - A23(IELEM)
01103      &             - A24(IELEM) - A26(IELEM)
01104       A36(IELEM) = - A31(IELEM) - A32(IELEM) - A33(IELEM)
01105      &             - A34(IELEM) - A35(IELEM)
01106       A43(IELEM) = - A41(IELEM) - A42(IELEM) - A44(IELEM)
01107      &             - A45(IELEM) - A46(IELEM)
01108       A55(IELEM) = - A51(IELEM) - A52(IELEM) - A53(IELEM)
01109      &             - A54(IELEM) - A56(IELEM)
01110       A66(IELEM) = - A61(IELEM) - A62(IELEM) - A63(IELEM)
01111      &             - A64(IELEM) - A65(IELEM)
01112 !
01113       ENDDO ! IELEM
01114 !
01115 !     OTHER TYPES OF FUNCTIONS F AND G
01116 !
01117 !-----------------------------------------------------------------------
01118 !
01119       ELSE
01120 !
01121         IF (LNG.EQ.1) WRITE(LU,100) IELMF,SF%NAME
01122         IF (LNG.EQ.1) WRITE(LU,110) IELMG,SG%NAME
01123         IF (LNG.EQ.1) WRITE(LU,200) IELMU,SU%NAME
01124         IF (LNG.EQ.1) WRITE(LU,300)
01125         IF (LNG.EQ.1) WRITE(LU,101) IELMF,SF%NAME
01126         IF (LNG.EQ.1) WRITE(LU,111) IELMG,SG%NAME
01127         IF (LNG.EQ.1) WRITE(LU,201) IELMU,SU%NAME
01128         IF (LNG.EQ.1) WRITE(LU,301)
01129 100     FORMAT(1X,'MT03CC (BIEF) :',/,
01130      &         1X,'DISCRETISATION DE F : ',1I6,
01131      &         1X,'NOM REEL : ',A6)
01132 110     FORMAT(1X,'DISCRETISATION DE G : ',1I6,
01133      &         1X,'NOM REEL : ',A6)
01134 200     FORMAT(1X,'DISCRETISATION DE U : ',1I6,
01135      &         1X,'NOM REEL : ',A6)
01136 300     FORMAT(1X,'CAS NON PREVU')
01137 101     FORMAT(1X,'MT03CC (BIEF) :',/,
01138      &         1X,'DISCRETIZATION OF F:',1I6,
01139      &         1X,'REAL NAME: ',A6)
01140 111     FORMAT(1X,'DISCRETIZATION OF G:',1I6,
01141      &         1X,'REAL NAME: ',A6)
01142 201     FORMAT(1X,'DISCRETIZATION OF U:',1I6,
01143      &         1X,'REAL NAME: ',A6)
01144 301     FORMAT(1X,'CASE NOT IMPLEMENTED')
01145         CALL PLANTE(1)
01146         STOP
01147 !
01148       ENDIF
01149 !
01150 !-----------------------------------------------------------------------
01151 !
01152       RETURN
01153       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0