point_waqtel.f

Go to the documentation of this file.

!                   ***********************
                    SUBROUTINE point_waqtel
!                   ***********************
!
     &(mesh2d,ielm1,mesh3d,ielm3)
!
!***********************************************************************
! WAQTEL   V8P2
!***********************************************************************
!
!brief    Memory allocation of structures, aliases, blocks...
!
!history  R ATA (LNHE)
!+        13/05/2015
!+        V7P1
!+
!+        Creation of the file
!
!history  S.E. BOURBAN (HRW)
!+        21/09/2017
!+        V7P3
!+        WAQPROCESS is now a prime number, so that multiple processes
!+        can be called by multiplication of the prime numbers.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_special
      USE declarations_telemac
      USE declarations_waqtel,ONLY:k2,rayeff,waqprocess,rayaed2,
     &                             deja_sw,deja_ca
      USE interface_waqtel, ex_point_waqtel => point_waqtel
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER,         INTENT(IN   ) :: IELM1
      TYPE(bief_mesh), INTENT(INOUT) :: MESH2D
      TYPE(bief_mesh), INTENT(INOUT),OPTIONAL :: MESH3D
      INTEGER,         INTENT(IN   ),OPTIONAL :: IELM3
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!
      deja_sw = .false.
      deja_ca = .false.
!
      IF( 2*int(waqprocess/2).EQ.waqprocess .OR.
     &    5*int(waqprocess/5).EQ.waqprocess ) THEN
        CALL bief_allvec(1,k2   ,'K2    ',ielm1,1,1,mesh2d)
      ELSE
        CALL bief_allvec(1,k2    ,'K2    ',    0,1,0,mesh2d)
      ENDIF
!
!     SUN RAY EFFECTS CAN BE 2D OR 3D
!
      IF( 3*int(waqprocess/3).EQ.waqprocess .OR.
     &    5*int(waqprocess/5).EQ.waqprocess ) THEN
        IF(PRESENT(mesh3d))THEN
          CALL bief_allvec(1,rayeff,'RAYEFF',ielm3,1,1,mesh3d)
        ELSE
          CALL bief_allvec(1,rayeff,'RAYEFF',ielm1,1,1,mesh2d)
        ENDIF
      ELSE
        IF(PRESENT(mesh3d))THEN
          CALL bief_allvec(1,rayeff,'RAYEFF',    0,1,0,mesh3d)
        ELSE
          CALL bief_allvec(1,rayeff,'RAYEFF',    0,1,0,mesh2d)
        ENDIF
      ENDIF
!
!     WINDX, WINDY, PATMOS, TAIR NOW ALLOCATED IN METEO_TELEMAC
!
!     WIND GIVEN IN P1
!
!     IF(.NOT.VENT) THEN
!       CALL BIEF_ALLVEC(1,WINDX,'WINDX ',IELM1,1,1,MESH2D)
!       CALL BIEF_ALLVEC(1,WINDY,'WINDY ',IELM1,1,1,MESH2D)
!     ENDIF
!
!     ATMOSPHERIC PRESSURE GIVEN IN P1
!
!     IF(.NOT.ATMOS) THEN
!       CALL BIEF_ALLVEC(1,PATMOS,'PATMOS',IELM1,1,1,MESH2D)
!     ENDIF
!
!     AIR TEMPERATURE GIVEN IN P1
!
!     CALL BIEF_ALLVEC(1,TAIR  ,'TAIR  ',IELM1,1,1,MESH2D)
!
!     SOLAR RADIATION FOR AED2
!
      IF(13*int(waqprocess/13).EQ.waqprocess) THEN
        CALL bief_allvec(1, rayaed2,  'RAYAED', ielm1,1,1,mesh2d)
      ELSE
        CALL bief_allvec(1, rayaed2,  'RAYAED', 0,1,0,mesh2d)
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END