Home

Context Navigation

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

Last change on this file since 208 was 206, checked in by raasch, 16 years ago
ocean-atmosphere coupling realized with MPI-1, adjustments in mrun, mbuild, subjob for lcxt4
Property svn:keywords set to `Id`
File size: 21.5 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |