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

Last change on this file since 231 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.