Changeset 262 for palm/trunk/SOURCE/data_output_dvrp.f90

Timestamp:

Mar 18, 2009 8:32:37 AM (15 years ago)

Author:

raasch

Message:

further updates for dvr output, bugfix in advec_particles concerning particle boundary condition

File:

• palm/trunk/SOURCE/data_output_dvrp.f90 (modified) (4 diffs)

palm/trunk/SOURCE/data_output_dvrp.f90

@@ -32 +32 @@
 ! Current revisions:
 ! -----------------
-! Output of messages replaced by message handling routine.
+! Clipping of dvr-output implemented,
+! output of messages replaced by message handling routine.
 ! TEST: different colours for isosurfaces
 !
@@ -167 +168 @@
 !--       number of particles and tails to be plotted; otherwise, all
 !--       particles/tails are plotted. dvrp_mask is used to mark the partikles.
-          IF ( .NOT. use_particle_tails )  THEN
+!          IF ( .NOT. use_particle_tails )  THEN
              ALLOCATE( dvrp_mask(number_of_particles) )
-          ELSE
-             ALLOCATE( dvrp_mask(number_of_tails*maximum_number_of_tailpoints) )
-          ENDIF
+!          ELSE
+!             ALLOCATE( dvrp_mask(number_of_tails*maximum_number_of_tailpoints) )
+!          ENDIF
           dvrp_mask = .TRUE.
           IF ( dvrp_total_overlap )  THEN
@@ -202 +203 @@
                             dvrp_not = dvrp_not + 1
                          ELSE
-                            dvrp_mask(k) = .FALSE.
+                            dvrp_mask(n) = .FALSE.
                          ENDIF
                       ENDIF
@@ -532 +533 @@
              CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, &
                              cyclic_dvrp, cyclic_dvrp, cyclic_dvrp )
-!             CALL DVRP_SLICER( m-1, section_mode, slicer_position )
-             tmp_pos = 1.0
-             CALL DVRP_SLICER( m-1, 2, tmp_pos )
+             tmp_pos = slicer_position
+             CALL DVRP_SLICER( m-1, section_mode, tmp_pos )
+!             tmp_pos = 1.0
+!             CALL DVRP_SLICER( m-1, 2, tmp_pos )
              WRITE (9,*) 'nx_dvrp=', nx_dvrp
              WRITE (9,*) 'ny_dvrp=', ny_dvrp