mt04cc.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\mt04cc.f
00002 !
00065                      SUBROUTINE MT04CC
00066 !                    *****************
00067 !
00068      &( A11 , A12 , A13 , A14 , A15 , A16 ,
00069      &        A22 , A23 , A24 , A25 , A26 ,
00070      &              A33 , A34 , A35 , A36 ,
00071      &                    A44 , A45 , A46 ,
00072      &                          A55 , A56 ,
00073      &                                A66 ,
00074      &  XMUL,SU,SV,U,V,XEL,YEL,IKLE1,IKLE2,IKLE3,
00075      &  IKLE4,IKLE5,IKLE6,NELEM,NELMAX)
00076 !
00077 !***********************************************************************
00078 ! BIEF   V6P1                                   21/08/2010
00079 !***********************************************************************
00080 !
00081 !
00082 !
00083 !
00084 !
00085 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00086 !| A11            |<--| ELEMENTS OF MATRIX
00087 !| ...            |<--| ELEMENTS OF MATRIX
00088 !| A66            |<--| ELEMENTS OF MATRIX
00089 !| IKLE1          |-->| FIRST POINTS OF TRIANGLES
00090 !| IKLE2          |-->| SECOND POINTS OF TRIANGLES
00091 !| IKLE3          |-->| THIRD POINTS OF TRIANGLES
00092 !| IKLE4          |-->| FOURTH POINTS OF TRIANGLES (QUADRATIC)
00093 !| IKLE5          |-->| FIFTH POINTS OF TRIANGLES (QUADRATIC)
00094 !| IKLE6          |-->| SIXTH POINTS OF TRIANGLES (QUADRATIC)
00095 !| NELEM          |-->| NUMBER OF ELEMENTS
00096 !| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
00097 !| SU             |-->| BIEF_OBJ STRUCTURE OF U
00098 !| SURFAC         |-->| AREA OF TRIANGLES
00099 !| SV             |-->| BIEF_OBJ STRUCTURE OF V
00100 !| U              |-->| FUNCTION U USED IN THE FORMULA
00101 !| V              |-->| FUNCTION V USED IN THE FORMULA
00102 !| XEL            |-->| ABSCISSAE OF POINTS IN THE MESH, PER ELEMENT
00103 !| YEL            |-->| ORDINATES OF POINTS IN THE MESH, PER ELEMENT
00104 !| XMUL           |-->| MULTIPLICATION FACTOR
00105 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00106 !
00107       USE BIEF!, EX_MT04CC => MT04CC
00108 !
00109       IMPLICIT NONE
00110       INTEGER LNG,LU
00111       COMMON/INFO/LNG,LU
00112 !
00113 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00114 !
00115       INTEGER, INTENT(IN) :: NELEM,NELMAX
00116       INTEGER, INTENT(IN) :: IKLE1(NELMAX),IKLE2(NELMAX)
00117       INTEGER, INTENT(IN) :: IKLE3(NELMAX),IKLE4(NELMAX)
00118       INTEGER, INTENT(IN) :: IKLE5(NELMAX),IKLE6(NELMAX)
00119 !
00120       DOUBLE PRECISION, INTENT(INOUT) :: A11(*),A12(*),A13(*)
00121       DOUBLE PRECISION, INTENT(INOUT) :: A14(*),A15(*),A16(*)
00122       DOUBLE PRECISION, INTENT(INOUT) :: A22(*),A23(*),A24(*)
00123       DOUBLE PRECISION, INTENT(INOUT) :: A25(*),A26(*),A33(*)
00124       DOUBLE PRECISION, INTENT(INOUT) :: A34(*),A35(*),A36(*)
00125       DOUBLE PRECISION, INTENT(INOUT) :: A44(*),A45(*),A46(*)
00126       DOUBLE PRECISION, INTENT(INOUT) :: A55(*),A56(*),A66(*)
00127 !
00128       DOUBLE PRECISION, INTENT(IN) :: XMUL
00129       DOUBLE PRECISION, INTENT(IN) :: U(*),V(*)
00130 !
00131 !     STRUCTURES OF      U, V
00132       TYPE(BIEF_OBJ), INTENT(IN) :: SU,SV
00133 !
00134       DOUBLE PRECISION, INTENT(IN) :: XEL(NELMAX,3),YEL(NELMAX,3)
00135 !
00136 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00137 !
00138 !     DECLARATIONS SPECIFIC TO THIS SUBROUTINE
00139 !
00140       INTEGER IELMU,IELMV,IELEM
00141 !
00142       DOUBLE PRECISION X2,X3,Y2,Y3,ANS1,ANS2,ANS3,ANS4
00143       DOUBLE PRECISION U1,U2,U3,U4,U5,U6
00144       DOUBLE PRECISION V1,V2,V3,V4,V5,V6
00145       DOUBLE PRECISION UNSU2, AUX360, AUX180, AUX45, AUX720, AUX1260
00146       DOUBLE PRECISION AUX2520, AUX630
00147 !
00148 !=======================================================================
00149 !
00150 !     EXTRACTS THE TYPE OF ELEMENT FOR VELOCITY
00151 !
00152       IELMU = SU%ELM
00153       IELMV = SV%ELM
00154 !
00155 !-----------------------------------------------------------------------
00156 !
00157       IF(IELMU.EQ.11.AND.IELMV.EQ.11) THEN
00158 !
00159 !-----------------------------------------------------------------------
00160 !
00161 !  P1 DISCRETISATION OF THE VELOCITY:
00162 !
00163       DO IELEM = 1 , NELEM
00164 !
00165 !   INITIALISES THE GEOMETRICAL VARIABLES
00166 !
00167         X2  =  XEL(IELEM,2)
00168         X3  =  XEL(IELEM,3)
00169 !
00170         Y2  =  YEL(IELEM,2)
00171         Y3  =  YEL(IELEM,3)
00172 !
00173         U1 = U(IKLE1(IELEM))
00174         U2 = U(IKLE2(IELEM))
00175         U3 = U(IKLE3(IELEM))
00176         V1 = V(IKLE1(IELEM))
00177         V2 = V(IKLE2(IELEM))
00178         V3 = V(IKLE3(IELEM))
00179 !
00180         UNSU2 = XMUL/(X2*Y3-Y2*X3)
00181 !
00182         AUX720 = UNSU2/360.D0
00183         AUX360 = UNSU2/180.D0
00184         AUX180 = UNSU2/90.D0
00185         AUX45 = 2.D0*UNSU2/45.D0
00186 !
00187 !
00188 !   INITIALISES THE INTERMEDIATE VARIABLES
00189 !
00190 !
00191 !  COMPUTES 15 OF THE 36 TERMS (SELECTED AMONG THE LEAST COMPLEX)
00192 !
00193 !
00194       ANS1 = 7.D0*U3**2*Y2**2-14.D0*U3*Y2*V3*X2-78.D0*U1*Y2*V1*X2+
00195      &       7.D0*V3**2*X2**2+39.D0*V1**2*X2**2-7.D0*U2*Y2*V3*X2-
00196      &       15.D0*U2*Y2*V1*X2+39.D0*U1**2*Y2**2-15.D0*V2*X2*U1*Y2-
00197      &       7.D0*V2*X2*U3*Y2-15.D0*U3*Y2*V1*X2-15.D0*U1*Y2*V3*X2+
00198      &       7.D0*V2**2*X2**2+7.D0*U2**2*Y2**2-14.D0*U2*Y2*V2*X2+
00199      &       7.D0*V3**2*X3**2-30.D0*U2*Y3*U1*Y2+14*U3*Y3*V3*X2+
00200      &       15.D0*U3*Y2*V1*X3-30.D0*V1*X3*V3*X2-14.D0*U3*Y3*V3*X3-
00201      &       78.D0*V1**2*X3*X2-15.D0*U3*Y3*V1*X3+15.D0*U1*Y2*V3*X3+
00202      &       15.D0*U3*Y3*V1*X2+7.D0*U2*Y3**2.D0*U3+14.D0*U3*Y2*V3*X3+
00203      &       78.D0*U1*Y2*V1*X3-14.D0*U3**2*Y3*Y2+15.D0*U1*Y3**2*U3+
00204      &       15.D0*U2*Y3**2.D0*U1-14.D0*V3**2*X3*X2+15.D0*V1*X3**2*V3-
00205      &       15.D0*U1*Y3*V3*X3-78.D0*U1**2*Y3*Y2+78.D0*U1*Y3*V1*X2+
00206      &       15.D0*U1*Y3*V3*X2-30.D0*U1*Y3*U3*Y2+15.D0*V2*X3**2*V1+
00207      &       7.D0*V2*X3**2*V3-14.D0*U2**2*Y3*Y2+7.D0*U3**2*Y3**2+
00208      &       39.D0*U1**2*Y3**2+14.D0*U2*Y2*V2*X3-78.D0*U1*Y3*V1*X3-
00209      &       14.D0*V2**2*X3*X2+15.D0*V1*X2**2*V3-14.D0*U2*Y3*V2*X3-
00210      &       15.D0*V2*X3*U1*Y3+15.D0*V2*X3*U1*Y2+7.D0*V2**2*X3**2+
00211      &       7.D0*U2**2*Y3**2+14.D0*U2*Y3*V2*X2+7.D0*V2*X2*U3*Y3
00212       A11(IELEM) = (7.D0*V2*X2**2*V3-14.D0*V2*X3*V3*X2+
00213      &              15.D0*V2*X2*U1*Y3+15.D0*U2*Y2*V1*X3+
00214      &              7.D0*V2*X3*U3*Y2-30.D0*V2*X3*V1*X2-
00215      &              7.D0*V2*X3*U3*Y3+39.D0*V1**2*X3**2+
00216      &              7.D0*U2*Y2**2*U3+15.D0*V2*X2**2*V1-
00217      &              14.D0*U2*Y3*U3*Y2-15.D0*U2*Y3*V1*X3-
00218      &              7.D0*U2*Y3*V3*X3+7.D0*U2*Y2*V3*X3+
00219      &              15.D0*U2*Y3*V1*X2+7.D0*U2*Y3*V3*X2+
00220      &              15.D0*U2*Y2**2*U1+15.D0*U1*Y2**2.D0*U3+
00221      &              ANS1) * AUX360
00222 !
00223       A12(IELEM)= (-18.D0*U1**2*Y3**2-18.D0*V2**2*X3**2-
00224      &             18.D0*V1**2*X3**2+18.D0*U1**2*Y3*Y2+
00225      &             18.D0*V1**2*X3*X2-10.D0*V1*X3**2*V3+
00226      &             18.D0*V2**2*X3*X2+18.D0*U2**2*Y3*Y2+2.D0*V3**2*X3*X2+
00227      &             2.D0*U3**2*Y3*Y2-5.D0*V1*X2*U3*Y3-V2*X3*U1*Y2-
00228      &             5.D0*V2*X3*U3*Y2+2.D0*V2*X3*V1*X2+10.D0*V2*X3*V3*X2-
00229      &             V2*X2*U1*Y3-5.D0*V2*X2*U3*Y3-U2*Y2*V1*X3-
00230      &             5.D0*U2*Y2*V3*X3-5.D0*V1*X3*U3*Y2+10.D0*V1*X3*V3*X2+
00231      &             10.D0*U1*Y3*U3*Y2-18.D0*U1*Y3*V1*X2-5.D0*U1*Y3*V3*X2-
00232      &             18.D0*U1*Y2*V1*X3-5.D0*U1*Y2*V3*X3+2.D0*U2*Y3*U1*Y2+
00233      &             10.D0*U2*Y3*U3*Y2-U2*Y3*V1*X2-5.D0*U2*Y3*V3*X2-
00234      &             18.D0*U2*Y3*V2*X2-18.D0*U2*Y2*V2*X3-2.D0*U3*Y3*V3*X2-
00235      &             2.D0*U3*Y2*V3*X3-18.D0*U2**2*Y3**2+2.D0*V2*X3*U1*Y3+
00236      &             10.D0*V2*X3*U3*Y3+10.D0*V1*X3*U3*Y3+
00237      &             36.D0*U1*Y3*V1*X3+10.D0*U1*Y3*V3*X3-
00238      &             2.D0*U3**2*Y3**2-2.D0*V3**2*X3**2-10.D0*V2*X3**2*V3-
00239      &             10.D0*U1*Y3**2*U3-2.D0*V2*X3**2*V1-10.D0*U2*Y3**2*U3-
00240      &             2.D0*U2*Y3**2*U1+10.D0*U2*Y3*V3*X3+36.D0*U2*Y3*V2*X3+
00241      &             4.D0*U3*Y3*V3*X3+2.D0*U2*Y3*V1*X3)*(-AUX720)
00242 !
00243       A13(IELEM) =(2.D0*U2**2*Y2**2+18.D0*U1**2*Y2**2+
00244      &            18.D0*V1**2*X2**2+2.D0*V2**2*X2**2+18.D0*V3**2*X2**2+
00245      &            10.D0*U2*Y2**2*U3-18.D0*U1**2*Y3*Y2-18.D0*V1**2*X3*X2+
00246      &            10.D0*U2*Y2**2*U1-2.D0*V2**2*X3*X2-2.D0*U2**2*Y3*Y2-
00247      &            18.D0*V3**2*X3*X2-18.D0*U3**2*Y3*Y2+V1*X2*U3*Y3-
00248      &            2.D0*V1*X2*U3*Y2+5.D0*V2*X3*U1*Y2+5.D0*V2*X3*U3*Y2-
00249      &            10.D0*V2*X3*V1*X2-10.D0*V2*X3*V3*X2+5.D0*V2*X2*U1*Y3-
00250      &            10.D0*V2*X2*U1*Y2+5.D0*V2*X2*U3*Y3-10.D0*V2*X2*U3*Y2+
00251      &            5.D0*U2*Y2*V1*X3-10.D0*U2*Y2*V1*X2+5.D0*U2*Y2*V3*X3-
00252      &            10.D0*U2*Y2*V3*X2+V1*X3*U3*Y2-2.D0*V1*X3*V3*X2-
00253      &            2.D0*U1*Y3*U3*Y2+18.D0*U1*Y3*V1*X2+U1*Y3*V3*X2+
00254      &            18.D0*U1*Y2*V1*X3-36.D0*U1*Y2*V1*X2+U1*Y2*V3*X3-
00255      &            2.D0*U1*Y2*V3*X2-10.D0*U2*Y3*U1*Y2-10.D0*U2*Y3*U3*Y2+
00256      &            5.D0*U2*Y3*V1*X2+5.D0*U2*Y3*V3*X2+2.D0*U2*Y3*V2*X2+
00257      &            2.D0*U2*Y2*V2*X3-4.D0*U2*Y2*V2*X2+18.D0*U3*Y3*V3*X2+
00258      &            18.D0*U3*Y2*V3*X3-36.D0*U3*Y2*V3*X2+18.D0*U3**2*Y2**2+
00259      &            10.D0*V2*X2**2*V3+10.D0*V2*X2**2*V1+2.D0*V1*X2**2*V3+
00260      &            2.D0*U1*Y2**2*U3) * AUX720
00261 !
00262       ANS1 = 2.D0*U3**2*Y2**2-4.D0*U3*Y2*V3*X2+12.D0*U1*Y2*V1*X2+
00263      &       2.D0*V3**2*X2**2-6.D0*V1**2*X2**2-4.D0*U2*Y2*V3*X2+
00264      &       4.D0*U2*Y2*V1*X2-6.D0*U1**2*Y2**2+4.D0*V2*X2*U1*Y2-
00265      &       4.D0*V2*X2*U3*Y2+2.D0*U3*Y2*V1*X2+2.D0*U1*Y2*V3*X2+
00266      &       6.D0*V2**2*X2**2+6.D0*U2**2*Y2**2-12.D0*U2*Y2*V2*X2+
00267      &       4.D0*V3**2*X3**2-4.D0*U2*Y3*U1*Y2+6.D0*U3*Y3*V3*X2+
00268      &       4.D0*U3*Y2*V1*X3-8.D0*V1*X3*V3*X2-8.D0*U3*Y3*V3*X3-
00269      &       18.D0*V1**2*X3*X2-10.D0*U3*Y3*V1*X3+4.D0*U1*Y2*V3*X3+
00270      &       4.D0*U3*Y3*V1*X2+6.D0*U2*Y3**2*U3+6.D0*U3*Y2*V3*X3+
00271      &       18.D0*U1*Y2*V1*X3-6.D0*U3**2*Y3*Y2+10.D0*U1*Y3**2*U3+
00272      &       8.D0*U2*Y3**2*U1-6.D0*V3**2*X3*X2+10.D0*V1*X3**2*V3-
00273      &       10.D0*U1*Y3*V3*X3-18.D0*U1**2*Y3*Y2+18.D0*U1*Y3*V1*X2+
00274      &       4.D0*U1*Y3*V3*X2-8.D0*U1*Y3*U3*Y2+8.D0*V2*X3**2*V1+
00275      &       6.D0*V2*X3**2*V3-14.D0*U2**2*Y3*Y2+4.D0*U3**2*Y3**2+
00276      &       24.D0*U1**2*Y3**2+14.D0*U2*Y2*V2*X3-48.D0*U1*Y3*V1*X3-
00277      &       14.D0*V2**2*X3*X2-2*V1*X2**2*V3-16.D0*U2*Y3*V2*X3-
00278      &       8.D0*V2*X3*U1*Y3+2.D0*V2*X3*U1*Y2+8.D0*V2**2*X3**2+
00279      &       8.D0*U2**2*Y3**2+14.D0*U2*Y3*V2*X2+5.D0*V2*X2*U3*Y3+
00280      &       4.D0*V2*X2**2*V3-10.D0*V2*X3*V3*X2+2.D0*V2*X2*U1*Y3
00281       A14(IELEM) = (2.D0*U2*Y2*V1*X3+5.D0*V2*X3*U3*Y2-4.D0*V2*X3*V1*X2-
00282      &             6.D0*V2*X3*U3*Y3+24.D0*V1**2*X3**2+4.D0*U2*Y2**2*U3-
00283      &             4.D0*V2*X2**2*V1-10.D0*U2*Y3*U3*Y2-8.D0*U2*Y3*V1*X3-
00284      &             6.D0*U2*Y3*V3*X3+5.D0*U2*Y2*V3*X3+2.D0*U2*Y3*V1*X2+
00285      &             5.D0*U2*Y3*V3*X2-4.D0*U2*Y2**2*U1-2.D0*U1*Y2**2*U3+
00286      &             ANS1)*(- AUX180)
00287 !
00288       ANS1 = 2.D0*U3**2*Y2**2-4.D0*U3*Y2*V3*X2+12.D0*U1*Y2*V1*X2+
00289      &       2.D0*V3**2*X2**2-6.D0*V1**2*X2**2-4.D0*U2*Y2*V3*X2+
00290      &       4.D0*U2*Y2*V1*X2-6.D0*U1**2*Y2**2+4.D0*V2*X2*U1*Y2-
00291      &       4.D0*V2*X2*U3*Y2+2.D0*U3*Y2*V1*X2+2.D0*U1*Y2*V3*X2+
00292      &       6.D0*V2**2*X2**2+6.D0*U2**2*Y2**2-12.D0*U2*Y2*V2*X2+
00293      &       6.D0*V3**2*X3**2+6.D0*U2*Y3*U1*Y2+8.D0*U3*Y3*V3*X2-
00294      &       3.D0*U3*Y2*V1*X3+6.D0*V1*X3*V3*X2-12.D0*U3*Y3*V3*X3+
00295      &       12.D0*V1**2*X3*X2+4.D0*U3*Y3*V1*X3-3.D0*U1*Y2*V3*X3-
00296      &       3.D0*U3*Y3*V1*X2+4.D0*U2*Y3**2*U3+8.D0*U3*Y2*V3*X3-
00297      &       12.D0*U1*Y2*V1*X3-8.D0*U3**2*Y3*Y2-4.D0*U1*Y3**2*U3-
00298      &       2.D0*U2*Y3**2.D0*U1-8.D0*V3**2*X3*X2-4.D0*V1*X3**2*V3+
00299      &       4.D0*U1*Y3*V3*X3+12.D0*U1**2*Y3*Y2-12.D0*U1*Y3*V1*X2-
00300      &       3.D0*U1*Y3*V3*X2+6.D0*U1*Y3*U3*Y2-2.D0*V2*X3**2*V1+
00301      &       4.D0*V2*X3**2*V3-8.D0*U2**2*Y3*Y2+6.D0*U3**2*Y3**2-
00302      &       6.D0*U1**2*Y3**2+8.D0*U2*Y2*V2*X3+12.D0*U1*Y3*V1*X3-
00303      &       8.D0*V2**2*X3*X2-2.D0*V1*X2**2*V3-4.D0*U2*Y3*V2*X3+
00304      &       2.D0*V2*X3*U1*Y3-3.D0*V2*X3*U1*Y2+2.D0*V2**2*X3**2
00305       A15(IELEM) = (2.D0*U2**2*Y3**2+8.D0*U2*Y3*V2*X2+4.D0*V2*X2*U3*Y3+
00306      &              4.D0*V2*X2**2*V3-8.D0*V2*X3*V3*X2-3.D0*V2*X2*U1*Y3-
00307      &              3.D0*U2*Y2*V1*X3+4.D0*V2*X3*U3*Y2+6.D0*V2*X3*V1*X2-
00308      &              4.D0*V2*X3*U3*Y3-6.D0*V1**2*X3**2+4.D0*U2*Y2**2*U3-
00309      &              4.D0*V2*X2**2*V1-8.D0*U2*Y3*U3*Y2+2.D0*U2*Y3*V1*X3-
00310      &              4.D0*U2*Y3*V3*X3+4.D0*U2*Y2*V3*X3-3.D0*U2*Y3*V1*X2+
00311      &              4.D0*U2*Y3*V3*X2-4.D0*U2*Y2**2*U1-2.D0*U1*Y2**2*U3+
00312      &              ANS1)* AUX180
00313 !
00314       ANS1 = 8.D0*U3**2*Y2**2-16.D0*U3*Y2*V3*X2-48.D0*U1*Y2*V1*X2+
00315      &       8.D0*V3**2*X2**2+24.D0*V1**2*X2**2-6.D0*U2*Y2*V3*X2-
00316      &       10.D0*U2*Y2*V1*X2+24.D0*U1**2*Y2**2-10.D0*V2*X2*U1*Y2-
00317      &       6.D0*V2*X2*U3*Y2-8.D0*U3*Y2*V1*X2-8.D0*U1*Y2*V3*X2+
00318      &       4.D0*V2**2*X2**2+4.D0*U2**2*Y2**2-8.D0*U2*Y2*V2*X2+
00319      &       6.D0*V3**2*X3**2-8.D0*U2*Y3*U1*Y2+14.D0*U3*Y3*V3*X2+
00320      &       2.D0*U3*Y2*V1*X3-4.D0*V1*X3*V3*X2-12.D0*U3*Y3*V3*X3-
00321      &       18.D0*V1**2*X3*X2+4.D0*U3*Y3*V1*X3+2.D0*U1*Y2*V3*X3+
00322      &       2.D0*U3*Y3*V1*X2+4.D0*U2*Y3**2*U3+14.D0*U3*Y2*V3*X3+
00323      &       18.D0*U1*Y2*V1*X3-14.D0*U3**2*Y3*Y2-4.D0*U1*Y3**2*U3-
00324      &       2.D0*U2*Y3**2*U1-14.D0*V3**2*X3*X2-4.D0*V1*X3**2*V3+
00325      &       4.D0*U1*Y3*V3*X3-18.D0*U1**2*Y3*Y2+18.D0*U1*Y3*V1*X2+
00326      &       2.D0*U1*Y3*V3*X2-4.D0*U1*Y3*U3*Y2-2.D0*V2*X3**2*V1+
00327      &       4.D0*V2*X3**2*V3-6.D0*U2**2*Y3*Y2+6.D0*U3**2*Y3**2-
00328      &       6.D0*U1**2*Y3**2+6.D0*U2*Y2*V2*X3+12.D0*U1*Y3*V1*X3-
00329      &       6.D0*V2**2*X3*X2+8.D0*V1*X2**2*V3-4.D0*U2*Y3*V2*X3+
00330      &       2.D0*V2*X3*U1*Y3+4.D0*V2*X3*U1*Y2+2.D0*V2**2*X3**2+
00331      &       2.D0*U2**2*Y3**2+6.D0*U2*Y3*V2*X2+5.D0*V2*X2*U3*Y3+
00332      &       6.D0*V2*X2**2*V3-10.D0*V2*X3*V3*X2+4.D0*V2*X2*U1*Y3+
00333      &       4.D0*U2*Y2*V1*X3+5.D0*V2*X3*U3*Y2-8.D0*V2*X3*V1*X2
00334       A16(IELEM) = (-4.D0*V2*X3*U3*Y3-6.D0*V1**2*X3**2+6.D0*U2*Y2**2*U3+
00335      &             10.D0*V2*X2**2*V1-10.D0*U2*Y3*U3*Y2+2.D0*U2*Y3*V1*X3-
00336      &             4.D0*U2*Y3*V3*X3+5.D0*U2*Y2*V3*X3+4.D0*U2*Y3*V1*X2+
00337      &             5.D0*U2*Y3*V3*X2+10.D0*U2*Y2**2*U1+8.D0*U1*Y2**2*U3+
00338      &             ANS1)*(-AUX180)
00339 !
00340       A22(IELEM)= ((7.D0*(V1*V3+V3**2+V1**2)+39.D0*V2**2+
00341      &             15.D0*(V2*V3+V2*V1))*X3**2+(-7.D0*U3*V1-
00342      &             14.D0*(U1*V1+U3*V3)-15.D0*(U2*V1+V2*U3+V2*U1+U2*V3)-
00343      &             78.D0*U2*V2-7.D0*U1*V3)*Y3*X3+(7.D0*(U1**2+U3**2+
00344      &             U1*U3)+39.D0*U2**2+15.D0*(U2*U3+U2*U1))*Y3**2)
00345      &             *AUX360
00346 !
00347       A23(IELEM)=(((10.D0*(V1*V3+V2*V1)+18.D0*(V3**2+V2**2)+2.D0*(V2*V3+
00348      &            V1**2))*X3+(-U2*V3-V2*U3-5.D0*(U1*V3+U3*V1+V2*U1+
00349      &            U2*V1)-18.D0*(U2*V2+U3*V3)-2.D0*U1*V1)*Y3)*X2+(-U2*V3-
00350      &            V2*U3-5.D0*(U1*V3+U3*V1+U2*V1+V2*U1)-18.D0*(U2*V2+
00351      &            U3*V3)-2.D0*U1*V1)*Y2*X3+(18.D0*(U2**2+U3**2)+
00352      &            2.D0*(U2*U3+U1**2)+10.D0*(U2*U1+U1*U3))*Y3*Y2)*AUX720
00353 !
00354       A24(IELEM) = (-8.D0*U1**2*Y3**2-24.D0*V2**2*X3**2-
00355      &             8.D0*V1**2*X3**2+2.D0*U1**2*Y3*Y2+2.D0*V1**2*X3*X2-
00356      &             6.D0*V1*X3**2*V3+30.D0*V2**2*X3*X2+
00357      &             30.D0*U2**2*Y3*Y2+2.D0*V3**2*X3*X2+2.D0*U3**2*Y3*Y2-
00358      &             V1*X2*U3*Y3-6.D0*V2*X3*U1*Y2-6.D0*V2*X3*U3*Y2+
00359      &             12.D0*V2*X3*V1*X2+12.D0*V2*X3*V3*X2-6.D0*V2*X2*U1*Y3-
00360      &             6.D0*V2*X2*U3*Y3-6.D0*U2*Y2*V1*X3-6.D0*U2*Y2*V3*X3-
00361      &             V1*X3*U3*Y2+2.D0*V1*X3*V3*X2+2.D0*U1*Y3*U3*Y2-
00362      &             2.D0*U1*Y3*V1*X2-U1*Y3*V3*X2-2.D0*U1*Y2*V1*X3-
00363      &             U1*Y2*V3*X3+12.D0*U2*Y3*U1*Y2+12.D0*U2*Y3*U3*Y2-
00364      &             6.D0*U2*Y3*V1*X2-6.D0*U2*Y3*V3*X2-30.D0*U2*Y3*V2*X2-
00365      &             30.D0*U2*Y2*V2*X3-2.D0*U3*Y3*V3*X2-2.D0*U3*Y2*V3*X3-
00366      &             24.D0*U2**2*Y3**2+8.D0*V2*X3*U1*Y3+10.D0*V2*X3*U3*Y3+
00367      &             6.D0*V1*X3*U3*Y3+16.D0*U1*Y3*V1*X3+6.D0*U1*Y3*V3*X3-
00368      &             4.D0*U3**2*Y3**2-4.D0*V3**2*X3**2-10.D0*V2*X3**2*V3-
00369      &             6.D0*U1*Y3**2*U3-8.D0*V2*X3**2*V1-10.D0*U2*Y3**2*U3-
00370      &             8.D0*U2*Y3**2*U1+10.D0*U2*Y3*V3*X3+48.D0*U2*Y3*V2*X3+
00371      &             8.D0*U3*Y3*V3*X3+8*U2*Y3*V1*X3)*AUX180
00372 !
00373       A25(IELEM) = (6.D0*V3**2*X3**2+12.D0*U2*Y3*U1*Y2-2.D0*U3*Y3*V3*X2-
00374      &             U3*Y2*V1*X3+2.D0*V1*X3*V3*X2-12.D0*U3*Y3*V3*X3+
00375      &             2.D0*V1**2*X3*X2-4.D0*U3*Y3*V1*X3-U1*Y2*V3*X3-
00376      &             U3*Y3*V1*X2-4.D0*U2*Y3**2*U3-2.D0*U3*Y2*V3*X3-
00377      &             2.D0*U1*Y2*V1*X3+2.D0*U3**2*Y3*Y2+4.D0*U1*Y3**2*U3-
00378      &             2.D0*U2*Y3**2*U1+2.D0*V3**2*X3*X2+4.D0*V1*X3**2*V3-
00379      &             4.D0*U1*Y3*V3*X3+2.D0*U1**2*Y3*Y2-2.D0*U1*Y3*V1*X2-
00380      &             U1*Y3*V3*X2+2.D0*U1*Y3*U3*Y2-2.D0*V2*X3**2.D0*V1-
00381      &             4.D0*V2*X3**2*V3+30.D0*U2**2*Y3*Y2+6.D0*U3**2*Y3**2+
00382      &             2.D0*U1**2*Y3**2-30.D0*U2*Y2*V2*X3-4.D0*U1*Y3*V1*X3+
00383      &             30.D0*V2**2*X3*X2+12.D0*U2*Y3*V2*X3+2.D0*V2*X3*U1*Y3-
00384      &             6.D0*V2*X3*U1*Y2-6.D0*V2**2*X3**2-6.D0*U2**2*Y3**2-
00385      &             30.D0*U2*Y3*V2*X2-6.D0*V2*X2*U3*Y3+12.D0*V2*X3*V3*X2-
00386      &             6.D0*V2*X2*U1*Y3-6.D0*U2*Y2*V1*X3-6.D0*V2*X3*U3*Y2+
00387      &             12.D0*V2*X3*V1*X2+4.D0*V2*X3*U3*Y3+2.D0*V1**2*X3**2+
00388      &             12.D0*U2*Y3*U3*Y2+2.D0*U2*Y3*V1*X3+4.D0*U2*Y3*V3*X3-
00389      &             6.D0*U2*Y2*V3*X3-6.D0*U2*Y3*V1*X2-6.D0*U2*Y3*V3*X2)*
00390      &             (-AUX180)
00391 !
00392       A26(IELEM) = (2.D0*U1**2*Y3**2-6.D0*V2**2*X3**2+2.D0*V1**2*X3**2+
00393      &             4.D0*U1**2*Y3*Y2+4.D0*V1**2*X3*X2+4.D0*V1*X3**2*V3-
00394      &             4.D0*V3**2*X3*X2-4.D0*U3**2*Y3*Y2+V2*X3*U1*Y2-
00395      &             V2*X3*U3*Y2-2.D0*V2*X3*V1*X2+2.D0*V2*X3*V3*X2+
00396      &             V2*X2*U1*Y3-V2*X2*U3*Y3+U2*Y2*V1*X3-U2*Y2*V3*X3-
00397      &             4.D0*U1*Y3*V1*X2-4.D0*U1*Y2*V1*X3-2.D0*U2*Y3*U1*Y2+
00398      &             2.D0*U2*Y3*U3*Y2+U2*Y3*V1*X2-U2*Y3*V3*X2+
00399      &             4.D0*U3*Y3*V3*X2+4.D0*U3*Y2*V3*X3-6.D0*U2**2*Y3**2+
00400      &             2.D0*V2*X3*U1*Y3+4.D0*V2*X3*U3*Y3-4.D0*V1*X3*U3*Y3-
00401      &             4.D0*U1*Y3*V1*X3-4.D0*U1*Y3*V3*X3+6.D0*U3**2*Y3**2+
00402      &             6.D0*V3**2*X3**2-4.D0*V2*X3**2*V3+4.D0*U1*Y3**2*U3-
00403      &             2.D0*V2*X3**2*V1-4.D0*U2*Y3**2*U3-2.D0*U2*Y3**2*U1+
00404      &             4.D0*U2*Y3*V3*X3+12.D0*U2*Y3*V2*X3-12.D0*U3*Y3*V3*X3+
00405      &             2.D0*U2*Y3*V1*X3)*AUX180
00406 !
00407       A33(IELEM) = (15.D0*U2*Y2**2*U3+7.D0*U2*Y2**2*U1+7.D0*U1**2*Y2**2+
00408      &             7.D0*V1**2*X2**2+7.D0*U2**2*Y2**2+7.D0*V2**2*X2**2+
00409      &             39.D0*V3**2*X2**2-15.D0*V1*X2*U3*Y2-7.D0*V2*X2*U1*Y2-
00410      &             15.D0*V2*X2*U3*Y2-7.D0*U2*Y2*V1*X2-15.D0*U2*Y2*V3*X2-
00411      &             14.D0*U1*Y2*V1*X2-15.D0*U1*Y2*V3*X2-
00412      &             14.D0*U2*Y2*V2*X2+15.D0*V2*X2**2*V3+
00413      &             39.D0*U3**2*Y2**2-78.D0*U3*Y2*V3*X2+
00414      &             7.D0*V2*X2**2*V1+15.D0*V1*X2**2*V3+
00415      &             15.D0*U1*Y2**2*U3)*AUX360
00416 !
00417       A34(IELEM) = (-(6.D0*(V3**2-V2**2)+4.D0*(V2*V3-V2*V1)-2.D0*(V1**2-
00418      &            V1*V3))*X2**2+(-(-2.D0*(V2*V3-V1*V3)+4.D0*(V2**2-
00419      &            V1**2))*X3-(-4.D0*(U2*V3-V2*U1+V2*U3-U1*V1-U2*V1)+
00420      &            12.D0*(U2*V2-U3*V3)-2.D0*(U1*V3+U3*V1))*Y2-
00421      &            (4.D0*(U1*V1-U2*V2)-U3*V1+V2*U3-U1*V3+U2*V3)*Y3)*X2-
00422      &            (4.D0*(U1*V1-U2*V2)-U3*V1+V2*U3-U1*V3+
00423      &            U2*V3)*Y2*X3-(-2.D0*(U1**2-U1*U3)-6.D0*(U2**2-U3**2)-
00424      &            4.D0*(U2*U1-U2*U3))*Y2**2-(4.D0*(U2**2-U1**2)+
00425      &            2.D0*(U1*U3-U2*U3))*Y3*Y2)*AUX180
00426 !
00427       A35(IELEM) = -(6.D0*U2**2*Y2**2+2.D0*U1**2*Y2**2+2.D0*V1**2*X2**2+
00428      &             6.D0*V2**2*X2**2-6.D0*V3**2*X2**2-4.D0*U2*Y2**2*U3+
00429      &             2.D0*U1**2*Y3*Y2+2.D0*V1**2*X3*X2+4.D0*U2*Y2**2*U1+
00430      &             2.D0*V2**2*X3*X2+2.D0*U2**2*Y3*Y2+30.D0*V3**2*X3*X2+
00431      &             30.D0*U3**2*Y3*Y2-6.D0*V1*X2*U3*Y3+2.D0*V1*X2*U3*Y2-
00432      &             V2*X3*U1*Y2-6.D0*V2*X3*U3*Y2+2.D0*V2*X3*V1*X2+
00433      &             12.D0*V2*X3*V3*X2-V2*X2*U1*Y3-4.D0*V2*X2*U1*Y2-
00434      &             6.D0*V2*X2*U3*Y3+4.D0*V2*X2*U3*Y2-U2*Y2*V1*X3-
00435      &             4.D0*U2*Y2*V1*X2-6.D0*U2*Y2*V3*X3+4.D0*U2*Y2*V3*X2-
00436      &             6.D0*V1*X3*U3*Y2+12.D0*V1*X3*V3*X2+12.D0*U1*Y3*U3*Y2-
00437      &             2.D0*U1*Y3*V1*X2-6.D0*U1*Y3*V3*X2-2.D0*U1*Y2*V1*X3-
00438      &             4.D0*U1*Y2*V1*X2-6.D0*U1*Y2*V3*X3+2.D0*U1*Y2*V3*X2+
00439      &             2.D0*U2*Y3*U1*Y2+12.D0*U2*Y3*U3*Y2-U2*Y3*V1*X2-
00440      &             6.D0*U2*Y3*V3*X2-2.D0*U2*Y3*V2*X2-2.D0*U2*Y2*V2*X3-
00441      &             12.D0*U2*Y2*V2*X2-30.D0*U3*Y3*V3*X2-
00442      &             30.D0*U3*Y2*V3*X3+12.D0*U3*Y2*V3*X2-6.D0*U3**2*Y2**2-
00443      &             4.D0*V2*X2**2*V3+4.D0*V2*X2**2*V1-2.D0*V1*X2**2*V3-
00444      &             2.D0*U1*Y2**2*U3)*AUX180
00445 !
00446       ANS1 = U3**2*Y2**2-2.D0*U3*Y2*V3*X2-2.D0*U1*Y2*V1*X2+V3**2*X2**2+
00447      &       V1**2*X2**2-3.D0*U2*Y2*V3*X2-3.D0*U2*Y2*V1*X2+U1**2*Y2**2-
00448      &       3.D0*V2*X2*U1*Y2-3.D0*V2*X2*U3*Y2-U3*Y2*V1*X2-U1*Y2*V3*X2+
00449      &       6.D0*V2**2*X2**2+6.D0*U2**2*Y2**2-12.D0*U2*Y2*V2*X2+
00450      &       V3**2*X3**2-2.D0*U2*Y3*U1*Y2+U3*Y3*V3*X2-2.D0*U3*Y3*V3*X3+
00451      &       V1**2*X3*X2-2.D0*U3*Y3*V1*X3+2.D0*U2*Y3**2*U3+U3*Y2*V3*X3-
00452      &       U1*Y2*V1*X3-U3**2*Y3*Y2+2.D0*U1*Y3**2*U3+2.D0*U2*Y3**2*U1-
00453      &       V3**2*X3*X2+2.D0*V1*X3**2*V3-2.D0*U1*Y3*V3*X3+U1**2*Y3*Y2-
00454      &       U1*Y3*V1*X2+2.D0*V2*X3**2*V1+2.D0*V2*X3**2*V3-
00455      &       9.D0*U2**2*Y3*Y2+U3**2*Y3**2+4.D0*U1**2*Y3**2+
00456      &       9.D0*U2*Y2*V2*X3-8.D0*U1*Y3*V1*X3-9.D0*V2**2*X3*X2+
00457      &       V1*X2**2*V3-8.D0*U2*Y3*V2*X3-2.D0*V2*X3*U1*Y3+
00458      &       V2*X3*U1*Y2+4.D0*V2**2*X3**2+4.D0*U2**2*Y3**2+
00459      &       9.D0*U2*Y3*V2*X2+2.D0*V2*X2*U3*Y3+3.D0*V2*X2**2*V3-
00460      &       4.D0*V2*X3*V3*X2+V2*X2*U1*Y3+U2*Y2*V1*X3+2.D0*V2*X3*U3*Y2-
00461      &       2.D0*V2*X3*V1*X2-2.D0*V2*X3*U3*Y3+4.D0*V1**2*X3**2+
00462      &       3.D0*U2*Y2**2.D0*U3+3.D0*V2*X2**2*V1-4.D0*U2*Y3*U3*Y2-
00463      &       2.D0*U2*Y3*V1*X3-2.D0*U2*Y3*V3*X3+2.D0*U2*Y2*V3*X3+
00464      &       U2*Y3*V1*X2+2.D0*U2*Y3*V3*X2+3.D0*U2*Y2**2*U1+U1*Y2**2*U3
00465       A44(IELEM) = ANS1*AUX45*2.D0
00466 !
00467       ANS1 = 2.D0*U3**2*Y2**2-4.D0*U3*Y2*V3*X2-4.D0*U1*Y2*V1*X2+
00468      &       2.D0*V3**2*X2**2+2.D0*V1**2*X2**2-6.D0*U2*Y2*V3*X2-
00469      &       6.D0*U2*Y2*V1*X2+2.D0*U1**2*Y2**2-6.D0*V2*X2*U1*Y2-
00470      &       6.D0*V2*X2*U3*Y2-2.D0*U3*Y2*V1*X2-2.D0*U1*Y2*V3*X2+
00471      &       12.D0*V2**2*X2**2+12.D0*U2**2*Y2**2-24.D0*U2*Y2*V2*X2-
00472      &       4.D0*U2*Y3*U1*Y2+4.D0*U3*Y3*V3*X2+U3*Y2*V1*X3-
00473      &       2.D0*V1*X3*V3*X2+2.D0*U3*Y3*V1*X3+U1*Y2*V3*X3+U3*Y3*V1*X2+
00474      &       2.D0*U2*Y3**2*U3+4.D0*U3*Y2*V3*X3-4.D0*U3**2*Y3*Y2-
00475      &       2.D0*U1*Y3**2*U3-4.D0*V3**2*X3*X2-2.D0*V1*X3**2*V3+
00476      &       2.D0*U1*Y3*V3*X3+U1*Y3*V3*X2-2.D0*U1*Y3*U3*Y2+
00477      &       2.D0*V2*X3**2*V3-12.D0*U2**2*Y3*Y2-2.D0*U1**2*Y3**2+
00478      &       12.D0*U2*Y2*V2*X3+4.D0*U1*Y3*V1*X3-12.D0*V2**2*X3*X2+
00479      &       2.D0*V1*X2**2*V3-4.D0*U2*Y3*V2*X3+2.D0*V2*X3*U1*Y2+
00480      &       2.D0*V2**2*X3**2+2.D0*U2**2*Y3**2+12.D0*U2*Y3*V2*X2+
00481      &       4.D0*V2*X2*U3*Y3+6.D0*V2*X2**2*V3-8.D0*V2*X3*V3*X2+
00482      &       2.D0*V2*X2*U1*Y3+2.D0*U2*Y2*V1*X3+4.D0*V2*X3*U3*Y2-
00483      &       4.D0*V2*X3*V1*X2-2.D0*V2*X3*U3*Y3-2.D0*V1**2*X3**2+
00484      &       6.D0*U2*Y2**2*U3+6.D0*V2*X2**2*V1-8.D0*U2*Y3*U3*Y2-
00485      &       2.D0*U2*Y3*V3*X3+4.D0*U2*Y2*V3*X3+2.D0*U2*Y3*V1*X2+
00486      &       4.D0*U2*Y3*V3*X2+6.D0*U2*Y2**2*U1+2.D0*U1*Y2**2*U3
00487       A45(IELEM) = ANS1*(-AUX45)
00488 !
00489       ANS1 = U3**2*Y2**2-2.D0*U3*Y2*V3*X2-2.D0*U1*Y2*V1*X2+V3**2*X2**2+
00490      &       V1**2*X2**2-3.D0*U2*Y2*V3*X2-3.D0*U2*Y2*V1*X2+U1**2*Y2**2-
00491      &       3.D0*V2*X2*U1*Y2-3.D0*V2*X2*U3*Y2-U3*Y2*V1*X2-U1*Y2*V3*X2+
00492      &       6.D0*V2**2*X2**2+6.D0*U2**2*Y2**2-12.D0*U2*Y2*V2*X2+
00493      &       6.D0*V3**2*X3**2-2.D0*U2*Y3*U1*Y2+3.D0*U3*Y3*V3*X2+
00494      &       U3*Y2*V1*X3-2.D0*V1*X3*V3*X2-12.D0*U3*Y3*V3*X3-V1**2*X3*X2-
00495      &       3.D0*U3*Y3*V1*X3+U1*Y2*V3*X3+U3*Y3*V1*X2+3.D0*U2*Y3**2*U3+
00496      &       3.D0*U3*Y2*V3*X3+U1*Y2*V1*X3-3.D0*U3**2*Y3*Y2+
00497      &       3.D0*U1*Y3**2*U3+U2*Y3**2*U1-3.D0*V3**2*X3*X2+
00498      &       3.D0*V1*X3**2*V3-3.D0*U1*Y3*V3*X3-U1**2*Y3*Y2+U1*Y3*V1*X2+
00499      &       U1*Y3*V3*X2-2.D0*U1*Y3*U3*Y2+V2*X3**2*V1+3.D0*V2*X3**2*V3-
00500      &       3.D0*U2**2*Y3*Y2+6.D0*U3**2*Y3**2+U1**2*Y3**2+
00501      &       3.D0*U2*Y2*V2*X3-2.D0*U1*Y3*V1*X3-3.D0*V2**2*X3*X2+
00502      &       V1*X2**2*V3-2.D0*U2*Y3*V2*X3-V2*X3*U1*Y3+V2*X3*U1*Y2+
00503      &       V2**2*X3**2+U2**2*Y3**2+3.D0*U2*Y3*V2*X2+2.D0*V2*X2*U3*Y3+
00504      &       3.D0*V2*X2**2*V3-4.D0*V2*X3*V3*X2+V2*X2*U1*Y3+U2*Y2*V1*X3+
00505      &       2.D0*V2*X3*U3*Y2-2.D0*V2*X3*V1*X2-3.D0*V2*X3*U3*Y3+
00506      &       V1**2*X3**2+3.D0*U2*Y2**2*U3+3.D0*V2*X2**2*V1-
00507      &       4.D0*U2*Y3*U3*Y2-U2*Y3*V1*X3-3.D0*U2*Y3*V3*X3+
00508      &       2.D0*U2*Y2*V3*X3+U2*Y3*V1*X2+2.D0*U2*Y3*V3*X2
00509       A55(IELEM) = (3.D0*U2*Y2**2*U1+U1*Y2**2*U3+ANS1)*2.D0*AUX45
00510 !
00511 !
00512 ! USES HERE THE 'MAGIC SQUARE' PROPERTIES
00513 ! (SUM OF EACH LINE = SUM OF EACH COLUMN = 0)
00514 !
00515           A16(IELEM) = - A11(IELEM) - A12(IELEM) - A13(IELEM)
00516      &                 - A14(IELEM) - A15(IELEM)
00517 !
00518           A26(IELEM) = - A12(IELEM) - A22(IELEM) - A23(IELEM)
00519      &                 - A24(IELEM) - A25(IELEM)
00520 !
00521           A36(IELEM) = - A13(IELEM) - A23(IELEM) - A33(IELEM)
00522      &                 - A34(IELEM) - A35(IELEM)
00523 !
00524           A46(IELEM) = - A14(IELEM) - A24(IELEM) - A34(IELEM)
00525      &                 - A44(IELEM) - A45(IELEM)
00526 !
00527           A56(IELEM) = - A15(IELEM) - A25(IELEM) - A35(IELEM)
00528      &                 - A45(IELEM) - A55(IELEM)
00529 !
00530           A66(IELEM) = - A16(IELEM) - A26(IELEM) - A36(IELEM)
00531      &                 - A46(IELEM) - A56(IELEM)
00532 !
00533         ENDDO ! IELEM
00534 !
00535 !-----------------------------------------------------------------------
00536 !
00537       ELSEIF(IELMU.EQ.13.AND.IELMV.EQ.13) THEN
00538 !
00539 !-----------------------------------------------------------------------
00540 !
00541 !  P2 DISCRETISATION OF THE VELOCITY:
00542 !
00543       DO IELEM = 1 , NELEM
00544 !
00545 !   INITIALISES THE GEOMETRICAL VARIABLES
00546 !
00547         X2  =  XEL(IELEM,2)
00548         X3  =  XEL(IELEM,3)
00549 !
00550         Y2  =  YEL(IELEM,2)
00551         Y3  =  YEL(IELEM,3)
00552 !
00553         U1 = U(IKLE1(IELEM))
00554         U2 = U(IKLE2(IELEM))
00555         U3 = U(IKLE3(IELEM))
00556         U4 = U(IKLE4(IELEM))
00557         U5 = U(IKLE5(IELEM))
00558         U6 = U(IKLE6(IELEM))
00559         V1 = V(IKLE1(IELEM))
00560         V2 = V(IKLE2(IELEM))
00561         V3 = V(IKLE3(IELEM))
00562         V4 = V(IKLE4(IELEM))
00563         V5 = V(IKLE5(IELEM))
00564         V6 = V(IKLE6(IELEM))
00565 !
00566 !   INITIALISES THE INTERMEDIATE VARIABLES
00567 !
00568         UNSU2 = XMUL/(X2*Y3-Y2*X3)
00569         AUX630  = UNSU2/630.D0
00570         AUX1260 = UNSU2/1260.D0
00571         AUX2520 = UNSU2/2520.D0
00572 !
00573 !  COMPUTES 15 OF THE 36 TERMS (SELECTED AMONG THE LEAST COMPLEX)
00574 !
00575       ANS1 = 112.D0*V4**2*X2**2-96.D0*U5*Y2*V5*X2+23.D0*V1*X2*U3*Y2+
00576      &       224.D0*U4*Y2*V4*X3+48.D0*U5*Y3*V4*X2+48.D0*V5*X2*U4*Y3+
00577      &       48.D0*U5**2*Y2**2-224.D0*U4*Y2*V4*X2-96.D0*U1*Y3*V6*X3+
00578      &       96.D0*U5*Y3*V5*X2+48.D0*V5*X3*U4*Y2-48.D0*U5*Y3*V4*X3+
00579      &       13.D0*V2**2*X3**2+48.D0*V5**2*X3**2+96.D0*U5*Y2*V5*X3-
00580      &       48.D0*V5*X3*U4*Y3+48.D0*U5*Y2*V4*X3-96.D0*V5*X3*V4*X2+
00581      &       224.D0*V4*X2*U4*Y3-224.D0*V4*X3*U4*Y3+13.D0*U2**2*Y3**2-
00582      &       32.D0*V2*X2*U4*Y3+112.D0*U4**2*Y3**2-192.D0*V1*X3*V6*X2-
00583      &       186.D0*V1*X3*U1*Y3+23.D0*U1*Y3*V3*X3-28.D0*U2*Y2**2*U6+
00584      &       13.D0*V2**2*X2**2-96.D0*U4*Y2*V1*X2-112.D0*U4*Y2*V6*X2+
00585      &       96.D0*U1*Y3*V6*X2-112.D0*V4*X3*U6*Y3+112.D0*V4**2*X3**2-
00586      &       112.D0*U4*Y3*V6*X3-96.D0*V4*X2*U1*Y2-112.D0*V4*X2*U6*Y2+
00587      &       112.D0*U6**2*Y3**2-32.D0*U2*Y2*V4*X3+32.D0*U2*Y2*V4*X2+
00588      &       28.D0*U2*Y3*V6*X3+13.D0*U2**2*Y2**2+13.D0*U3**2*Y3**2-
00589      &       96.D0*V5*X3*U5*Y3+48.D0*U5*Y3**2*U6-96.D0*U5*Y2*U4*Y3+
00590      &       93.D0*U1**2*Y3**2+93.D0*V1**2*X3**2+13.D0*V3**2*X3**2+
00591      &       112.D0*V6**2*X3**2+48.D0*U5*Y2**2*U6-224.D0*V6**2*X3*X2-
00592      &       26.D0*V3**2*X3*X2-26.D0*V2**2*X2*X3-224.D0*U6**2*Y2*Y3+
00593      &       48.D0*U5**2*Y3**2+48.D0*V5*X2**2*V4-32.D0*V2*X2**2*V4+
00594      &       4.D0*U5*Y2**2*U1+32.D0*V2*X2*U4*Y2+112.D0*U6**2*Y2**2
00595       ANS2 = -32.D0*V2*X3**2*V4-28.D0*V2*X2**2*V6-
00596      &       32.D0*U2*Y2**2*U4+V2*X2**2*V3-28.D0*V4*X2**2*V3-
00597      &       26.D0*V2*X2*U2*Y2-32.D0*U3*Y3**2*U6+93.D0*V1**2*X2**2-
00598      &       186.D0*U1**2*Y3*Y2+96.D0*U1*Y2**2*U6+112.D0*V6**2*X2**2+
00599      &       4.D0*U5*Y3**2*U1+V2*X3**2*V3-26.D0*U2**2*Y2*Y3+
00600      &       96.D0*V4*X2**2*V1+4.D0*V5*X2**2*V1+96.D0*V4*X3**2*V1+
00601      &       96.D0*V1*X2**2*V6-U2*Y2*V3*X2+48.D0*V5**2*X2**2-
00602      &       23.D0*V1*X3**2*V3-V2*X3*U3*Y3+13.D0*V3**2*X2**2-
00603      &       28.D0*U4*Y2**2*U3-U2*Y3*V3*X3-28.D0*U2*Y3**2*U6-
00604      &       23.D0*U2*Y3**2*U1-23.D0*U1*Y3**2*U3+112.D0*U4*Y3**2*U6-
00605      &       186.D0*V1*X2*U1*Y2-96.D0*V1*X3*U6*Y3-224.D0*U4**2*Y2*Y3-
00606      &       224.D0*V4**2*X3*X2-23.D0*U1*Y2**2*U3-32.D0*V6*X3**2*V3-
00607      &       48.D0*V5*X2*U4*Y2-48.D0*U5*Y2*V4*X2-26.D0*U3**2*Y3*Y2-
00608      &       186.D0*V1**2*X3*X2-23.D0*V2*X3**2*V1-V2*X2*U3*Y2+
00609      &       96.D0*U4*Y2**2*U1-23.D0*V1*X2**2*V3+48.D0*V5*X2**2*V6+
00610      &       U2*Y2**2*U3-28.D0*U4*Y3**2*U3-32.D0*U2*Y3**2*U4+
00611      &       4.D0*V5*X3**2*V1-96.D0*U5**2*Y2*Y3-96.D0*V5**2*X3*X2+
00612      &       48.D0*V5*X3**2*V4+48.D0*U5*Y3**2*U4+13.D0*U3**2*Y2**2+
00613      &       224.D0*V6*X2*U6*Y3-224.D0*V6*X2*U6*Y2+32.D0*V6*X2*U3*Y2-
00614      &       32.D0*V3*X2**2*V6-32.D0*U2*Y3*V4*X2+32.D0*U2*Y3*V4*X3
00615       ANS3 = 48.D0*V5*X3**2*V6-26.D0*V3*X2*U3*Y2+32.D0*V3*X2*U6*Y2-
00616      &       32.D0*V2*X3*U4*Y2+93.D0*U1**2*Y2**2+23.D0*U1*Y2*V3*X2+
00617      &       32.D0*V2*X3*U4*Y3+224.D0*V6*X3*U6*Y2-96.D0*U1*Y2*V6*X2+
00618      &       28.D0*V4*X2*U3*Y2+26.D0*V2*X2*U2*Y3+23.D0*V1*X3*U3*Y3+
00619      &       26.D0*U2*Y2*V2*X3-26.D0*U2*Y3*V2*X3+112.D0*V4*X3**2*V6+
00620      &       4.D0*U5*Y3*V1*X2-4.D0*U5*Y3*V1*X3+U2*Y3**2*U3-
00621      &       192.D0*U1*Y3*U6*Y2+112.D0*U4*Y2**2*U6+96.D0*V1*X3*U6*Y2+
00622      &       112.D0*U4**2*Y2**2-96.D0*U5*Y2*U6*Y3+46.D0*U1*Y3*U3*Y2+
00623      &       96.D0*V1*X3**2*V6+28.D0*U4*Y2*V3*X2+48.D0*V5*X3*U6*Y2-
00624      &       23.D0*U1*Y3*V3*X2-4.D0*V5*X3*U1*Y3-8.D0*V5*X3*V1*X2-
00625      &       48.D0*V5*X2*U6*Y2+96.D0*U1*Y3**2*U6-4.D0*V5*X2*U1*Y2-
00626      &       8.D0*U5*Y2*U1*Y3+48.D0*U5*Y2**2*U4+26.D0*V3*X3*U3*Y2+
00627      &       26.D0*V3*X2*U3*Y3-32.D0*V3*X2*U6*Y3+64.D0*U3*Y3*U6*Y2-
00628      &       96.D0*V5*X3*V6*X2-32.D0*U3*Y3*V6*X2+186.D0*V1*X2*U1*Y3+
00629      &       96.D0*V1*X2*U6*Y3+4.D0*V5*X3*U1*Y2+46.D0*V1*X3*V3*X2+
00630      &       96.D0*U4*Y3**2*U1+64.D0*V6*X3*V3*X2-26.D0*V3*X3*U3*Y3-
00631      &       224.D0*V6*X3*U6*Y3+4.D0*U5*Y2*V1*X3+112.D0*V4*X2**2*V6+
00632      &       32.D0*V3*X3*U6*Y3-32.D0*V3*X3*U6*Y2+48.D0*U5*Y2*V6*X3+
00633      &       46.D0*U2*Y3*U1*Y2+112.D0*U4*Y3*V6*X2+U2*Y3*V3*X2-
00634      &       28.D0*U4*Y3*V3*X2+96.D0*U4*Y3*V1*X2-96.D0*U4*Y3*V1*X3
00635       ANS4 = -23.D0*V1*X3*U3*Y2-23.D0*U1*Y2*V3*X3-23.D0*U2*Y3*V1*X2+
00636      &       96.D0*U1*Y2*V6*X3-32.D0*V6*X3*U3*Y2+32.D0*V6*X3*U3*Y3-
00637      &       28.D0*V2*X3*U6*Y2-192.D0*V4*X3*V1*X2-28.D0*V4*X2*U3*Y3+
00638      &       V2*X3*U3*Y2+96.D0*V4*X2*U1*Y3+56.D0*V2*X3*V6*X2-
00639      &       48.D0*U5*Y2*V6*X2+23.D0*U2*Y3*V1*X3+56.D0*U2*Y3*U6*Y2+
00640      &       112.D0*V4*X2*U6*Y3+28.D0*U4*Y3*V3*X3-2.D0*U2*Y3*U3*Y2-
00641      &       28.D0*U2*Y3*V6*X2-23.D0*V2*X3*U1*Y2-96.D0*V1*X2*U6*Y2+
00642      &       112.D0*V4*X3*U6*Y2-28.D0*V4*X3*U3*Y2-2.D0*V2*X3*V3*X2-
00643      &       224.D0*V4*X3*V6*X2+96.D0*V4*X3*U1*Y2+28.D0*V4*X3*U3*Y3+
00644      &       46.D0*V2*X3*V1*X2-28.D0*V2*X3**2*V6-4.D0*U5*Y2*V1*X2+
00645      &       23.D0*V2*X3*U1*Y3+28.D0*V2*X3*U6*Y3-28.D0*V2*X2*U6*Y3+
00646      &       112.D0*U4*Y2*V6*X3+V2*X2*U3*Y3-32.D0*U3*Y2**2*U6+
00647      &       48.D0*V5*X2*U6*Y3-23.D0*V1*X2*U3*Y3+56.D0*V4*X3*V3*X2-
00648      &       28.D0*V4*X3**2*V3-96.D0*V4*X3*U1*Y3-224.D0*U4*Y2*U6*Y3+
00649      &       23.D0*U2*Y2*V1*X2-23.D0*V2*X2*U1*Y3+56.D0*U4*Y2*U3*Y3+
00650      &       96.D0*U4*Y2*V1*X3-28.D0*U4*Y2*V3*X3+U2*Y2*V3*X3-
00651      &       48.D0*U5*Y3*V6*X3-48.D0*V5*X3*U6*Y3+4.D0*V5*X2*U1*Y3-
00652      &       28.D0*U2*Y2*V6*X3+28.D0*U2*Y2*V6*X2+23.D0*V2*X2*U1*Y2-
00653      &       192.D0*U4*Y2*U1*Y3+48.D0*U5*Y3*V6*X2+28.D0*V2*X2*U6*Y2-
00654      &       23.D0*U2*Y2*V1*X3-23.D0*V2*X2**2*V1+186.D0*V1*X3*U1*Y2
00655       A11(IELEM) = (64.D0*U2*Y3*U4*Y2+64.D0*V2*X3*V4*X2+
00656      &              ANS1+ANS2+ANS3+ANS4-23.D0*U2*Y2**2*U1)*AUX1260
00657 !
00658       ANS1 =(8.D0*U2*U6+8.D0*U5*U1-32.D0*U5*U6+8.D0*U4*U3-32.D0*U4*U6+
00659      &      10.D0*U3*U1-32.D0*U4*U1+32.D0*U6*U3-54.D0*U2**2-
00660      &      64.D0*U5*U2+32.D0*U3*U5-96.D0*U5**2+18.D0*U2*U1+
00661      &      10.D0*U3**2+10.D0*U2*U3-64.D0*U6*U1-32.D0*U4*U5-
00662      &      32.D0*U2*U4-54.D0*U1**2+32.D0*U4**2-96.D0*U6**2)*Y3*Y2+
00663      &      (-32.D0*U4**2-32.D0*U3*U5+32.D0*U5*U6-8.D0*U2*U6+
00664      &      32.D0*U2*U4+54.D0*U2**2-18.D0*U2*U1+32.D0*U4*U6-
00665      &      8.D0*U5*U1+32.D0*U4*U1-32.D0*U6*U3-10.D0*U2*U3+64.D0*U5*U2-
00666      &      10.D0*U3**2+32.D0*U4*U5+96.D0*U5**2+96.D0*U6**2-8.D0*U4*U3-
00667      &      10.D0*U3*U1+64.D0*U6*U1+54.D0*U1**2)*Y3**2+((-96.D0*V5**2+
00668      &      10.D0*V2*V3+8.D0*V4*V3+8.D0*V2*V6-64.D0*V5*V2-54.D0*V1**2+
00669      &      32.D0*V3*V5+32.D0*V4**2+32.D0*V3*V6-96.D0*V6**2+
00670      &      10.D0*V3**2-32.D0*V5*V6-32.D0*V4*V6-32.D0*V4*V5-
00671      &      32.D0*V4*V1-64.D0*V1*V6+8.D0*V5*V1+18.D0*V2*V1-
00672      &      32.D0*V2*V4+10.D0*V1*V3-54.D0*V2**2)*X3+(96.D0*U5*V5+
00673      &      54.D0*V1*U1-5.D0*U2*V3+16.D0*U4*V1-16.D0*U6*V3-4.D0*V2*U6-
00674      &      4.D0*V4*U3-5.D0*V1*U3-5.D0*U1*V3+54.D0*U2*V2+16.D0*V4*U5-
00675      &      4.D0*U5*V1+16.D0*U4*V5-4.D0*U2*V6-4.D0*V5*U1-4.D0*U4*V3+
00676      &      96.D0*U6*V6+16.D0*V4*U6+16.D0*U4*V6+16.D0*U5*V6)*Y3)*X2
00677       ANS2 =(16.D0*V5*U6+32.D0*U5*V2-16.D0*U5*V3-5.D0*V2*U3+
00678      &      16.D0*U2*V4+32.D0*U1*V6-9.D0*V2*U1-16.D0*U3*V6-9.D0*U2*V1+
00679      &      16.D0*V2*U4+16.D0*V4*U1-10.D0*U3*V3+32.D0*U2*V5-16.D0*U3*V5-
00680      &      32.D0*U4*V4+32.D0*V1*U6)*Y3*X2+X3*Y2*(96.D0*U5*V5+
00681      &      54.D0*V1*U1-5.D0*U2*V3+16.D0*U4*V1-16.D0*U6*V3-4.D0*V2*U6-
00682      &      4.D0*V4*U3-5.D0*V1*U3-5.D0*U1*V3+54.D0*U2*V2+16.D0*V4*U5-
00683      &      4.D0*U5*V1+16.D0*U4*V5-4.D0*U2*V6-4.D0*V5*U1-4.D0*U4*V3+
00684      &      96.D0*U6*V6+16.D0*V4*U6+16.D0*U4*V6+16.D0*U5*V6+
00685      &      16.D0*V5*U6+32.D0*U5*V2-16.D0*U5*V3-5.D0*V2*U3+16.D0*U2*V4+
00686      &      32.D0*U1*V6-9.D0*V2*U1-16.D0*U3*V6-9.D0*U2*V1+16.D0*V2*U4+
00687      &      16.D0*V4*U1-10.D0*U3*V3+32.D0*U2*V5-16.D0*U3*V5-
00688      &      32.D0*U4*V4+32.D0*V1*U6)+(32.D0*V4*V5-18.D0*V2*V1+
00689      &      32.D0*V4*V1-10.D0*V1*V3+32.D0*V5*V6- 8.D0*V5*V1+
00690      &      32.D0*V2*V4-10.D0*V2*V3-8.D0*V4*V3+32.D0*V4*V6+
00691      &      54.D0*V2**2-8.D0*V2*V6+96.D0*V5**2+54.D0*V1**2-32.D0*V3*V6-
00692      &      32.D0*V4**2+96.D0*V6**2+64.D0*V1*V6-10.D0*V3**2+
00693      &      64.D0*V5*V2-32.D0*V3*V5)*X3**2+(-192.D0*U5*V5-
00694      &      108.D0*V1*U1+10.D0*U2*V3-32.D0*U4*V1+32.D0*U6*V3+
00695      &      8.D0*V2*U6+8.D0*V4*U3+10.D0*V1*U3+10.D0*U1*V3-
00696      &      108.D0*U2*V2-32.D0*V4*U5+8.D0*U5*V1-32.D0 *U4*V5)*Y3*X3
00697       A12(IELEM) = ((8.D0*U2*V6+8.D0*V5*U1+8.D0*U4*V3-192.D0*U6*V6-
00698      &             32.D0*V4*U6-32.D0*U4*V6-32.D0*U5*V6-32.D0*V5*U6-
00699      &             64.D0*U5*V2+32.D0*U5*V3+10.D0*V2*U3-32.D0*U2*V4-
00700      &             64.D0*U1*V6+18.D0*V2*U1+32.D0*U3*V6+18.D0*U2*V1-
00701      &             32.D0*V2*U4-32.D0*V4*U1+20.D0*U3*V3-64.D0*U2*V5+
00702      &             32.D0*U3*V5+64.D0*U4*V4-64.D0*V1*U6)*Y3*X3+
00703      &             ANS1+ANS2)*AUX2520
00704 !
00705       ANS1 = -96.D0*V4**2*X2**2+192.D0*U5*Y2*V5*X2-18.D0*V1*X2*U3*Y2-
00706      &       96.D0*U4*Y2*V4*X3-16.D0*U5*Y3*V4*X2-16.D0*V5*X2*U4*Y3-
00707      &       96.D0*U5**2*Y2**2+192.D0*U4*Y2*V4*X2-96.D0*U5*Y3*V5*X2-
00708      &       16.D0*V5*X3*U4*Y2-96.D0*U5*Y2*V5*X3-16.D0*U5*Y2*V4*X3+
00709      &       32.D0*V5*X3*V4*X2-96.D0*V4*X2*U4*Y3+16.D0*V2*X2*U4*Y3+
00710      &       32.D0*V1*X3*V6*X2+8.D0*U2*Y2**2*U6+10.D0*V2**2*X2**2+
00711      &       16.D0*U2*Y2*V5*X3+16.D0*U2*Y3*V5*X2+64.D0*U4*Y2*V1*X2+
00712      &       32.D0*U4*Y2*V6*X2-16.D0*U1*Y3*V6*X2-32.D0*U2*Y3*U5*Y2+
00713      &       64.D0*V4*X2*U1*Y2+32.D0*V4*X2*U6*Y2+16.D0*U2*Y2*V4*X3-
00714      &       32.D0*U2*Y2*V4*X2+10.D0*U2**2*Y2**2+32.D0*U5*Y2*U4*Y3+
00715      &       32.D0*V5*X2**2*V2-32.D0*U5*Y2**2*U6-32.D0*V6**2*X3*X2+
00716      &       54.D0*V3**2*X3*X2-10.D0*V2**2*X2*X3-32.D0*U6**2*Y2*Y3-
00717      &       32.D0*V5*X2**2*V4+32.D0*V2*X2**2*V4+8.D0*U5*Y2**2*U1-
00718      &       32.D0*V2*X2*U4*Y2+32.D0*U6**2*Y2**2+8.D0*V2*X2**2*V6+
00719      &       32.D0*U2*Y2**2*U4+10.D0*V2*X2**2*V3+8.D0*V4*X2**2*V3-
00720      &       20.D0*V2*X2*U2*Y2-54.D0*V1**2*X2**2+54.D0*U1**2*Y3*Y2-
00721      &       32.D0*U1*Y2**2*U6+32.D0*V6**2*X2**2-10.D0*U2**2*Y2*Y3-
00722      &       64.D0*V4*X2**2*V1+8.D0*V5*X2**2*V1-32.D0*V5*X2*U2*Y2-
00723      &       32.D0*U5*Y2*V2*X2-64.D0*V3*X2**2*V5-32.D0*V1*X2**2*V6
00724       ANS2 = -10.D0*U2*Y2*V3*X2-96.D0*V5**2*X2**2-54.D0*V3**2*X2**2+
00725      &        8.D0*U4*Y2**2*U3+64.D0*V3*X2*U5*Y2-64.D0*U3*Y2**2*U5+
00726      &        108.D0*V1*X2*U1*Y2+96.D0*U4**2*Y2*Y3+96.D0*V4**2*X3*X2+
00727      &        18.D0*U1*Y2**2*U3+64.D0*U3*Y2*V5*X2+32.D0*V5*X2*U4*Y2+
00728      &        32.D0*U5*Y2*V4*X2+54.D0*U3**2*Y3*Y2+32.D0*U5*Y2**2*U2+
00729      &        54.D0*V1**2*X3*X2-10.D0*V2*X2*U3*Y2-64.D0*U4*Y2**2*U1+
00730      &        18.D0*V1*X2**2*V3-32.D0*V5*X2**2*V6+10.D0*U2*Y2**2*U3+
00731      &        96.D0*U5**2*Y2*Y3+96.D0*V5**2*X3*X2-54.D0*U3**2*Y2**2+
00732      &        32.D0*V6*X2*U6*Y3-64.D0*V6*X2*U6*Y2+32.D0*V6*X2*U3*Y2+
00733      &        16.D0*V2*X3*U5*Y2-32.D0*V3*X2**2*V6+16.D0*U2*Y3*V4*X2-
00734      &        32.D0*U2*Y3*U4*Y2-32.D0*V2*X3*V4*X2+108.D0*V3*X2*U3*Y2+
00735      &        32.D0*V3*X2*U6*Y2+16.D0*V2*X3*U4*Y2-32.D0*V2*X3*V5*X2-
00736      &        54.D0*U1**2*Y2**2-18.D0*U1*Y2*V3*X2+32.D0*V6*X3*U6*Y2+
00737      &        32.D0*U1*Y2*V6*X2-8.D0*V4*X2*U3*Y2+10.D0*V2*X2*U2*Y3+
00738      &        10.D0*U2*Y2*V2*X3+16.D0*V2*X2*U5*Y3-32.D0*V5*X3*U3*Y2+
00739      &        4.D0*U5*Y3*V1*X2+32.D0*U1*Y3*U6*Y2-32.D0*U4*Y2**2*U6-
00740      &        16.D0*V1*X3*U6*Y2-96.D0*U4**2*Y2**2+32.D0*U5*Y2*U6*Y3-
00741      &        18.D0*U1*Y3*U3*Y2+64.D0*U5*Y2*U3*Y3-8.D0*U4*Y2*V3*X2-
00742      &        16.D0*V5*X3*U6*Y2+9.D0*U1*Y3*V3*X2+64.D0*V5*X3*V3*X2-
00743      &        8.D0*V5*X3*V1*X2+32.D0*V5*X2*U6*Y2-8.D0*V5*X2*U1*Y2
00744       ANS3 = -8.D0*U5*Y2*U1*Y3-32.D0*U5*Y2**2*U4-54.D0*V3*X3*U3*Y2-
00745      &       54.D0*V3*X2*U3*Y3-16.D0*V3*X2*U6*Y3+32.D0*U3*Y3*U6*Y2+
00746      &       32.D0*V5*X3*V6*X2-16.D0*U3*Y3*V6*X2-54.D0*V1*X2*U1*Y3-
00747      &       16.D0*V1*X2*U6*Y3+4.D0*V5*X3*U1*Y2-18.D0*V1*X3*V3*X2+
00748      &       32.D0*V6*X3*V3*X2+4.D0*U5*Y2*V1*X3-32.D0*U5*Y2*V3*X3-
00749      &       32.D0*V4*X2**2*V6-16.D0*V3*X3*U6*Y2-16.D0*U5*Y2*V6*X3-
00750      &       10.D0*U2*Y3*U1*Y2-16.D0*U4*Y3*V6*X2+5.D0*U2*Y3*V3*X2+
00751      &       4.D0*U4*Y3*V3*X2-32.D0*U4*Y3*V1*X2+9.D0*V1*X3*U3*Y2+
00752      &       9.D0*U1*Y2*V3*X3+5.D0*U2*Y3*V1*X2-16.D0*U1*Y2*V6*X3-
00753      &       16.D0*V6*X3*U3*Y2+4.D0*V2*X3*U6*Y2+64.D0*V4*X3*V1*X2+
00754      &       4.D0*V4*X2*U3*Y3+5.D0*V2*X3*U3*Y2-32.D0*V4*X2*U1*Y3-
00755      &       8.D0*V2*X3*V6*X2+32.D0*U5*Y2*V6*X2-8.D0*U2*Y3*U6*Y2-
00756      &       16.D0*V4*X2*U6*Y3-10.D0*U2*Y3*U3*Y2+4.D0*U2*Y3*V6*X2+
00757      &       5.D0*V2*X3*U1*Y2+32.D0*V1*X2*U6*Y2-16.D0*V4*X3*U6*Y2+
00758      &       4.D0*V4*X3*U3*Y2-10.D0*V2*X3*V3*X2+32.D0*V4*X3*V6*X2-
00759      &       32.D0*V4*X3*U1*Y2-10.D0*V2*X3*V1*X2-8.D0*U5*Y2*V1*X2+
00760      &       4.D0*V2*X2*U6*Y3-16.D0*U4*Y2*V6*X3+5.D0*V2*X2*U3*Y3-
00761      &       32.D0*U3*Y2**2*U6-16.D0*V5*X2*U6*Y3+9.D0*V1*X2*U3*Y3-
00762      &       8.D0*V4*X3*V3*X2+32.D0*U4*Y2*U6*Y3-10.D0*U2*Y2*V1*X2
00763       A13(IELEM) = (5.D0*V2*X2*U1*Y3-8.D0*U4*Y2*U3*Y3-32.D0*U4*Y2*V1*X3+
00764      &             4.D0*U4*Y2*V3*X3+5.D0*U2*Y2*V3*X3-32.D0*V5*X2*U3*Y3+
00765      &             4.D0*V5*X2*U1*Y3+4.D0*U2*Y2*V6*X3-8.D0*U2*Y2*V6*X2-
00766      &             10.D0*V2*X2*U1*Y2+64.D0*U4*Y2*U1*Y3-
00767      &             16.D0*U5*Y3*V6*X2-8.D0*V2*X2*U6*Y2-
00768      &             32.D0*U5*Y3*V3*X2+5.D0*U2*Y2*V1*X3+
00769      &             10.D0*V2*X2**2*V1-54.D0*V1*X3*U1*Y2+
00770      &             10.D0*U2*Y2**2*U1+ANS1+ANS2+ANS3)*(-AUX2520)
00771 !
00772       ANS1 = -48.D0*V4**2*X2**2-96.D0*U5*Y2*V5*X2-2.D0*V1*X2*U3*Y2+
00773      &       16.D0*U5*Y3*V4*X2+16.D0*V5*X2*U4*Y3+48.D0*U5**2*Y2**2+
00774      &       96.D0*U4*Y2*V4*X2-64.D0*U1*Y3*V6*X3+96.D0*U5*Y3*V5*X2+
00775      &       16.D0*V5*X3*U4*Y2-32.D0*U5*Y3*V4*X3+20.D0*V2**2*X3**2+
00776      &       48.D0*V5**2*X3**2+96.D0*U5*Y2*V5*X3-32.D0*V5*X3*U4*Y3+
00777      &       16.D0*U5*Y2*V4*X3-32.D0*V5*X3*V4*X2-96.D0*V4*X3*U4*Y3+
00778      &       20.D0*U2**2*Y3**2+8.D0*V2*X2*U4*Y3+48.D0*U4**2*Y3**2-
00779      &       56.D0*V1*X3*V6*X2-120.D0*V1*X3*U1*Y3+14.D0*U1*Y3*V3*X3+
00780      &       8.D0*U2*Y2**2*U6+18.D0*V2**2*X2**2+20.D0*U2*Y2*V5*X3-
00781      &       16.D0*U2*Y3*V5*X3+20.D0*U2*Y3*V5*X2+16.D0*U4*Y2*V1*X2+
00782      &       32.D0*U4*Y2*V6*X2+28.D0*U1*Y3*V6*X2-40.D0*U2*Y3*U5*Y2-
00783      &       64.D0*V4*X3*U6*Y3+48.D0*V4**2*X3**2-64.D0*U4*Y3*V6*X3+
00784      &       16.D0*V4*X2*U1*Y2+32.D0*V4*X2*U6*Y2+80.D0*U6**2*Y3**2+
00785      &       8.D0*U2*Y2*V4*X3-24.D0*U2*Y2*V4*X2+16.D0*U2*Y3*V6*X3+
00786      &       18.D0*U2**2*Y2**2+4.D0*U3**2*Y3**2-96.D0*V5*X3*U5*Y3+
00787      &       32.D0*U5*Y3**2*U6-32.D0*U5*Y2*U4*Y3+24.D0*V5*X2**2*V2+
00788      &       60.D0*U1**2*Y3**2+60.D0*V1**2*X3**2+4.D0*V3**2*X3**2+
00789      &       80.D0*V6**2*X3**2-64.D0*V6**2*X3*X2-6.D0*V3**2*X3*X2-
00790      &       8.D0*U5*Y3**2*U3-38.D0*V2**2*X2*X3-64.D0*U6**2*Y2*Y3+
00791      &       48.D0*U5**2*Y3**2+24.D0*V2*X2**2*V4-8.D0*U5*Y2**2*U1
00792       ANS2 = -24.D0*V2*X2*U4*Y2-16.D0*U6**2*Y2**2-8.D0*V2*X3**2*V4+
00793      &       8.D0*V2*X2**2*V6+24.D0*U2*Y2**2*U4-8.D0*V5*X3**2*V3-
00794      &       8.D0*V2*X2**2*V3+8.D0*V4*X2**2*V3-36.D0*V2*X2*U2*Y2-
00795      &       24.D0*U3*Y3**2*U6-6.D0*V1**2*X2**2-54.D0*U1**2*Y3*Y2-
00796      &       8.D0*U1*Y2**2*U6-16.D0*V6**2*X2**2-2.D0*V2*X3**2*V3-
00797      &       38.D0*U2**2*Y2*Y3-16.D0*V4*X2**2*V1-8.D0*V5*X2**2*V1-
00798      &       24.D0*V5*X2*U2*Y2-24.D0*U5*Y2*V2*X2+16.D0*V2*X3**2*V5+
00799      &       56.D0*V4*X3**2*V1-8.D0*V1*X2**2*V6+8.D0*U2*Y2*V3*X2+
00800      &       48.D0*V5**2*X2**2-14.D0*V1*X3**2*V3+2.D0*V2*X3*U3*Y3+
00801      &       2.D0*V3**2*X2**2+8.D0*U4*Y2**2*U3+2.D0*U2*Y3*V3*X3-
00802      &       16.D0*U2*Y3**2*U6-16.D0*U2*Y3**2*U1-14.D0*U1*Y3**2*U3+
00803      &       64.D0*U4*Y3**2*U6+12.D0*V1*X2*U1*Y2-64.D0*V1*X3*U6*Y3+
00804      &       2.D0*U1*Y2**2*U3-24.D0*V6*X3**2*V3-6.D0*U3**2*Y3*Y2+
00805      &       24.D0*U5*Y2**2*U2-54.D0*V1**2*X3*X2-16.D0*V2*X3**2*V1+
00806      &       8.D0*V2*X2*U3*Y2-16.D0*U4*Y2**2*U1+2.D0*V1*X2**2*V3-
00807      &       8.D0*U2*Y2**2*U3+16.D0*U2*Y3**2*U5-16.D0*U4*Y3**2*U3-
00808      &       8.D0*U2*Y3**2*U4-96.D0*U5**2*Y2*Y3-96.D0*V5**2*X3*X2+
00809      &       32.D0*V5*X3**2*V4+32.D0*U5*Y3**2*U4+2.D0*U3**2*Y2**2-
00810      &       16.D0*V2*X3*U5*Y3+64.D0*V6*X2*U6*Y3+32.D0*V6*X2*U6*Y2-
00811      &       8.D0*V6*X2*U3*Y2+20.D0*V2*X3*U5*Y2+8.D0*V3*X2**2*V6
00812       ANS3 = 8.D0*U2*Y3*V4*X2+8.D0*U2*Y3*V4*X3-16.D0*U2*Y3*U4*Y2-
00813      &       16.D0*V2*X3*V4*X2+32.D0*V5*X3**2*V6-4.D0*V3*X2*U3*Y2-
00814      &       8.D0*V3*X2*U6*Y2+8.D0*V2*X3*U4*Y2-40.D0*V2*X3*V5*X2-
00815      &       6.D0*U1**2*Y2**2-2.D0*U1*Y2*V3*X2+8.D0*V2*X3*U4*Y3+
00816      &       64.D0*V6*X3*U6*Y2+8.D0*U1*Y2*V6*X2-8.D0*V4*X2*U3*Y2+
00817      &       38.D0*V2*X2*U2*Y3+14.D0*V1*X3*U3*Y3+38.D0*U2*Y2*V2*X3-
00818      &       40.D0*U2*Y3*V2*X3+20.D0*V2*X2*U5*Y3-4.D0*V5*X3*U3*Y2+
00819      &       64.D0*V4*X3**2*V6-4.D0*U5*Y3*V1*X2-2.D0*U2*Y3**2*U3-
00820      &       56.D0*U1*Y3*U6*Y2-32.D0*U4*Y2**2*U6+8.D0*U5*Y3*V3*X3+
00821      &       28.D0*V1*X3*U6*Y2-48.D0*U4**2*Y2**2-32.D0*U5*Y2*U6*Y3+
00822      &       12.D0*U1*Y3*U3*Y2+8.D0*U5*Y2*U3*Y3+64.D0*V1*X3**2*V6-
00823      &       8.D0*U4*Y2*V3*X2+16.D0*V5*X3*U6*Y2-6.D0*U1*Y3*V3*X2+
00824      &       8.D0*V5*X3*V3*X2+8.D0*V5*X3*V1*X2+64.D0*U1*Y3**2*U6+
00825      &       8.D0*V5*X2*U1*Y2+8.D0*U5*Y2*U1*Y3+6.D0*V3*X3*U3*Y2+
00826      &       6.D0*V3*X2*U3*Y3-8.D0*V3*X2*U6*Y3+16.D0*U3*Y3*U6*Y2-
00827      &       32.D0*V5*X3*V6*X2-8.D0*U3*Y3*V6*X2+54.D0*V1*X2*U1*Y3+
00828      &       28.D0*V1*X2*U6*Y3-4.D0*V5*X3*U1*Y2+8.D0*V5*X3*U3*Y3+
00829      &       12.D0*V1*X3*V3*X2+56.D0*U4*Y3**2*U1+16.D0*V6*X3*V3*X2-
00830      &       8.D0*V3*X3*U3*Y3-160.D0*V6*X3*U6*Y3-4.D0*U5*Y2*V1*X3-
00831      &       4.D0*U5*Y2*V3*X3-32.D0*V4*X2**2*V6+24.D0*V3*X3*U6*Y3
00832       ANS4 = -8.D0*V3*X3*U6*Y2+16.D0*U5*Y2*V6*X3+16.D0*U2*Y3*U1*Y2+
00833      &       16.D0*U4*Y3*V6*X2-5.D0*U2*Y3*V3*X2-4.D0*U4*Y3*V3*X2+
00834      &       20.D0*U4*Y3*V1*X2-56.D0*U4*Y3*V1*X3-6.D0*V1*X3*U3*Y2-
00835      &       6.D0*U1*Y2*V3*X3-8.D0*U2*Y3*V1*X2+28.D0*U1*Y2*V6*X3-
00836      &       8.D0*V6*X3*U3*Y2+24.D0*V6*X3*U3*Y3-4.D0*V2*X3*U6*Y2-
00837      &       40.D0*V4*X3*V1*X2-4.D0*V4*X2*U3*Y3-5.D0*V2*X3*U3*Y2+
00838      &       20.D0*V4*X2*U1*Y3+8.D0*V2*X3*V6*X2+16.D0*U2*Y3*V1*X3+
00839      &       8.D0*U2*Y3*U6*Y2+16.D0*V4*X2*U6*Y3+16.D0*U4*Y3*V3*X3+
00840      &       10.D0*U2*Y3*U3*Y2-4.D0*U2*Y3*V6*X2-8.D0*V2*X3*U1*Y2+
00841      &       8.D0*V1*X2*U6*Y2+16.D0*V4*X3*U6*Y2-4.D0*V4*X3*U3*Y2+
00842      &       10.D0*V2*X3*V3*X2-32.D0*V4*X3*V6*X2+20.D0*V4*X3*U1*Y2+
00843      &       16.D0*V4*X3*U3*Y3+16.D0*V2*X3*V1*X2-16.D0*V2*X3**2*V6+
00844      &       8.D0*U5*Y2*V1*X2+16.D0*V2*X3*U1*Y3+16.D0*V2*X3*U6*Y3-
00845      &       4.D0*V2*X2*U6*Y3+16.D0*U4*Y2*V6*X3-5.D0*V2*X2*U3*Y3+
00846      &       8.D0*U3*Y2**2*U6+16.D0*V5*X2*U6*Y3-6.D0*V1*X2*U3*Y3+
00847      &       8.D0*V4*X3*V3*X2-16.D0*V4*X3**2*V3-56.D0*V4*X3*U1*Y3-
00848      &       32.D0*U4*Y2*U6*Y3-8.D0*V2*X2*U1*Y3+8.D0*U4*Y2*U3*Y3+
00849      &       20.D0*U4*Y2*V1*X3-4.D0*U4*Y2*V3*X3-5.D0*U2*Y2*V3*X3-
00850      &       4.D0*V5*X2*U3*Y3-32.D0*U5*Y3*V6*X3-32.D0*V5*X3*U6*Y3
00851       A14(IELEM) = (-4.D0*V5*X2*U1*Y3-4.D0*U2*Y2*V6*X3-8.D0*U2*Y2*V6*X2-
00852      &       40.D0*U4*Y2*U1*Y3+16.D0*U5*Y3*V6*X2-8.D0*V2*X2*U6*Y2-
00853      &       4.D0*U5*Y3*V3*X2-8.D0*U2*Y2*V1*X3+54.D0*V1*X3*U1*Y2+
00854      &       ANS1+ANS2+ANS3+ANS4) * (-AUX630)
00855 !
00856       ANS1 = (-96.D0*U5*V5+12.D0*V1*U1+8.D0*U2*V3+12.D0*U4*V1-
00857      &       16.D0*U6*V3-8.D0*V2*U6-8.D0*V4*U3-V1*U3-U1*V3-20.D0*U2*V2+
00858      &       8.D0*U5*V1-8.D0*U2*V6+8.D0*V5*U1-8.D0*U4*V3+64.D0*U6*V6+
00859      &       32.D0*V4*U6+32.D0*U4*V6-12.D0*U5*V2-12.D0*U5*V3+
00860      &       8.D0*V2*U3-16.D0*U2*V4+12.D0*U1*V6-V2*U1-16.D0*U3*V6-
00861      &       U2*V1-16.D0*V2*U4+12.D0*V4*U1-20.D0*U3*V3-12.D0*U2*V5-
00862      &       12.D0*U3*V5+64.D0*U4*V4+12.D0*V1*U6)*Y3*X2+((24.D0*V3*V5+
00863      &       24.D0*V5*V2+20.D0*V3**2+32.D0*V3*V6-24.D0*V1*V6-
00864      &       16.D0*V2*V3+2.D0*V2*V1+16.D0*V4*V3-64.D0*V4**2+
00865      &       96.D0*V5**2-64.D0*V6**2+20.D0*V2**2-16.D0*V5*V1-
00866      &       64.D0*V4*V6+32.D0*V2*V4-12.D0*V1**2-24.D0*V4*V1+
00867      &       2.D0*V1*V3+16.D0*V2*V6)*X3+(96.D0*U5*V5-12.D0*V1*U1-
00868      &       8.D0*U2*V3-16.D0*U4*V1+8.D0*U6*V3+8.D0*V2*U6+8.D0*V4*U3+
00869      &       2.D0*V1*U3+2.D0*U1*V3+36.D0*U2*V2-8.D0*U5*V1+8.D0*U2*V6-
00870      &       8.D0*V5*U1+8.D0*U4*V3-32.D0*U6*V6-32.D0*V4*U6-32.D0*U4*V6+
00871      &       24.D0*U5*V2-8.D0*V2*U3+24.D0*U2*V4-8.D0*U1*V6+8.D0*U3*V6+
00872      &       24.D0*V2*U4-16.D0*V4*U1+4.D0*U3*V3+24.D0*U2*V5-
00873      &       96.D0*U4*V4-8.D0*V1*U6)*Y2)*X2+(-96.D0*U5*V5+12.D0*V1*U1+
00874      &       8.D0*U2*V3+12.D0*U4*V1-16.D0*U6*V3-8.D0*V2*U6)*Y2*X3
00875       ANS2 = (-8.D0*V4*U3-V1*U3-U1*V3-20.D0*U2*V2+8.D0*U5*V1-8.D0*U2*V6+
00876      &       8.D0*V5*U1-8.D0*U4*V3+64.D0*U6*V6+32.D0*V4*U6+32.D0*U4*V6-
00877      &       12.D0*U5*V2-12.D0*U5*V3+8.D0*V2*U3-16.D0*U2*V4+12.D0*U1*V6-
00878      &       V2*U1-16.D0*U3*V6-U2*V1-16.D0*V2*U4+12.D0*V4*U1-
00879      &       20.D0*U3*V3-12.D0*U2*V5-12.D0*U3*V5+64.D0*U4*V4+
00880      &       12.D0*V1*U6)*Y2*X3+(-8.D0*U6*U3-24.D0*U2*U4+48.D0*U4**2-
00881      &       2.D0*U3*U1+32.D0*U4*U6+8.D0*U6*U1-48.D0*U5**2+8.D0*U5*U1-
00882      &       18.D0*U2**2-2.D0*U3**2+6.D0*U1**2+16.D0*U6**2+16.D0*U4*U1-
00883      &       8.D0*U4*U3+8.D0*U2*U3-8.D0*U2*U6-24.D0*U5*U2)*Y2**2+
00884      &       (-64.D0*U4**2-24.D0*U6*U1+24.D0*U5*U2+24.D0*U3*U5-
00885      &       64.D0*U4*U6+16.D0*U2*U6-16.D0*U2*U3-16.D0*U5*U1-
00886      &       12.D0*U1**2+20.D0*U2**2+20.D0*U3**2+32.D0*U6*U3+
00887      &       96.D0*U5**2+2.D0*U2*U1+32.D0*U2*U4+2.D0*U3*U1-
00888      &       24.D0*U4*U1-64.D0*U6**2+16.D0*U4*U3)*Y3*Y2+(96.D0*U5*V5-
00889      &       12.D0*V1*U1-8.D0*U2*V3-8.D0*U4*V1+24.D0*U6*V3+
00890      &       8.D0*V2 *U6+8.D0*V4*U3+4.D0*U2*V2-8.D0*U5*V1+8.D0*U2*V6-
00891      &       8.D0*V5*U1+8.D0*U4*V3-96.D0*U6*V6-32.D0*V4*U6-32.D0*U4*V6+
00892      &       24.D0*U5*V3-8.D0*V2*U3+8.D0*U2*V4-16.D0*U1*V6+2.D0*V2*U1+
00893      &       24.D0*U3*V6+2.D0*U2*V1+8.D0*V2*U4-8.D0*V4*U1+36.D0*U3*V3+
00894      &       24.D0*U3*V5-32.D0*U4*V4-16.D0*V1*U6)*Y3*X3
00895       A15(IELEM) = ((16.D0*V1*V6-24.D0*V3*V6-18.D0*V3**2+48.D0*V6**2-
00896      &             8.D0*V4*V3+ 6.D0*V1**2+8.D0*V5*V1+8.D0*V2*V3+
00897      &             16.D0*V4**2-8.D0*V2*V4-48.D0*V5**2+32.D0*V4*V6+
00898      &             8.D0*V4*V1-2.D0*V2**2-2.D0*V2*V1-8.D0*V2*V6-
00899      &             24.D0*V3*V5)*X3**2+(-24.D0*V2*V4+6.D0*V1**2+
00900      &             32.D0*V4*V6+16.D0*V4*V1-8.D0*V4*V3+8.D0*V5*V1+
00901      &             48.D0*V4**2+8.D0*V2*V3+16.D0*V6**2-8.D0*V2*V6+
00902      &             8.D0*V1*V6-18.D0*V2**2-48.D0*V5**2-2.D0*V3**2-
00903      &             24.D0*V5*V2-2.D0*V1*V3-8.D0*V3*V6)*X2**2+
00904      &             (-48.D0*U5**2-24.D0*U3*U5-8.D0*U2*U4-2.D0*U2*U1+
00905      &             6.D0*U1**2-2.D0*U2**2-18.D0*U3**2+16.D0*U4**2+
00906      &             32.D0*U4*U6+8.D0*U2*U3+8.D0*U5*U1+48.D0*U6**2-
00907      &             24.D0*U6*U3-8.D0*U2*U6-8.D0*U4*U3+8.D0*U4*U1+
00908      &             16.D0*U6*U1)*Y3**2+ANS1+ANS2) * (-AUX630)
00909 !
00910       ANS1 = -186.D0*U2*Y3*V2*X3-23.D0*U2*Y3**2*U1+48.D0*V4*X3**2*V6-
00911      &       96.D0*U4*Y3*V2*X3-23.D0*V2*X3**2*V3+112.D0*V4*X3**2*V5-
00912      &       32.D0*V4*X3**2*V1+4.D0*V2*X3**2*V6+96.D0*V4*X3**2*V2-
00913      &       28.D0*V5*X3**2*V1+48.D0*U4*Y3**2*U6+4.D0*U2*Y3**2*U6+
00914      &       96.D0*U2*Y3**2*U4+48.D0*U5*Y3**2*U6-32.D0*V3*X3**2*V5-
00915      &       23.D0*U2*Y3**2*U3-28.D0*U5*Y3**2*U1-28.D0*U4*Y3**2*U3-
00916      &       23.D0*V2*X3**2*V1+112.D0*U4**2*Y3**2+V3*X3**2*V1-
00917      &       32.D0*U4*Y3**2*U1+112.D0*U5**2*Y3**2+13.D0*U3**2*Y3**2+
00918      &       112.D0*V5**2*X3**2+13.D0*V3**2*X3**2+112.D0*U5*Y3**2*U4+
00919      &       48.D0*V6**2*X3**2-26.D0*V1*X3*U1*Y3+13.D0*U1**2*Y3**2-
00920      &       V3*X3*U1*Y3-26.D0*V3*X3*U3*Y3+U1*Y3**2*U3+
00921      &       48.D0*V5*X3**2*V6+112.D0*V4**2*X3**2-224.D0*V4*X3*U4*Y3-
00922      &       224.D0*V5*X3*U5*Y3-112.D0*V5*X3*U4*Y3-112.D0*V4*X3*U5*Y3-
00923      &       4.D0*V2*X3*U6*Y3+23.D0*U2*Y3*V3*X3+96.D0*V2*X3**2*V5+
00924      &       23.D0*V2*X3*U1*Y3+23.D0*V2*X3*U3*Y3+23.D0*U2*Y3*V1*X3-
00925      &       4.D0*U2*Y3*V6*X3+28.D0*U4*Y3*V3*X3+32.D0*U4*Y3*V1*X3-
00926      &       48.D0*U4*Y3*V6*X3+28.D0*U5*Y3*V1*X3-48.D0*V4*X3*U6*Y3+
00927      &       32.D0*V4*X3*U1*Y3+28.D0*V4*X3*U3*Y3-32.D0*U3*Y3**2*U5
00928       A22(IELEM) = (ANS1-48.D0*U5*Y3*V6*X3-48.D0*V5*X3*U6*Y3+
00929      &             28.D0*V5*X3*U1*Y3+96.D0*U2*Y3**2*U5+
00930      &             93.D0*U2**2*Y3**2-96.D0*U2*Y3*V5*X3-
00931      &             96.D0*V2*X3*U5*Y3+32.D0*V3*X3*U5*Y3+
00932      &             32.D0*U3*Y3*V5*X3+93.D0*V2**2*X3**2-
00933      &             96.D0*U2*Y3*V4*X3+48.D0*U6**2*Y3**2+
00934      &             13.D0*V1**2*X3**2-28.D0*V4*X3**2*V3-
00935      &             V1*X3*U3*Y3-96.D0*V6*X3*U6*Y3) * (AUX1260)
00936 !
00937       ANS1 = (-32.D0*V5**2-18.D0*V2*V3-8.D0*V4*V3-8.D0*V2*V6+
00938      &       32.D0*V5*V2-10.D0*V1**2+32.D0*V3*V5+96.D0*V4**2+
00939      &       64.D0*V3*V6+96.D0*V6**2+54.D0*V3**2+32.D0*V5*V6+
00940      &       32.D0*V4*V6+32.D0*V4*V5-32.D0*V4*V1-32.D0*V1*V6-
00941      &       8.D0*V5*V1-10.D0*V2*V1+64.D0*V2*V4-10.D0*V1*V3+
00942      &       54.D0*V2**2)*X3+(32.D0*U5*V5+10.D0*V1*U1+9.D0*U2*V3+
00943      &       16.D0*U4*V1-32.D0*U6*V3+4.D0*V2*U6+4.D0*V4*U3+5.D0*V1*U3+
00944      &       5.D0*U1*V3-54.D0*U2*V2-16.D0*V4*U5+4.D0*U5*V1-16.D0*U4*V5+
00945      &       4.D0*U2*V6+4.D0*V5*U1+4.D0*U4*V3-96.D0*U6*V6-16.D0*V4*U6-
00946      &       16.D0*U4*V6-16.D0*U5*V6-16.D0*V5*U6-16.D0*U5*V2-
00947      &       16.D0*U5*V3+9.D0*V2*U3-32.D0*U2*V4+16.D0*U1*V6+5.D0*V2*U1-
00948      &       32.D0*U3*V6+5.D0*U2*V1-32.D0*V2*U4+16.D0*V4*U1-
00949      &       54.D0*U3*V3-16.D0*U2*V5-16.D0*U3*V5-96.D0*U4*V4+
00950      &       16.D0*V1*U6)*Y3
00951       A23(IELEM) = (ANS1*X2+(32.D0*U5*V5+10.D0*V1*U1+9.D0*U2*V3+
00952      &             16.D0*U4*V1-32.D0*U6*V3+ 4.D0*V2*U6+4.D0*V4*U3+
00953      &             5.D0*V1*U3+5.D0*U1*V3-54.D0*U2*V2-16.D0*V4*U5+
00954      &             4.D0*U5*V1-16.D0*U4*V5+4.D0*U2*V6+4.D0*V5*U1+
00955      &             4.D0*U4*V3-96.D0*U6*V6-16.D0*V4*U6-16.D0*U4*V6-
00956      &             16.D0*U5*V6-16.D0*V5*U6-16.D0*U5*V2-16.D0*U5*V3+
00957      &             9.D0*V2*U3-32.D0*U2*V4+16.D0*U1*V6+5.D0*V2*U1-
00958      &             32.D0*U3*V6+5.D0*U2*V1-32.D0*V2*U4+16.D0*V4*U1-
00959      &             54.D0*U3*V3-16.D0*U2*V5-16.D0*U3*V5-96.D0*U4*V4+
00960      &             16.D0*V1*U6)*Y2*X3+(96.D0*U4**2-10.D0*U1**2+
00961      &             96.D0*U6**2+54.D0*U2**2+54.D0*U3**2-32.D0*U5**2+
00962      &             64.D0*U6*U3+32.D0*U5*U6-8.D0*U5*U1-10.D0*U2*U1+
00963      &             64.D0*U2*U4-10.D0*U3*U1+32.D0*U4*U5-32.D0*U6*U1-
00964      &             8.D0*U2*U6-18.D0*U2*U3-8.D0*U4*U3+32.D0*U4*U6+
00965      &             32.D0*U3*U5-32.D0*U4*U1+32.D0*U5*U2)*Y3*Y2) *
00966      &             AUX2520
00967 !
00968       ANS1 = ((8.D0*V2*V6+96.D0*V4**2+66.D0*V2**2+32.D0*V4*V6-
00969      &       16.D0*V2*V1-32.D0*V4*V1+32.D0*V5* V6-24.D0*V5*V1-
00970      &       8.D0*V3*V6+72.D0*V5*V2+ 2.D0*V1**2-32.D0*V3*V5+
00971      &       96.D0*V4*V5+96.D0*V5**2-16.D0*V2*V3+2.D0*V3**2-
00972      &       8.D0*V1*V6+72.D0*V2*V4+ 6.D0*V1*V3-24.D0*V4*V3)*X3+
00973      &       (-96.D0*U5*V5-2.D0* V1*U1+8.D0*U2*V3+16.D0*U4*V1+
00974      &       4.D0*U6*V3-4.D0*V2*U6+12.D0*V4*U3-3.D0*V1*U3-
00975      &       3.D0*U1*V3-66.D0*U2*V2-48.D0*V4*U5+12.D0*U5*V1-
00976      &       48.D0*U4*V5-4.D0*U2*V6+12.D0*V5*U1+ 12.D0*U4*V3-
00977      &       16.D0*V4*U6-16.D0*U4*V6-16.D0*U5*V6-16.D0*V5*U6-
00978      &       36.D0*U5*V2+16.D0*U5*V3+8.D0*V2*U3-36.D0*U2*V4+
00979      &       4.D0*U1*V6+8.D0*V2*U1+4.D0*U3*V6+8.D0*U2*V1-
00980      &       36.D0*V2*U4+16.D0*V4*U1-2.D0*U3*V3-36.D0*U2*V5+
00981      &       16.D0*U3*V5-96.D0*U4*V4+4.D0*V1*U6)*Y3)*X2+
00982      &       (-96.D0*U5*V5-2.D0*V1*U1+8.D0*U2*V3+16.D0*U4*V1+
00983      &       4.D0*U6*V3-4.D0*V2*U6+12.D0*V4*U3-3.D0*V1*U3-
00984      &       3.D0*U1*V3-66.D0*U2*V2-48.D0*V4*U5+12.D0*U5*V1-
00985      &       48.D0*U4*V5-4.D0*U2*V6+12.D0*V5*U1+12.D0*U4*V3-
00986      &       16.D0*V4*U6-16.D0*U4*V6-16.D0*U5*V6-16.D0*V5*U6-
00987      &       36.D0*U5*V2+16.D0*U5*V3+8.D0*V2*U3-36.D0*U2*V4)*Y2*X3
00988       ANS2 = (4.D0*U1*V6+8.D0*V2*U1+4.D0*U3*V6+8.D0*U2*V1-
00989      &        36.D0*V2*U4+16.D0*V4*U1-2.D0*U3*V3-36.D0*U2*V5+
00990      &        16.D0*U3*V5-96.D0*U4*V4+4.D0*V1*U6)*Y2*X3+
00991      &        (-8.D0*U6*U3+32.D0*U5*U6-32.D0*U3*U5-32.D0*U4*U1-
00992      &        24.D0*U4*U3-24.D0*U5*U1+96.D0*U4**2+66.D0*U2**2+
00993      &        2.D0*U3**2+8.D0*U2*U6+96.D0*U5**2-16.D0*U2*U1+
00994      &        72.D0*U2*U4-16.D0*U2*U3+32.D0*U4*U6+6.D0*U3*U1-
00995      &        8.D0*U6*U1+96.D0*U4*U5+2.D0*U1**2+72.D0*U5*U2)*Y3*Y2+
00996      &        (-64.D0*U4*U5+24.D0*U3*U5-32.D0*U5*U6-16.D0*U6*U1-
00997      &        56.D0*U2*U4-60.D0*U2**2+16.D0*U2*U1+16.D0*U5*U1+
00998      &        8.D0*U4*U1+8.D0*U6*U3-32.D0*U4*U6+16.D0*U4*U3-
00999      &        64.D0*U5*U2+14.D0*U2*U3+2.D0*U3*U1-20.D0*U1**2-
01000      &        4.D0*U3**2-80.D0*U5**2-48.D0*U6**2-48.D0*U4**2)*Y3**2
01001       A24(IELEM) = (ANS1+ANS2+(-56.D0*V2*V4+14.D0*V2*V3+16.D0*V5*V1+
01002      &             16.D0*V2*V1-64.D0*V4*V5-32.D0*V5*V6-16.D0*V1*V6-
01003      &             80.D0*V5**2+16.D0*V4*V3-60.D0*V2**2-4.D0*V3**2-
01004      &             48.D0*V6**2-20.D0*V1**2+2.D0*V1*V3-48.D0*V4**2-
01005      &             32.D0*V4*V6+8.D0*V4*V1+8.D0*V3*V6-64.D0*V5*V2+
01006      &             24.D0*V3*V5)*X3**2+(160.D0*U5*V5+40.D0*V1*U1-
01007      &             14.D0*U2*V3-8.D0*U4*V1-8.D0*U6*V3-16.D0*V4*U3-
01008      &             2.D0*V1*U3-2.D0*U1*V3+120.D0*U2*V2+64.D0*V4*U5-
01009      &             16.D0*U5*V1+64.D0*U4*V5-16.D0*V5*U1-16.D0*U4*V3+
01010      &             96.D0*U6*V6+32.D0*V4*U6+32.D0*U4*V6+32.D0*U5*V6+
01011      &             32.D0*V5*U6+64.D0*U5*V2-24.D0*U5*V3-14.D0*V2*U3+
01012      &             56.D0*U2*V4+16.D0*U1*V6-16.D0*V2*U1-8.D0*U3*V6-
01013      &             16.D0*U2*V1+ 56.D0*V2*U4-8.D0*V4*U1+8.D0*U3*V3+
01014      &             64.D0*U2*V5-24.D0*U3*V5+96.D0*U4*V4+
01015      &             16.D0*V1*U6)*Y3*X3)*AUX630
01016 !
01017       ANS1 = ((-8.D0*V5*V2+16.D0*V3*V5+24.D0*V3*V6+16.D0*V3**2-
01018      &       2.D0*V2*V3+2.D0*V2*V1+8.D0*V2*V4-24.D0*V1*V6-32.D0*V5**2-
01019      &       16.D0*V1**2-16.D0*V4*V1+32.D0*V4**2)*X3+(12.D0*V1*U6+
01020      &       4.D0*U5*V2-32.D0*U4*V4-8.D0*U3*V5+4.D0*U2*V5-8.D0*U5*V3+
01021      &       V2*U3-4.D0*U2*V4+12.D0*U1*V6-V2*U1-12.D0*U3*V6-U2*V1-
01022      &       4.D0*V2*U4+8.D0*V4*U1-16.D0*U3*V3+32.D0*U5*V5+16.D0*V1*U1+
01023      &       U2*V3+8.D0*U4*V1-12.D0*U6*V3)*Y3)*X2+((12.D0*V1*U6+
01024      &       4.D0*U5*V2-32.D0*U4*V4-8.D0*U3*V5+4.D0*U2*V5-8.D0*U5*V3+
01025      &       V2*U3-4.D0*U2*V4+12.D0*U1*V6-V2*U1-12.D0*U3*V6-U2*V1-
01026      &       4.D0*V2*U4+8.D0*V4*U1-16.D0*U3*V3+32.D0*U5*V5+
01027      &       16.D0*V1*U1+U2*V3+8.D0*U4*V1-12.D0*U6*V3)*Y2+(-96.D0*U5*V5+
01028      &       4.D0*V1*U1+8.D0*U4*V1+24.D0*U6*V3-8.D0*V2*U6+8.D0*V4*U3-
01029      &       8.D0*V1*U3-8.D0*U1*V3-12.D0*U2*V2-32.D0*V4*U5+8.D0*U5*V1-
01030      &       32.D0*U4*V5-8.D0*U2*V6+8.D0*V5*U1+8.D0*U4*V3+96.D0*U6*V6-
01031      &       16.D0*U5*V2+24.D0*U5*V3-8.D0*U2*V4+2.D0*V2*U1+24.D0*U3*V6+
01032      &       2.D0*U2*V1-8.D0*V2*U4+8.D0*V4*U1+36.D0*U3*V3-16.D0*U2*V5+
01033      &       24.D0*U3*V5-32.D0*U4*V4)*Y3)*X3
01034       A26(IELEM)=(ANS1+(32.D0*V4*V5+8.D0*V2*V4-18.D0*V3**2+8.D0*V1*V3-
01035      &           8.D0*V5*V1+48.D0*V5**2-24.D0*V3*V6+6.D0*V2**2-
01036      &           8.D0*V4*V1-48.D0*V6**2-2.D0*V1**2-8.D0*V4*V3+
01037      &           16.D0*V4**2-2.D0*V2*V1+8.D0*V2*V6+16.D0*V5*V2-
01038      &           24.D0*V3*V5)*X3**2+(-16.D0*U4*U1-24.D0*U6*U1
01039      &           -8.D0*U5*U2+16.D0*U3*U5-16.D0*U1**2+8.D0*U2*U4+
01040      &           16.D0*U3**2-32.D0*U5**2+2.D0*U2*U1+32.D0*U4**2-
01041      &           2.D0*U2*U3+24.D0*U6*U3)*Y3*Y2+(-8.D0*U5*U1-24.D0*U3*U5-
01042      &           24.D0*U6*U3+6.D0*U2**2+8.D0*U2*U6+32.D0*U4*U5+
01043      &           16.D0*U5*U2-2.D0*U2*U1-8.D0*U4*U3+8.D0*U2*U4+
01044      &           8.D0*U3*U1+16.D0*U4**2+48.D0*U5**2-2.D0*U1**2-
01045      &           8.D0*U4*U1-18.D0*U3**2-48.D0*U6**2)*Y3**2)*(-AUX630)
01046 !
01047       A33(IELEM) = ((-224.D0*U5*V5-26.D0*V1*U1+23.D0*U2*V3-96.D0*U6*V3+
01048      &             28.D0*V2*U6-4.D0*V4*U3+23.D0*V1*U3+23.D0*U1*V3-
01049      &             26.D0*U2*V2-48.D0*V4*U5+28.D0*U5*V1-48.D0*U4*V5+
01050      &             28.D0*U2*V6+28.D0*V5*U1-4.D0*U4*V3-224.D0*U6*V6-
01051      &             48.D0*V4*U6-48.D0*U4*V6-112.D0*U5*V6-112.D0*V5*U6+
01052      &             32.D0*U5*V2-96.D0*U5*V3+23.D0*V2*U3+32.D0*U1*V6-
01053      &             V2*U1-96.D0*U3*V6-U2*V1-186.D0*U3*V3+32.D0*U2*V5-
01054      &             96.D0*U3*V5-96.D0*U4*V4+32.D0*V1*U6)*Y2*X2+
01055      &             (4.D0*U4*U3+U2*U1+112.D0*U6**2+48.D0*U4**2+
01056      &             93.D0*U3**2-32.D0*U6*U1-23.D0*U3*U1+13.D0*U2**2-
01057      &             28.D0*U5*U1-28.D0*U2*U6+112.D0*U5**2-32.D0*U5*U2+
01058      &             13.D0*U1**2+48.D0*U4*U6+96.D0*U3*U5-23.D0*U2*U3+
01059      &             112.D0*U5*U6+48.D0*U4*U5+96.D0*U6*U3)*Y2**2+
01060      &             (-32.D0*V1*V6+48.D0*V4**2+96.D0*V3*V6+112.D0*V5*V6+
01061      &             48.D0*V4*V5+V2*V1+93.D0*V3**2-23.D0*V2*V3-
01062      &             23.D0*V1*V3-28.D0*V2*V6+112.D0*V5**2+4.D0*V4*V3+
01063      &             48.D0*V4*V6-28.D0*V5*V1+13.D0*V2**2+112.D0*V6**2+
01064      &             13.D0*V1**2+96.D0*V3*V5-32.D0*V5*V2)*X2**2)
01065      &             *AUX1260
01066 !
01067       ANS1 = (-8.D0*V3*V6+32.D0*V5**2-24.D0*V2*V4-16.D0*V2**2+
01068      &       2.D0*V2*V3-32.D0*V6**2+16.D0*V1**2+16.D0*V1*V6+
01069      &       24.D0*V4*V1-16.D0*V5*V2-2.D0*V1*V3+8.D0*V3*V5)*X3*X2+
01070      &       (96.D0*U5*V5-4.D0*V1*U1+8.D0*U6*V3-8.D0*V2*U6+
01071      &       8.D0*V4*U3-2.D0*V1*U3-2.D0*U1*V3-36.D0*U2*V2-8.D0*U5*V1-
01072      &       8.D0*U2*V6-8.D0*V5*U1+8.D0*U4*V3+32.D0*U6*V6+32.D0*U5*V6+
01073      &       32.D0*V5*U6-24.D0*U5*V2+16.D0*U5*V3-24.D0*U2*V4-
01074      &       8.D0*U1*V6+8.D0*V2*U1+8.D0*U3*V6+8.D0*U2*V1-24.D0*V2*U4+
01075      &       12.D0*U3*V3-24.D0*U2*V5+16.D0*U3*V5-96.D0*U4*V4-
01076      &       8.D0*V1*U6)*Y2*X2+(16.D0*U2*V2+U1*V3-8.D0*V1*U6+
01077      &       12.D0*V2*U4+8.D0*U5*V2-4.D0*U3*V5+8.D0*U2*V5-4.D0*U5*V3-
01078      &       8.D0*U1*V6+12.D0*U2*V4-U2*V3-V2*U3-16.D0*V1*U1+4.D0*U6*V3+
01079      &       4.D0*U3*V6-12.D0*U4*V1+32.D0*U6*V6+V1*U3-12.D0*V4*U1-
01080      &       32.D0*U5*V5)*Y2*X3+(2.D0*U3*U1-8.D0*U4*U3-8.D0*U2*U1-
01081      &       32.D0*U5*U6+8.D0*U5*U1+24.D0*U2*U4+8.D0*U6*U1-8.D0*U6*U3+
01082      &       48.D0*U4**2-48.D0*U5**2+18.D0*U2**2+2.D0*U1**2-
01083      &       16.D0*U6**2-6.D0*U3**2+8.D0*U2*U6-16.D0*U3*U5+
01084      &       24.D0*U5*U2)*Y2**2
01085       A34(IELEM) = (ANS1+(8.D0*V1*V6-8.D0*V4*V3+18.D0*V2**2-16.D0*V6**2+
01086      &             2.D0*V1**2-6.D0*V3**2+2.D0*V1*V3+8.D0*V2*V6+
01087      &             8.D0*V5*V1-32.D0*V5*V6+24.D0*V2*V4-48.D0*V5**2-
01088      &             8.D0*V3*V6-8.D0*V2*V1+48.D0*V4**2-16.D0*V3*V5+
01089      &             24.D0*V5*V2)*X2**2+(16.D0*U2*V2+U1*V3-8.D0*V1*U6+
01090      &             12.D0*V2*U4+8.D0*U5*V2-4.D0*U3*V5+8.D0*U2*V5-
01091      &             4.D0*U5*V3-8.D0*U1*V6+12.D0*U2*V4-U2*V3-V2*U3-
01092      &             16.D0*V1*U1+4.D0*U6*V3+4.D0*U3*V6-12.D0*U4*V1+
01093      &             32.D0*U6*V6+V1*U3-12.D0*V4*U1-32.D0*U5*V5)*Y3*X2+
01094      &             (-8.D0*U6*U3-16.D0*U5*U2+24.D0*U4*U1-2.D0*U3*U1+
01095      &             16.D0*U6*U1+32.D0*U5**2+8.D0*U3*U5+16.D0*U1**2-
01096      &             32.D0*U6**2-24.D0*U2*U4+2.D0*U2*U3-
01097      &             16.D0*U2**2)*Y3*Y2)*AUX630
01098 !
01099       ANS1 = ((96.D0*U5*V5-4.D0*V1*U1+8.D0*U6*V3-8.D0*V2*U6+8.D0*V4*U3-
01100      &        2.D0*V1*U3-2.D0*U1*V3-36.D0*U2*V2-8.D0*U5*V1-8.D0*U2*V6-
01101      &        8.D0*V5*U1+8.D0*U4*V3+32.D0*U6*V6+32.D0*U5*V6+
01102      &        32.D0*V5*U6-24.D0*U5*V2+16.D0*U5*V3-24.D0*U2*V4-
01103      &        8.D0*U1*V6+8.D0*V2*U1+8.D0*U3*V6+8.D0*U2*V1-24.D0*V2*U4+
01104      &        12.D0*U3*V3-24.D0*U2*V5+16.D0*U3*V5-96.D0*U4*V4-
01105      &        8.D0*V1*U6)*Y2+(-96.D0*U5*V5-2.D0*V1*U1+8.D0*U2*V3+
01106      &        4.D0*U4*V1-36.D0*U6*V3+12.D0*V2*U6-4.D0*V4*U3+8.D0*V1*U3+
01107      &        8.D0*U1*V3-2.D0*U2*V2-16.D0*V4*U5+12.D0*U5*V1-16.D0*U4*V5+
01108      &        12.D0*U2*V6+12.D0*V5*U1-4.D0*U4*V3-96.D0*U6*V6-
01109      &        16.D0*V4*U6-16.D0*U4*V6-48.D0*U5*V6-48.D0*V5*U6+
01110      &        16.D0*U5*V2-36.D0*U5*V3+8.D0*V2*U3+4.D0*U2*V4+16.D0*U1*V6-
01111      &        3.D0*V2*U1-36.D0*U3*V6-3.D0*U2*V1+4.D0*V2*U4+4.D0*V4*U1-
01112      &        66.D0*U3*V3+16.D0*U2*V5-36.D0*U3*V5+16.D0*V1*U6)*Y3)*X2+
01113      &        (-96.D0*U5*V5-2.D0*V1*U1+8.D0*U2*V3+4.D0*U4*V1-
01114      &        36.D0*U6*V3+12.D0*V2*U6-4.D0*V4*U3+8.D0*V1*U3+8.D0*U1*V3-
01115      &        2.D0*U2*V2-16.D0*V4*U5+12.D0*U5*V1-16.D0*U4*V5+
01116      &        12.D0*U2*V6+12.D0*V5*U1-4.D0*U4*V3-96.D0*U6*V6-
01117      &        16.D0*V4*U6-16.D0*U4*V6-48.D0*U5*V6-48.D0*V5*U6+
01118      &        16.D0*U5*V2-36.D0*U5*V3+8.D0*V2*U3+4.D0*U2*V4)*Y2*X3
01119       ANS2 = (16.D0*U1*V6-3.D0*V2*U1-36.D0*U3*V6-3.D0*U2*V1+
01120      &        4.D0*V2*U4+4.D0*V4*U1-66.D0*U3*V3+16.D0*U2*V5-
01121      &        36.D0*U3*V5+16.D0*V1*U6)*Y2*X3+(72.D0*V3*V6+
01122      &        2.D0*V1**2+66.D0*V3**2+72.D0*V3*V5-32.D0*V1*V6-
01123      &        16.D0*V2*V3-32.D0*V5*V2+6.D0*V2*V1+96.D0*V5**2+
01124      &        96.D0*V6**2+ 32.D0*V4*V5+2.D0*V2**2+
01125      &        32.D0*V4*V6-8.D0*V2*V4+96.D0*V5*V6-24.D0*V5*V1+
01126      &        8.D0*V4*V3-24.D0*V2*V6-16.D0*V1*V3-8.D0*V4*V1)*X3*X2+
01127      &        (8.D0*V1*V6-8.D0*V4*V3+18.D0*V2**2-16.D0*V6**2+
01128      &        2.D0*V1**2-6.D0*V3**2+2.D0*V1*V3+8.D0*V2*V6+
01129      &        8.D0*V5*V1-32.D0*V5*V6+24.D0*V2*V4-48.D0*V5**2-
01130      &        8.D0*V3*V6-8.D0*V2*V1+48.D0*V4**2-16.D0*V3*V5+
01131      &        24.D0*V5*V2)*X2**2+(2.D0*U3*U1-8.D0*U4*U3-8.D0*U2*U1-
01132      &        32.D0*U5*U6+8.D0*U5*U1+24.D0*U2*U4+8.D0*U6*U1-
01133      &        8.D0*U6*U3+48.D0*U4**2-48.D0*U5**2+18.D0*U2**2+
01134      &        2.D0*U1**2-16.D0*U6**2-6.D0*U3**2+8.D0*U2*U6-
01135      &        16.D0*U3*U5+24.D0*U5*U2)*Y2**2
01136       A35(IELEM) = ((72.D0*U6*U3+96.D0*U5*U6- 24.D0*U5*U1-
01137      &             16.D0*U3*U1-8.D0*U4*U1+8.D0*U4*U3+32.D0*U4*U6+
01138      &             72.D0*U3*U5-32.D0*U6*U1+32.D0*U4*U5-32.D0*U5*U2+
01139      &             66.D0*U3**2-24.D0*U2*U6-8.D0*U2*U4-16.D0*U2*U3+
01140      &             96.D0*U5**2+2.D0*U1**2+96.D0*U6**2+6.D0*U2*U1+
01141      &             2.D0*U2**2)*Y3*Y2+ANS1+ANS2)*(-AUX630)
01142 !
01143       ANS1 = ((4.D0*V1*V6+4.D0*V6**2+8.D0*V3*V5+4.D0*V2*V1+6.D0*V2*V3-
01144      &        12.D0*V4**2-36.D0*V5**2-8.D0*V5*V6-21.D0*V2**2+
01145      &        8.D0*V5*V1-V3**2-24.D0*V4*V5+V1**2+12.D0*V4*V1+
01146      &        4.D0*V4*V3-4.D0*V2*V6-24.D0*V2*V4-2.D0*V1*V3-
01147      &        24.D0*V5*V2)*X3+(-48.D0*U5*V5-2.D0*V1*U1+3.D0*U2*V3+
01148      &        6.D0*U4*V1+2.D0*U6*V3+6.D0*V4*U3-V1*U3-U1*V3-
01149      &        24.D0*U2*V2-24.D0*V4*U5+6.D0*U5*V1-24.D0*U4*V5+6.D0*V5*U1+
01150      &        6.D0*U4*V3-8.D0*U6*V6-12.D0*V4*U6-12.D0*U4*V6-12.D0*U5*V6-
01151      &        12.D0*V5*U6-12.D0*U5*V2+6.D0*U5*V3+3.D0*V2*U3-12.D0*U2*V4+
01152      &        2.D0*U1*V6+3.D0*V2*U1+2.D0*U3*V6+3.D0*U2*V1-12.D0*V2*U4+
01153      &        6.D0*V4*U1-2.D0*U3*V3-12.D0*U2*V5+6.D0*U3*V5-48.D0*U4*V4+
01154      &        2.D0*V1*U6)*Y2)*X2+(12.D0*V4*V5+6.D0*V4*V1-4.D0*V4*V3+
01155      &        10.D0*V1*V6+10.D0*V2**2-2.D0*V2*V3+16.D0*V6**2+6.D0*V2*V4+
01156      &        16.D0*V5**2+12.D0*V4*V6-4.D0*V3*V6+12.D0*V4**2+
01157      &        10.D0*V1**2-4.D0*V2*V1-2.D0*V2*V6-2.D0*V1*V3+8.D0*V5*V6-
01158      &        2.D0*V5*V1+V3**2+10.D0*V5*V2-4.D0*V3*V5)*X3**2+
01159      &        (24.D0*V4**2+V1**2-6.D0*V5*V1+12.D0*V2**2+24.D0*V5**2-
01160      &        2.D0*V1*V6+12.D0*V4*V6-3*V2*V1-3.D0*V2*V3+4.D0*V6**2+
01161      &        V1*V3+24.D0*V4*V5-6.D0*V4*V1+12.D0*V5*V6+V3**2-6.D0*V4*V3-
01162      &        6.D0*V3*V5+12.D0*V5*V2+12.D0*V2*V4-2.D0*V3*V6)*X2**2
01163       ANS2 = (-32.D0*U5*V5-20.D0*V1*U1+2.D0*U2*V3-6.D0*U4*V1+4.D0*U6*V3+
01164      &       2.D0*V2*U6+4.D0*V4*U3+2.D0*V1*U3+2.D0*U1*V3-20.D0*U2*V2-
01165      &       12.D0*V4*U5+2.D0*U5*V1-12.D0*U4*V5+2.D0*U2*V6+2.D0*V5*U1+
01166      &       4.D0*U4*V3-32.D0*U6*V6-12.D0*V4*U6-12.D0*U4*V6-8.D0*U5*V6-
01167      &       8.D0*V5*U6-10.D0*U5*V2+4.D0*U5*V3+2.D0*V2*U3-6.D0*U2*V4-
01168      &       10.D0*U1*V6+4.D0*V2*U1+4.D0*U3*V6+4.D0*U2*V1-6.D0*V2*U4-
01169      &       6.D0*V4*U1-2.D0*U3*V3-10.D0*U2*V5+4.D0*U3*V5-24.D0*U4*V4-
01170      &       10.D0*V1*U6)*Y3*X3+(U1**2+8.D0*U5*U1-8.D0*U5*U6+4.D0*U4*U3+
01171      &       6.D0*U2*U3+4.D0*U6*U1-24.D0*U4*U5-4*U2*U6-24.D0*U2*U4-
01172      &       12.D0*U4**2+8.D0*U3*U5+4.D0*U6**2+12.D0*U4*U1-2.D0*U3*U1-
01173      &       21.D0*U2**2-U3**2-36.D0*U5**2+4.D0*U2*U1-
01174      &       24.D0*U5*U2)*Y3*Y2+(U3**2-4.D0*U3*U5+12.D0*U4**2-
01175      &       2.D0*U2*U3+16.D0*U6**2+10.D0*U2**2+16.D0*U5**2-4.D0*U2*U1+
01176      &       12.D0*U4*U6-2.D0*U5*U1+6.D0*U4*U1+10.D0*U5*U2-4.D0*U6*U3-
01177      &       4.D0*U4*U3+8.D0*U5*U6-2.D0*U3*U1-2.D0*U2*U6+6.D0*U2*U4+
01178      &       12.D0*U4*U5+10.D0*U6*U1+10.D0*U1**2)*Y3**2+
01179      &       (36.D0*U5*V5-V1*U1-3.D0*U2*V3-6.D0*U4*V1+2.D0*V2*U6)*Y2*X3
01180       A44(IELEM) = ((-2.D0*V4*U3+ V1*U3+U1*V3+21.D0*U2*V2+12.D0*V4*U5-
01181      &             4.D0*U5*V1+12.D0*U4*V5+2.D0*U2*V6-4.D0*V5*U1-
01182      &             2.D0*U4*V3-4.D0*U6*V6+4.D0*U5*V6+4.D0*V5*U6+
01183      &             12.D0*U5*V2-4.D0*U5*V3-3.D0*V2*U3+12.D0*U2*V4-
01184      &             2.D0*U1*V6-2.D0*V2*U1-2.D0*U2*V1+12.D0*V2*U4-
01185      &             6.D0*V4*U1+U3*V3+12.D0*U2*V5-4.D0*U3*V5+12.D0*U4*V4-
01186      &             2.D0*V1*U6)*Y2*X3+(36.D0*U5*V5- V1*U1-3.D0*U2*V3-
01187      &             6.D0*U4*V1+2.D0*V2*U6-2.D0*V4*U3+V1*U3+U1*V3+
01188      &             21.D0*U2*V2+12.D0*V4*U5-4.D0*U5*V1+12.D0*U4*V5+
01189      &             2.D0*U2*V6-4.D0*V5*U1-2.D0*U4*V3-4.D0*U6*V6+
01190      &             4.D0*U5*V6+4.D0*V5*U6+12.D0*U5*V2-4.D0*U5*V3-
01191      &             3.D0*V2*U3+12.D0*U2*V4-2.D0*U1*V6-2.D0*V2*U1-
01192      &             2.D0*U2*V1+12.D0*V2*U4-6.D0*V4*U1+U3*V3+12.D0*U2*V5-
01193      &             4.D0*U3*V5+12.D0*U4*V4-2.D0*V1*U6)*Y3*X2+(U3*U1-
01194      &             2.D0*U6*U1+12.D0*U4*U6-3.D0*U2*U1+ 24.D0*U4**2+
01195      &             24.D0*U5**2+12.D0*U2*U4-6.D0*U5*U1-2.D0*U6*U3+U1**2+
01196      &             4.D0*U6**2+12.D0*U2**2+U3**2+24.D0*U4*U5-6.D0*U4*U3-
01197      &             6.D0*U4*U1+12.D0*U5*U6-3.D0*U2*U3-6.D0*U3*U5+
01198      &             12.D0*U5*U2)*Y2**2+ANS1+ANS2)*8.D0*AUX630
01199 !
01200       ANS1 = (-2.D0*V2*U3-4.D0*V1*U1+4.D0*U6*V3-4.D0*V4*U1+4.D0*V4*U3+
01201      &       4.D0*U4*V3-4.D0*V5*U1-8.D0*V4*U6+4.D0*U3*V5-8.D0*U4*V6-
01202      &       4.D0*U5*V1-4.D0*U4*V1+8.D0*V5*U6+4.D0*U5*V3+4.D0*U3*V6+
01203      &       2.D0*U2*V1-2.D0*U2*V3+2.D0*V2*U1-48.D0*U4*V4-4.D0*U1*V6+
01204      &       8.D0*U5*V6-4.D0*V1*U6+48.D0*U5*V5+4.D0*U3*V3)*Y2*X2+
01205      &       (-56.D0*U5*V5-20.D0*V1*U1+U2*V3-8.D0*U4*V1+4.D0*U6*V3+
01206      &       4.D0*V2*U6+4.D0*V4*U3+2.D0*V1*U3+2.D0*U1*V3-4.D0*U2*V2-
01207      &       16.D0*V4*U5+4.D0*U5*V1-16.D0*U4*V5+4.D0*U2*V6+4.D0*V5*U1+
01208      &       4.D0*U4*V3-24.D0*U6*V6-16.D0*V4*U6-16.D0*U4*V6-
01209      &       16.D0*U5*V6-16.D0*V5*U6+V2*U3+4.D0*U2*V4-8.D0*U1*V6+
01210      &       2.D0*V2*U1+4.D0*U3*V6+2.D0*U2*V1+4.D0*V2*U4-8.D0*V4*U1-
01211      &       4.D0*U3*V3-24.D0*U4*V4-8.D0*V1*U6)*Y3*X2+((-56.D0*U5*V5-
01212      &       20.D0*V1*U1+U2*V3-8.D0*U4*V1+4.D0*U6*V3+4.D0*V2*U6+4*V4*U3+
01213      &       2.D0*V1*U3+2.D0*U1*V3-4.D0*U2*V2-16.D0*V4*U5+4.D0*U5*V1-
01214      &       16.D0*U4*V5+4.D0*U2*V6+4.D0*V5*U1+4.D0*U4*V3-24.D0*U6*V6-
01215      &       16.D0*V4*U6-16.D0*U4*V6-16.D0*U5*V6-16.D0*V5*U6+V2*U3+
01216      &       4.D0*U2*V4-8.D0*U1*V6+2.D0*V2*U1+4.D0*U3*V6+2.D0*U2*V1+
01217      &       4.D0*V2*U4-8.D0*V4*U1-4.D0*U3*V3-24.D0*U4*V4-
01218      &       8.D0*V1*U6)*Y2)*X3
01219       ANS2 = (4.D0*U5*V2+4.D0*U2*V5+48.D0*U5*V5-8.D0*V4*U6+4.D0*U2*V4-
01220      &       8.D0*U4*V6+2.D0*U1*V3+8.D0*U4*V5-2.D0*U2*V3+8.D0*V4*U5-
01221      &       4.D0*V1*U6+4.D0*V2*U4+2.D0*V1*U3-4.D0*V4*U1-4.D0*V5*U1-
01222      &       4.D0*U5*V1-4.D0*U4*V1+4.D0*U2*V2-4.D0*U1*V6+4.D0*V2*U6-
01223      &       2.D0*V2*U3-4.D0*V1*U1-48.D0*U6*V6+4.D0*U2*V6)*Y3*X3+
01224      &       (-4.D0*U3*U1+32.D0*U4*U5-8.D0*U5*U1+16.D0*U6*U1-2.D0*U2*U3+
01225      &       32.D0*U5*U6+32.D0*U4*U6+24.D0*U4**2+24.D0*U6**2-8.D0*U4*U3+
01226      &       20.D0*U1**2-8.D0*U2*U4+4.D0*U3**2-8.D0*U6*U3+4.D0*U2**2-
01227      &       8.D0*U2*U6+56.D0*U5**2-4.D0*U2*U1+16.D0*U4*U1)*Y3*Y2+
01228      &       (24.D0*U6**2+2.D0*U1**2-4.D0*U2*U4-2.D0*U3*U1-2.D0*U2**2+
01229      &       4.D0*U5*U1+4.D0*U4*U1+8.D0*U4*U6-8.D0*U4*U5-4.D0*U2*U6+
01230      &       4.D0*U6*U1-4.D0*U5*U2-24.D0*U5**2+2.D0*U2*U3)*Y3**2+
01231      &       (-8.D0*V3*V6+16.D0*V1*V6+32.D0*V5*V6+4.D0*V3**2+4.D0*V2**2-
01232      &       4.D0*V2*V1+56.D0*V5**2-2.D0*V2*V3+20.D0*V1**2-8.D0*V2*V4+
01233      &       32.D0*V4*V6+24.D0*V6**2-8.D0*V5*V1+32.D0*V4*V5+16.D0*V4*V1-
01234      &       8.D0*V4*V3-8.D0*V2*V6-4.D0*V1*V3+24.D0*V4**2)*X3*X2
01235       A46(IELEM) = ((2.D0*V2*V3-4.D0*V2*V6-2.D0*V1*V3+4.D0*V1*V6+
01236      &       8.D0*V4*V6+4.D0*V4*V1-2.D0*V2**2-8.D0*V4*V5-24.D0*V5**2+
01237      &       24.D0*V6**2+2.D0*V1**2-4.D0*V2*V4+4.D0*V5*V1-
01238      &       4.D0*V5*V2)*X3**2+(-4.D0*V4*V3-8.D0*V5*V6-2.D0*V2*V1-
01239      &       2.D0*V3**2+4.D0*V1*V6-24.D0*V5**2+4.D0*V5*V1+2.D0*V2*V3+
01240      &       8.D0*V4*V6+24.D0*V4**2-4.D0*V3*V6+4.D0*V4*V1-4.D0*V3*V5+
01241      &       2.D0*V1**2)*X2**2+(-2.D0*U2*U1+4.D0*U6*U1+4.D0*U5*U1-
01242      &       4.D0*U4*U3-24.D0*U5**2+8.D0*U4*U6-8.D0*U5*U6+24.D0*U4**2-
01243      &       4.D0*U6*U3-2.D0*U3**2+4.D0*U4*U1+2.D0*U1**2+2.D0*U2*U3-
01244      &       4.D0*U3*U5)*Y2**2+ANS1+ANS2)*4.D0*(-AUX630)
01245 !
01246       ANS1 = 24.D0*V4**2*X2**2-48.D0*U5*Y2*V5*X2-V1*X2*U3*Y2+
01247      &       12.D0*U4*Y2*V4*X3+12.D0*U5*Y3*V4*X2+12.D0*V5*X2*U4*Y3+
01248      &       24.D0*U5**2*Y2**2-48.D0*U4*Y2*V4*X2+6.D0*U1*Y3*V6*X3+
01249      &       36.D0*U5*Y3*V5*X2+12.D0*V5*X3*U4*Y2-12.D0*U5*Y3*V4*X3+
01250      &       V2**2*X3**2+24.D0*V5**2*X3**2+36.D0*U5*Y2*V5*X3-
01251      &       12.D0*V5*X3*U4*Y3+12.D0*U5*Y2*V4*X3-24.D0*V5*X3*V4*X2+
01252      &       12.D0*V4*X2*U4*Y3-8.D0*V4*X3*U4*Y3+U2**2*Y3**2+
01253      &       4.D0*U4**2*Y3**2+4.D0*V1*X3*V6*X2-2.D0*V1*X3*U1*Y3+
01254      &       3.D0*U1*Y3*V3*X3+12.D0*V2**2*X2**2+6.D0*U2*Y3*V5*X3+
01255      &       6.D0*U4*Y2*V1*X2-12.D0*U4*Y2*V6*X2-2.D0*U1*Y3*V6*X2-
01256      &       12.D0*V4*X3*U6*Y3+4.D0*V4**2*X3**2-12.D0*U4*Y3*V6*X3+
01257      &       6.D0*V4*X2*U1*Y2-12.D0*V4*X2*U6*Y2+24.D0*U6**2*Y3**2-
01258      &       12.D0*U2*Y2*V4*X2+6.D0*U2*Y3*V6*X3+12.D0*U2**2*Y2**2+
01259      &       12.D0*U3**2*Y3**2-48.D0*V5*X3*U5*Y3+24.D0*U5*Y3**2*U6-
01260      &       24.D0*U5*Y2*U4*Y3+12.D0*V5*X2**2*V2+U1**2*Y3**2+
01261      &       V1**2*X3**2+12.D0*V3**2*X3**2+24.D0*V6**2*X3**2+
01262      &       12.D0*U5*Y2**2*U6-12.D0*V6**2*X3*X2-3.D0*V3**2*X3*X2+
01263      &       12.D0*U5*Y3**2*U3-3.D0*V2**2*X2*X3-12.D0*U6**2*Y2*Y3+
01264      &       24.D0*U5**2*Y3**2+24.D0*V5*X2**2*V4+12.D0*V2*X2**2*V4-
01265      &       6.D0*U5*Y2**2*U1-12.D0*V2*X2*U4*Y2+4.D0*U6**2*Y2**2
01266       ANS2 = -2.D0*V2*X3**2*V4+12.D0*U2*Y2**2*U4+12.D0*V5*X3**2*V3-
01267      &        3.D0*V2*X2**2*V3-6.D0*V4*X2**2*V3-24.D0*V2*X2*U2*Y2+
01268      &        12.D0*U3*Y3**2*U6+V1**2*X2**2-U1**2*Y3*Y2-
01269      &        2.D0*U1*Y2**2*U6+4.D0*V6**2*X2**2-6.D0*U5*Y3**2*U1-
01270      &        3.D0*V2*X3**2*V3-3.D0*U2**2*Y2*Y3-6.D0*V4*X2**2*V1-
01271      &        6.D0*V5*X2**2*V1-12.D0*V5*X2*U2*Y2-12.D0*U5*Y2*V2*X2-
01272      &        6.D0*V2*X3**2*V5-2.D0*V4*X3**2*V1-6.D0*V3*X2**2*V5-
01273      &        2.D0*V1*X2**2*V6+3.D0*U2*Y2*V3*X2+24.D0*V5**2*X2**2-
01274      &        3.D0*V1*X3**2*V3+3.D0*V2*X3*U3*Y3+V3**2*X2**2-
01275      &        6.D0*U4*Y2**2*U3+6.D0*V3*X2*U5*Y2+3.D0*U2*Y3*V3*X3-
01276      &        6.D0*U2*Y3**2*U6-6.D0*U3*Y2**2*U5+U2*Y3**2*U1-
01277      &        3.D0*U1*Y3**2*U3+12.D0*U4*Y3**2*U6-2.D0*V1*X2*U1*Y2+
01278      &        6.D0*V1*X3*U6*Y3-12.D0*U4**2*Y2*Y3-12.D0*V4**2*X3*X2+
01279      &        U1*Y2**2*U3+12.D0*V6*X3**2*V3+6.D0*U3*Y2*V5*X2-
01280      &        24.D0*V5*X2*U4*Y2-24.D0*U5*Y2*V4*X2-3.D0*U3**2*Y3*Y2+
01281      &        12.D0*U5*Y2**2*U2-V1**2*X3*X2+V2*X3**2*V1+
01282      &        3.D0*V2*X2*U3*Y2-6*U4*Y2**2*U1+V1*X2**2*V3+
01283      &        12.D0*V5*X2**2*V6-3.D0*U2*Y2**2*U3
01284       ANS3 = -6.D0*U2*Y3**2*U5-2.D0*U2*Y3**2*U4-6.D0*V5*X3**2*V1-
01285      &       36.D0*U5**2*Y2*Y3-36.D0*V5**2*X3*X2+12.D0*V5*X3**2*V4+
01286      &       12.D0*U5*Y3**2*U4+U3**2*Y2**2+6.D0*V2*X3*U5*Y3+
01287      &       12.D0*V6*X2*U6*Y3-8.D0*V6*X2*U6*Y2+2.D0*V6*X2*U3*Y2-
01288      &       2.D0*V3*X2**2*V6+2*U2*Y3*V4*X3+24.D0*V5*X3**2*V6-
01289      &       2.D0*V3*X2*U3*Y2+2.D0*V3*X2*U6*Y2+U1**2*Y2**2-
01290      &       U1*Y2*V3*X2+2.D0*V2*X3*U4*Y3+12.D0*V6*X3*U6*Y2+
01291      &       2.D0*U1*Y2*V6*X2+6.D0*V4*X2*U3*Y2+3.D0*V2*X2*U2*Y3+
01292      &       3.D0*V1*X3*U3*Y3+3.D0*U2*Y2*V2*X3-2.D0*U2*Y3*V2*X3+
01293      &       12.D0*V4*X3**2*V6-4.D0*U5*Y3*V1*X2+6.D0*U5*Y3*V1*X3-
01294      &       3.D0*U2*Y3**2*U3+4.D0*U1*Y3*U6*Y2+12.D0*U4*Y2**2*U6-
01295      &       12.D0*U5*Y3*V3*X3-2.D0*V1*X3*U6*Y2+24.D0*U4**2*Y2**2-
01296      &       24.D0*U5*Y2*U6*Y3-6.D0*V1*X3**2*V6+6.D0*U4*Y2*V3*X2+
01297      &       12.D0*V5*X3*U6*Y2+6.D0*V5*X3*U1*Y3+8.D0*V5*X3*V1*X2-
01298      &       12.D0*V5*X2*U6*Y2-6.D0*U1*Y3**2*U6+6.D0*V5*X2*U1*Y2+
01299      &       8.D0*U5*Y2*U1*Y3+24.D0*U5*Y2**2*U4+3.D0*V3*X3*U3*Y2+
01300      &       3.D0*V3*X2*U3*Y3-24.D0*V5*X3*V6*X2+V1*X2*U1*Y3-
01301      &       2.D0*V1*X2*U6*Y3-4.D0*V5*X3*U1*Y2-12.D0*V5*X3*U3*Y3-
01302      &       2.D0*U4*Y3**2*U1-24*V3*X3*U3*Y3
01303       ANS4 = -48.D0*V6*X3*U6*Y3-4.D0*U5*Y2*V1*X3+12.D0*V4*X2**2*V6-
01304      &       12.D0*V3*X3*U6*Y3+12.D0*U5*Y2*V6*X3+8.D0*U4*Y3*V6*X2-
01305      &       U2*Y3*V3*X2-2.D0*U4*Y3*V3*X2-2.D0*U4*Y3*V1*X2+
01306      &       2.D0*U4*Y3*V1*X3-2.D0*U1*Y2*V6*X3-12.D0*V6*X3*U3*Y3-
01307      &       2.D0*V2*X3*U6*Y2+4.D0*V4*X3*V1*X2-2.D0*V4*X2*U3*Y3-
01308      &       V2*X3*U3*Y2-2.D0*V4*X2*U1*Y3+4.D0*V2*X3*V6*X2-
01309      &       12.D0*U5*Y2*V6*X2-U2*Y3*V1*X3+4.D0*U2*Y3*U6*Y2+
01310      &       8.D0*V4*X2*U6*Y3+2.D0*U2*Y3*U3*Y2-2.D0*U2*Y3*V6*X2+
01311      &       2.D0*V1*X2*U6*Y2+8.D0*V4*X3*U6*Y2-2.D0*V4*X3*U3*Y2+
01312      &       2.D0*V2*X3*V3*X2-16.D0*V4*X3*V6*X2-2.D0*V4*X3*U1*Y2-
01313      &       6.D0*V2*X3**2*V6+6.D0*U5*Y2*V1*X2-V2*X3*U1*Y3+
01314      &       6.D0*V2*X3*U6*Y3-2.D0*V2*X2*U6*Y3+8.D0*U4*Y2*V6*X3-
01315      &       V2*X2*U3*Y3-2.D0*U3*Y2**2*U6+12.D0*V5*X2*U6*Y3+
01316      &       4.D0*V4*X3*V3*X2+2.D0*V4*X3*U1*Y3-16.D0*U4*Y2*U6*Y3+
01317      &       3.D0*U2*Y2*V1*X2+4.D0*U4*Y2*U3*Y3-2.D0*U4*Y2*V1*X3-
01318      &       2.D0*U4*Y2*V3*X3-U2*Y2*V3*X3-24.D0*U5*Y3*V6*X3-
01319      &       24.D0*V5*X3*U6*Y3-4.D0*V5*X2*U1*Y3-2.D0*U2*Y2*V6*X3+
01320      &       3.D0*V2*X2*U1*Y2+4.D0*U4*Y2*U1*Y3+12.D0*U5*Y3*V6*X2-
01321      &       3.D0*V2*X2**2*V1+V1*X3*U1*Y2-3.D0*U2*Y2**2*U1
01322         A55(IELEM) = (ANS1+ANS2+ANS3+ANS4)*8.D0*AUX630
01323 !
01324 !
01325 ! USES HERE THE 'MAGIC SQUARE' PROPERTIES TO GET THE REMAINING TERMS
01326 ! (SUM OF EACH LINE = SUM OF EACH COLUMN = 0)
01327 !
01328         A16(IELEM) = - A11(IELEM) - A12(IELEM) - A13(IELEM)
01329      &               - A14(IELEM) - A15(IELEM)
01330 !
01331         A25(IELEM) = - A12(IELEM) - A22(IELEM) - A23(IELEM)
01332      &               - A24(IELEM) - A26(IELEM)
01333 !
01334         A36(IELEM) = - A13(IELEM) - A23(IELEM) - A33(IELEM)
01335      &               - A34(IELEM) - A35(IELEM)
01336 !
01337         A45(IELEM) = - A14(IELEM) - A24(IELEM) - A34(IELEM)
01338      &               - A44(IELEM) - A46(IELEM)
01339 !
01340         A56(IELEM) = - A15(IELEM) - A25(IELEM) - A35(IELEM)
01341      &               - A45(IELEM) - A55(IELEM)
01342 !
01343         A66(IELEM) = - A16(IELEM) - A26(IELEM) - A36(IELEM)
01344      &                - A46(IELEM) - A56(IELEM)
01345 !
01346       ENDDO ! IELEM
01347 !
01348 !-----------------------------------------------------------------------
01349 !
01350       ELSE
01351         IF(IELMU.EQ.IELMV) THEN
01352         IF (LNG.EQ.1) WRITE(LU,100) IELMU
01353         IF (LNG.EQ.2) WRITE(LU,101) IELMU
01354 100     FORMAT(1X,'MT04CC (BIEF) :',/,
01355      &         1X,'DISCRETISATION DE U ET V : ',1I6,' NON PREVUE')
01356 101     FORMAT(1X,'MT04CC (BIEF) :',/,
01357      &         1X,'DISCRETIZATION OF U AND V : ',1I6,' NOT AVAILABLE')
01358         ELSE
01359         IF (LNG.EQ.1) WRITE(LU,200) IELMU,IELMV
01360         IF (LNG.EQ.2) WRITE(LU,201) IELMU,IELMV
01361 200     FORMAT(1X,'MT04CC (BIEF) :',/,
01362      &         1X,'U ET V DE DISCRETISATIONS DIFFERENTES :',1I6,3X,1I6)
01363 201     FORMAT(1X,'MT04CC (BIEF) :',/,
01364      &         1X,'U AND V OF A DIFFERENT DISCRETISATION:',1I6,3X,1I6)
01365         ENDIF
01366 !
01367         CALL PLANTE(1)
01368         STOP
01369 !
01370       ENDIF
01371 !
01372 !-----------------------------------------------------------------------
01373 !
01374       RETURN
01375       END

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