point_postel3d.f

Go to the documentation of this file.

!                       *************************
                        SUBROUTINE point_postel3d
!                       *************************
!
!***********************************************************************
! POSTEL3D VERSION 5.1   01/09/99   T. DENOT (LNH) 01 30 87 74 89
! FORTRAN90
!
!history Y AUDOUIN (LNHE)
!+       25/05/2015
!+       V7P0
!+       Modification to comply with the hermes module
!
!***********************************************************************
!
      USE bief
      USE declarations_telemac
      USE declarations_postel3d
!
      USE declarations_special
      IMPLICIT NONE
!
!
!***********************************************************************
!
      INTEGER CFG(2),NGEO
      ngeo=pos_files(posgeo)%LU
!
!-----------------------------------------------------------------------
!
      WRITE(lu,21)
 21   FORMAT(1x,/,1x,'POINT_TELEMAC3D: MEMORY ALLOCATION',/)
!
!-----------------------------------------------------------------------
! discretisation types are declared here
!
      ielm0 = 10*(ielmh/10) ! FOR TELEMAC2D
      ielm1 = ielm0 + 1     ! FOR TELEMAC2D
!
! Telemac3D discretisation types: 3D, 2D horizontal boundary,
! 2D vertical boundary
!
      ielm3  = 41     ! TELEMAC3D PRISMS
      ielm2h = 11     ! PRISM TRIANGULAR BOTTOM AND SURFACE
      ielm2v = 21     ! PRISM QUADRILATERAL LATERAL BOUNDARIES
!
      ielmx=max(ielmu,ielm2h,ielmh) ! IT WILL BE MAX. DISCR. IN 2D
!
      cfg(1) = 1
      cfg(2) = 1
!
!=======================================================================
!
!                     *********************
!                     *  MESH - GEOMETRY  *
!                     *********************
!
! TWO meshes are allocated: (1) 2D base mesh, (2) 3D sigma-mesh
!
! allocation of the 2D mesh structure for Telemac2D
! discretisation ielmh given in lecdon
! ielmx = ielmu if quasi-bubble element required, otherwise ielmh
!
      equa = 'NO_EQUATION_IS_GIVEN'
!
      CALL almesh(mesh2d,'MESH2D',ielmx,spheri,cfg,
     &            pos_files(posgeo)%FMT,ngeo,
     &            equa,0,nplan=1)
!
! aliases for certain components of the 2D mesh structure
!
      x2      => mesh2d%X
      y2      => mesh2d%Y
      z2      => mesh2d%Z
      xnebor2 => mesh2d%XNEBOR
      ynebor2 => mesh2d%YNEBOR
      xsgbor2 => mesh2d%XSGBOR
      ysgbor2 => mesh2d%YSGBOR
      ikle2   => mesh2d%IKLE
      nbor2   => mesh2d%NBOR   ! PREV. SIMPLY NBOR
!
!
      nelem2  => mesh2d%NELEM
      nelmax2 => mesh2d%NELMAX  ! PREVIOUSLY NELMA2 (ADAPTIVITY OUTLOOK)
      nptfr2  => mesh2d%NPTFR   ! PREVIOUSLY SIMPLY NPTFR
      nptfrx2 => mesh2d%NPTFRX
      dim2    => mesh2d%DIM1
      typelm2 => mesh2d%TYPELM
      npoin2  => mesh2d%NPOIN
      npmax2  => mesh2d%NPMAX
      mxptvs2 => mesh2d%MXPTVS
      mxelvs2 => mesh2d%MXELVS
!
!-----------------------------------------------------------------------
! ALLOCATION OF THE 3D MESH STRUCTURE (EQUA=EMPTY) (READ AGAIN?)
!
      equa = 'NO_EQUATION_IS_GIVEN'
!
      CALL almesh(mesh3d,'MESH3D',ielm3,spheri,cfg,
     &            pos_files(posgeo)%FMT,ngeo,
     &            equa,0,nplan=nplan)
!
! alias for certain components of the 3D mesh structure
! they are defined in declarations
!
      x       => mesh3d%X%R  ! REAL VALUE!!!
      y       => mesh3d%Y%R
      z       => mesh3d%Z%R
      x3      => mesh3d%X    ! POINTERS
      y3      => mesh3d%Y
      z3      => mesh3d%Z
      xnebor3 => mesh3d%XNEBOR
      ynebor3 => mesh3d%YNEBOR
      znebor3 => mesh3d%ZNEBOR
      xsgbor3 => mesh3d%XSGBOR
      ysgbor3 => mesh3d%YSGBOR
      zsgbor3 => mesh3d%ZSGBOR
      ikle3   => mesh3d%IKLE
      nbor3   => mesh3d%NBOR
!
      nelem3  => mesh3d%NELEM
      nelmax3 => mesh3d%NELMAX   ! PREVIOUSLY NELMA3 (ADAPTIVITY?)
      neleb   => mesh3d%NELEB
      nelebx  => mesh3d%NELEBX
      nptfr3  => mesh3d%NPTFR
      nptfrx3 => mesh3d%NPTFRX
      dim3    => mesh3d%DIM1
      typelm3 => mesh3d%TYPELM
      npoin3  => mesh3d%NPOIN
      npmax3  => mesh3d%NPMAX
      mxptvs3 => mesh3d%MXPTVS
      mxelvs3 => mesh3d%MXELVS
!
!
      WRITE(lu,32)
     &           typelm2,npoin2,nelem2,nptfr2,typelm3,npoin3,nelem3,
     &           nplan,neleb,nptfr3+2*npoin2,nptfr3,npoin2,npoin2
!
 32   FORMAT(/,' 2D MESH',/,
     &         ' -------',//,
     &         ' 2D ELEMENT TYPE                : ',i8,/,
     &         ' NUMBER OF 2D NODES             : ',i8,/,
     &         ' NUMBER OF 2D ELEMENTS          : ',i8,/,
     &         ' NUMBER OF 2D BOUNDARY NODES    : ',i8,///,
     &         ' 3D MESH',/,
     &         ' -------',//,
     &         ' 3D ELEMENT TYPE                : ',i8,/,
     &         ' NUMBER OF 3D NODES             : ',i8,/,
     &         ' NUMBER OF 3D ELEMENTS          : ',i8,/,
     &         ' NUMBER OF LEVELS               : ',i8,/,
     &         ' NUMBER OF BOUNDARY ELEMENTS    : ',i8,/,
     &         ' TOTAL NUMBER OF BOUNDARY NODES : ',i8,/,
     &         ' INCLUDING   LATERAL BOUNDARIES : ',i8,/,
     &         '                        SURFACE : ',i8,/,
     &         '                         BOTTOM : ',i8,/)
!
!
! DEFINITION DES POINTEURS
!
      CALL bief_allvec(1, u,      'U     ', ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, v,      'V     ', ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, w,      'W     ', ielm3,  1,1,mesh3d)
!
      CALL allblo(tab,'TAB   ')
      IF (nva3.GT.4) THEN
        CALL bief_allvec_in_block(tab,nva3-4,1,'TAB   ',ielm3,
     &                            1,1,mesh3d)
      ENDIF
!
!=======================================================================
!
! IMPRESSIONS :
!
      WRITE(lu,23)
23    FORMAT(1x,///,21x,'*************************************',/,
     &21x,              '*    END OF MEMORY ORGANIZATION:    *',/,
     &21x,              '*************************************',/)
!
!-----------------------------------------------------------------------
!
      RETURN
      END