om3181.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\om3181.f
00002 !
00089                      SUBROUTINE OM3181
00090 !                    *****************
00091 !
00092      &(OP ,  DM,TYPDIM,XM,TYPEXM,   DN,TYPDIN,XN,TYPEXN,   C,
00093      & NULONE,NELBOR,NBOR,NELMAX,SIZDN,NELEB,SZMXN)
00094 !
00095 !***********************************************************************
00096 ! BIEF   V6P1                                   21/08/2010
00097 !***********************************************************************
00098 !
00099 !
00100 !
00101 !
00102 !
00103 !
00104 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00105 !| C              |-->| A GIVEN CONSTANT USED IN OPERATION OP
00106 !| DM             |<->| DIAGONAL OF M
00107 !| DN             |-->| DIAGONAL OF N
00108 !| NBOR           |-->| GLOBAL NUMBER OF BOUNDARY POINTS
00109 !| NELBOR         |-->| FOR THE KTH BOUNDARY EDGE, GIVES THE CORRESPONDING
00110 !|                |   | ELEMENT.
00111 !| NELEB          |-->| NUMBER OF BOUNDARY ELEMENTS
00112 !| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
00113 !| NULONE         |-->| GOES WITH ARRAY NELBOR. NELBOR GIVES THE
00114 !|                |   | ADJACENT ELEMENT, NULONE GIVES THE LOCAL
00115 !|                |   | NUMBER OF THE FIRST NODE OF THE BOUNDARY EDGE
00116 !|                |   | I.E. 1, 2 OR 3 FOR TRIANGLES.
00117 !| OP             |-->| OPERATION TO BE DONE (SEE ABOVE)
00118 !| SIZDN          |-->| SIZE OF DIAGONAL DN
00119 !| SZMXN          |-->| SIZE OF OFF-DIAGONAL TERMS XN
00120 !| TYPDIM         |<->| TYPE OF DIAGONAL OF M:
00121 !|                |   | TYPDIM = 'Q' : ANY VALUE
00122 !|                |   | TYPDIM = 'I' : IDENTITY
00123 !|                |   | TYPDIM = '0' : ZERO
00124 !| TYPDIN         |<->| TYPE OF DIAGONAL OF N:
00125 !|                |   | TYPDIN = 'Q' : ANY VALUE
00126 !|                |   | TYPDIN = 'I' : IDENTITY
00127 !|                |   | TYPDIN = '0' : ZERO
00128 !| TYPEXM         |-->| TYPE OF OFF-DIAGONAL TERMS OF M:
00129 !|                |   | TYPEXM = 'Q' : ANY VALUE
00130 !|                |   | TYPEXM = 'S' : SYMMETRIC
00131 !|                |   | TYPEXM = '0' : ZERO
00132 !| TYPEXN         |-->| TYPE OF OFF-DIAGONAL TERMS OF N:
00133 !|                |   | TYPEXN = 'Q' : ANY VALUE
00134 !|                |   | TYPEXN = 'S' : SYMMETRIC
00135 !|                |   | TYPEXN = '0' : ZERO
00136 !| XM             |-->| OFF-DIAGONAL TERMS OF M
00137 !| XN             |-->| OFF-DIAGONAL TERMS OF N
00138 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00139 !
00140       USE BIEF
00141 !
00142       IMPLICIT NONE
00143       INTEGER LNG,LU
00144       COMMON/INFO/LNG,LU
00145 !
00146 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00147 !
00148       INTEGER, INTENT(IN)             :: NELMAX,SIZDN,NELEB,SZMXN
00149       CHARACTER(LEN=8), INTENT(IN)    :: OP
00150       INTEGER, INTENT(IN)             :: NULONE(3*NELEB)
00151       INTEGER, INTENT(IN)             :: NELBOR(*),NBOR(*)
00152       DOUBLE PRECISION, INTENT(IN)    :: DN(*),XN(NELEB,*)
00153       DOUBLE PRECISION, INTENT(INOUT) :: DM(*),XM(NELMAX,*)
00154       CHARACTER(LEN=1), INTENT(INOUT) :: TYPDIM,TYPEXM,TYPDIN,TYPEXN
00155       DOUBLE PRECISION, INTENT(IN)    :: C
00156 !
00157 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00158 !
00159       INTEGER K,IEL,NUL1,NUL2,NUL3
00160 !
00161       DOUBLE PRECISION Z(1)
00162 !
00163       INTEGER CONVSYM(4,4),CONVNSY(4,4)
00164 !
00165 !-----------------------------------------------------------------------
00166 !
00167 !     BEWARE: IN FORTRAN THE FIRST INDEX VARIES MOST QUICKLY
00168 !     123456789 SHOULD NOT BE USED
00169 !
00170       DATA CONVNSY/ 123456789 ,     7     ,     8     ,    9      ,
00171      &                  1     , 123456789 ,    10     ,   11      ,
00172      &                  2     ,     4     , 123456789 ,   12      ,
00173      &                  3     ,     5     ,     6     , 123456789  /
00174 !
00175       DATA CONVSYM/ 123456789 ,     1     ,     2     ,    3      ,
00176      &                  1     , 123456789 ,     4     ,    5      ,
00177      &                  2     ,     4     , 123456789 ,    6      ,
00178      &                  3     ,     5     ,     6     , 123456789  /
00179 !
00180 !***********************************************************************
00181 !
00182       IF(OP(1:8).EQ.'M=M+N   ') THEN
00183         IF(TYPDIM.EQ.'Q'.AND.TYPDIN.EQ.'Q') THEN
00184           CALL OVDB( 'X=X+Y   ' , DM , DN , Z , C , NBOR , SIZDN )
00185         ELSE
00186           IF (LNG.EQ.1) WRITE(LU,198) TYPDIM(1:1),OP(1:8),TYPDIN(1:1)
00187           IF (LNG.EQ.2) WRITE(LU,199) TYPDIM(1:1),OP(1:8),TYPDIN(1:1)
00188 198       FORMAT(1X,'OM5161 (BIEF) : TYPDIM = ',A1,' NON PROGRAMME',
00189      &      /,1X,'POUR L''OPERATION : ',A8,' AVEC TYPDIN = ',A1)
00190 199       FORMAT(1X,'OM5161 (BIEF) : TYPDIM = ',A1,' NOT IMPLEMENTED',
00191      &      /,1X,'FOR THE OPERATION : ',A8,' WITH TYPDIN = ',A1)
00192           CALL PLANTE(1)
00193           STOP
00194         ENDIF
00195 !
00196         IF(TYPEXM(1:1).EQ.'Q'.AND.TYPEXN(1:1).EQ.'Q') THEN
00197 !
00198 !          CASE WHERE BOTH MATRICES ARE NONSYMMETRICAL
00199 !
00200            IF(NCSIZE.GT.1) THEN
00201            DO K = 1 , NELEB
00202              IEL = NELBOR(K)
00203              IF(IEL.GT.0) THEN
00204              NUL1=NULONE(K)
00205              NUL2=NULONE(K+NELEB)
00206              NUL3=NULONE(K+2*NELEB)
00207              XM( IEL , CONVNSY(NUL1,NUL2) ) =
00208      &       XM( IEL , CONVNSY(NUL1,NUL2) ) + XN(K, 1)
00209              XM( IEL , CONVNSY(NUL1,NUL3) ) =
00210      &       XM( IEL , CONVNSY(NUL1,NUL3) ) + XN(K, 2)
00211              XM( IEL , CONVNSY(NUL2,NUL3) ) =
00212      &       XM( IEL , CONVNSY(NUL2,NUL3) ) + XN(K, 3)
00213              XM( IEL , CONVNSY(NUL2,NUL1) ) =
00214      &       XM( IEL , CONVNSY(NUL2,NUL1) ) + XN(K, 4)
00215              XM( IEL , CONVNSY(NUL3,NUL1) ) =
00216      &       XM( IEL , CONVNSY(NUL3,NUL1) ) + XN(K, 5)
00217              XM( IEL , CONVNSY(NUL3,NUL2) ) =
00218      &       XM( IEL , CONVNSY(NUL3,NUL2) ) + XN(K, 6)
00219              ENDIF
00220            ENDDO
00221            ELSE
00222            DO K = 1 , NELEB
00223              IEL = NELBOR(K)
00224              NUL1=NULONE(K)
00225              NUL2=NULONE(K+NELEB)
00226              NUL3=NULONE(K+2*NELEB)
00227              XM( IEL , CONVNSY(NUL1,NUL2) ) =
00228      &       XM( IEL , CONVNSY(NUL1,NUL2) ) + XN(K, 1)
00229              XM( IEL , CONVNSY(NUL1,NUL3) ) =
00230      &       XM( IEL , CONVNSY(NUL1,NUL3) ) + XN(K, 2)
00231              XM( IEL , CONVNSY(NUL2,NUL3) ) =
00232      &       XM( IEL , CONVNSY(NUL2,NUL3) ) + XN(K, 3)
00233              XM( IEL , CONVNSY(NUL2,NUL1) ) =
00234      &       XM( IEL , CONVNSY(NUL2,NUL1) ) + XN(K, 4)
00235              XM( IEL , CONVNSY(NUL3,NUL1) ) =
00236      &       XM( IEL , CONVNSY(NUL3,NUL1) ) + XN(K, 5)
00237              XM( IEL , CONVNSY(NUL3,NUL2) ) =
00238      &       XM( IEL , CONVNSY(NUL3,NUL2) ) + XN(K, 6)
00239            ENDDO
00240            ENDIF
00241 !
00242         ELSEIF(TYPEXM(1:1).EQ.'Q'.AND.TYPEXN(1:1).EQ.'S') THEN
00243 !
00244 !          CASE WHERE M CAN BE ANYTHING AND N IS SYMMETRICAL
00245 !
00246            IF(NCSIZE.GT.1) THEN
00247            DO K = 1 , NELEB
00248              IEL = NELBOR(K)
00249              IF(IEL.GT.0) THEN
00250              NUL1=NULONE(K)
00251              NUL2=NULONE(K+NELEB)
00252              NUL3=NULONE(K+2*NELEB)
00253              XM( IEL , CONVNSY(NUL1,NUL2) ) =
00254      &       XM( IEL , CONVNSY(NUL1,NUL2) ) + XN(K, 1)
00255              XM( IEL , CONVNSY(NUL1,NUL3) ) =
00256      &       XM( IEL , CONVNSY(NUL1,NUL3) ) + XN(K, 2)
00257              XM( IEL , CONVNSY(NUL2,NUL3) ) =
00258      &       XM( IEL , CONVNSY(NUL2,NUL3) ) + XN(K, 3)
00259              XM( IEL , CONVNSY(NUL2,NUL1) ) =
00260      &       XM( IEL , CONVNSY(NUL2,NUL1) ) + XN(K, 1)
00261              XM( IEL , CONVNSY(NUL3,NUL1) ) =
00262      &       XM( IEL , CONVNSY(NUL3,NUL1) ) + XN(K, 2)
00263              XM( IEL , CONVNSY(NUL3,NUL2) ) =
00264      &       XM( IEL , CONVNSY(NUL3,NUL2) ) + XN(K, 3)
00265              ENDIF
00266            ENDDO
00267            ELSE
00268            DO K = 1 , NELEB
00269              IEL = NELBOR(K)
00270              NUL1=NULONE(K)
00271              NUL2=NULONE(K+NELEB)
00272              NUL3=NULONE(K+2*NELEB)
00273              XM( IEL , CONVNSY(NUL1,NUL2) ) =
00274      &       XM( IEL , CONVNSY(NUL1,NUL2) ) + XN(K, 1)
00275              XM( IEL , CONVNSY(NUL1,NUL3) ) =
00276      &       XM( IEL , CONVNSY(NUL1,NUL3) ) + XN(K, 2)
00277              XM( IEL , CONVNSY(NUL2,NUL3) ) =
00278      &       XM( IEL , CONVNSY(NUL2,NUL3) ) + XN(K, 3)
00279              XM( IEL , CONVNSY(NUL2,NUL1) ) =
00280      &       XM( IEL , CONVNSY(NUL2,NUL1) ) + XN(K, 1)
00281              XM( IEL , CONVNSY(NUL3,NUL1) ) =
00282      &       XM( IEL , CONVNSY(NUL3,NUL1) ) + XN(K, 2)
00283              XM( IEL , CONVNSY(NUL3,NUL2) ) =
00284      &       XM( IEL , CONVNSY(NUL3,NUL2) ) + XN(K, 3)
00285            ENDDO
00286            ENDIF
00287 !
00288         ELSEIF(TYPEXM(1:1).EQ.'S'.AND.TYPEXN(1:1).EQ.'S') THEN
00289 !
00290 !          CASE WHERE BOTH MATRICES ARE SYMMETRICAL
00291 !
00292            IF(NCSIZE.GT.1) THEN
00293            DO K = 1 , NELEB
00294              IEL = NELBOR(K)
00295              IF(IEL.GT.0) THEN
00296              NUL1=NULONE(K)
00297              NUL2=NULONE(K+NELEB)
00298              NUL3=NULONE(K+2*NELEB)
00299              XM( IEL , CONVSYM(NUL1,NUL2) ) =
00300      &       XM( IEL , CONVSYM(NUL1,NUL2) ) + XN(K, 1)
00301              XM( IEL , CONVSYM(NUL1,NUL3) ) =
00302      &       XM( IEL , CONVSYM(NUL1,NUL3) ) + XN(K, 2)
00303              XM( IEL , CONVSYM(NUL2,NUL3) ) =
00304      &       XM( IEL , CONVSYM(NUL2,NUL3) ) + XN(K, 3)
00305              ENDIF
00306            ENDDO
00307            ELSE
00308            DO K = 1 , NELEB
00309              IEL = NELBOR(K)
00310              NUL1=NULONE(K)
00311              NUL2=NULONE(K+NELEB)
00312              NUL3=NULONE(K+2*NELEB)
00313              XM( IEL , CONVSYM(NUL1,NUL2) ) =
00314      &       XM( IEL , CONVSYM(NUL1,NUL2) ) + XN(K, 1)
00315              XM( IEL , CONVSYM(NUL1,NUL3) ) =
00316      &       XM( IEL , CONVSYM(NUL1,NUL3) ) + XN(K, 2)
00317              XM( IEL , CONVSYM(NUL2,NUL3) ) =
00318      &       XM( IEL , CONVSYM(NUL2,NUL3) ) + XN(K, 3)
00319            ENDDO
00320            ENDIF
00321 !
00322         ELSE
00323            IF (LNG.EQ.1) WRITE(LU,98) TYPEXM(1:1),OP(1:8),TYPEXN(1:1)
00324            IF (LNG.EQ.2) WRITE(LU,99) TYPEXM(1:1),OP(1:8),TYPEXN(1:1)
00325 98         FORMAT(1X,'OM3181 (BIEF) : TYPEXM = ',A1,' NE CONVIENT PAS',
00326      &       /,1X,'POUR L''OPERATION : ',A8,' AVEC TYPEXN = ',A1)
00327 99         FORMAT(1X,'OM3181 (BIEF) : TYPEXM = ',A1,' DOES NOT GO',
00328      &       /,1X,'FOR THE OPERATION : ',A8,' WITH TYPEXN = ',A1)
00329            CALL PLANTE(1)
00330            STOP
00331         ENDIF
00332 !
00333 !-----------------------------------------------------------------------
00334 !
00335       ELSEIF(OP(1:8).EQ.'M=M+TN  ') THEN
00336 !
00337         CALL OVDB( 'X=X+Y   ' , DM , DN , Z , C , NBOR , NELEB )
00338 !
00339         IF(TYPEXM(1:1).EQ.'Q'.AND.TYPEXN(1:1).EQ.'Q') THEN
00340 !
00341 !          CASE WHERE BOTH MATRICES ARE NONSYMMETRICAL
00342 !
00343            IF(NCSIZE.GT.1) THEN
00344            DO K = 1 , NELEB
00345              IEL = NELBOR(K)
00346              IF(IEL.GT.0) THEN
00347              NUL1=NULONE(K)
00348              NUL2=NULONE(K+NELEB)
00349              NUL3=NULONE(K+2*NELEB)
00350              XM( IEL , CONVNSY(NUL1,NUL2) ) =
00351      &       XM( IEL , CONVNSY(NUL1,NUL2) ) + XN(K, 4)
00352              XM( IEL , CONVNSY(NUL1,NUL3) ) =
00353      &       XM( IEL , CONVNSY(NUL1,NUL3) ) + XN(K, 5)
00354              XM( IEL , CONVNSY(NUL2,NUL3) ) =
00355      &       XM( IEL , CONVNSY(NUL2,NUL3) ) + XN(K, 6)
00356              XM( IEL , CONVNSY(NUL2,NUL1) ) =
00357      &       XM( IEL , CONVNSY(NUL2,NUL1) ) + XN(K, 1)
00358              XM( IEL , CONVNSY(NUL3,NUL1) ) =
00359      &       XM( IEL , CONVNSY(NUL3,NUL1) ) + XN(K, 2)
00360              XM( IEL , CONVNSY(NUL3,NUL2) ) =
00361      &       XM( IEL , CONVNSY(NUL3,NUL2) ) + XN(K, 3)
00362              ENDIF
00363            ENDDO
00364            ELSE
00365            DO K = 1 , NELEB
00366              IEL = NELBOR(K)
00367              NUL1=NULONE(K)
00368              NUL2=NULONE(K+NELEB)
00369              NUL3=NULONE(K+2*NELEB)
00370              XM( IEL , CONVNSY(NUL1,NUL2) ) =
00371      &       XM( IEL , CONVNSY(NUL1,NUL2) ) + XN(K, 4)
00372              XM( IEL , CONVNSY(NUL1,NUL3) ) =
00373      &       XM( IEL , CONVNSY(NUL1,NUL3) ) + XN(K, 5)
00374              XM( IEL , CONVNSY(NUL2,NUL3) ) =
00375      &       XM( IEL , CONVNSY(NUL2,NUL3) ) + XN(K, 6)
00376              XM( IEL , CONVNSY(NUL2,NUL1) ) =
00377      &       XM( IEL , CONVNSY(NUL2,NUL1) ) + XN(K, 1)
00378              XM( IEL , CONVNSY(NUL3,NUL1) ) =
00379      &       XM( IEL , CONVNSY(NUL3,NUL1) ) + XN(K, 2)
00380              XM( IEL , CONVNSY(NUL3,NUL2) ) =
00381      &       XM( IEL , CONVNSY(NUL3,NUL2) ) + XN(K, 3)
00382            ENDDO
00383            ENDIF
00384 !
00385         ELSEIF(TYPEXM(1:1).EQ.'Q'.AND.TYPEXN(1:1).EQ.'S') THEN
00386 !
00387 !          CASE WHERE M CAN BE ANYTHING AND N IS SYMMETRICAL
00388 !
00389            IF(NCSIZE.GT.1) THEN
00390            DO K = 1 , NELEB
00391              IEL = NELBOR(K)
00392              IF(IEL.GT.0) THEN
00393              NUL1=NULONE(K)
00394              NUL2=NULONE(K+NELEB)
00395              NUL3=NULONE(K+2*NELEB)
00396              XM( IEL , CONVNSY(NUL1,NUL2) ) =
00397      &       XM( IEL , CONVNSY(NUL1,NUL2) ) + XN(K, 1)
00398              XM( IEL , CONVNSY(NUL1,NUL3) ) =
00399      &       XM( IEL , CONVNSY(NUL1,NUL3) ) + XN(K, 2)
00400              XM( IEL , CONVNSY(NUL2,NUL3) ) =
00401      &       XM( IEL , CONVNSY(NUL2,NUL3) ) + XN(K, 3)
00402              XM( IEL , CONVNSY(NUL2,NUL1) ) =
00403      &       XM( IEL , CONVNSY(NUL2,NUL1) ) + XN(K, 1)
00404              XM( IEL , CONVNSY(NUL3,NUL1) ) =
00405      &       XM( IEL , CONVNSY(NUL3,NUL1) ) + XN(K, 2)
00406              XM( IEL , CONVNSY(NUL3,NUL2) ) =
00407      &       XM( IEL , CONVNSY(NUL3,NUL2) ) + XN(K, 3)
00408              ENDIF
00409            ENDDO
00410            ELSE
00411            DO K = 1 , NELEB
00412              IEL = NELBOR(K)
00413              NUL1=NULONE(K)
00414              NUL2=NULONE(K+NELEB)
00415              NUL3=NULONE(K+2*NELEB)
00416              XM( IEL , CONVNSY(NUL1,NUL2) ) =
00417      &       XM( IEL , CONVNSY(NUL1,NUL2) ) + XN(K, 1)
00418              XM( IEL , CONVNSY(NUL1,NUL3) ) =
00419      &       XM( IEL , CONVNSY(NUL1,NUL3) ) + XN(K, 2)
00420              XM( IEL , CONVNSY(NUL2,NUL3) ) =
00421      &       XM( IEL , CONVNSY(NUL2,NUL3) ) + XN(K, 3)
00422              XM( IEL , CONVNSY(NUL2,NUL1) ) =
00423      &       XM( IEL , CONVNSY(NUL2,NUL1) ) + XN(K, 1)
00424              XM( IEL , CONVNSY(NUL3,NUL1) ) =
00425      &       XM( IEL , CONVNSY(NUL3,NUL1) ) + XN(K, 2)
00426              XM( IEL , CONVNSY(NUL3,NUL2) ) =
00427      &       XM( IEL , CONVNSY(NUL3,NUL2) ) + XN(K, 3)
00428            ENDDO
00429            ENDIF
00430 !
00431         ELSEIF(TYPEXM(1:1).EQ.'S'.AND.TYPEXN(1:1).EQ.'S') THEN
00432 !
00433 !          CASE WHERE BOTH MATRICES ARE SYMMETRICAL
00434 !
00435            IF(NCSIZE.GT.1) THEN
00436            DO K = 1 , NELEB
00437              IEL = NELBOR(K)
00438              IF(IEL.GT.0) THEN
00439              NUL1=NULONE(K)
00440              NUL2=NULONE(K+NELEB)
00441              NUL3=NULONE(K+2*NELEB)
00442              XM( IEL , CONVSYM(NUL1,NUL2) ) =
00443      &       XM( IEL , CONVSYM(NUL1,NUL2) ) + XN(K, 1)
00444              XM( IEL , CONVSYM(NUL1,NUL3) ) =
00445      &       XM( IEL , CONVSYM(NUL1,NUL3) ) + XN(K, 2)
00446              XM( IEL , CONVSYM(NUL2,NUL3) ) =
00447      &       XM( IEL , CONVSYM(NUL2,NUL3) ) + XN(K, 3)
00448              ENDIF
00449            ENDDO
00450            ELSE
00451            DO K = 1 , NELEB
00452              IEL = NELBOR(K)
00453              NUL1=NULONE(K)
00454              NUL2=NULONE(K+NELEB)
00455              NUL3=NULONE(K+2*NELEB)
00456              XM( IEL , CONVSYM(NUL1,NUL2) ) =
00457      &       XM( IEL , CONVSYM(NUL1,NUL2) ) + XN(K, 1)
00458              XM( IEL , CONVSYM(NUL1,NUL3) ) =
00459      &       XM( IEL , CONVSYM(NUL1,NUL3) ) + XN(K, 2)
00460              XM( IEL , CONVSYM(NUL2,NUL3) ) =
00461      &       XM( IEL , CONVSYM(NUL2,NUL3) ) + XN(K, 3)
00462            ENDDO
00463            ENDIF
00464 !
00465         ELSE
00466            IF (LNG.EQ.1) WRITE(LU,98) TYPEXM(1:1),OP(1:8),TYPEXN(1:1)
00467            IF (LNG.EQ.2) WRITE(LU,99) TYPEXM(1:1),OP(1:8),TYPEXN(1:1)
00468            CALL PLANTE(1)
00469            STOP
00470         ENDIF
00471 !
00472 !-----------------------------------------------------------------------
00473 !
00474       ELSE
00475 !
00476         IF (LNG.EQ.1) WRITE(LU,70) OP
00477         IF (LNG.EQ.2) WRITE(LU,71) OP
00478 70      FORMAT(1X,'OM3181 (BIEF) : OPERATION INCONNUE : ',A8)
00479 71      FORMAT(1X,'OM3181 (BIEF) : UNKNOWN OPERATION : ',A8)
00480         CALL PLANTE(1)
00481         STOP
00482 !
00483       ENDIF
00484 !
00485 !-----------------------------------------------------------------------
00486 !
00487       RETURN
00488       END

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