Home

Context Navigation

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

Last change on this file since 70 was 39, checked in by raasch, 18 years ago
comments prepared for 3.1c
Property svn:keywords set to `Id`
File size: 23.5 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |