Home

Context Navigation

source: palm/trunk/SOURCE/global_min_max.f90 @ 4216

Last change on this file since 4216 was 4182, checked in by scharf, 5 years ago
corrected "Former revisions" section minor formatting in "Former revisions" section added "Author" section
Property svn:keywords set to `Id`
File size: 10.2 KB

Line
1	!> @file global_min_max.f90
2	!------------------------------------------------------------------------------!
3	! This file is part of the PALM model system.
4	!
5	! PALM is free software: you can redistribute it and/or modify it under the
6	! terms of the GNU General Public License as published by the Free Software
7	! Foundation, either version 3 of the License, or (at your option) any later
8	! version.
9	!
10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
13	!
14	! You should have received a copy of the GNU General Public License along with
15	! PALM. If not, see <http://www.gnu.org/licenses/>.
16	!
17	! Copyright 1997-2019 Leibniz Universitaet Hannover
18	!------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! ------------------
22	!
23	!
24	! Former revisions:
25	! -----------------
26	! $Id: global_min_max.f90 4182 2019-08-22 15:20:23Z suehring $
27	! Corrected "Former revisions" section
28	!
29	! 3655 2019-01-07 16:51:22Z knoop
30	! Corrected "Former revisions" section
31	!
32	! Revision 1.1 1997/07/24 11:14:03 raasch
33	! Initial revision
34	!
35	!
36	! Description:
37	! ------------
38	!> Determine the array minimum/maximum and the corresponding indices.
39	!------------------------------------------------------------------------------!
40	SUBROUTINE global_min_max( i1, i2, j1, j2, k1, k2, ar, mode, offset, value, &
41	value_ijk, value1, value1_ijk )
42
43
44	USE indices, &
45	ONLY: nbgp, ny, nx
46
47	USE kinds
48
49	USE pegrid
50
51	IMPLICIT NONE
52
53	CHARACTER (LEN=*) :: mode !<
54
55	INTEGER(iwp) :: i !<
56	INTEGER(iwp) :: i1 !<
57	INTEGER(iwp) :: i2 !<
58	INTEGER(iwp) :: id_fmax !<
59	INTEGER(iwp) :: id_fmin !<
60	INTEGER(iwp) :: j !<
61	INTEGER(iwp) :: j1 !<
62	INTEGER(iwp) :: j2 !<
63	INTEGER(iwp) :: k !<
64	INTEGER(iwp) :: k1 !<
65	INTEGER(iwp) :: k2 !<
66	INTEGER(iwp) :: fmax_ijk(3) !<
67	INTEGER(iwp) :: fmax_ijk_l(3) !<
68	INTEGER(iwp) :: fmin_ijk(3) !<
69	INTEGER(iwp) :: fmin_ijk_l(3) !<
70	INTEGER(iwp) :: value_ijk(3) !<
71
72	INTEGER(iwp), OPTIONAL :: value1_ijk(3) !<
73
74	REAL(wp) :: offset !<
75	REAL(wp) :: value !<
76	REAL(wp) :: ar(i1:i2,j1:j2,k1:k2) !<
77
78	#if defined( __ibm )
79	REAL(sp) :: fmax(2) !<
80	REAL(sp) :: fmax_l(2) !<
81	REAL(sp) :: fmin(2) !<
82	REAL(sp) :: fmin_l(2) !<
83	! on 32bit-machines MPI_2REAL must not be replaced
84	! by MPI_2DOUBLE_PRECISION
85	#else
86	REAL(wp) :: fmax(2) !<
87	REAL(wp) :: fmax_l(2) !<
88	REAL(wp) :: fmin(2) !<
89	REAL(wp) :: fmin_l(2) !<
90	#endif
91	REAL(wp), OPTIONAL :: value1 !<
92
93
94	!
95	!-- Determine array minimum
96	IF ( mode == 'min' .OR. mode == 'minmax' ) THEN
97
98	!
99	!-- Determine the local minimum
100	fmin_ijk_l = MINLOC( ar )
101	fmin_ijk_l(1) = i1 + fmin_ijk_l(1) - 1 ! MINLOC assumes lowerbound = 1
102	fmin_ijk_l(2) = j1 + fmin_ijk_l(2) - nbgp
103	fmin_ijk_l(3) = k1 + fmin_ijk_l(3) - nbgp
104	fmin_l(1) = ar(fmin_ijk_l(1),fmin_ijk_l(2),fmin_ijk_l(3))
105
106	#if defined( __parallel )
107	fmin_l(2) = myid
108	IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr )
109	CALL MPI_ALLREDUCE( fmin_l, fmin, 1, MPI_2REAL, MPI_MINLOC, comm2d, &
110	ierr )
111
112	!
113	!-- Determine the global minimum. Result stored on PE0.
114	id_fmin = fmin(2)
115	IF ( id_fmin /= 0 ) THEN
116	IF ( myid == 0 ) THEN
117	CALL MPI_RECV( fmin_ijk, 3, MPI_INTEGER, id_fmin, 0, comm2d, &
118	status, ierr )
119	ELSEIF ( myid == id_fmin ) THEN
120	CALL MPI_SEND( fmin_ijk_l, 3, MPI_INTEGER, 0, 0, comm2d, ierr )
121	ENDIF
122	ELSE
123	fmin_ijk = fmin_ijk_l
124	ENDIF
125	!
126	!-- Send the indices of the just determined array minimum to other PEs
127	CALL MPI_BCAST( fmin_ijk, 3, MPI_INTEGER, 0, comm2d, ierr )
128	#else
129	fmin(1) = fmin_l(1)
130	fmin_ijk = fmin_ijk_l
131	#endif
132
133	ENDIF
134
135	!
136	!-- Determine array maximum
137	IF ( mode == 'max' .OR. mode == 'minmax' ) THEN
138
139	!
140	!-- Determine the local maximum
141	fmax_ijk_l = MAXLOC( ar )
142	fmax_ijk_l(1) = i1 + fmax_ijk_l(1) - 1 ! MAXLOC assumes lowerbound = 1
143	fmax_ijk_l(2) = j1 + fmax_ijk_l(2) - nbgp
144	fmax_ijk_l(3) = k1 + fmax_ijk_l(3) - nbgp
145	fmax_l(1) = ar(fmax_ijk_l(1),fmax_ijk_l(2),fmax_ijk_l(3))
146
147	#if defined( __parallel )
148	fmax_l(2) = myid
149	IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr )
150	CALL MPI_ALLREDUCE( fmax_l, fmax, 1, MPI_2REAL, MPI_MAXLOC, comm2d, &
151	ierr )
152
153	!
154	!-- Determine the global maximum. Result stored on PE0.
155	id_fmax = fmax(2)
156	IF ( id_fmax /= 0 ) THEN
157	IF ( myid == 0 ) THEN
158	CALL MPI_RECV( fmax_ijk, 3, MPI_INTEGER, id_fmax, 0, comm2d, &
159	status, ierr )
160	ELSEIF ( myid == id_fmax ) THEN
161	CALL MPI_SEND( fmax_ijk_l, 3, MPI_INTEGER, 0, 0, comm2d, ierr )
162	ENDIF
163	ELSE
164	fmax_ijk = fmax_ijk_l
165	ENDIF
166	!
167	!-- send the indices of the just determined array maximum to other PEs
168	CALL MPI_BCAST( fmax_ijk, 3, MPI_INTEGER, 0, comm2d, ierr )
169	#else
170	fmax(1) = fmax_l(1)
171	fmax_ijk = fmax_ijk_l
172	#endif
173
174	ENDIF
175
176	!
177	!-- Determine absolute array maximum
178	IF ( mode == 'abs' ) THEN
179
180	!
181	!-- Determine the local absolut maximum
182	fmax_l(1) = 0.0_wp
183	fmax_ijk_l(1) = i1
184	fmax_ijk_l(2) = j1
185	fmax_ijk_l(3) = k1
186	DO k = k1, k2
187	DO j = j1, j2
188	DO i = i1, i2
189	IF ( ABS( ar(i,j,k) ) > fmax_l(1) ) THEN
190	fmax_l(1) = ABS( ar(i,j,k) )
191	fmax_ijk_l(1) = i
192	fmax_ijk_l(2) = j
193	fmax_ijk_l(3) = k
194	ENDIF
195	ENDDO
196	ENDDO
197	ENDDO
198
199	!
200	!-- Set a flag in case that the determined value is negative.
201	!-- A constant offset has to be subtracted in order to handle the special
202	!-- case i=0 correctly
203	IF ( ar(fmax_ijk_l(1),fmax_ijk_l(2),fmax_ijk_l(3)) < 0.0_wp ) THEN
204	fmax_ijk_l(1) = -fmax_ijk_l(1) - 10
205	ENDIF
206
207	#if defined( __parallel )
208	fmax_l(2) = myid
209	IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr )
210	CALL MPI_ALLREDUCE( fmax_l, fmax, 1, MPI_2REAL, MPI_MAXLOC, comm2d, &
211	ierr )
212
213	!
214	!-- Determine the global absolut maximum. Result stored on PE0.
215	id_fmax = fmax(2)
216	IF ( id_fmax /= 0 ) THEN
217	IF ( myid == 0 ) THEN
218	CALL MPI_RECV( fmax_ijk, 3, MPI_INTEGER, id_fmax, 0, comm2d, &
219	status, ierr )
220	ELSEIF ( myid == id_fmax ) THEN
221	CALL MPI_SEND( fmax_ijk_l, 3, MPI_INTEGER, 0, 0, comm2d, ierr )
222	ENDIF
223	ELSE
224	fmax_ijk = fmax_ijk_l
225	ENDIF
226	!
227	!-- Send the indices of the just determined absolut maximum to other PEs
228	CALL MPI_BCAST( fmax_ijk, 3, MPI_INTEGER, 0, comm2d, ierr )
229	#else
230	fmax(1) = fmax_l(1)
231	fmax_ijk = fmax_ijk_l
232	#endif
233
234	ENDIF
235
236	!
237	!-- Determine absolute maximum of ( array - offset )
238	IF ( mode == 'absoff' ) THEN
239
240	!
241	!-- Determine the local absolut maximum
242	fmax_l(1) = 0.0_wp
243	fmax_ijk_l(1) = i1
244	fmax_ijk_l(2) = j1
245	fmax_ijk_l(3) = k1
246	DO k = k1, k2
247	DO j = j1, j2
248	!
249	!-- Attention: the lowest gridpoint is excluded here, because there
250	!-- --------- is no advection at nzb=0 and mode 'absoff' is only
251	!-- used for calculating u,v extrema for CFL-criteria
252	DO i = i1+1, i2
253	IF ( ABS( ar(i,j,k) - offset ) > fmax_l(1) ) THEN
254	fmax_l(1) = ABS( ar(i,j,k) - offset )
255	fmax_ijk_l(1) = i
256	fmax_ijk_l(2) = j
257	fmax_ijk_l(3) = k
258	ENDIF
259	ENDDO
260	ENDDO
261	ENDDO
262
263	!
264	!-- Set a flag in case that the determined value is negative.
265	!-- A constant offset has to be subtracted in order to handle the special
266	!-- case i=0 correctly
267	IF ( ar(fmax_ijk_l(1),fmax_ijk_l(2),fmax_ijk_l(3)) < 0.0_wp ) THEN
268	fmax_ijk_l(1) = -fmax_ijk_l(1) - 10
269	ENDIF
270
271	#if defined( __parallel )
272	fmax_l(2) = myid
273	IF ( collective_wait ) CALL MPI_BARRIER( comm2d, ierr )
274	CALL MPI_ALLREDUCE( fmax_l, fmax, 1, MPI_2REAL, MPI_MAXLOC, comm2d, &
275	ierr )
276
277	!
278	!-- Determine the global absolut maximum. Result stored on PE0.
279	id_fmax = fmax(2)
280	IF ( id_fmax /= 0 ) THEN
281	IF ( myid == 0 ) THEN
282	CALL MPI_RECV( fmax_ijk, 3, MPI_INTEGER, id_fmax, 0, comm2d, &
283	status, ierr )
284	ELSEIF ( myid == id_fmax ) THEN
285	CALL MPI_SEND( fmax_ijk_l, 3, MPI_INTEGER, 0, 0, comm2d, ierr )
286	ENDIF
287	ELSE
288	fmax_ijk = fmax_ijk_l
289	ENDIF
290	!
291	!-- Send the indices of the just determined absolut maximum to other PEs
292	CALL MPI_BCAST( fmax_ijk, 3, MPI_INTEGER, 0, comm2d, ierr )
293	#else
294	fmax(1) = fmax_l(1)
295	fmax_ijk = fmax_ijk_l
296	#endif
297
298	ENDIF
299
300	!
301	!-- Determine output parameters
302	SELECT CASE( mode )
303
304	CASE( 'min' )
305
306	value = fmin(1)
307	value_ijk = fmin_ijk
308
309	CASE( 'max' )
310
311	value = fmax(1)
312	value_ijk = fmax_ijk
313
314	CASE( 'minmax' )
315
316	value = fmin(1)
317	value_ijk = fmin_ijk
318	value1 = fmax(1)
319	value1_ijk = fmax_ijk
320
321	CASE( 'abs', 'absoff' )
322
323	value = fmax(1)
324	value_ijk = fmax_ijk
325	IF ( fmax_ijk(1) < 0 ) THEN
326	value = -value
327	value_ijk(1) = -value_ijk(1) - 10 !???
328	ENDIF
329
330	END SELECT
331
332	!
333	!-- Limit index values to the range 0..nx, 0..ny
334	IF ( value_ijk(3) < 0 ) value_ijk(3) = nx +1 + value_ijk(3)
335	IF ( value_ijk(3) > nx ) value_ijk(3) = value_ijk(3) - (nx+1)
336	IF ( value_ijk(2) < 0 ) value_ijk(2) = ny +1 + value_ijk(2)
337	IF ( value_ijk(2) > ny ) value_ijk(2) = value_ijk(2) - (ny+1)
338
339
340	END SUBROUTINE global_min_max

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |