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

Last change on this file since 133 was 130, checked in by letzel, 16 years ago

DVRP output modifications:

  • The user can now visualize user-defined quantities using dvrp.

data_output_dvrp calls the new user_interface subroutine
user_data_output_dvrp in case of unknown variables (CASE DEFAULT).

  • Two instead of one digit are allowed to specify isosurface and slicer

variables with the parameter mode_dvrp.

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