poros.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\poros.f
00002 !
00073                      SUBROUTINE POROS
00074 !                    ****************
00075 !
00076      &(TETA,ZF,HN,MESH)
00077 !
00078 !***********************************************************************
00079 ! TELEMAC2D   V6P1                                   21/08/2010
00080 !***********************************************************************
00081 !
00082 !
00083 !
00084 !
00085 !
00086 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00087 !| HN             |-->| WATER DEPTH AT TIME T(N)
00088 !| MESH           |-->| MESH STRUCTURE
00089 !| TETA           |<--| POROSITY PER ELEMENT
00090 !| ZF             |-->| BOTTOM TOPOGRAPHY
00091 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00092 !
00093       USE BIEF
00094       USE INTERFACE_TELEMAC2D, EX_POROS => POROS
00095 !
00096       IMPLICIT NONE
00097       INTEGER LNG,LU
00098       COMMON/INFO/LNG,LU
00099 !
00100 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00101 !
00102       TYPE(BIEF_OBJ), INTENT(IN)    :: ZF,HN
00103       TYPE(BIEF_OBJ), INTENT(INOUT) :: TETA
00104       TYPE(BIEF_MESH), INTENT(IN)   :: MESH
00105 !
00106 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00107 !
00108       INTEGER IELMZ,IELMH
00109 !
00110       IELMZ=ZF%ELM
00111       IELMH=HN%ELM
00112 !
00113 !-----------------------------------------------------------------------
00114 !
00115 !     1) COMPUTES POROSITY ON TIDAL FLATS
00116 !
00117       IF(IELMZ.EQ.11.AND.IELMH.EQ.11) THEN
00118 !
00119         CALL PORO11(TETA%R,ZF%R,
00120      &              HN%R,MESH%IKLE%I,MESH%NELEM,MESH%NELMAX)
00121 !
00122       ELSE
00123 !
00124         IF(LNG.EQ.1) WRITE(LU,10) IELMH,IELMZ
00125         IF(LNG.EQ.2) WRITE(LU,11) IELMH,IELMZ
00126 10      FORMAT(1X,'POROS : DISCRETISATION NON PREVUE :'    ,I6,' ',I6)
00127 11      FORMAT(1X,'POROS : DISCRETIZATION NOT IMPLEMENTED:',I6,' ',I6)
00128         CALL PLANTE(1)
00129         STOP
00130 !
00131       ENDIF
00132 !
00133 !     2) CORRECTED BY A USER SUBROUTINE
00134 !
00135       CALL CORPOR(TETA)
00136 !
00137 !-----------------------------------------------------------------------
00138 !
00139       RETURN
00140       END

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