cpik13.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\cpik13.f
00002 !
00072                      SUBROUTINE CPIK13
00073 !                    *****************
00074 !
00075      &(IKLE,IKLBOR,ELTSEG,NBOR,NELBOR,NULONE,NELEM,NELMAX,NPOIN,NPTFR,
00076      & NELEB,NELEBX)
00077 !
00078 !***********************************************************************
00079 ! BIEF   V7P0                                     28/03/2014
00080 !***********************************************************************
00081 !
00082 !
00083 !
00084 !
00085 !
00086 !
00087 !
00088 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00089 !| ELTSEG         |-->| SEGMENTS NUMBERS IN EVERY ELEMENT
00090 !| IKLBOR         |-->| CONNECTIVITY TABLE OF BOUNDARY ELEMENTS
00091 !| IKLE           |<->| CONNECTIVITY TABLE
00092 !| NBOR           |-->| GLOBAL NUMBERS OF BOUNDARY POINTS
00093 !| NELBOR         |-->| BOUNDARY ELEMENT THAT CONTAINS SEGMENT K
00094 !| NELEB          |-->| NUMBER OF BOUNDARY SEGMENTS
00095 !| NELEBX         |-->| MAXIMUM NUMBER OF BOUNDARY SEGMENTS
00096 !| NULONE         |-->| LOCAL NUMBER OF K IN ELEMENT NELBOR
00097 !| NELEM          |-->| NUMBER OF ELEMENTS
00098 !| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
00099 !| NPOIN          |-->| NUMBER OF POINTS IN THE MESH
00100 !| NPTFR          |-->| NUMBER OF BOUNDARY POINTS
00101 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00102 !
00103       USE BIEF, ONLY : NCSIZE
00104 !
00105       IMPLICIT NONE
00106       INTEGER LNG,LU
00107       COMMON/INFO/LNG,LU
00108 !
00109 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00110 !
00111       INTEGER, INTENT(IN)    :: NELEM,NELMAX,NPOIN,NPTFR,NELEB,NELEBX
00112       INTEGER, INTENT(IN)    :: ELTSEG(NELMAX,3)
00113       INTEGER, INTENT(IN)    :: NELBOR(NELEBX),NULONE(NELEBX)
00114       INTEGER, INTENT(INOUT) :: IKLE(NELMAX,6),IKLBOR(NELEBX,3),NBOR(*)
00115 !
00116 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00117 !
00118       INTEGER IELEM,K,IELEB
00119 !
00120 !-----------------------------------------------------------------------
00121 !
00122 !     CONNECTIVITY TABLE OF QUADRATIC GLOBAL POINTS
00123 !
00124       DO IELEM = 1 , NELEM
00125 !
00126 !       NUMBER=NPOIN+NUMBER OF THE SEGMENT CONTAINING THE POINT
00127 !
00128         IKLE(IELEM,4) = NPOIN + ELTSEG(IELEM,1)
00129         IKLE(IELEM,5) = NPOIN + ELTSEG(IELEM,2)
00130         IKLE(IELEM,6) = NPOIN + ELTSEG(IELEM,3)
00131 !
00132       ENDDO
00133 !
00134 !-----------------------------------------------------------------------
00135 !
00136 !     CONNECTIVITY TABLE OF QUADRATIC BOUNDARY POINTS
00137 !     GLOBAL NUMBERS OF BOUNDARY QUADRATIC POINTS
00138 !
00139       DO IELEB=1,NELEB
00140         K=IKLBOR(IELEB,1)
00141 !       WE DECIDE HERE THE BOUNDARY NUMBERING OF QUADRATIC POINTS
00142         IKLBOR(IELEB,3)=K+NPTFR
00143         IELEM=NELBOR(IELEB)
00144         NBOR(IKLBOR(IELEB,3))=IKLE(IELEM,NULONE(IELEB)+3)
00145       ENDDO
00146 !
00147 !-----------------------------------------------------------------------
00148 !
00149 !     SECURITY CHECK IN SCALAR MODE
00150 !
00151       IF(NCSIZE.LE.1) THEN
00152         DO IELEB=1,NELEB
00153           K=IKLBOR(IELEB,1)
00154           IF(NBOR(IKLBOR(IELEB,3)).NE.K+NPOIN) THEN
00155             WRITE(LU,*) 'CPIK13: PROBLEM OF NUMBERING FOR K=',K
00156             CALL PLANTE(1)
00157             STOP
00158           ENDIF
00159         ENDDO
00160       ENDIF
00161 !
00162 !-----------------------------------------------------------------------
00163 !
00164       RETURN
00165       END

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