Home

Context Navigation

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

Last change on this file since 235 was 226, checked in by raasch, 16 years ago
preparations for the next release
Property svn:keywords set to `Id`
File size: 23.5 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |