corpor.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\corpor.f
00002 !
00082                      SUBROUTINE CORPOR
00083 !                    *****************
00084 !
00085      &(POROS)
00086 !
00087 !***********************************************************************
00088 ! TELEMAC2D   V6P1                                   21/08/2010
00089 !***********************************************************************
00090 !
00091 !
00092 !
00093 !
00094 !
00095 !
00096 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00097 !| POROS          |<->| POROSITY TO BE MODIFIED.
00098 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00099 !
00100       USE BIEF
00101       USE DECLARATIONS_TELEMAC2D
00102 !
00103       IMPLICIT NONE
00104       INTEGER LNG,LU
00105       COMMON/INFO/LNG,LU
00106 !
00107 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00108 !
00109       TYPE(BIEF_OBJ), INTENT(INOUT) :: POROS
00110 !
00111 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00112 !
00113 !     DOUBLE PRECISION XSOM(4),YSOM(4),XX1,YY1
00114 !     INTEGER NSOM,IELEM
00115 !
00116 !-----------------------------------------------------------------------
00117 !
00118 !     EXAMPLE : POROSITY IS SET TO 0.5 IN A QUADRILATERAL
00119 !
00120 !     NSOM = 4
00121 !     XSOM(1) = 8020.88D0
00122 !     XSOM(2) = 7761.81D0
00123 !     XSOM(3) = 8679.17D0
00124 !     XSOM(4) = 8988.75D0
00125 !     YSOM(1) =-1547.11D0
00126 !     YSOM(2) =-2415.26D0
00127 !     YSOM(3) =-2604.16D0
00128 !     YSOM(4) =-1543.75D0
00129 !
00130 !-----------------------------------------------------------------------
00131 !
00132 !     CALL OS( 'X=C     ' , POROS , POROS , POROS , 1.D0 )
00133 !
00134 !--------------------------------------------------------------
00135 !
00136 !     DO IELEM = 1 , NELEM
00137 !
00138 !       XX1 = (  X(IKLE%I(IELEM)          )+
00139 !    *           X(IKLE%I(IELEM+NELMAX)   )+
00140 !    *           X(IKLE%I(IELEM+2*NELMAX) ))/3.D0
00141 !       YY1 = (  Y(IKLE%I(IELEM)          )+
00142 !    *           Y(IKLE%I(IELEM+NELMAX)   )+
00143 !    *           Y(IKLE%I(IELEM+2*NELMAX) ))/3.D0
00144 !
00145 !       IF(INPOLY(XX1,YY1,XSOM,YSOM,NSOM)) THEN
00146 !         POROS%R(IELEM) = 0.5D0
00147 !       ENDIF
00148 !
00149 !      ENDDO ! IELEM
00150 !
00151 !-----------------------------------------------------------------------
00152 !
00153       RETURN
00154       END

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