mt11ac.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\mt11ac.f
00002 !
00075                      SUBROUTINE MT11AC
00076 !                    *****************
00077 !
00078      &(  A11 , A12 , A13 , A14 , A15, A16,
00079      &   A21 , A22 , A23 , A24 , A25, A26,
00080      &   A31 , A32 , A33 , A34 , A35, A36,
00081      &   XMUL,SF,F,XEL,YEL,IKLE1,IKLE2,IKLE3,
00082      &   IKLE4,IKLE5,IKLE6,
00083      &   NELEM,NELMAX,ICOORD)
00084 !
00085 !***********************************************************************
00086 ! BIEF   V6P1                                   21/08/2010
00087 !***********************************************************************
00088 !
00089 !
00090 !
00091 !
00092 !
00093 !
00094 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00095 !| A11            |<--| ELEMENTS OF MATRIX
00096 !| ...            |<--| ELEMENTS OF MATRIX
00097 !| A36            |<--| ELEMENTS OF MATRIX
00098 !| F              |-->| FUNCTION USED IN THE FORMULA
00099 !| ICOORD         |-->| 1: DERIVATIVE ALONG X, 2: ALONG Y
00100 !| IKLE1          |-->| FIRST POINTS OF TRIANGLES
00101 !| IKLE2          |-->| SECOND POINTS OF TRIANGLES
00102 !| IKLE3          |-->| THIRD POINTS OF TRIANGLES
00103 !| IKLE4          |-->| FOURTH POINTS OF TRIANGLES (QUADRATIC)
00104 !| IKLE5          |-->| FIFTH POINTS OF TRIANGLES (QUADRATIC)
00105 !| IKLE6          |-->| SIXTH POINTS OF TRIANGLES (QUADRATIC)
00106 !| NELEM          |-->| NUMBER OF ELEMENTS
00107 !| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
00108 !| SF             |-->| STRUCTURE OF FUNCTIONS F
00109 !| XEL            |-->| ABSCISSAE OF POINTS IN THE MESH, PER ELEMENT
00110 !| YEL            |-->| ORDINATES OF POINTS IN THE MESH, PER ELEMENT
00111 !| XMUL           |-->| MULTIPLICATION FACTOR
00112 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00113 !
00114       USE BIEF!, EX_MT11AC => MT11AC
00115 !
00116       IMPLICIT NONE
00117       INTEGER LNG,LU
00118       COMMON/INFO/LNG,LU
00119 !
00120 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00121 !
00122       INTEGER, INTENT(IN) :: NELEM,NELMAX,ICOORD
00123       INTEGER, INTENT(IN) :: IKLE1(NELMAX),IKLE2(NELMAX)
00124       INTEGER, INTENT(IN) :: IKLE3(NELMAX),IKLE4(NELMAX)
00125       INTEGER, INTENT(IN) :: IKLE5(NELMAX),IKLE6(NELMAX)
00126 !
00127       DOUBLE PRECISION, INTENT(INOUT) :: A11(*),A12(*),A13(*)
00128       DOUBLE PRECISION, INTENT(INOUT) :: A14(*),A15(*),A16(*)
00129       DOUBLE PRECISION, INTENT(INOUT) :: A21(*),A22(*),A23(*)
00130       DOUBLE PRECISION, INTENT(INOUT) :: A24(*),A25(*),A26(*)
00131       DOUBLE PRECISION, INTENT(INOUT) :: A31(*),A32(*),A33(*)
00132       DOUBLE PRECISION, INTENT(INOUT) :: A34(*),A35(*),A36(*)
00133 !
00134       DOUBLE PRECISION, INTENT(IN) :: XMUL
00135       DOUBLE PRECISION, INTENT(IN) :: F(*)
00136 !
00137 !     STRUCTURE OF F
00138       TYPE(BIEF_OBJ), INTENT(IN) :: SF
00139 !
00140       DOUBLE PRECISION, INTENT(IN) :: XEL(NELMAX,3),YEL(NELMAX,3)
00141 !
00142 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00143 !
00144       INTEGER IELEM,IELMF
00145       DOUBLE PRECISION X2,X3,Y2,Y3,F1,F2,F3,F4,F5,F6
00146       DOUBLE PRECISION XSUR30,XSUR120
00147       DOUBLE PRECISION XSUR90,XSUR360
00148 !
00149 !-----------------------------------------------------------------------
00150 !
00151 !
00152 !-----------------------------------------------------------------------
00153 !
00154       XSUR30  = XMUL/30.D0
00155       XSUR120 = XMUL/120.D0
00156       XSUR90  = XMUL/90.D0
00157       XSUR360 = XMUL/360.D0
00158 !
00159       IELMF=SF%ELM
00160 !
00161 !-----------------------------------------------------------------------
00162 !  CASE WHERE F IS OF TYPE P1
00163 !-----------------------------------------------------------------------
00164 !
00165       IF(IELMF.EQ.11) THEN
00166 !
00167 !================================
00168 !  DERIVATIVE WRT X =
00169 !================================
00170 !
00171         IF(ICOORD.EQ.1) THEN
00172 !
00173 !   LOOP ON THE ELEMENTS
00174 !
00175         DO IELEM = 1 , NELEM
00176 !
00177 !   INITIALISES THE GEOMETRICAL VARIABLES
00178 !
00179         Y2 = YEL(IELEM,2)
00180         Y3 = YEL(IELEM,3)
00181 !
00182         F1  =  F(IKLE1(IELEM))
00183         F2  =  F(IKLE2(IELEM))
00184         F3  =  F(IKLE3(IELEM))
00185 !
00186 !   EXTRA-DIAGONAL TERMS
00187 !
00188         A12(IELEM) = ( 2.D0*(F1-F2)*Y2 + (3.D0*F2-2.D0*F1-F3)*Y3)
00189      &             * XSUR120
00190         A13(IELEM) = ( 2.D0*(F3-F1)*Y3 + (2.D0*F1-3.D0*F3+F2)*Y2)
00191      &             * XSUR120
00192         A14(IELEM) = (( 4.D0*F1+F3)*Y3 + (F3-4.D0*F1-2.D0*F2)*Y2)
00193      &             * XSUR30
00194         A15(IELEM) = ((-2.D0*(F1+F2)-F3)*Y2  +
00195      &                ( 2.D0*(F1+F3)+F2)*Y3) * XSUR30
00196         A16(IELEM) = ((-4.D0*F1-F2)*Y2 + (4.D0*F1-F2+2.D0*F3)*Y3)
00197      &             * XSUR30
00198         A21(IELEM) = (      (F1-F3)*Y2 + (F3-3.D0*F1+2.D0*F2)*Y3)
00199      &             * XSUR120
00200         A23(IELEM) = (      2.D0*(F2-F3)*Y3 +        (F1-F3) *Y2)
00201      &             * XSUR120
00202         A24(IELEM) = ( 2.D0*(F3-F1)*Y2 - (        F3+4.D0*F2)*Y3)
00203      &             * XSUR30
00204         A25(IELEM) = ( 2.D0*(F3-F1)*Y2 + (F1-2.D0*F3-4.D0*F2)*Y3)
00205      &             * XSUR30
00206         A26(IELEM) = (      (F3-F1)*Y2 - (F1+2.D0*F2+2.D0*F3)*Y3)
00207      &             * XSUR30
00208         A31(IELEM) = (      (F2-F1)*Y3 + (3.D0*F1-F2-2.D0*F3)*Y2)
00209      &             * XSUR120
00210         A32(IELEM) = ( 2.D0*(F2-F3)*Y2 +             (-F1+F2)*Y3)
00211      &             * XSUR120
00212         A34(IELEM) = (      (F1-F2)*Y3 + (   F1+2.D0*(F2+F3))*Y2)
00213      &             * XSUR30
00214         A35(IELEM) = ( 2.D0*(F1-F2)*Y3 + (4.D0*F3-F1+2.D0*F2)*Y2)
00215      &             * XSUR30
00216         A36(IELEM) = (( 4.D0*F3+F2)*Y2 +         2.D0*(F1-F2)*Y3)
00217      &             * XSUR30
00218 !
00219 !   DIAGONAL TERMS
00220 !
00221         A11(IELEM) = - A21(IELEM) - A31(IELEM)
00222         A22(IELEM) = - A12(IELEM) - A32(IELEM)
00223         A33(IELEM) = - A13(IELEM) - A23(IELEM)
00224 !
00225       ENDDO ! IELEM
00226 !
00227         ELSEIF(ICOORD.EQ.2) THEN
00228 !
00229 !================================
00230 !  DERIVATIVE WRT Y  =
00231 !================================
00232 !
00233         DO IELEM = 1 , NELEM
00234 !
00235 !   INITIALISES THE GEOMETRICAL VARIABLES
00236 !
00237         X2  =  XEL(IELEM,2)
00238         X3  =  XEL(IELEM,3)
00239 !
00240         F1  =  F(IKLE1(IELEM))
00241         F2  =  F(IKLE2(IELEM))
00242         F3  =  F(IKLE3(IELEM))
00243 !
00244 !   EXTRA-DIAGONAL TERMS
00245 !
00246         A12(IELEM) = (  2.D0*(F2-F1)*X2 + (2.D0*F1-3.D0*F2+F3)*X3)
00247      &             * XSUR120
00248         A13(IELEM) = (  2.D0*(F1-F3)*X3 - (2.D0*F1-3.D0*F3+F2)*X2)
00249      &             * XSUR120
00250         A14(IELEM) = ((-4.D0*F1-F3 )*X3 + (4.D0*F1-F3+2.D0*F2)*X2)
00251      &             * XSUR30
00252         A15(IELEM) = (( 2.D0*(F1+F2)+F3)*X2  -
00253      &                ( 2.D0*(F1+F3)+F2)*X3) * XSUR30
00254         A16(IELEM) = ((4.D0*F1+F2  )*X2 + (F2-4.D0*F1-2.D0*F3)*X3)
00255      &             * XSUR30
00256         A21(IELEM) = (       (F3-F1)*X2 + (3.D0*F1-2.D0*F2-F3)*X3)
00257      &             * XSUR120
00258         A23(IELEM) = (       (F3-F1)*X2 +         2.D0*(F3-F2)*X3)
00259      &             * XSUR120
00260         A24(IELEM) = (  2.D0*(F1-F3)*X2 + (        F3+4.D0*F2)*X3)
00261      &             * XSUR30
00262         A25(IELEM) = (  2.D0*(F1-F3)*X2 + (2.D0*F3-F1+4.D0*F2)*X3)
00263      &             * XSUR30
00264         A26(IELEM) = (       (F1-F3)*X2 + (   F1+2.D0*(F2+F3))*X3)
00265      &             * XSUR30
00266         A31(IELEM) = (       (F1-F2)*X3 + (F2+2.D0*F3-3.D0*F1)*X2)
00267      &             * XSUR120
00268         A32(IELEM) = (  2.D0*(F3-F2)*X2 + (             F1-F2)*X3)
00269      &             * XSUR120
00270         A34(IELEM) = (       (F2-F1)*X3 - (F1+2.D0*F2+2.D0*F3)*X2)
00271      &             * XSUR30
00272         A35(IELEM) = (  2.D0*(F2-F1)*X3 - (4.D0*F3-F1+2.D0*F2)*X2)
00273      &             * XSUR30
00274         A36(IELEM) = ((-4.D0*F3-F2 )*X2 +         2.D0*(F2-F1)*X3)
00275      &             * XSUR30
00276 !
00277 !   DIAGONAL TERMS
00278 !
00279         A11(IELEM) = - A21(IELEM) - A31(IELEM)
00280         A22(IELEM) = - A12(IELEM) - A32(IELEM)
00281         A33(IELEM) = - A13(IELEM) - A23(IELEM)
00282 !
00283         ENDDO ! IELEM
00284 !
00285         ELSE
00286 !
00287           IF (LNG.EQ.1) WRITE(LU,200) ICOORD
00288           IF (LNG.EQ.2) WRITE(LU,201) ICOORD
00289           CALL PLANTE(1)
00290           STOP
00291         ENDIF
00292 !
00293 !-----------------------------------------------------------------------
00294 !  CASE WHERE F IS OF TYPE P2
00295 !-----------------------------------------------------------------------
00296 !
00297       ELSEIF(IELMF.EQ.13) THEN
00298 !
00299 !================================
00300 !  DERIVATIVE WRT X  =
00301 !================================
00302 !
00303         IF(ICOORD.EQ.1) THEN
00304 !
00305 !   LOOP ON THE ELEMENTS
00306 !
00307           DO IELEM = 1 , NELEM
00308 !
00309 !   INITIALISES THE GEOMETRICAL VARIABLES
00310 !
00311             Y2 = YEL(IELEM,2)
00312             Y3 = YEL(IELEM,3)
00313 !
00314             F1  =  F(IKLE1(IELEM))
00315             F2  =  F(IKLE2(IELEM))
00316             F3  =  F(IKLE3(IELEM))
00317             F4  =  F(IKLE4(IELEM))
00318             F5  =  F(IKLE5(IELEM))
00319             F6  =  F(IKLE6(IELEM))
00320 !
00321 !   EXTRADIAGONAL TERMS
00322 !
00323             A12(IELEM) = ((3.D0*F2-6.D0*F1-8.D0*(F6-F4)-F3+4.D0*F5)*Y3
00324      &                 +   6.D0*(F1-F2)*Y2 ) * XSUR360
00325             A13(IELEM) = ((8.D0*(F4-F6)-4.D0*F5+F2+6.D0*F1-3.D0*F3)*Y2
00326      &                 +   6.D0*(F3-F1)*Y3 ) * XSUR360
00327             A14(IELEM) = ((-16.D0*F4+4.D0*(F5+F6)-6.D0*F1-F3      )*Y2
00328      &                 +  (6.D0*F1-2.D0*F2+4.D0*F4+8.D0*F6-F3)*Y3)
00329      &                  * XSUR90
00330             A15(IELEM) = ((F3-8.D0*F4-4.D0*(F5+F6))*Y2
00331      &                 +  (4.D0*(F4+F5)+8.D0*F6-F2)*Y3) * XSUR90
00332             A16(IELEM) = ((F2+6.D0*F1+16.D0*F6-4.D0*(F4+F5)       )*Y3
00333      &                 -  (8.D0*F4-F2+6.D0*F1-2.D0*F3+4.D0*F6)*Y2)
00334      &                  * XSUR90
00335             A21(IELEM) = ((8.D0*(F4-F5)-3.D0*F1-F3+4.D0*F6        )*Y2
00336      &                 - (3.D0*F1-6.D0*F2+8.D0*(F4-F5)+4.D0*F6-F3 )*Y3)
00337      &                 * XSUR360
00338             A23(IELEM) = ((8.D0*(F4-F5)+F1+3.D0*F3-4.D0*F6        )*Y2
00339      &                 +  6.D0*(F2-F3)*Y3 ) * XSUR360
00340             A24(IELEM) = ((12.D0*(F5-F4)-2.D0*(F1+F3)+4.D0*F6     )*Y2
00341      &                 + (2.D0*F1-4.D0*F4+F3-6.D0*F2-8.D0*F5)*Y3 )
00342      &                 * XSUR90
00343             A25(IELEM) = ((12.D0*(F5-F4)+2.D0*(F1+F3)-4.D0*F6     )*Y2
00344      &                 +  (-F1+4.D0*(F4+F6)-6.D0*F2-16.D0*F5)*Y3 )
00345      &                 * XSUR90
00346             A26(IELEM) = ((F1-8.D0*F5-4.D0*(F6+F4)                )*Y3
00347      &                 - (4.D0*(F4-F5)+F1-F3)*Y2                 )
00348      &                 * XSUR90
00349             A31(IELEM) = ((4.D0*F4-8.D0*(F5-F6)+3.D0*F1-F2-6.D0*F3)*Y2
00350      &                 +  (3.D0*F1-8.D0*(F6-F5)-4.D0*F4+F2)*Y3 )
00351      &                 * XSUR360
00352             A32(IELEM) = ((4.D0*F4-F1-8.D0*(F6-F5)-3.D0*F2)*Y3
00353      &                 +  6.D0*(F2-F3)*Y2) * XSUR360
00354             A34(IELEM) = ((4.D0*(F4+F6)-F1+8.D0*F5)*Y2
00355      &                 +  (F1-4.D0*(F5-F6)-F2)*Y3 ) * XSUR90
00356             A35(IELEM) = ((16.D0*F5-4.D0*(F4+F6)+F1+6.D0*F3)*Y2
00357      &                 +  (-2.D0*(F1+F2)+4.D0*F4+12.D0*(F6-F5))*Y3)
00358      &                 * XSUR90
00359             A36(IELEM) = ((8.D0*F5-F2-2.D0*F1+6.D0*F3+4.D0*F6)*Y2
00360      &                 +  (2.D0*(F2+F1)+12.D0*(F6-F5)-4.D0*F4)*Y3)
00361      &                 * XSUR90
00362 !
00363 !   DIAGONAL TERMS
00364 !
00365             A11(IELEM) = ((24.D0*(F6-F1)-5.D0*F3+4.D0*F5+F2)*Y2
00366      &                 +  (24.D0*(F1-F4)+5.D0*F2-4.D0*F5-F3)*Y3)*XSUR360
00367             A22(IELEM) = ((6.D0*(F1-F3)+24.D0*(F5-F4)      )*Y2
00368      &                 +  (24.D0*(F4-F2)+4.D0*F6+F3-5.D0*F1)*Y3)*XSUR360
00369             A33(IELEM) = ((24.D0*(F3-F6)+5.D0*F1-4.D0*F4-F2)*Y2
00370      &                 +  (6.D0*(F2-F1)+24.D0*(F6-F5)      )*Y3)*XSUR360
00371 !
00372           ENDDO ! IELEM
00373 !
00374         ELSEIF(ICOORD.EQ.2) THEN
00375 !
00376 !================================
00377 !  DERIVATIVE WRT Y  =
00378 !================================
00379 !
00380           DO IELEM = 1 , NELEM
00381 !
00382 !   INITIALISES THE GEOMETRICAL VARIABLES
00383 !
00384             X2  =  XEL(IELEM,2)
00385             X3  =  XEL(IELEM,3)
00386 !
00387             F1  =  F(IKLE1(IELEM))
00388             F2  =  F(IKLE2(IELEM))
00389             F3  =  F(IKLE3(IELEM))
00390             F4  =  F(IKLE4(IELEM))
00391             F5  =  F(IKLE5(IELEM))
00392             F6  =  F(IKLE6(IELEM))
00393 !
00394 !   EXTRADIAGONAL TERMS
00395 !
00396             A12(IELEM) = ((F3+8.D0*(F6-F4)+6.D0*F1-3.D0*F2-4.D0*F5)*X3
00397      &                 +   6.D0*(F2-F1)*X2 ) * XSUR360
00398             A13(IELEM) = ((8.D0*(F6-F4)-6.D0*F1+3.D0*F3+4.D0*F5-F2)*X2
00399      &                 -   6.D0*(F3-F1)*X3 ) * XSUR360
00400             A14(IELEM) = ((6.D0*F1+F3-4.D0*(F6+F5)+16.D0*F4  )*X2
00401      &                 +  (F3-8.D0*F6-6.D0*F1+2.D0*F2-4.D0*F4)*X3)
00402      &                 *  XSUR90
00403             A15(IELEM) = ((4.D0*(F6+F5)+8.D0*F4-F3)*X2
00404      &                 +  (F2-8.D0*F6-4.D0*(F4+F5))*X3 ) * XSUR90
00405             A16(IELEM) = ((6.D0*F1-2.D0*F3+4.D0*F6+8.D0*F4-F2)*X2
00406      &                 +  (-16.D0*F6-6.D0*F1-F2+4.D0*(F5+F4) )*X3)
00407      &                  * XSUR90
00408             A21(IELEM) = ((3.D0*F1+F3-4.D0*F6-8.D0*(F4-F5))*X2
00409      &                 +  (-F3+4.D0*F6+3.D0*F1-6.D0*F2+8.D0*(F4-F5))*X3)
00410      &                 *  XSUR360
00411             A23(IELEM) = ((4.D0*F6-F1-3.D0*F3-8.D0*(F4-F5))*X2
00412      &                 -   6.D0*(F2-F3)*X3 ) * XSUR360
00413             A24(IELEM) = ((2.D0*(F1+F3)-4.D0*F6+12.D0*(F4-F5) )*X2
00414      &                 +  (-F3-2.D0*F1+6.D0*F2+4.D0*F4+8.D0*F5)*X3)
00415      &                 * XSUR90
00416             A25(IELEM) = ((4.D0*F6-2.D0*(F1+F3)+12.D0*(F4-F5) )*X2
00417      &                 -  (4.D0*F6-F1-6.D0*F2+4.D0*F4-16.D0*F5)*X3)
00418      &                 * XSUR90
00419             A26(IELEM) = ((F1-F3+4.D0*(F4-F5)     )*X2
00420      &                 +  (4.D0*(F6+F4)-F1+8.D0*F5)*X3 )
00421      &                 * XSUR90
00422             A31(IELEM) = ((F2-3.D0*F1+6.D0*F3-8.D0*(F6-F5)-4.D0*F4)*X2
00423      &                 -  (-8.D0*F6+3.D0*F1+F2-4.D0*F4+8.D0*F5    )*X3)
00424      &                 * XSUR360
00425             A32(IELEM) = ((F1-8.D0*(F5-F6)+3.D0*F2-4.D0*F4)*X3
00426      &                 +   6.D0*(F3-F2)*X2 ) * XSUR360
00427             A34(IELEM) = ((F1-4.D0*(F6+F4)-8.D0*F5)*X2
00428      &                 +  (-4.D0*(F6-F5)-F1+F2    )*X3 ) * XSUR90
00429             A35(IELEM) = ((4.D0*(F6+F4)-F1-6.D0*F3-16.D0*F5  )*X2
00430      &                 -  (12.D0*(F6-F5)-2.D0*(F1+F2)+4.D0*F4)*X3)
00431      &                 * XSUR90
00432             A36(IELEM) = ((2.D0*F1-6.D0*F3-4.D0*F6-8.D0*F5+F2)*X2
00433      &                 +  (12.D0*(F5-F6)-2.D0*(F1+F2)+4.D0*F4)*X3)
00434      &                 * XSUR90
00435 !
00436 !   DIAGONAL TERMS
00437 !
00438             A11(IELEM) = ((24.D0*(F1-F6)-4.D0*F5+5.D0*F3-F2)*X2
00439      &                 +  (24.D0*(F4-F1)-5.D0*F2+4.D0*F5+F3)*X3)*XSUR360
00440             A22(IELEM) = ((24.D0*(F4-F5)+6.D0*(F3-F1)      )*X2
00441      &                 +  (5.D0*F1+24.D0*(F2-F4)-4.D0*F6-F3)*X3)*XSUR360
00442             A33(IELEM) = ((4.D0*F4-5.D0*F1+24.D0*(F6-F3)+F2)*X2
00443      &                 +  (6.D0*(F1-F2)+24.D0*(F5-F6)      )*X3)*XSUR360
00444 !
00445           ENDDO ! IELEM
00446 !
00447         ELSE
00448 !
00449           IF (LNG.EQ.1) WRITE(LU,200) ICOORD
00450           IF (LNG.EQ.2) WRITE(LU,201) ICOORD
00451           CALL PLANTE(1)
00452           STOP
00453         ENDIF
00454 !
00455 !-----------------------------------------------------------------------
00456 !
00457       ELSE
00458         IF (LNG.EQ.1) WRITE(LU,100) IELMF
00459         IF (LNG.EQ.2) WRITE(LU,101) IELMF
00460 100     FORMAT(1X,'MT11AC (BIEF) :',/,
00461      &         1X,'DISCRETISATION DE F : ',1I6,' NON PREVUE')
00462 101     FORMAT(1X,'MT11AC (BIEF) :',/,
00463      &         1X,'DISCRETIZATION OF F : ',1I6,' NOT AVAILABLE')
00464         CALL PLANTE(1)
00465         STOP
00466       ENDIF
00467 !
00468 200   FORMAT(1X,'MT11AC (BIEF) : COMPOSANTE IMPOSSIBLE ',
00469      &          1I6,' VERIFIER ICOORD')
00470 201   FORMAT(1X,'MT11AC (BIEF) : IMPOSSIBLE COMPONENT ',
00471      &          1I6,' CHECK ICOORD')
00472 !
00473 !-----------------------------------------------------------------------
00474 !
00475       RETURN
00476       END

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