kepicl.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\kepicl.f
00002 !
00065                      SUBROUTINE KEPICL
00066 !                    *****************
00067 !
00068      & (LIKBOF,LIEBOF,LIUBOF,LIKBOL,LIEBOL,LIUBOL,LIKBOS,LIEBOS,LIUBOS,
00069      &  NPTFR,NPLAN,NPOIN2,KENT,KSORT,KADH,KLOG)
00070 !
00071 !***********************************************************************
00072 ! TELEMAC3D   V6P1                                   21/08/2010
00073 !***********************************************************************
00074 !
00075 !
00076 !
00077 !
00078 !
00079 !
00080 !
00081 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00082 !| KADH           |-->| CONVENTION FOR NO SLIP BOUNDARY CONDITION
00083 !| KENT           |-->| CONVENTION FOR LIQUID INPUT WITH PRESCRIBED VALUE
00084 !| KLOG           |-->| CONVENTION FOR LOGARITHMIC SOLID BOUNDARY
00085 !| KSORT          |-->| CONVENTION FOR FREE OUTPUT
00086 !| LIEBOF         |<->| TYPE OF BOUNDARY CONDITIONS ON EPSILON AT THE BOTTOM
00087 !| LIEBOL         |<->| TYPE OF BOUNDARY CONDITIONS ON EPSILON ON THE LATERAL WALLS
00088 !| LIEBOS         |<->| TYPE OF BOUNDARY CONDITIONS ON EPSILON AT THE SURFACE
00089 !| LIKBOF         |<->| TYPE OF BOUNDARY CONDITIONS ON K AT THE BOTTOM
00090 !| LIKBOL         |<->| TYPE OF BOUNDARY CONDITIONS ON K ON THE LATERAL WALLS
00091 !| LIKBOS         |<->| TYPE OF BOUNDARY CONDITIONS ON K AT THE SURFACE
00092 !| LIUBOF         |-->| TYPE OF BOUNDARY CONDITIONS ON U AT THE BOTTOM
00093 !| LIUBOL         |-->| TYPE OF BOUNDARY CONDITIONS ON U ON THE LATERAL WALLS
00094 !| LIUBOS         |-->| TYPE OF BOUNDARY CONDITIONS ON U AT THE SURFACE
00095 !| NPLAN          |-->| NUMBER OF PLANES IN THE 3D MESH OF PRISMS
00096 !| NPOIN2         |-->| NUMBER OF POINTS IN 2D
00097 !| NPTFR          |-->| NUMBER OF BOUNDARY POINTS IN 2D
00098 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00099 !
00100       USE DECLARATIONS_TELEMAC3D, ONLY : LIMKF,LIMEF,LIMKS,LIMES
00101 !
00102       USE BIEF
00103 !
00104       IMPLICIT NONE
00105       INTEGER LNG,LU
00106       COMMON/INFO/LNG,LU
00107 !
00108 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00109 !
00110       INTEGER, INTENT(IN)    :: NPTFR, NPLAN, NPOIN2
00111       INTEGER, INTENT(IN)    :: KENT, KSORT, KADH, KLOG
00112       INTEGER, INTENT(IN)    :: LIUBOF(NPOIN2), LIUBOS(NPOIN2)
00113       INTEGER, INTENT(IN)    :: LIUBOL(NPTFR*NPLAN*2)
00114       INTEGER, INTENT(INOUT) :: LIKBOF(NPOIN2), LIKBOS(NPOIN2)
00115       INTEGER, INTENT(INOUT) :: LIKBOL(NPTFR*NPLAN*2)
00116       INTEGER, INTENT(INOUT) :: LIEBOF(NPOIN2), LIEBOS(NPOIN2)
00117       INTEGER, INTENT(INOUT) :: LIEBOL(NPTFR*NPLAN*2)
00118 !
00119 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00120 !
00121       INTEGER IPTFR, IPOIN2,NPTFR3
00122 !
00123 !-----------------------------------------------------------------------
00124 !
00125 !     LATERAL BOUNDARIES :
00126 !
00127 !-----------------------------------------------------------------------
00128 !
00129       NPTFR3=NPLAN*NPTFR
00130 !
00131       DO IPTFR=1,NPTFR3
00132         IF(LIUBOL(IPTFR).EQ.KENT) THEN
00133           LIKBOL(IPTFR) = KENT
00134           LIEBOL(IPTFR) = KENT
00135         ELSE
00136           LIKBOL(IPTFR) = KSORT
00137           LIEBOL(IPTFR) = KSORT
00138         ENDIF
00139 !       SAVING VALUES IN THE SECOND DIMENSION (SEE POINT_TELEMAC3D)
00140         LIKBOL(IPTFR+NPTFR3) = KSORT
00141         LIEBOL(IPTFR+NPTFR3) = KSORT
00142       ENDDO
00143 !
00144 !-----------------------------------------------------------------------
00145 !
00146 !     BOTTOM
00147 !
00148 !-----------------------------------------------------------------------
00149 !
00150       IF(LIMKF.EQ.2) THEN
00151         DO IPOIN2 = 1,NPOIN2
00152           LIKBOF(IPOIN2) = KENT
00153         ENDDO
00154       ELSE
00155         DO IPOIN2 = 1,NPOIN2
00156           IF(LIUBOF(IPOIN2).EQ.KSORT) THEN
00157             LIKBOF(IPOIN2) = KSORT
00158           ELSE
00159             LIKBOF(IPOIN2) = KENT
00160           ENDIF
00161         ENDDO
00162       ENDIF
00163 !
00164       IF(LIMEF.EQ.2) THEN
00165         DO IPOIN2 = 1,NPOIN2
00166           LIEBOF(IPOIN2) = KENT
00167         ENDDO
00168       ELSE
00169         DO IPOIN2 = 1,NPOIN2
00170           IF(LIUBOF(IPOIN2).EQ.KSORT) THEN
00171             LIEBOF(IPOIN2) = KSORT
00172           ELSE
00173             LIEBOF(IPOIN2) = KENT
00174           ENDIF
00175         ENDDO
00176       ENDIF
00177 !
00178 !-----------------------------------------------------------------------
00179 !
00180 !     FREE SURFACE
00181 !
00182 !-----------------------------------------------------------------------
00183 !
00184       IF(LIMES.EQ.2) THEN
00185         DO IPOIN2 = 1,NPOIN2
00186           LIKBOS(IPOIN2) = KENT
00187         ENDDO
00188       ELSE
00189         DO IPOIN2 = 1,NPOIN2
00190           LIKBOS(IPOIN2) = KSORT
00191         ENDDO
00192       ENDIF
00193 !
00194       IF(LIMKS.EQ.2) THEN
00195         DO IPOIN2 = 1,NPOIN2
00196           LIEBOS(IPOIN2) = KENT
00197         ENDDO
00198       ELSE
00199         DO IPOIN2 = 1,NPOIN2
00200           LIEBOS(IPOIN2) = KSORT
00201         ENDDO
00202       ENDIF
00203 !
00204 !-----------------------------------------------------------------------
00205 !
00206       RETURN
00207       END

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