mv0606.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\mv0606.f
00002 !
00083                      SUBROUTINE MV0606
00084 !                    *****************
00085 !
00086      &(OP, X , DA,TYPDIA,XA,TYPEXT, Y,C,
00087      & IKLE1,IKLE2,IKLE3,IKLE4,IKLE5,IKLE6,
00088      & NPOIN,NELEM,NELMAX,W1,W2,W3,W4,W5,W6)
00089 !
00090 !***********************************************************************
00091 ! BIEF   V6P1                                   21/08/2010
00092 !***********************************************************************
00093 !
00094 !
00095 !
00096 !
00097 !
00098 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00099 !| C              |-->| A GIVEN CONSTANT
00100 !| DA             |-->| MATRIX DIAGONAL
00101 !| IKLE1          |-->| FIRST POINTS OF ELEMENTS
00102 !| IKLE2          |-->| SECOND POINTS OF ELEMENTS
00103 !| IKLE3          |-->| THIRD POINTS OF ELEMENTS
00104 !| IKLE4          |-->| FOURTH POINTS OF ELEMENTS
00105 !| IKLE5          |-->| FIFTH POINTS OF ELEMENTS
00106 !| IKLE6          |-->| SIXTH POINTS OF ELEMENTS
00107 !| NELEM          |-->| NUMBER OF ELEMENTS
00108 !| NPOIN          |-->| NUMBER OF LINEAR POINTS
00109 !| NPT2           |-->| NUMBER OF QUADRATIC POINTS
00110 !| OP             |-->| OPERATION TO BE DONE (SEE ABOVE)
00111 !| TYPDIA         |-->| TYPE OF DIAGONAL:
00112 !|                |   | TYPDIA = 'Q' : ANY VALUE
00113 !|                |   | TYPDIA = 'I' : IDENTITY
00114 !|                |   | TYPDIA = '0' : ZERO
00115 !| TYPEXT         |-->| TYPE OF OFF-DIAGONAL TERMS
00116 !|                |   | TYPEXT = 'Q' : ANY VALUE
00117 !|                |   | TYPEXT = 'S' : SYMMETRIC
00118 !|                |   | TYPEXT = '0' : ZERO
00119 !| W1             |<->| RESULT IN NON ASSEMBLED FORM
00120 !| W2             |<->| RESULT IN NON ASSEMBLED FORM
00121 !| W3             |<->| RESULT IN NON ASSEMBLED FORM
00122 !| W4             |<->| RESULT IN NON ASSEMBLED FORM
00123 !| W5             |<->| RESULT IN NON ASSEMBLED FORM
00124 !| W6             |<->| RESULT IN NON ASSEMBLED FORM
00125 !| XA             |<->| RESULT IN ASSEMBLED FORM
00126 !| Y              |-->| VECTOR USED IN THE OPERATION
00127 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00128 !
00129       USE BIEF, EX_MV0606 => MV0606
00130 !
00131       IMPLICIT NONE
00132       INTEGER LNG,LU
00133       COMMON/INFO/LNG,LU
00134 !
00135 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00136 !
00137       INTEGER, INTENT(IN) :: NELEM,NELMAX,NPOIN
00138       INTEGER, INTENT(IN) :: IKLE1(*),IKLE2(*),IKLE3(*)
00139       INTEGER, INTENT(IN) :: IKLE4(*),IKLE5(*),IKLE6(*)
00140 !
00141       DOUBLE PRECISION, INTENT(INOUT) :: W1(*),W2(*),W3(*)
00142       DOUBLE PRECISION, INTENT(INOUT) :: W4(*),W5(*),W6(*)
00143       DOUBLE PRECISION, INTENT(IN) :: Y(*),DA(*)
00144       DOUBLE PRECISION, INTENT(INOUT) :: X(*)
00145       DOUBLE PRECISION, INTENT(IN) :: XA(NELMAX,*)
00146       DOUBLE PRECISION, INTENT(IN) :: C
00147 !
00148       CHARACTER(LEN=8), INTENT(IN) :: OP
00149       CHARACTER(LEN=1), INTENT(IN) :: TYPDIA,TYPEXT
00150 !
00151 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00152 !
00153       INTEGER IELEM,I1,I2,I3,I4,I5,I6
00154       DOUBLE PRECISION Z(1)
00155 !
00156 !-----------------------------------------------------------------------
00157 !
00158       IF(OP(1:8).EQ.'X=AY    ') THEN
00159 !
00160 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
00161 !
00162         IF(TYPEXT(1:1).EQ.'S') THEN
00163 !
00164           DO IELEM = 1 , NELEM
00165 !
00166           I1 = IKLE1(IELEM)
00167           I2 = IKLE2(IELEM)
00168           I3 = IKLE3(IELEM)
00169           I4 = IKLE4(IELEM)
00170           I5 = IKLE5(IELEM)
00171           I6 = IKLE6(IELEM)
00172 !
00173           W1(IELEM) =
00174      &                  + XA(IELEM,1) * Y(I2)
00175      &                  + XA(IELEM,2) * Y(I3)
00176      &                  + XA(IELEM,3) * Y(I4)
00177      &                  + XA(IELEM,4) * Y(I5)
00178      &                  + XA(IELEM,5) * Y(I6)
00179 !
00180           W2(IELEM) =
00181      &                  + XA(IELEM,1) * Y(I1)
00182      &                  + XA(IELEM,6) * Y(I3)
00183      &                  + XA(IELEM,7) * Y(I4)
00184      &                  + XA(IELEM,8) * Y(I5)
00185      &                  + XA(IELEM,9) * Y(I6)
00186 !
00187           W3(IELEM) =
00188      &                  + XA(IELEM,2)  * Y(I1)
00189      &                  + XA(IELEM,6)  * Y(I2)
00190      &                  + XA(IELEM,10) * Y(I4)
00191      &                  + XA(IELEM,11) * Y(I5)
00192      &                  + XA(IELEM,12) * Y(I6)
00193 !
00194           W4(IELEM) =
00195      &                  + XA(IELEM,3)  * Y(I1)
00196      &                  + XA(IELEM,7)  * Y(I2)
00197      &                  + XA(IELEM,10) * Y(I3)
00198      &                  + XA(IELEM,13) * Y(I5)
00199      &                  + XA(IELEM,14) * Y(I6)
00200 !
00201           W5(IELEM) =
00202      &                  + XA(IELEM,4)  * Y(I1)
00203      &                  + XA(IELEM,8)  * Y(I2)
00204      &                  + XA(IELEM,11) * Y(I3)
00205      &                  + XA(IELEM,13) * Y(I4)
00206      &                  + XA(IELEM,15) * Y(I6)
00207 !
00208           W6(IELEM) =
00209      &                  + XA(IELEM,5)  * Y(I1)
00210      &                  + XA(IELEM,9)  * Y(I2)
00211      &                  + XA(IELEM,12) * Y(I3)
00212      &                  + XA(IELEM,14) * Y(I4)
00213      &                  + XA(IELEM,15) * Y(I5)
00214 !
00215           ENDDO ! IELEM
00216 !
00217         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
00218 !
00219           DO IELEM = 1 , NELEM
00220 !
00221           I1 = IKLE1(IELEM)
00222           I2 = IKLE2(IELEM)
00223           I3 = IKLE3(IELEM)
00224           I4 = IKLE4(IELEM)
00225           I5 = IKLE5(IELEM)
00226           I6 = IKLE6(IELEM)
00227 !
00228           W1(IELEM) =
00229      &                  + XA(IELEM, 1) * Y(I2)
00230      &                  + XA(IELEM, 2) * Y(I3)
00231      &                  + XA(IELEM, 3) * Y(I4)
00232      &                  + XA(IELEM, 4) * Y(I5)
00233      &                  + XA(IELEM, 5) * Y(I6)
00234 !
00235           W2(IELEM) =
00236      &                  + XA(IELEM,16) * Y(I1)
00237      &                  + XA(IELEM, 6) * Y(I3)
00238      &                  + XA(IELEM, 7) * Y(I4)
00239      &                  + XA(IELEM, 8) * Y(I5)
00240      &                  + XA(IELEM, 9) * Y(I6)
00241 !
00242           W3(IELEM) =
00243      &                  + XA(IELEM,17) * Y(I1)
00244      &                  + XA(IELEM,21) * Y(I2)
00245      &                  + XA(IELEM,10) * Y(I4)
00246      &                  + XA(IELEM,11) * Y(I5)
00247      &                  + XA(IELEM,12) * Y(I6)
00248 !
00249           W4(IELEM) =
00250      &                  + XA(IELEM,18) * Y(I1)
00251      &                  + XA(IELEM,22) * Y(I2)
00252      &                  + XA(IELEM,25) * Y(I3)
00253      &                  + XA(IELEM,13) * Y(I5)
00254      &                  + XA(IELEM,14) * Y(I6)
00255 !
00256           W5(IELEM) =
00257      &                  + XA(IELEM,19) * Y(I1)
00258      &                  + XA(IELEM,23) * Y(I2)
00259      &                  + XA(IELEM,26) * Y(I3)
00260      &                  + XA(IELEM,28) * Y(I4)
00261      &                  + XA(IELEM,15) * Y(I6)
00262 !
00263           W6(IELEM) =
00264      &                  + XA(IELEM,20) * Y(I1)
00265      &                  + XA(IELEM,24) * Y(I2)
00266      &                  + XA(IELEM,27) * Y(I3)
00267      &                  + XA(IELEM,29) * Y(I4)
00268      &                  + XA(IELEM,30) * Y(I5)
00269 !
00270           ENDDO ! IELEM
00271 !
00272         ELSEIF(TYPEXT(1:1).EQ.'0') THEN
00273 !
00274           CALL OV ('X=C     ', W1 , Y , Z , 0.D0 , NELEM )
00275           CALL OV ('X=C     ', W2 , Y , Z , 0.D0 , NELEM )
00276           CALL OV ('X=C     ', W3 , Y , Z , 0.D0 , NELEM )
00277           CALL OV ('X=C     ', W4 , Y , Z , 0.D0 , NELEM )
00278           CALL OV ('X=C     ', W5 , Y , Z , 0.D0 , NELEM )
00279           CALL OV ('X=C     ', W6 , Y , Z , 0.D0 , NELEM )
00280 !
00281         ELSE
00282 !
00283           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
00284           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
00285           CALL PLANTE(0)
00286           STOP
00287 !
00288         ENDIF
00289 !
00290 !   CONTRIBUTION OF THE DIAGONAL:
00291 !
00292         IF(TYPDIA(1:1).EQ.'Q') THEN
00293           CALL OV ('X=YZ    ', X , Y , DA , C  , NPOIN )
00294         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
00295           CALL OV ('X=Y     ', X , Y , Z  , C  , NPOIN )
00296         ELSEIF(TYPDIA(1:1).EQ.'0') THEN
00297           CALL OV ('X=C     ', X , Y , DA , 0.D0 , NPOIN )
00298         ELSE
00299           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
00300           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
00301           CALL PLANTE(0)
00302           STOP
00303         ENDIF
00304 !
00305 !-----------------------------------------------------------------------
00306 !
00307       ELSEIF(OP(1:8).EQ.'X=-AY   ') THEN
00308 !
00309 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
00310 !
00311         IF(TYPEXT(1:1).EQ.'S') THEN
00312 !
00313           DO IELEM = 1 , NELEM
00314 !
00315           I1 = IKLE1(IELEM)
00316           I2 = IKLE2(IELEM)
00317           I3 = IKLE3(IELEM)
00318           I4 = IKLE4(IELEM)
00319           I5 = IKLE5(IELEM)
00320           I6 = IKLE6(IELEM)
00321 !
00322           W1(IELEM) =
00323      &                  - XA(IELEM,1) * Y(I2)
00324      &                  - XA(IELEM,2) * Y(I3)
00325      &                  - XA(IELEM,3) * Y(I4)
00326      &                  - XA(IELEM,4) * Y(I5)
00327      &                  - XA(IELEM,5) * Y(I6)
00328 !
00329           W2(IELEM) =
00330      &                  - XA(IELEM,1) * Y(I1)
00331      &                  - XA(IELEM,6) * Y(I3)
00332      &                  - XA(IELEM,7) * Y(I4)
00333      &                  - XA(IELEM,8) * Y(I5)
00334      &                  - XA(IELEM,9) * Y(I6)
00335 !
00336           W3(IELEM) =
00337      &                  - XA(IELEM,2)  * Y(I1)
00338      &                  - XA(IELEM,6)  * Y(I2)
00339      &                  - XA(IELEM,10) * Y(I4)
00340      &                  - XA(IELEM,11) * Y(I5)
00341      &                  - XA(IELEM,12) * Y(I6)
00342 !
00343           W4(IELEM) =
00344      &                  - XA(IELEM,3)  * Y(I1)
00345      &                  - XA(IELEM,7)  * Y(I2)
00346      &                  - XA(IELEM,10) * Y(I3)
00347      &                  - XA(IELEM,13) * Y(I5)
00348      &                  - XA(IELEM,14) * Y(I6)
00349 !
00350           W5(IELEM) =
00351      &                  - XA(IELEM,4)  * Y(I1)
00352      &                  - XA(IELEM,8)  * Y(I2)
00353      &                  - XA(IELEM,11) * Y(I3)
00354      &                  - XA(IELEM,13) * Y(I4)
00355      &                  - XA(IELEM,15) * Y(I6)
00356 !
00357           W6(IELEM) =
00358      &                  - XA(IELEM,5)  * Y(I1)
00359      &                  - XA(IELEM,9)  * Y(I2)
00360      &                  - XA(IELEM,12) * Y(I3)
00361      &                  - XA(IELEM,14) * Y(I4)
00362      &                  - XA(IELEM,15) * Y(I5)
00363 !
00364           ENDDO ! IELEM
00365 !
00366         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
00367 !
00368           DO IELEM = 1 , NELEM
00369 !
00370           I1 = IKLE1(IELEM)
00371           I2 = IKLE2(IELEM)
00372           I3 = IKLE3(IELEM)
00373           I4 = IKLE4(IELEM)
00374           I5 = IKLE5(IELEM)
00375           I6 = IKLE6(IELEM)
00376 !
00377           W1(IELEM) =
00378      &                  - XA(IELEM, 1) * Y(I2)
00379      &                  - XA(IELEM, 2) * Y(I3)
00380      &                  - XA(IELEM, 3) * Y(I4)
00381      &                  - XA(IELEM, 4) * Y(I5)
00382      &                  - XA(IELEM, 5) * Y(I6)
00383 !
00384           W2(IELEM) =
00385      &                  - XA(IELEM,16) * Y(I1)
00386      &                  - XA(IELEM, 6) * Y(I3)
00387      &                  - XA(IELEM, 7) * Y(I4)
00388      &                  - XA(IELEM, 8) * Y(I5)
00389      &                  - XA(IELEM, 9) * Y(I6)
00390 !
00391           W3(IELEM) =
00392      &                  - XA(IELEM,17) * Y(I1)
00393      &                  - XA(IELEM,21) * Y(I2)
00394      &                  - XA(IELEM,10) * Y(I4)
00395      &                  - XA(IELEM,11) * Y(I5)
00396      &                  - XA(IELEM,12) * Y(I6)
00397 !
00398           W4(IELEM) =
00399      &                  - XA(IELEM,18) * Y(I1)
00400      &                  - XA(IELEM,22) * Y(I2)
00401      &                  - XA(IELEM,25) * Y(I3)
00402      &                  - XA(IELEM,13) * Y(I5)
00403      &                  - XA(IELEM,14) * Y(I6)
00404 !
00405           W5(IELEM) =
00406      &                  - XA(IELEM,19) * Y(I1)
00407      &                  - XA(IELEM,23) * Y(I2)
00408      &                  - XA(IELEM,26) * Y(I3)
00409      &                  - XA(IELEM,28) * Y(I4)
00410      &                  - XA(IELEM,15) * Y(I6)
00411 !
00412           W6(IELEM) =
00413      &                  - XA(IELEM,20) * Y(I1)
00414      &                  - XA(IELEM,24) * Y(I2)
00415      &                  - XA(IELEM,27) * Y(I3)
00416      &                  - XA(IELEM,29) * Y(I4)
00417      &                  - XA(IELEM,30) * Y(I5)
00418 !
00419           ENDDO ! IELEM
00420 !
00421         ELSEIF(TYPEXT(1:1).EQ.'0') THEN
00422 !
00423           CALL OV ('X=C     ', W1 , Y , Z , 0.D0 , NELEM )
00424           CALL OV ('X=C     ', W2 , Y , Z , 0.D0 , NELEM )
00425           CALL OV ('X=C     ', W3 , Y , Z , 0.D0 , NELEM )
00426           CALL OV ('X=C     ', W4 , Y , Z , 0.D0 , NELEM )
00427           CALL OV ('X=C     ', W5 , Y , Z , 0.D0 , NELEM )
00428           CALL OV ('X=C     ', W6 , Y , Z , 0.D0 , NELEM )
00429 !
00430         ELSE
00431 !
00432           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
00433           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
00434           CALL PLANTE(0)
00435           STOP
00436 !
00437         ENDIF
00438 !
00439 !   CONTRIBUTION OF THE DIAGONAL:
00440 !
00441         IF(TYPDIA(1:1).EQ.'Q') THEN
00442           CALL OV ('X=-YZ   ', X , Y , DA , C  , NPOIN )
00443         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
00444           CALL OV ('X=-Y    ', X , Y , Z  , C  , NPOIN )
00445         ELSEIF(TYPDIA(1:1).EQ.'0') THEN
00446           CALL OV ('X=C     ', X , Y , DA , 0.D0 , NPOIN )
00447         ELSE
00448           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
00449           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
00450           CALL PLANTE(0)
00451           STOP
00452         ENDIF
00453 !
00454 !-----------------------------------------------------------------------
00455 !
00456       ELSEIF(OP(1:8).EQ.'X=X+AY  ') THEN
00457 !
00458 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
00459 !
00460         IF(TYPEXT(1:1).EQ.'S') THEN
00461 !
00462           DO IELEM = 1 , NELEM
00463 !
00464           I1 = IKLE1(IELEM)
00465           I2 = IKLE2(IELEM)
00466           I3 = IKLE3(IELEM)
00467           I4 = IKLE4(IELEM)
00468           I5 = IKLE5(IELEM)
00469           I6 = IKLE6(IELEM)
00470 !
00471           W1(IELEM) = W1(IELEM)
00472      &                  + XA(IELEM,1) * Y(I2)
00473      &                  + XA(IELEM,2) * Y(I3)
00474      &                  + XA(IELEM,3) * Y(I4)
00475      &                  + XA(IELEM,4) * Y(I5)
00476      &                  + XA(IELEM,5) * Y(I6)
00477 !
00478           W2(IELEM) = W2(IELEM)
00479      &                  + XA(IELEM,1) * Y(I1)
00480      &                  + XA(IELEM,6) * Y(I3)
00481      &                  + XA(IELEM,7) * Y(I4)
00482      &                  + XA(IELEM,8) * Y(I5)
00483      &                  + XA(IELEM,9) * Y(I6)
00484 !
00485           W3(IELEM) = W3(IELEM)
00486      &                  + XA(IELEM,2)  * Y(I1)
00487      &                  + XA(IELEM,6)  * Y(I2)
00488      &                  + XA(IELEM,10) * Y(I4)
00489      &                  + XA(IELEM,11) * Y(I5)
00490      &                  + XA(IELEM,12) * Y(I6)
00491 !
00492           W4(IELEM) = W4(IELEM)
00493      &                  + XA(IELEM,3)  * Y(I1)
00494      &                  + XA(IELEM,7)  * Y(I2)
00495      &                  + XA(IELEM,10) * Y(I3)
00496      &                  + XA(IELEM,13) * Y(I5)
00497      &                  + XA(IELEM,14) * Y(I6)
00498 !
00499           W5(IELEM) = W5(IELEM)
00500      &                  + XA(IELEM,4)  * Y(I1)
00501      &                  + XA(IELEM,8)  * Y(I2)
00502      &                  + XA(IELEM,11) * Y(I3)
00503      &                  + XA(IELEM,13) * Y(I4)
00504      &                  + XA(IELEM,15) * Y(I6)
00505 !
00506           W6(IELEM) = W6(IELEM)
00507      &                  + XA(IELEM,5)  * Y(I1)
00508      &                  + XA(IELEM,9)  * Y(I2)
00509      &                  + XA(IELEM,12) * Y(I3)
00510      &                  + XA(IELEM,14) * Y(I4)
00511      &                  + XA(IELEM,15) * Y(I5)
00512 !
00513           ENDDO ! IELEM
00514 !
00515         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
00516 !
00517           DO IELEM = 1 , NELEM
00518 !
00519           I1 = IKLE1(IELEM)
00520           I2 = IKLE2(IELEM)
00521           I3 = IKLE3(IELEM)
00522           I4 = IKLE4(IELEM)
00523           I5 = IKLE5(IELEM)
00524           I6 = IKLE6(IELEM)
00525 !
00526           W1(IELEM) = W1(IELEM)
00527      &                  + XA(IELEM, 1) * Y(I2)
00528      &                  + XA(IELEM, 2) * Y(I3)
00529      &                  + XA(IELEM, 3) * Y(I4)
00530      &                  + XA(IELEM, 4) * Y(I5)
00531      &                  + XA(IELEM, 5) * Y(I6)
00532 !
00533           W2(IELEM) = W2(IELEM)
00534      &                  + XA(IELEM,16) * Y(I1)
00535      &                  + XA(IELEM, 6) * Y(I3)
00536      &                  + XA(IELEM, 7) * Y(I4)
00537      &                  + XA(IELEM, 8) * Y(I5)
00538      &                  + XA(IELEM, 9) * Y(I6)
00539 !
00540           W3(IELEM) = W3(IELEM)
00541      &                  + XA(IELEM,17) * Y(I1)
00542      &                  + XA(IELEM,21) * Y(I2)
00543      &                  + XA(IELEM,10) * Y(I4)
00544      &                  + XA(IELEM,11) * Y(I5)
00545      &                  + XA(IELEM,12) * Y(I6)
00546 !
00547           W4(IELEM) = W4(IELEM)
00548      &                  + XA(IELEM,18) * Y(I1)
00549      &                  + XA(IELEM,22) * Y(I2)
00550      &                  + XA(IELEM,25) * Y(I3)
00551      &                  + XA(IELEM,13) * Y(I5)
00552      &                  + XA(IELEM,14) * Y(I6)
00553 !
00554           W5(IELEM) = W5(IELEM)
00555      &                  + XA(IELEM,19) * Y(I1)
00556      &                  + XA(IELEM,23) * Y(I2)
00557      &                  + XA(IELEM,26) * Y(I3)
00558      &                  + XA(IELEM,28) * Y(I4)
00559      &                  + XA(IELEM,15) * Y(I6)
00560 !
00561           W6(IELEM) = W6(IELEM)
00562      &                  + XA(IELEM,20) * Y(I1)
00563      &                  + XA(IELEM,24) * Y(I2)
00564      &                  + XA(IELEM,27) * Y(I3)
00565      &                  + XA(IELEM,29) * Y(I4)
00566      &                  + XA(IELEM,30) * Y(I5)
00567 !
00568           ENDDO ! IELEM
00569 !
00570         ELSEIF(TYPEXT(1:1).NE.'0') THEN
00571 !
00572           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
00573           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
00574           CALL PLANTE(0)
00575           STOP
00576 !
00577         ENDIF
00578 !
00579 !   CONTRIBUTION OF THE DIAGONAL:
00580 !
00581         IF(TYPDIA(1:1).EQ.'Q') THEN
00582           CALL OV ('X=X+YZ  ', X , Y , DA , C , NPOIN )
00583         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
00584           CALL OV ('X=X+Y   ', X , Y , Z  , C  , NPOIN )
00585         ELSEIF(TYPDIA(1:1).NE.'0') THEN
00586           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
00587           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
00588           CALL PLANTE(0)
00589           STOP
00590         ENDIF
00591 !
00592 !-----------------------------------------------------------------------
00593 !
00594       ELSEIF(OP(1:8).EQ.'X=X-AY  ') THEN
00595 !
00596 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
00597 !
00598         IF(TYPEXT(1:1).EQ.'S') THEN
00599 !
00600           DO IELEM = 1 , NELEM
00601 !
00602           I1 = IKLE1(IELEM)
00603           I2 = IKLE2(IELEM)
00604           I3 = IKLE3(IELEM)
00605           I4 = IKLE4(IELEM)
00606           I5 = IKLE5(IELEM)
00607           I6 = IKLE6(IELEM)
00608 !
00609           W1(IELEM) = W1(IELEM)
00610      &                  - XA(IELEM,1) * Y(I2)
00611      &                  - XA(IELEM,2) * Y(I3)
00612      &                  - XA(IELEM,3) * Y(I4)
00613      &                  - XA(IELEM,4) * Y(I5)
00614      &                  - XA(IELEM,5) * Y(I6)
00615 !
00616           W2(IELEM) = W2(IELEM)
00617      &                  - XA(IELEM,1) * Y(I1)
00618      &                  - XA(IELEM,6) * Y(I3)
00619      &                  - XA(IELEM,7) * Y(I4)
00620      &                  - XA(IELEM,8) * Y(I5)
00621      &                  - XA(IELEM,9) * Y(I6)
00622 !
00623           W3(IELEM) = W3(IELEM)
00624      &                  - XA(IELEM,2)  * Y(I1)
00625      &                  - XA(IELEM,6)  * Y(I2)
00626      &                  - XA(IELEM,10) * Y(I4)
00627      &                  - XA(IELEM,11) * Y(I5)
00628      &                  - XA(IELEM,12) * Y(I6)
00629 !
00630           W4(IELEM) = W4(IELEM)
00631      &                  - XA(IELEM,3)  * Y(I1)
00632      &                  - XA(IELEM,7)  * Y(I2)
00633      &                  - XA(IELEM,10) * Y(I3)
00634      &                  - XA(IELEM,13) * Y(I5)
00635      &                  - XA(IELEM,14) * Y(I6)
00636 !
00637           W5(IELEM) = W5(IELEM)
00638      &                  - XA(IELEM,4)  * Y(I1)
00639      &                  - XA(IELEM,8)  * Y(I2)
00640      &                  - XA(IELEM,11) * Y(I3)
00641      &                  - XA(IELEM,13) * Y(I4)
00642      &                  - XA(IELEM,15) * Y(I6)
00643 !
00644           W6(IELEM) = W6(IELEM)
00645      &                  - XA(IELEM,5)  * Y(I1)
00646      &                  - XA(IELEM,9)  * Y(I2)
00647      &                  - XA(IELEM,12) * Y(I3)
00648      &                  - XA(IELEM,14) * Y(I4)
00649      &                  - XA(IELEM,15) * Y(I5)
00650 !
00651           ENDDO ! IELEM
00652 !
00653         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
00654 !
00655           DO IELEM = 1 , NELEM
00656 !
00657           I1 = IKLE1(IELEM)
00658           I2 = IKLE2(IELEM)
00659           I3 = IKLE3(IELEM)
00660           I4 = IKLE4(IELEM)
00661           I5 = IKLE5(IELEM)
00662           I6 = IKLE6(IELEM)
00663 !
00664           W1(IELEM) = W1(IELEM)
00665      &                  - XA(IELEM, 1) * Y(I2)
00666      &                  - XA(IELEM, 2) * Y(I3)
00667      &                  - XA(IELEM, 3) * Y(I4)
00668      &                  - XA(IELEM, 4) * Y(I5)
00669      &                  - XA(IELEM, 5) * Y(I6)
00670 !
00671           W2(IELEM) = W2(IELEM)
00672      &                  - XA(IELEM,16) * Y(I1)
00673      &                  - XA(IELEM, 6) * Y(I3)
00674      &                  - XA(IELEM, 7) * Y(I4)
00675      &                  - XA(IELEM, 8) * Y(I5)
00676      &                  - XA(IELEM, 9) * Y(I6)
00677 !
00678           W3(IELEM) = W3(IELEM)
00679      &                  - XA(IELEM,17) * Y(I1)
00680      &                  - XA(IELEM,21) * Y(I2)
00681      &                  - XA(IELEM,10) * Y(I4)
00682      &                  - XA(IELEM,11) * Y(I5)
00683      &                  - XA(IELEM,12) * Y(I6)
00684 !
00685           W4(IELEM) = W4(IELEM)
00686      &                  - XA(IELEM,18) * Y(I1)
00687      &                  - XA(IELEM,22) * Y(I2)
00688      &                  - XA(IELEM,25) * Y(I3)
00689      &                  - XA(IELEM,13) * Y(I5)
00690      &                  - XA(IELEM,14) * Y(I6)
00691 !
00692           W5(IELEM) = W5(IELEM)
00693      &                  - XA(IELEM,19) * Y(I1)
00694      &                  - XA(IELEM,23) * Y(I2)
00695      &                  - XA(IELEM,26) * Y(I3)
00696      &                  - XA(IELEM,28) * Y(I4)
00697      &                  - XA(IELEM,15) * Y(I6)
00698 !
00699           W6(IELEM) = W6(IELEM)
00700      &                  - XA(IELEM,20) * Y(I1)
00701      &                  - XA(IELEM,24) * Y(I2)
00702      &                  - XA(IELEM,27) * Y(I3)
00703      &                  - XA(IELEM,29) * Y(I4)
00704      &                  - XA(IELEM,30) * Y(I5)
00705 !
00706           ENDDO ! IELEM
00707 !
00708         ELSEIF(TYPEXT(1:1).NE.'0') THEN
00709 !
00710           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
00711           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
00712           CALL PLANTE(0)
00713           STOP
00714 !
00715         ENDIF
00716 !
00717 !   CONTRIBUTION OF THE DIAGONAL:
00718 !
00719         IF(TYPDIA(1:1).EQ.'Q') THEN
00720           CALL OV ('X=X-YZ  ', X , Y , DA , C , NPOIN )
00721         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
00722           CALL OV ('X=X-Y   ', X , Y , Z  , C  , NPOIN )
00723         ELSEIF(TYPDIA(1:1).NE.'0') THEN
00724           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
00725           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
00726           CALL PLANTE(0)
00727           STOP
00728         ENDIF
00729 !
00730 !-----------------------------------------------------------------------
00731 !
00732       ELSEIF(OP(1:8).EQ.'X=X+CAY ') THEN
00733 !
00734 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
00735 !
00736         IF(TYPEXT(1:1).EQ.'S') THEN
00737 !
00738           DO IELEM = 1 , NELEM
00739 !
00740           I1 = IKLE1(IELEM)
00741           I2 = IKLE2(IELEM)
00742           I3 = IKLE3(IELEM)
00743           I4 = IKLE4(IELEM)
00744           I5 = IKLE5(IELEM)
00745           I6 = IKLE6(IELEM)
00746 !
00747           W1(IELEM) = W1(IELEM) + C * (
00748      &                  + XA(IELEM,1) * Y(I2)
00749      &                  + XA(IELEM,2) * Y(I3)
00750      &                  + XA(IELEM,3) * Y(I4)
00751      &                  + XA(IELEM,4) * Y(I5)
00752      &                  + XA(IELEM,5) * Y(I6)  )
00753 !
00754           W2(IELEM) = W2(IELEM) + C * (
00755      &                  + XA(IELEM,1) * Y(I1)
00756      &                  + XA(IELEM,6) * Y(I3)
00757      &                  + XA(IELEM,7) * Y(I4)
00758      &                  + XA(IELEM,8) * Y(I5)
00759      &                  + XA(IELEM,9) * Y(I6)  )
00760 !
00761           W3(IELEM) = W3(IELEM) + C * (
00762      &                  + XA(IELEM,2)  * Y(I1)
00763      &                  + XA(IELEM,6)  * Y(I2)
00764      &                  + XA(IELEM,10) * Y(I4)
00765      &                  + XA(IELEM,11) * Y(I5)
00766      &                  + XA(IELEM,12)*  Y(I6) )
00767 !
00768           W4(IELEM) = W4(IELEM) + C * (
00769      &                  + XA(IELEM,3)  * Y(I1)
00770      &                  + XA(IELEM,7)  * Y(I2)
00771      &                  + XA(IELEM,10) * Y(I3)
00772      &                  + XA(IELEM,13) * Y(I5)
00773      &                  + XA(IELEM,14) * Y(I6) )
00774 !
00775           W5(IELEM) = W5(IELEM) + C * (
00776      &                  + XA(IELEM,4)  * Y(I1)
00777      &                  + XA(IELEM,8)  * Y(I2)
00778      &                  + XA(IELEM,11) * Y(I3)
00779      &                  + XA(IELEM,13) * Y(I4)
00780      &                  + XA(IELEM,15) * Y(I6) )
00781 !
00782           W6(IELEM) = W6(IELEM) + C * (
00783      &                  + XA(IELEM,5)  * Y(I1)
00784      &                  + XA(IELEM,9)  * Y(I2)
00785      &                  + XA(IELEM,12) * Y(I3)
00786      &                  + XA(IELEM,14) * Y(I4)
00787      &                  + XA(IELEM,15) * Y(I5) )
00788 !
00789           ENDDO ! IELEM
00790 !
00791         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
00792 !
00793           DO IELEM = 1 , NELEM
00794 !
00795           I1 = IKLE1(IELEM)
00796           I2 = IKLE2(IELEM)
00797           I3 = IKLE3(IELEM)
00798           I4 = IKLE4(IELEM)
00799           I5 = IKLE5(IELEM)
00800           I6 = IKLE6(IELEM)
00801 !
00802           W1(IELEM) = W1(IELEM) + C * (
00803      &                  + XA(IELEM, 1) * Y(I2)
00804      &                  + XA(IELEM, 2) * Y(I3)
00805      &                  + XA(IELEM, 3) * Y(I4)
00806      &                  + XA(IELEM, 4) * Y(I5)
00807      &                  + XA(IELEM, 5) * Y(I6) )
00808 !
00809           W2(IELEM) = W2(IELEM) + C * (
00810      &                  + XA(IELEM,16) * Y(I1)
00811      &                  + XA(IELEM, 6) * Y(I3)
00812      &                  + XA(IELEM, 7) * Y(I4)
00813      &                  + XA(IELEM, 8) * Y(I5)
00814      &                  + XA(IELEM, 9) * Y(I6) )
00815 !
00816           W3(IELEM) = W3(IELEM) + C * (
00817      &                  + XA(IELEM,17) * Y(I1)
00818      &                  + XA(IELEM,21) * Y(I2)
00819      &                  + XA(IELEM,10) * Y(I4)
00820      &                  + XA(IELEM,11) * Y(I5)
00821      &                  + XA(IELEM,12) * Y(I6) )
00822 !
00823           W4(IELEM) = W4(IELEM) + C * (
00824      &                  + XA(IELEM,18) * Y(I1)
00825      &                  + XA(IELEM,22) * Y(I2)
00826      &                  + XA(IELEM,25) * Y(I3)
00827      &                  + XA(IELEM,13) * Y(I5)
00828      &                  + XA(IELEM,14) * Y(I6) )
00829 !
00830           W5(IELEM) = W5(IELEM) + C * (
00831      &                  + XA(IELEM,19) * Y(I1)
00832      &                  + XA(IELEM,23) * Y(I2)
00833      &                  + XA(IELEM,26) * Y(I3)
00834      &                  + XA(IELEM,28) * Y(I4)
00835      &                  + XA(IELEM,15) * Y(I6) )
00836 !
00837           W6(IELEM) = W6(IELEM) + C * (
00838      &                  + XA(IELEM,20) * Y(I1)
00839      &                  + XA(IELEM,24) * Y(I2)
00840      &                  + XA(IELEM,27) * Y(I3)
00841      &                  + XA(IELEM,29) * Y(I4)
00842      &                  + XA(IELEM,30) * Y(I5) )
00843 !
00844           ENDDO ! IELEM
00845 !
00846         ELSEIF(TYPEXT(1:1).NE.'0') THEN
00847 !
00848           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
00849           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
00850           CALL PLANTE(0)
00851           STOP
00852 !
00853         ENDIF
00854 !
00855 !   CONTRIBUTION OF THE DIAGONAL:
00856 !
00857         IF(TYPDIA(1:1).EQ.'Q') THEN
00858           CALL OV ('X=X+CYZ ', X , Y , DA , C , NPOIN )
00859         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
00860           CALL OV ('X=X+CY  ', X , Y , Z  , C  , NPOIN )
00861         ELSEIF(TYPDIA(1:1).NE.'0') THEN
00862           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
00863           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
00864           CALL PLANTE(0)
00865           STOP
00866         ENDIF
00867 !
00868 !-----------------------------------------------------------------------
00869 !
00870       ELSEIF(OP(1:8).EQ.'X=TAY   ') THEN
00871 !
00872 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
00873 !
00874         IF(TYPEXT(1:1).EQ.'S') THEN
00875 !
00876           DO IELEM = 1 , NELEM
00877 !
00878           I1 = IKLE1(IELEM)
00879           I2 = IKLE2(IELEM)
00880           I3 = IKLE3(IELEM)
00881           I4 = IKLE4(IELEM)
00882           I5 = IKLE5(IELEM)
00883           I6 = IKLE6(IELEM)
00884 !
00885           W1(IELEM) =
00886      &                  + XA(IELEM,1) * Y(I2)
00887      &                  + XA(IELEM,2) * Y(I3)
00888      &                  + XA(IELEM,3) * Y(I4)
00889      &                  + XA(IELEM,4) * Y(I5)
00890      &                  + XA(IELEM,5) * Y(I6)
00891 !
00892           W2(IELEM) =
00893      &                  + XA(IELEM,1) * Y(I1)
00894      &                  + XA(IELEM,6) * Y(I3)
00895      &                  + XA(IELEM,7) * Y(I4)
00896      &                  + XA(IELEM,8) * Y(I5)
00897      &                  + XA(IELEM,9) * Y(I6)
00898 !
00899           W3(IELEM) =
00900      &                  + XA(IELEM,2)  * Y(I1)
00901      &                  + XA(IELEM,6)  * Y(I2)
00902      &                  + XA(IELEM,10) * Y(I4)
00903      &                  + XA(IELEM,11) * Y(I5)
00904      &                  + XA(IELEM,12) * Y(I6)
00905 !
00906           W4(IELEM) =
00907      &                  + XA(IELEM,3)  * Y(I1)
00908      &                  + XA(IELEM,7)  * Y(I2)
00909      &                  + XA(IELEM,10) * Y(I3)
00910      &                  + XA(IELEM,13) * Y(I5)
00911      &                  + XA(IELEM,14) * Y(I6)
00912 !
00913           W5(IELEM) =
00914      &                  + XA(IELEM,4)  * Y(I1)
00915      &                  + XA(IELEM,8)  * Y(I2)
00916      &                  + XA(IELEM,11) * Y(I3)
00917      &                  + XA(IELEM,13) * Y(I4)
00918      &                  + XA(IELEM,15) * Y(I6)
00919 !
00920           W6(IELEM) =
00921      &                  + XA(IELEM,5)  * Y(I1)
00922      &                  + XA(IELEM,9)  * Y(I2)
00923      &                  + XA(IELEM,12) * Y(I3)
00924      &                  + XA(IELEM,14) * Y(I4)
00925      &                  + XA(IELEM,15) * Y(I5)
00926 !
00927           ENDDO ! IELEM
00928 !
00929         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
00930 !
00931           DO IELEM = 1 , NELEM
00932 !
00933           I1 = IKLE1(IELEM)
00934           I2 = IKLE2(IELEM)
00935           I3 = IKLE3(IELEM)
00936           I4 = IKLE4(IELEM)
00937           I5 = IKLE5(IELEM)
00938           I6 = IKLE6(IELEM)
00939 !
00940           W1(IELEM) =
00941      &                  + XA(IELEM,16) * Y(I2)
00942      &                  + XA(IELEM,17) * Y(I3)
00943      &                  + XA(IELEM,18) * Y(I4)
00944      &                  + XA(IELEM,19) * Y(I5)
00945      &                  + XA(IELEM,20) * Y(I6)
00946 !
00947           W2(IELEM) =
00948      &                  + XA(IELEM, 1) * Y(I1)
00949      &                  + XA(IELEM,21) * Y(I3)
00950      &                  + XA(IELEM,22) * Y(I4)
00951      &                  + XA(IELEM,23) * Y(I5)
00952      &                  + XA(IELEM,24) * Y(I6)
00953 !
00954           W3(IELEM) =
00955      &                  + XA(IELEM, 2) * Y(I1)
00956      &                  + XA(IELEM, 6) * Y(I2)
00957      &                  + XA(IELEM,25) * Y(I4)
00958      &                  + XA(IELEM,26) * Y(I5)
00959      &                  + XA(IELEM,27) * Y(I6)
00960 !
00961           W4(IELEM) =
00962      &                  + XA(IELEM, 3) * Y(I1)
00963      &                  + XA(IELEM, 7) * Y(I2)
00964      &                  + XA(IELEM,10) * Y(I3)
00965      &                  + XA(IELEM,28) * Y(I5)
00966      &                  + XA(IELEM,29) * Y(I6)
00967 !
00968           W5(IELEM) =
00969      &                  + XA(IELEM, 4) * Y(I1)
00970      &                  + XA(IELEM, 8) * Y(I2)
00971      &                  + XA(IELEM,11) * Y(I3)
00972      &                  + XA(IELEM,13) * Y(I4)
00973      &                  + XA(IELEM,30) * Y(I6)
00974 !
00975           W6(IELEM) =
00976      &                  + XA(IELEM, 5) * Y(I1)
00977      &                  + XA(IELEM, 9) * Y(I2)
00978      &                  + XA(IELEM,12) * Y(I3)
00979      &                  + XA(IELEM,14) * Y(I4)
00980      &                  + XA(IELEM,15) * Y(I5)
00981 !
00982           ENDDO ! IELEM
00983 !
00984         ELSEIF(TYPEXT(1:1).EQ.'0') THEN
00985 !
00986           CALL OV ('X=C     ', W1 , Y , Z , 0.D0 , NELEM )
00987           CALL OV ('X=C     ', W2 , Y , Z , 0.D0 , NELEM )
00988           CALL OV ('X=C     ', W3 , Y , Z , 0.D0 , NELEM )
00989           CALL OV ('X=C     ', W4 , Y , Z , 0.D0 , NELEM )
00990           CALL OV ('X=C     ', W5 , Y , Z , 0.D0 , NELEM )
00991           CALL OV ('X=C     ', W6 , Y , Z , 0.D0 , NELEM )
00992 !
00993         ELSE
00994 !
00995           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
00996           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
00997           CALL PLANTE(0)
00998           STOP
00999 !
01000         ENDIF
01001 !
01002 !   CONTRIBUTION OF THE DIAGONAL
01003 !
01004         IF(TYPDIA(1:1).EQ.'Q') THEN
01005           CALL OV ('X=YZ    ', X , Y , DA , C  , NPOIN )
01006         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
01007           CALL OV ('X=Y     ', X , Y , Z  , C  , NPOIN )
01008         ELSEIF(TYPDIA(1:1).EQ.'0') THEN
01009           CALL OV ('X=C     ', X , Y , DA , 0.D0 , NPOIN )
01010         ELSE
01011           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
01012           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
01013           CALL PLANTE(0)
01014           STOP
01015         ENDIF
01016 !
01017 !-----------------------------------------------------------------------
01018 !
01019       ELSEIF(OP(1:8).EQ.'X=-TAY  ') THEN
01020 !
01021 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
01022 !
01023         IF(TYPEXT(1:1).EQ.'S') THEN
01024 !
01025           DO IELEM = 1 , NELEM
01026 !
01027           I1 = IKLE1(IELEM)
01028           I2 = IKLE2(IELEM)
01029           I3 = IKLE3(IELEM)
01030           I4 = IKLE4(IELEM)
01031           I5 = IKLE5(IELEM)
01032           I6 = IKLE6(IELEM)
01033 !
01034           W1(IELEM) =
01035      &                  - XA(IELEM,1) * Y(I2)
01036      &                  - XA(IELEM,2) * Y(I3)
01037      &                  - XA(IELEM,3) * Y(I4)
01038      &                  - XA(IELEM,4) * Y(I5)
01039      &                  - XA(IELEM,5) * Y(I6)
01040 !
01041           W2(IELEM) =
01042      &                  - XA(IELEM,1) * Y(I1)
01043      &                  - XA(IELEM,6) * Y(I3)
01044      &                  - XA(IELEM,7) * Y(I4)
01045      &                  - XA(IELEM,8) * Y(I5)
01046      &                  - XA(IELEM,9) * Y(I6)
01047 !
01048           W3(IELEM) =
01049      &                  - XA(IELEM,2)  * Y(I1)
01050      &                  - XA(IELEM,6)  * Y(I2)
01051      &                  - XA(IELEM,10) * Y(I4)
01052      &                  - XA(IELEM,11) * Y(I5)
01053      &                  - XA(IELEM,12) * Y(I6)
01054 !
01055           W4(IELEM) =
01056      &                  - XA(IELEM,3)  * Y(I1)
01057      &                  - XA(IELEM,7)  * Y(I2)
01058      &                  - XA(IELEM,10) * Y(I3)
01059      &                  - XA(IELEM,13) * Y(I5)
01060      &                  - XA(IELEM,14) * Y(I6)
01061 !
01062           W5(IELEM) =
01063      &                  - XA(IELEM,4)  * Y(I1)
01064      &                  - XA(IELEM,8)  * Y(I2)
01065      &                  - XA(IELEM,11) * Y(I3)
01066      &                  - XA(IELEM,13) * Y(I4)
01067      &                  - XA(IELEM,15) * Y(I6)
01068 !
01069           W6(IELEM) =
01070      &                  - XA(IELEM,5)  * Y(I1)
01071      &                  - XA(IELEM,9)  * Y(I2)
01072      &                  - XA(IELEM,12) * Y(I3)
01073      &                  - XA(IELEM,14) * Y(I4)
01074      &                  - XA(IELEM,15) * Y(I5)
01075 !
01076           ENDDO ! IELEM
01077 !
01078         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
01079 !
01080           DO IELEM = 1 , NELEM
01081 !
01082           I1 = IKLE1(IELEM)
01083           I2 = IKLE2(IELEM)
01084           I3 = IKLE3(IELEM)
01085           I4 = IKLE4(IELEM)
01086           I5 = IKLE5(IELEM)
01087           I6 = IKLE6(IELEM)
01088 !
01089           W1(IELEM) =
01090      &                  - XA(IELEM,16) * Y(I2)
01091      &                  - XA(IELEM,17) * Y(I3)
01092      &                  - XA(IELEM,18) * Y(I4)
01093      &                  - XA(IELEM,19) * Y(I5)
01094      &                  - XA(IELEM,20) * Y(I6)
01095 !
01096           W2(IELEM) =
01097      &                  - XA(IELEM, 1) * Y(I1)
01098      &                  - XA(IELEM,21) * Y(I3)
01099      &                  - XA(IELEM,22) * Y(I4)
01100      &                  - XA(IELEM,23) * Y(I5)
01101      &                  - XA(IELEM,24) * Y(I6)
01102 !
01103           W3(IELEM) =
01104      &                  - XA(IELEM, 2) * Y(I1)
01105      &                  - XA(IELEM, 6) * Y(I2)
01106      &                  - XA(IELEM,25) * Y(I4)
01107      &                  - XA(IELEM,26) * Y(I5)
01108      &                  - XA(IELEM,27) * Y(I6)
01109 !
01110           W4(IELEM) =
01111      &                  - XA(IELEM, 3) * Y(I1)
01112      &                  - XA(IELEM, 7) * Y(I2)
01113      &                  - XA(IELEM,10) * Y(I3)
01114      &                  - XA(IELEM,28) * Y(I5)
01115      &                  - XA(IELEM,29) * Y(I6)
01116 !
01117           W5(IELEM) =
01118      &                  - XA(IELEM, 4) * Y(I1)
01119      &                  - XA(IELEM, 8) * Y(I2)
01120      &                  - XA(IELEM,11) * Y(I3)
01121      &                  - XA(IELEM,13) * Y(I4)
01122      &                  - XA(IELEM,30) * Y(I6)
01123 !
01124           W6(IELEM) =
01125      &                  - XA(IELEM, 5) * Y(I1)
01126      &                  - XA(IELEM, 9) * Y(I2)
01127      &                  - XA(IELEM,12) * Y(I3)
01128      &                  - XA(IELEM,14) * Y(I4)
01129      &                  - XA(IELEM,15) * Y(I5)
01130 !
01131           ENDDO ! IELEM
01132 !
01133         ELSEIF(TYPEXT(1:1).EQ.'0') THEN
01134 !
01135           CALL OV ('X=C     ', W1 , Y , Z , 0.D0 , NELEM )
01136           CALL OV ('X=C     ', W2 , Y , Z , 0.D0 , NELEM )
01137           CALL OV ('X=C     ', W3 , Y , Z , 0.D0 , NELEM )
01138           CALL OV ('X=C     ', W4 , Y , Z , 0.D0 , NELEM )
01139           CALL OV ('X=C     ', W5 , Y , Z , 0.D0 , NELEM )
01140           CALL OV ('X=C     ', W6 , Y , Z , 0.D0 , NELEM )
01141 !
01142         ELSE
01143 !
01144           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
01145           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
01146           CALL PLANTE(0)
01147           STOP
01148 !
01149         ENDIF
01150 !
01151 !   CONTRIBUTION OF THE DIAGONAL
01152 !
01153         IF(TYPDIA(1:1).EQ.'Q') THEN
01154           CALL OV ('X=-YZ   ', X , Y , DA , C  , NPOIN )
01155         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
01156           CALL OV ('X=-Y    ', X , Y , Z  , C  , NPOIN )
01157         ELSEIF(TYPDIA(1:1).EQ.'0') THEN
01158           CALL OV ('X=C     ', X , Y , DA , 0.D0 , NPOIN )
01159         ELSE
01160           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
01161           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
01162           CALL PLANTE(0)
01163           STOP
01164         ENDIF
01165 !
01166 !-----------------------------------------------------------------------
01167 !
01168       ELSEIF(OP(1:8).EQ.'X=X+TAY ') THEN
01169 !
01170 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
01171 !
01172         IF(TYPEXT(1:1).EQ.'S') THEN
01173 !
01174           DO IELEM = 1 , NELEM
01175 !
01176           I1 = IKLE1(IELEM)
01177           I2 = IKLE2(IELEM)
01178           I3 = IKLE3(IELEM)
01179           I4 = IKLE4(IELEM)
01180           I5 = IKLE5(IELEM)
01181           I6 = IKLE6(IELEM)
01182 !
01183           W1(IELEM) = W1(IELEM)
01184      &                  + XA(IELEM,1) * Y(I2)
01185      &                  + XA(IELEM,2) * Y(I3)
01186      &                  + XA(IELEM,3) * Y(I4)
01187      &                  + XA(IELEM,4) * Y(I5)
01188      &                  + XA(IELEM,5) * Y(I6)
01189 !
01190           W2(IELEM) = W2(IELEM)
01191      &                  + XA(IELEM,1) * Y(I1)
01192      &                  + XA(IELEM,6) * Y(I3)
01193      &                  + XA(IELEM,7) * Y(I4)
01194      &                  + XA(IELEM,8) * Y(I5)
01195      &                  + XA(IELEM,9) * Y(I6)
01196 !
01197           W3(IELEM) = W3(IELEM)
01198      &                  + XA(IELEM,2)  * Y(I1)
01199      &                  + XA(IELEM,6)  * Y(I2)
01200      &                  + XA(IELEM,10) * Y(I4)
01201      &                  + XA(IELEM,11) * Y(I5)
01202      &                  + XA(IELEM,12) * Y(I6)
01203 !
01204           W4(IELEM) = W4(IELEM)
01205      &                  + XA(IELEM,3)  * Y(I1)
01206      &                  + XA(IELEM,7)  * Y(I2)
01207      &                  + XA(IELEM,10) * Y(I3)
01208      &                  + XA(IELEM,13) * Y(I5)
01209      &                  + XA(IELEM,14) * Y(I6)
01210 !
01211           W5(IELEM) = W5(IELEM)
01212      &                  + XA(IELEM,4)  * Y(I1)
01213      &                  + XA(IELEM,8)  * Y(I2)
01214      &                  + XA(IELEM,11) * Y(I3)
01215      &                  + XA(IELEM,13) * Y(I4)
01216      &                  + XA(IELEM,15) * Y(I6)
01217 !
01218           W6(IELEM) = W6(IELEM)
01219      &                  + XA(IELEM,5)  * Y(I1)
01220      &                  + XA(IELEM,9)  * Y(I2)
01221      &                  + XA(IELEM,12) * Y(I3)
01222      &                  + XA(IELEM,14) * Y(I4)
01223      &                  + XA(IELEM,15) * Y(I5)
01224 !
01225         ENDDO ! IELEM
01226 !
01227         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
01228 !
01229           DO IELEM = 1 , NELEM
01230 !
01231           I1 = IKLE1(IELEM)
01232           I2 = IKLE2(IELEM)
01233           I3 = IKLE3(IELEM)
01234           I4 = IKLE4(IELEM)
01235           I5 = IKLE5(IELEM)
01236           I6 = IKLE6(IELEM)
01237 !
01238           W1(IELEM) = W1(IELEM)
01239      &                  + XA(IELEM,16) * Y(I2)
01240      &                  + XA(IELEM,17) * Y(I3)
01241      &                  + XA(IELEM,18) * Y(I4)
01242      &                  + XA(IELEM,19) * Y(I5)
01243      &                  + XA(IELEM,20) * Y(I6)
01244 !
01245           W2(IELEM) = W2(IELEM)
01246      &                  + XA(IELEM, 1) * Y(I1)
01247      &                  + XA(IELEM,21) * Y(I3)
01248      &                  + XA(IELEM,22) * Y(I4)
01249      &                  + XA(IELEM,23) * Y(I5)
01250      &                  + XA(IELEM,24) * Y(I6)
01251 !
01252           W3(IELEM) = W3(IELEM)
01253      &                  + XA(IELEM, 2) * Y(I1)
01254      &                  + XA(IELEM, 6) * Y(I2)
01255      &                  + XA(IELEM,25) * Y(I4)
01256      &                  + XA(IELEM,26) * Y(I5)
01257      &                  + XA(IELEM,27) * Y(I6)
01258 !
01259           W4(IELEM) = W4(IELEM)
01260      &                  + XA(IELEM, 3) * Y(I1)
01261      &                  + XA(IELEM, 7) * Y(I2)
01262      &                  + XA(IELEM,10) * Y(I3)
01263      &                  + XA(IELEM,28) * Y(I5)
01264      &                  + XA(IELEM,29) * Y(I6)
01265 !
01266           W5(IELEM) = W5(IELEM)
01267      &                  + XA(IELEM, 4) * Y(I1)
01268      &                  + XA(IELEM, 8) * Y(I2)
01269      &                  + XA(IELEM,11) * Y(I3)
01270      &                  + XA(IELEM,13) * Y(I4)
01271      &                  + XA(IELEM,30) * Y(I6)
01272 !
01273           W6(IELEM) = W6(IELEM)
01274      &                  + XA(IELEM, 5) * Y(I1)
01275      &                  + XA(IELEM, 9) * Y(I2)
01276      &                  + XA(IELEM,12) * Y(I3)
01277      &                  + XA(IELEM,14) * Y(I4)
01278      &                  + XA(IELEM,15) * Y(I5)
01279 !
01280           ENDDO ! IELEM
01281 !
01282         ELSEIF(TYPEXT(1:1).NE.'0') THEN
01283 !
01284           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
01285           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
01286           CALL PLANTE(0)
01287           STOP
01288 !
01289         ENDIF
01290 !
01291 !   CONTRIBUTION OF THE DIAGONAL
01292 !
01293         IF(TYPDIA(1:1).EQ.'Q') THEN
01294           CALL OV ('X=X+YZ  ', X , Y , DA , C , NPOIN )
01295         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
01296           CALL OV ('X=X+Y   ', X , Y , Z  , C  , NPOIN )
01297         ELSEIF(TYPDIA(1:1).NE.'0') THEN
01298           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
01299           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
01300           CALL PLANTE(0)
01301           STOP
01302         ENDIF
01303 !
01304 !-----------------------------------------------------------------------
01305 !
01306       ELSEIF(OP(1:8).EQ.'X=X-TAY ') THEN
01307 !
01308 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
01309 !
01310         IF(TYPEXT(1:1).EQ.'S') THEN
01311 !
01312           DO IELEM = 1 , NELEM
01313 !
01314           I1 = IKLE1(IELEM)
01315           I2 = IKLE2(IELEM)
01316           I3 = IKLE3(IELEM)
01317           I4 = IKLE4(IELEM)
01318           I5 = IKLE5(IELEM)
01319           I6 = IKLE6(IELEM)
01320 !
01321           W1(IELEM) = W1(IELEM)
01322      &                  - XA(IELEM,1) * Y(I2)
01323      &                  - XA(IELEM,2) * Y(I3)
01324      &                  - XA(IELEM,3) * Y(I4)
01325      &                  - XA(IELEM,4) * Y(I5)
01326      &                  - XA(IELEM,5) * Y(I6)
01327 !
01328           W2(IELEM) = W2(IELEM)
01329      &                  - XA(IELEM,1) * Y(I1)
01330      &                  - XA(IELEM,6) * Y(I3)
01331      &                  - XA(IELEM,7) * Y(I4)
01332      &                  - XA(IELEM,8) * Y(I5)
01333      &                  - XA(IELEM,9) * Y(I6)
01334 !
01335           W3(IELEM) = W3(IELEM)
01336      &                  - XA(IELEM,2)  * Y(I1)
01337      &                  - XA(IELEM,6)  * Y(I2)
01338      &                  - XA(IELEM,10) * Y(I4)
01339      &                  - XA(IELEM,11) * Y(I5)
01340      &                  - XA(IELEM,12) * Y(I6)
01341 !
01342           W4(IELEM) = W4(IELEM)
01343      &                  - XA(IELEM,3)  * Y(I1)
01344      &                  - XA(IELEM,7)  * Y(I2)
01345      &                  - XA(IELEM,10) * Y(I3)
01346      &                  - XA(IELEM,13) * Y(I5)
01347      &                  - XA(IELEM,14) * Y(I6)
01348 !
01349           W5(IELEM) = W5(IELEM)
01350      &                  - XA(IELEM,4)  * Y(I1)
01351      &                  - XA(IELEM,8)  * Y(I2)
01352      &                  - XA(IELEM,11) * Y(I3)
01353      &                  - XA(IELEM,13) * Y(I4)
01354      &                  - XA(IELEM,15) * Y(I6)
01355 !
01356           W6(IELEM) = W6(IELEM)
01357      &                  - XA(IELEM,5)  * Y(I1)
01358      &                  - XA(IELEM,9)  * Y(I2)
01359      &                  - XA(IELEM,12) * Y(I3)
01360      &                  - XA(IELEM,14) * Y(I4)
01361      &                  - XA(IELEM,15) * Y(I5)
01362 !
01363         ENDDO ! IELEM
01364 !
01365         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
01366 !
01367           DO IELEM = 1 , NELEM
01368 !
01369           I1 = IKLE1(IELEM)
01370           I2 = IKLE2(IELEM)
01371           I3 = IKLE3(IELEM)
01372           I4 = IKLE4(IELEM)
01373           I5 = IKLE5(IELEM)
01374           I6 = IKLE6(IELEM)
01375 !
01376           W1(IELEM) = W1(IELEM)
01377      &                  - XA(IELEM,16) * Y(I2)
01378      &                  - XA(IELEM,17) * Y(I3)
01379      &                  - XA(IELEM,18) * Y(I4)
01380      &                  - XA(IELEM,19) * Y(I5)
01381      &                  - XA(IELEM,20) * Y(I6)
01382 !
01383           W2(IELEM) = W2(IELEM)
01384      &                  - XA(IELEM, 1) * Y(I1)
01385      &                  - XA(IELEM,21) * Y(I3)
01386      &                  - XA(IELEM,22) * Y(I4)
01387      &                  - XA(IELEM,23) * Y(I5)
01388      &                  - XA(IELEM,24) * Y(I6)
01389 !
01390           W3(IELEM) = W3(IELEM)
01391      &                  - XA(IELEM, 2) * Y(I1)
01392      &                  - XA(IELEM, 6) * Y(I2)
01393      &                  - XA(IELEM,25) * Y(I4)
01394      &                  - XA(IELEM,26) * Y(I5)
01395      &                  - XA(IELEM,27) * Y(I6)
01396 !
01397           W4(IELEM) = W4(IELEM)
01398      &                  - XA(IELEM, 3) * Y(I1)
01399      &                  - XA(IELEM, 7) * Y(I2)
01400      &                  - XA(IELEM,10) * Y(I3)
01401      &                  - XA(IELEM,28) * Y(I5)
01402      &                  - XA(IELEM,29) * Y(I6)
01403 !
01404           W5(IELEM) = W5(IELEM)
01405      &                  - XA(IELEM, 4) * Y(I1)
01406      &                  - XA(IELEM, 8) * Y(I2)
01407      &                  - XA(IELEM,11) * Y(I3)
01408      &                  - XA(IELEM,13) * Y(I4)
01409      &                  - XA(IELEM,30) * Y(I6)
01410 !
01411           W6(IELEM) = W6(IELEM)
01412      &                  - XA(IELEM, 5) * Y(I1)
01413      &                  - XA(IELEM, 9) * Y(I2)
01414      &                  - XA(IELEM,12) * Y(I3)
01415      &                  - XA(IELEM,14) * Y(I4)
01416      &                  - XA(IELEM,15) * Y(I5)
01417 !
01418           ENDDO ! IELEM
01419 !
01420         ELSEIF(TYPEXT(1:1).NE.'0') THEN
01421 !
01422           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
01423           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
01424           CALL PLANTE(0)
01425           STOP
01426 !
01427         ENDIF
01428 !
01429 !   CONTRIBUTION OF THE DIAGONAL
01430 !
01431         IF(TYPDIA(1:1).EQ.'Q') THEN
01432           CALL OV ('X=X-YZ  ', X , Y , DA , C , NPOIN )
01433         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
01434           CALL OV ('X=X-Y   ', X , Y , Z  , C  , NPOIN )
01435         ELSEIF(TYPDIA(1:1).NE.'0') THEN
01436           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
01437           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
01438           CALL PLANTE(0)
01439           STOP
01440         ENDIF
01441 !
01442 !-----------------------------------------------------------------------
01443 !
01444       ELSEIF(OP(1:8).EQ.'X=X+CTAY') THEN
01445 !
01446 !   CONTRIBUTION OF EXTRADIAGONAL TERMS:
01447 !
01448         IF(TYPEXT(1:1).EQ.'S') THEN
01449 !
01450           DO IELEM = 1 , NELEM
01451 !
01452           I1 = IKLE1(IELEM)
01453           I2 = IKLE2(IELEM)
01454           I3 = IKLE3(IELEM)
01455           I4 = IKLE4(IELEM)
01456           I5 = IKLE5(IELEM)
01457           I6 = IKLE6(IELEM)
01458 !
01459           W1(IELEM) = W1(IELEM) + C * (
01460      &                  + XA(IELEM,1) * Y(I2)
01461      &                  + XA(IELEM,2) * Y(I3)
01462      &                  + XA(IELEM,3) * Y(I4)
01463      &                  + XA(IELEM,4) * Y(I5)
01464      &                  + XA(IELEM,5) * Y(I6)  )
01465 !
01466           W2(IELEM) = W2(IELEM) + C * (
01467      &                  + XA(IELEM,1) * Y(I1)
01468      &                  + XA(IELEM,6) * Y(I3)
01469      &                  + XA(IELEM,7) * Y(I4)
01470      &                  + XA(IELEM,8) * Y(I5)
01471      &                  + XA(IELEM,9) * Y(I6)  )
01472 !
01473           W3(IELEM) = W3(IELEM) + C * (
01474      &                  + XA(IELEM,2)  * Y(I1)
01475      &                  + XA(IELEM,6)  * Y(I2)
01476      &                  + XA(IELEM,10) * Y(I4)
01477      &                  + XA(IELEM,11) * Y(I5)
01478      &                  + XA(IELEM,12)*  Y(I6) )
01479 !
01480           W4(IELEM) = W4(IELEM) + C * (
01481      &                  + XA(IELEM,3)  * Y(I1)
01482      &                  + XA(IELEM,7)  * Y(I2)
01483      &                  + XA(IELEM,10) * Y(I3)
01484      &                  + XA(IELEM,13) * Y(I5)
01485      &                  + XA(IELEM,14) * Y(I6) )
01486 !
01487           W5(IELEM) = W5(IELEM) + C * (
01488      &                  + XA(IELEM,4)  * Y(I1)
01489      &                  + XA(IELEM,8)  * Y(I2)
01490      &                  + XA(IELEM,11) * Y(I3)
01491      &                  + XA(IELEM,13) * Y(I4)
01492      &                  + XA(IELEM,15) * Y(I6) )
01493 !
01494           W6(IELEM) = W6(IELEM) + C * (
01495      &                  + XA(IELEM,5)  * Y(I1)
01496      &                  + XA(IELEM,9)  * Y(I2)
01497      &                  + XA(IELEM,12) * Y(I3)
01498      &                  + XA(IELEM,14) * Y(I4)
01499      &                  + XA(IELEM,15) * Y(I5) )
01500 !
01501         ENDDO ! IELEM
01502 !
01503         ELSEIF(TYPEXT(1:1).EQ.'Q') THEN
01504 !
01505           DO IELEM = 1 , NELEM
01506 !
01507           I1 = IKLE1(IELEM)
01508           I2 = IKLE2(IELEM)
01509           I3 = IKLE3(IELEM)
01510           I4 = IKLE4(IELEM)
01511           I5 = IKLE5(IELEM)
01512           I6 = IKLE6(IELEM)
01513 !
01514           W1(IELEM) = W1(IELEM) + C * (
01515      &                  + XA(IELEM,16) * Y(I2)
01516      &                  + XA(IELEM,17) * Y(I3)
01517      &                  + XA(IELEM,18) * Y(I4)
01518      &                  + XA(IELEM,19) * Y(I5)
01519      &                  + XA(IELEM,20) * Y(I6) )
01520 !
01521           W2(IELEM) = W2(IELEM) + C * (
01522      &                  + XA(IELEM, 1) * Y(I1)
01523      &                  + XA(IELEM,21) * Y(I3)
01524      &                  + XA(IELEM,22) * Y(I4)
01525      &                  + XA(IELEM,23) * Y(I5)
01526      &                  + XA(IELEM,24) * Y(I6) )
01527 !
01528           W3(IELEM) = W3(IELEM) + C * (
01529      &                  + XA(IELEM, 2) * Y(I1)
01530      &                  + XA(IELEM, 6) * Y(I2)
01531      &                  + XA(IELEM,25) * Y(I4)
01532      &                  + XA(IELEM,26) * Y(I5)
01533      &                  + XA(IELEM,27) * Y(I6) )
01534 !
01535           W4(IELEM) = W4(IELEM) + C * (
01536      &                  + XA(IELEM, 3) * Y(I1)
01537      &                  + XA(IELEM, 7) * Y(I2)
01538      &                  + XA(IELEM,10) * Y(I3)
01539      &                  + XA(IELEM,28) * Y(I5)
01540      &                  + XA(IELEM,29) * Y(I6) )
01541 !
01542           W5(IELEM) = W5(IELEM) + C * (
01543      &                  + XA(IELEM, 4) * Y(I1)
01544      &                  + XA(IELEM, 8) * Y(I2)
01545      &                  + XA(IELEM,11) * Y(I3)
01546      &                  + XA(IELEM,13) * Y(I4)
01547      &                  + XA(IELEM,30) * Y(I6) )
01548 !
01549           W6(IELEM) = W6(IELEM) + C * (
01550      &                  + XA(IELEM, 5) * Y(I1)
01551      &                  + XA(IELEM, 9) * Y(I2)
01552      &                  + XA(IELEM,12) * Y(I3)
01553      &                  + XA(IELEM,14) * Y(I4)
01554      &                  + XA(IELEM,15) * Y(I5) )
01555 !
01556           ENDDO ! IELEM
01557 !
01558         ELSEIF(TYPEXT(1:1).NE.'0') THEN
01559 !
01560           IF (LNG.EQ.1) WRITE(LU,1000) TYPEXT
01561           IF (LNG.EQ.2) WRITE(LU,1001) TYPEXT
01562           CALL PLANTE(0)
01563           STOP
01564 !
01565         ENDIF
01566 !
01567 !   CONTRIBUTION OF THE DIAGONAL
01568 !
01569         IF(TYPDIA(1:1).EQ.'Q') THEN
01570           CALL OV ('X=X+CYZ ', X , Y , DA , C , NPOIN )
01571         ELSEIF(TYPDIA(1:1).EQ.'I') THEN
01572           CALL OV ('X=X+CY  ', X , Y , Z  , C  , NPOIN )
01573         ELSEIF(TYPDIA(1:1).NE.'0') THEN
01574           IF (LNG.EQ.1) WRITE(LU,2000) TYPDIA
01575           IF (LNG.EQ.2) WRITE(LU,2001) TYPDIA
01576           CALL PLANTE(0)
01577           STOP
01578         ENDIF
01579 !
01580 !-----------------------------------------------------------------------
01581 !
01582       ELSE
01583 !
01584         IF (LNG.EQ.1) WRITE(LU,3000) OP
01585         IF (LNG.EQ.2) WRITE(LU,3001) OP
01586         CALL PLANTE(0)
01587 !
01588 !-----------------------------------------------------------------------
01589 !
01590       ENDIF
01591 !
01592 !-----------------------------------------------------------------------
01593 !
01594       RETURN
01595 !
01596 1000  FORMAT(1X,'MV0606 (BIEF) : TERMES EXTRADIAG. TYPE INCONNU: ',A1)
01597 1001  FORMAT(1X,'MV0606 (BIEF) : EXTRADIAG. TERMS  UNKNOWN TYPE : ',A1)
01598 2000  FORMAT(1X,'MV0606 (BIEF) : DIAGONALE : TYPE INCONNU: ',A1)
01599 2001  FORMAT(1X,'MV0606 (BIEF) : DIAGONAL : UNKNOWN TYPE : ',A1)
01600 3000  FORMAT(1X,'MV0606 (BIEF) : OPERATION INCONNUE : ',A8)
01601 3001  FORMAT(1X,'MV0606 (BIEF) : UNKNOWN OPERATION : ',A8)
01602 !
01603       END

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