The TELEMAC-MASCARET system  trunk
point_waqtel.f
Go to the documentation of this file.
1 ! ***********************
2  SUBROUTINE point_waqtel
3 ! ***********************
4 !
5  &(mesh2d,ielm1,mesh3d,ielm3)
6 !
7 !***********************************************************************
8 ! WAQTEL V8P2
9 !***********************************************************************
10 !
11 !brief Memory allocation of structures, aliases, blocks...
12 !
13 !history R ATA (LNHE)
14 !+ 13/05/2015
15 !+ V7P1
16 !+
17 !+ Creation of the file
18 !
19 !history S.E. BOURBAN (HRW)
20 !+ 21/09/2017
21 !+ V7P3
22 !+ WAQPROCESS is now a prime number, so that multiple processes
23 !+ can be called by multiplication of the prime numbers.
24 !
25 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
26 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
27 !
28  USE bief
32  & deja_sw,deja_ca
33  USE interface_waqtel, ex_point_waqtel => point_waqtel
34 !
35  IMPLICIT NONE
36 !
37 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
38 !
39  INTEGER, INTENT(IN ) :: IELM1
40  TYPE(bief_mesh), INTENT(INOUT) :: MESH2D
41  TYPE(bief_mesh), INTENT(INOUT),OPTIONAL :: MESH3D
42  INTEGER, INTENT(IN ),OPTIONAL :: IELM3
43 !
44 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
45 !
46 !
47  deja_sw = .false.
48  deja_ca = .false.
49 !
50  IF( 2*int(waqprocess/2).EQ.waqprocess .OR.
51  & 5*int(waqprocess/5).EQ.waqprocess ) THEN
52  CALL bief_allvec(1,k2 ,'K2 ',ielm1,1,1,mesh2d)
53  ELSE
54  CALL bief_allvec(1,k2 ,'K2 ', 0,1,0,mesh2d)
55  ENDIF
56 !
57 ! SUN RAY EFFECTS CAN BE 2D OR 3D
58 !
59  IF( 3*int(waqprocess/3).EQ.waqprocess .OR.
60  & 5*int(waqprocess/5).EQ.waqprocess ) THEN
61  IF(PRESENT(mesh3d))THEN
62  CALL bief_allvec(1,rayeff,'RAYEFF',ielm3,1,1,mesh3d)
63  ELSE
64  CALL bief_allvec(1,rayeff,'RAYEFF',ielm1,1,1,mesh2d)
65  ENDIF
66  ELSE
67  IF(PRESENT(mesh3d))THEN
68  CALL bief_allvec(1,rayeff,'RAYEFF', 0,1,0,mesh3d)
69  ELSE
70  CALL bief_allvec(1,rayeff,'RAYEFF', 0,1,0,mesh2d)
71  ENDIF
72  ENDIF
73 !
74 ! WINDX, WINDY, PATMOS, TAIR NOW ALLOCATED IN METEO_TELEMAC
75 !
76 ! WIND GIVEN IN P1
77 !
78 ! IF(.NOT.VENT) THEN
79 ! CALL BIEF_ALLVEC(1,WINDX,'WINDX ',IELM1,1,1,MESH2D)
80 ! CALL BIEF_ALLVEC(1,WINDY,'WINDY ',IELM1,1,1,MESH2D)
81 ! ENDIF
82 !
83 ! ATMOSPHERIC PRESSURE GIVEN IN P1
84 !
85 ! IF(.NOT.ATMOS) THEN
86 ! CALL BIEF_ALLVEC(1,PATMOS,'PATMOS',IELM1,1,1,MESH2D)
87 ! ENDIF
88 !
89 ! AIR TEMPERATURE GIVEN IN P1
90 !
91 ! CALL BIEF_ALLVEC(1,TAIR ,'TAIR ',IELM1,1,1,MESH2D)
92 !
93 ! SOLAR RADIATION FOR AED2
94 !
95  IF(13*int(waqprocess/13).EQ.waqprocess) THEN
96  CALL bief_allvec(1, rayaed2, 'RAYAED', ielm1,1,1,mesh2d)
97  ELSE
98  CALL bief_allvec(1, rayaed2, 'RAYAED', 0,1,0,mesh2d)
99  ENDIF
100 !
101 !-----------------------------------------------------------------------
102 !
103  RETURN
104  END
subroutine bief_allvec(NAT, VEC, NOM, IELM, DIM2, STATUT, MESH)
Definition: bief_allvec.f:7
type(bief_obj), target k2
subroutine point_waqtel(MESH2D, IELM1, MESH3D, IELM3)
Definition: point_waqtel.f:7
type(bief_obj), target rayaed2
type(bief_obj), target rayeff
Definition: bief.f:3