Home

Context Navigation

source: palm/trunk/SOURCE/init_dvrp.f90 @ 164

Last change on this file since 164 was 155, checked in by letzel, 17 years ago
Introduce prefix_chr to ensure unique dvrp_file path.
Property svn:keywords set to `Id`
File size: 21.1 KB

Line
1	SUBROUTINE init_dvrp
2
3	!------------------------------------------------------------------------------!
4	! Actual revisions:
5	! -----------------
6	!
7	! introduce prefix_chr to ensure unique dvrp_file path
8	! TEST: print* statements
9	! ToDo: checking of mode_dvrp for legal values is not correct
10	!
11	! Former revisions:
12	! -----------------
13	! $Id: init_dvrp.f90 155 2008-03-28 10:56:30Z raasch $
14	!
15	! 130 2007-11-13 14:08:40Z letzel
16	! allow two instead of one digit to specify isosurface and slicer variables
17	! Test output of isosurface on camera file
18	!
19	! 82 2007-04-16 15:40:52Z raasch
20	! Preprocessor strings for different linux clusters changed to "lc",
21	! routine local_flush is used for buffer flushing
22	!
23	! 17 2007-02-19 01:57:39Z raasch
24	! dvrp_output_local activated for all streams
25	!
26	! 13 2007-02-14 12:15:07Z raasch
27	! RCS Log replace by Id keyword, revision history cleaned up
28	!
29	! Revision 1.12 2006/02/23 12:30:22 raasch
30	! ebene renamed section, pl.. replaced by do..,
31	!
32	! Revision 1.1 2000/04/27 06:24:39 raasch
33	! Initial revision
34	!
35	!
36	! Description:
37	! ------------
38	! Initializing actions needed when using dvrp-software
39	!------------------------------------------------------------------------------!
40	#if defined( __dvrp_graphics )
41
42	USE arrays_3d
43	USE DVRP
44	USE dvrp_variables
45	USE grid_variables
46	USE indices
47	USE pegrid
48	USE control_parameters
49
50	IMPLICIT NONE
51
52	CHARACTER (LEN=2) :: section_chr
53	CHARACTER (LEN=3) :: prefix_chr
54	CHARACTER (LEN=80) :: dvrp_file_local
55	INTEGER :: i, j, k, l, m, pn, tv, vn
56	LOGICAL :: allocated
57	REAL :: center(3), distance
58
59	REAL, DIMENSION(:,:,:), ALLOCATABLE :: local_pf
60
61	TYPE(CSTRING), SAVE :: dvrp_directory_c, dvrp_file_c, &
62	dvrp_file_local_c,dvrp_host_c, &
63	dvrp_password_c, dvrp_username_c, name_c
64
65	!
66	!-- Set the maximum time the program can be suspended on user request (by
67	!-- dvrp steering). This variable is defined in module DVRP.
68	DVRP_MAX_SUSPEND_TIME = 7200
69
70	!
71	!-- Allocate array holding the names and limits of the steering variables
72	!-- (must have the same number of elements as array mode_dvrp!)
73	ALLOCATE( steering_dvrp(10) )
74
75	!
76	!-- Check, if output parameters are given and/or allowed
77	!-- and set default-values, where necessary
78	IF ( dvrp_username == ' ' ) THEN
79	IF ( myid == 0 ) THEN
80	PRINT*, '+++ init_dvrp: dvrp_username is undefined'
81	CALL local_stop
82	ENDIF
83	ENDIF
84
85	IF ( dvrp_output /= 'ftp' .AND. dvrp_output /= 'rtsp' .AND. &
86	dvrp_output /= 'local' ) THEN
87	IF ( myid == 0 ) THEN
88	PRINT*, '+++ init_dvrp: dvrp_output="', dvrp_output, '" not allowed'
89	CALL local_stop
90	ENDIF
91	ENDIF
92
93	IF ( dvrp_directory == 'default' ) THEN
94	dvrp_directory = TRIM( dvrp_username ) // '/' // TRIM( run_identifier )
95	ENDIF
96
97	IF ( dvrp_output /= 'local' ) THEN
98	IF ( dvrp_file /= 'default' .AND. dvrp_file /= '/dev/null' ) THEN
99	IF ( myid == 0 ) THEN
100	PRINT*, '+++ init_dvrp: dvrp_file="', dvrp_file, '" not allowed'
101	CALL local_stop
102	ENDIF
103	ENDIF
104	ENDIF
105
106	!
107	!-- Strings are assigned to strings of special type which have a CHAR( 0 )
108	!-- (C end-of-character symbol) at their end. This is needed when strings are
109	!-- passed to C routines.
110	dvrp_directory_c = dvrp_directory
111	dvrp_file_c = dvrp_file
112	dvrp_host_c = dvrp_host
113	dvrp_password_c = dvrp_password
114	dvrp_username_c = dvrp_username
115
116	!
117	!-- Loop over all output modes choosed
118	m = 1
119	allocated = .FALSE.
120	DO WHILE ( mode_dvrp(m) /= ' ' )
121
122	!
123	!-- Check, if mode is allowed
124	IF ( mode_dvrp(m)(1:10) /= 'isosurface' .AND. &
125	mode_dvrp(m)(1:6) /= 'slicer' .AND. &
126	mode_dvrp(m)(1:9) /= 'particles' ) THEN
127
128	IF ( myid == 0 ) THEN
129	PRINT*, '+++ init_dvrp: mode_dvrp="', mode_dvrp, '" not allowed'
130	ENDIF
131	CALL local_stop
132
133	ENDIF
134	!
135	!-- Determine prefix for dvrp_file
136	WRITE ( prefix_chr, '(I2.2,''_'')' ) m
137	!
138	!-- Camera position must be computed and written on file when no dvrp-output
139	!-- has been generated so far (in former runs)
140	! IF ( dvrp_filecount == 0 ) THEN
141	!
142	!-- Compute center of domain and distance of camera from center
143	center(1) = ( nx + 1.0 ) * dx * 0.5 * superelevation_x
144	center(2) = ( ny + 1.0 ) * dy * 0.5 * superelevation_y
145	center(3) = ( zu(nz_do3d) - zu(nzb) ) * 0.5 * superelevation
146	distance = 1.5 * MAX( ( nx + 1.0 ) * dx * superelevation_x, &
147	( ny + 1.0 ) * dy * superelevation_y, &
148	( zu(nz_do3d) - zu(nzb) ) * superelevation )
149
150	!
151	!-- Write camera position on file
152	CALL DVRP_INIT( m-1, 0 )
153
154	!
155	!-- Create filename for camera
156	IF ( dvrp_output == 'rtsp' ) THEN
157
158	WRITE ( 9, * ) '*** vor dvrp_output_rtsp'
159	CALL local_flush( 9 )
160
161	dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) // '/camera.dvr'
162	dvrp_file_c = dvrp_file
163	CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host_c, dvrp_username_c, &
164	dvrp_password_c, dvrp_directory_c, &
165	dvrp_file_c )
166	WRITE ( 9, * ) '*** nach dvrp_output_rtsp'
167	CALL local_flush( 9 )
168
169	ELSEIF ( dvrp_output == 'ftp' ) THEN
170
171	dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) // '.camera.dvr'
172	dvrp_file_c = dvrp_file
173	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host_c, dvrp_username_c, &
174	! dvrp_password_c, dvrp_directory_c, &
175	! dvrp_file_c )
176
177	ELSE
178
179	IF ( dvrp_file(1:9) /= '/dev/null' ) THEN
180	dvrp_file_local = prefix_chr // TRIM( mode_dvrp(m) ) &
181	// '.camera.dvr'
182	dvrp_file_local_c = dvrp_file_local
183	ELSE
184	dvrp_file_local_c = dvrp_file_c
185	ENDIF
186	CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )
187
188	ENDIF
189
190	CALL DVRP_CAMERA( m-1, center, distance )
191	WRITE ( 9, * ) '*** #1'
192	CALL local_flush( 9 )
193
194	!
195	!-- Define bounding box material and create a bounding box
196	CALL DVRP_MATERIAL_RGB( m-1, 1, 0.5, 0.5, 0.5, 0.0 )
197	CALL DVRP_BOUNDINGBOX( m-1, 1, 0.01, 0.0, 0.0, 0.0, &
198	(nx+1) * dx * superelevation_x, &
199	(ny+1) * dy * superelevation_y, &
200	zu(nz_do3d) * superelevation )
201
202	CALL DVRP_VISUALIZE( m-1, 0, 0 )
203	CALL DVRP_EXIT( m-1 )
204	WRITE ( 9, * ) '*** #2'
205	CALL local_flush( 9 )
206
207
208	!
209	!-- Write topography isosurface on file
210	CALL DVRP_INIT( m-1, 0 )
211
212	!
213	!-- Create filename for buildings
214	IF ( dvrp_output == 'rtsp' ) THEN
215
216	dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) &
217	// '/buildings.dvr'
218	dvrp_file_c = dvrp_file
219	CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host_c, dvrp_username_c, &
220	dvrp_password_c, dvrp_directory_c, &
221	dvrp_file_c )
222	WRITE ( 9, * ) '*** #3'
223	CALL local_flush( 9 )
224
225	ELSEIF ( dvrp_output == 'ftp' ) THEN
226
227	dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) &
228	// '.buildings.dvr'
229	dvrp_file_c = dvrp_file
230	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host_c, dvrp_username_c, &
231	! dvrp_password_c, dvrp_directory_c, &
232	! dvrp_file_c )
233
234	ELSE
235
236	IF ( dvrp_file(1:9) /= '/dev/null' ) THEN
237	dvrp_file_local = prefix_chr // TRIM( mode_dvrp(m) ) &
238	// '.buildings.dvr'
239	dvrp_file_local_c = dvrp_file_local
240	ELSE
241	dvrp_file_local_c = dvrp_file_c
242	ENDIF
243	CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )
244
245	ENDIF
246
247	!
248	!-- Determine local gridpoint coordinates
249	IF ( .NOT. allocated ) THEN
250	ALLOCATE( xcoor_dvrp(nxl:nxr+1), ycoor_dvrp(nys:nyn+1), &
251	zcoor_dvrp(nzb:nz_do3d) )
252	allocated = .TRUE.
253
254	DO i = nxl, nxr+1
255	xcoor_dvrp(i) = i * dx * superelevation_x
256	ENDDO
257	DO j = nys, nyn+1
258	ycoor_dvrp(j) = j * dy * superelevation_y
259	ENDDO
260	zcoor_dvrp = zu(nzb:nz_do3d) * superelevation
261	nx_dvrp = nxr+1 - nxl + 1
262	ny_dvrp = nyn+1 - nys + 1
263	nz_dvrp = nz_do3d - nzb + 1
264	ENDIF
265
266	!
267	!-- Define the grid used by dvrp
268	CALL DVRP_GRID( m-1, nx_dvrp, ny_dvrp, nz_dvrp, xcoor_dvrp, &
269	ycoor_dvrp, zcoor_dvrp )
270	CALL DVRP_MATERIAL_RGB( m-1, 1, 0.8, 0.7, 0.6, 0.0 )
271	WRITE ( 9, * ) '*** #4'
272	CALL local_flush( 9 )
273
274	!
275	!-- Compute and plot isosurface in dvr-format
276	ALLOCATE( local_pf(nxl:nxr+1,nys:nyn+1,nzb:nz_do3d) )
277	local_pf = 0.0
278	DO i = nxl, nxr+1
279	DO j = nys, nyn+1
280	IF ( nzb_s_inner(j,i) > 0 ) THEN
281	local_pf(i,j,nzb:nzb_s_inner(j,i)) = 1.0
282	ENDIF
283	ENDDO
284	ENDDO
285	WRITE ( 9, * ) '*** #4.1'
286	CALL local_flush( 9 )
287	CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, &
288	cyclic_dvrp, cyclic_dvrp, cyclic_dvrp )
289	WRITE ( 9, * ) '*** #4.2'
290	CALL local_flush( 9 )
291	CALL DVRP_THRESHOLD( m-1, 1.0 )
292	WRITE ( 9, * ) '*** #4.3'
293	CALL local_flush( 9 )
294	CALL DVRP_VISUALIZE( m-1, 1, 0 )
295	WRITE ( 9, * ) '*** #4.4'
296	CALL local_flush( 9 )
297
298	DEALLOCATE( local_pf )
299
300	CALL DVRP_EXIT( m-1 )
301	WRITE ( 9, * ) '*** #5'
302	CALL local_flush( 9 )
303
304	!
305	!-- Write the surface isosurface on file
306	CALL DVRP_INIT( m-1, 0 )
307
308	!
309	!-- Create filename for surface
310	IF ( dvrp_output == 'rtsp' ) THEN
311
312	dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) // '/surface.dvr'
313	dvrp_file_c = dvrp_file
314	CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host_c, dvrp_username_c, &
315	dvrp_password_c, dvrp_directory_c, &
316	dvrp_file_c )
317	WRITE ( 9, * ) '*** #6'
318	CALL local_flush( 9 )
319
320	ELSEIF ( dvrp_output == 'ftp' ) THEN
321
322	dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) // '.surface.dvr'
323	dvrp_file_c = dvrp_file
324	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host_c, dvrp_username_c, &
325	! dvrp_password_c, dvrp_directory_c, &
326	! dvrp_file_c )
327
328	ELSE
329
330	IF ( dvrp_file(1:9) /= '/dev/null' ) THEN
331	dvrp_file_local = prefix_chr // TRIM( mode_dvrp(m) ) &
332	// '.surface.dvr'
333	dvrp_file_local_c = dvrp_file_local
334	ELSE
335	dvrp_file_local_c = dvrp_file_c
336	ENDIF
337	CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )
338
339	ENDIF
340
341	!
342	!-- Determine local gridpoint coordinates
343	IF ( .NOT. allocated ) THEN
344	ALLOCATE( xcoor_dvrp(nxl:nxr+1), ycoor_dvrp(nys:nyn+1), &
345	zcoor_dvrp(nzb:nz_do3d) )
346	allocated = .TRUE.
347
348	DO i = nxl, nxr+1
349	xcoor_dvrp(i) = i * dx * superelevation_x
350	ENDDO
351	DO j = nys, nyn+1
352	ycoor_dvrp(j) = j * dy * superelevation_y
353	ENDDO
354	zcoor_dvrp = zu(nzb:nz_do3d) * superelevation
355	nx_dvrp = nxr+1 - nxl + 1
356	ny_dvrp = nyn+1 - nys + 1
357	nz_dvrp = nz_do3d - nzb + 1
358	ENDIF
359
360	!
361	!-- Define the grid used by dvrp
362	CALL DVRP_GRID( m-1, nx_dvrp, ny_dvrp, nz_dvrp, xcoor_dvrp, &
363	ycoor_dvrp, zcoor_dvrp )
364	CALL DVRP_MATERIAL_RGB( m-1, 1, 0.0, 0.6, 0.0, 0.0 )
365	WRITE ( 9, * ) '*** #7'
366	CALL local_flush( 9 )
367
368	!
369	!-- Compute and plot isosurface in dvr-format
370	ALLOCATE( local_pf(nxl:nxr+1,nys:nyn+1,nzb:nz_do3d) )
371	local_pf = 0.0
372	local_pf(:,:,0) = 1.0
373
374	CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, &
375	cyclic_dvrp, cyclic_dvrp, cyclic_dvrp )
376	CALL DVRP_THRESHOLD( m-1, 1.0 )
377	CALL DVRP_VISUALIZE( m-1, 1, 0 )
378
379	DEALLOCATE( local_pf )
380
381	CALL DVRP_EXIT( m-1 )
382	WRITE ( 9, * ) '*** #8'
383	CALL local_flush( 9 )
384
385
386	! ENDIF
387
388
389	!
390	!-- Initialize dvrp for all dvrp-calls during the run
391	CALL DVRP_INIT( m-1, 0 )
392
393	!
394	!-- Preliminary definition of filename for dvrp-output
395	IF ( dvrp_output == 'rtsp' ) THEN
396
397	!
398	!-- First initialize parameters for possible interactive steering.
399	!-- Every parameter has to be passed to the respective stream.
400	pn = 1
401	!
402	!-- Initialize threshold counter needed for initialization of the
403	!-- isosurface steering variables
404	tv = 0
405
406	DO WHILE ( mode_dvrp(pn) /= ' ' )
407
408	IF ( mode_dvrp(pn)(1:10) == 'isosurface' ) THEN
409
410	READ ( mode_dvrp(pn), '(10X,I2)' ) vn
411	steering_dvrp(pn)%name = do3d(0,vn)
412	tv = tv + 1
413
414	IF ( do3d(0,vn)(1:1) == 'w' ) THEN
415	steering_dvrp(pn)%min = -4.0
416	steering_dvrp(pn)%max = 5.0
417	ELSE
418	steering_dvrp(pn)%min = 288.0
419	steering_dvrp(pn)%max = 292.0
420	ENDIF
421
422	name_c = TRIM( do3d(0,vn) )
423	WRITE ( 9, * ) '*** #9'
424	CALL local_flush( 9 )
425	CALL DVRP_STEERING_INIT( m-1, name_c, steering_dvrp(pn)%min, &
426	steering_dvrp(pn)%max, threshold(tv) )
427	WRITE ( 9, * ) '*** #10'
428	CALL local_flush( 9 )
429
430	ELSEIF ( mode_dvrp(pn)(1:6) == 'slicer' ) THEN
431
432	READ ( mode_dvrp(pn), '(6X,I2)' ) vn
433	steering_dvrp(pn)%name = do2d(0,vn)
434	name_c = TRIM( do2d(0,vn) )
435
436	l = MAX( 2, LEN_TRIM( do2d(0,vn) ) )
437	section_chr = do2d(0,vn)(l-1:l)
438	SELECT CASE ( section_chr )
439	CASE ( 'xy' )
440	steering_dvrp(pn)%imin = 0
441	steering_dvrp(pn)%imax = nz_do3d
442	slicer_position_dvrp(pn) = section(1,1)
443	CALL DVRP_STEERING_INIT( m-1, name_c, &
444	steering_dvrp(pn)%imin, &
445	steering_dvrp(pn)%imax, &
446	slicer_position_dvrp(pn) )
447	CASE ( 'xz' )
448	steering_dvrp(pn)%imin = 0
449	steering_dvrp(pn)%imax = ny
450	slicer_position_dvrp(pn) = section(1,2)
451	CALL DVRP_STEERING_INIT( m-1, name_c, &
452	steering_dvrp(pn)%imin, &
453	steering_dvrp(pn)%imax, &
454	slicer_position_dvrp(pn) )
455	CASE ( 'yz' )
456	steering_dvrp(pn)%imin = 0
457	steering_dvrp(pn)%imax = nx
458	slicer_position_dvrp(pn) = section(1,3)
459	CALL DVRP_STEERING_INIT( m-1, name_c, &
460	steering_dvrp(pn)%imin, &
461	steering_dvrp(pn)%imax, &
462	slicer_position_dvrp(pn) )
463	END SELECT
464
465	ENDIF
466
467	pn = pn + 1
468
469	ENDDO
470
471	WRITE ( 9, * ) '*** #11'
472	CALL local_flush( 9 )
473
474	dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) // '/*****.dvr'
475	dvrp_file_c = dvrp_file
476	CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host_c, dvrp_username_c, &
477	dvrp_password_c, dvrp_directory_c, &
478	dvrp_file_c )
479	WRITE ( 9, * ) '*** #12'
480	CALL local_flush( 9 )
481
482	ELSEIF ( dvrp_output == 'ftp' ) THEN
483
484	dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) // '.%05d.dvr'
485	dvrp_file_c = dvrp_file
486	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host_c, dvrp_username_c, &
487	! dvrp_password_c, dvrp_directory_c, dvrp_file_c )
488
489	ELSE
490
491	IF ( dvrp_file(1:9) /= '/dev/null' ) THEN
492	dvrp_file_local = prefix_chr // TRIM( mode_dvrp(m) ) &
493	// '_%05d.dvr'
494	dvrp_file_local_c = dvrp_file_local
495	ELSE
496	dvrp_file_local_c = dvrp_file_c
497	ENDIF
498	CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )
499
500	ENDIF
501
502	! dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) // '.%05d.dvr' &
503	! // CHAR( 0 )
504	! dvrp_file = prefix_chr // TRIM( mode_dvrp(m) ) // '/*****.dvr' &
505	! // CHAR( 0 )
506	! dvrp_file = '/dev/null' // CHAR( 0 )
507	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host, dvrp_username, dvrp_password, &
508	! dvrp_directory, dvrp_file )
509	! CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host, dvrp_username, dvrp_password, &
510	! dvrp_directory, dvrp_file )
511	! CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file )
512
513	!
514	!-- Determine local gridpoint coordinates
515	IF ( .NOT. allocated ) THEN
516	ALLOCATE( xcoor_dvrp(nxl:nxr+1), ycoor_dvrp(nys:nyn+1), &
517	zcoor_dvrp(nzb:nz_do3d) )
518	allocated = .TRUE.
519
520	DO i = nxl, nxr+1
521	xcoor_dvrp(i) = i * dx * superelevation_x
522	ENDDO
523	DO j = nys, nyn+1
524	ycoor_dvrp(j) = j * dy * superelevation_y
525	ENDDO
526	zcoor_dvrp = zu(nzb:nz_do3d) * superelevation
527	nx_dvrp = nxr+1 - nxl + 1
528	ny_dvrp = nyn+1 - nys + 1
529	nz_dvrp = nz_do3d - nzb + 1
530	ENDIF
531
532	!
533	!-- Define the grid used by dvrp
534	WRITE ( 9, * ) '*** #13'
535	CALL local_flush( 9 )
536
537	CALL DVRP_GRID( m-1, nx_dvrp, ny_dvrp, nz_dvrp, xcoor_dvrp, ycoor_dvrp, &
538	zcoor_dvrp )
539	WRITE ( 9, * ) '*** #14'
540	CALL local_flush( 9 )
541
542
543	m = m + 1
544
545	ENDDO
546
547	#endif
548	END SUBROUTINE init_dvrp
549
550
551	SUBROUTINE init_dvrp_logging
552
553	!------------------------------------------------------------------------------!
554	! Description:
555	! ------------
556	! Initializes logging events for time measurement with dvrp software
557	! and splits one PE from the global communicator in case that dvrp output
558	! shall be done by one single PE.
559	!------------------------------------------------------------------------------!
560	#if defined( __dvrp_graphics )
561
562	USE dvrp_variables
563	USE pegrid
564
565	IMPLICIT NONE
566
567	CHARACTER (LEN=4) :: chr
568	INTEGER :: idummy
569
570	!
571	!-- Initialize logging of calls by DVRP graphic software
572	WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_INIT'
573	CALL local_flush( 9 )
574	CALL DVRP_LOG_INIT( 'DVRP_LOG' // CHAR( 0 ), 0 )
575	WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_INIT'
576	CALL local_flush( 9 )
577
578	!
579	!-- User-defined logging events: #1 (total time needed by PALM)
580	WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_SYMBOL'
581	CALL local_flush( 9 )
582	CALL DVRP_LOG_SYMBOL( 1, 'PALM_total' // CHAR( 0 ) )
583	WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_SYMBOL'
584	CALL local_flush( 9 )
585	CALL DVRP_LOG_SYMBOL( 2, 'PALM_timestep' // CHAR( 0 ) )
586	WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_EVENT'
587	CALL local_flush( 9 )
588	CALL DVRP_LOG_EVENT( 1, 1 )
589	WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_EVENT'
590	CALL local_flush( 9 )
591
592	#if defined( __parallel )
593	!
594	!-- Find out, if dvrp output shall be done by a dedicated PE
595	CALL local_getenv( 'use_seperate_pe_for_dvrp_output', 31, chr, idummy )
596	IF ( chr == 'true' ) THEN
597	use_seperate_pe_for_dvrp_output = .TRUE.
598	WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_SPLIT'
599	CALL local_flush( 9 )
600	CALL DVRP_SPLIT( MPI_COMM_WORLD, comm_palm )
601	WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_SPLIT'
602	CALL local_flush( 9 )
603	CALL MPI_COMM_SIZE( comm_palm, numprocs, ierr )
604	ENDIF
605	#endif
606
607	#endif
608	END SUBROUTINE init_dvrp_logging
609
610
611	SUBROUTINE close_dvrp
612
613	!------------------------------------------------------------------------------!
614	! Description:
615	! ------------
616	! Exit of dvrp software and finish dvrp logging
617	!------------------------------------------------------------------------------!
618	#if defined( __dvrp_graphics )
619
620	USE control_parameters
621	USE dvrp
622	USE dvrp_variables
623
624	INTEGER :: m
625
626	!
627	!-- If required, close dvrp-software and logging of dvrp-calls
628	IF ( dt_dvrp /= 9999999.9 ) THEN
629	m = 1
630	DO WHILE ( mode_dvrp(m) /= ' ' )
631	CALL DVRP_EXIT( m-1 )
632	m = m + 1
633	ENDDO
634	CALL DVRP_LOG_EVENT( -1, 1 ) ! Logging of total cpu-time used by PALM
635	IF ( use_seperate_pe_for_dvrp_output ) THEN
636	CALL DVRP_SPLIT_EXIT( 1 ) ! Argument 0: reduced output
637	ELSE
638	CALL DVRP_LOG_EXIT( 1 ) ! Argument 0: reduced output
639	ENDIF
640	ENDIF
641
642	#endif
643	END SUBROUTINE close_dvrp

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |