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

Timestamp:

Aug 16, 2007 2:30:26 PM (17 years ago)

Author:

raasch

Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:

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

palm/trunk/SOURCE/data_output_dvrp.f90

@@ -151 +151 @@
 !
 !--       Definition of characteristics of particle material
-          CALL DVRP_MATERIAL_RGB( m-1, 1, 0.1, 0.7, 0.1, 0.0 )
+!          CALL DVRP_MATERIAL_RGB( m-1, 1, 0.1, 0.7, 0.1, 0.0 )
+          CALL DVRP_MATERIAL_RGB( m-1, 1, 0.0, 0.0, 0.0, 0.0 )
 
 !
@@ -318 +319 @@
                    ENDDO
                 ENDDO
-                DO  k = nzb, nz_do3d
-                   DO  j = nys+1, nyn
-                      DO  i = nxl, nxr+1
-                         local_pf(i,j,k) = 0.25 * local_pf(i,j-1,k) + &
-                                           0.50 * local_pf(i,j,k)   + &
-                                           0.25 * local_pf(i,j+1,k)
-                      ENDDO
-                   ENDDO
-                ENDDO
+! Averaging for Langmuir circulation
+!                DO  k = nzb, nz_do3d
+!                   DO  j = nys+1, nyn
+!                      DO  i = nxl, nxr+1
+!                         local_pf(i,j,k) = 0.25 * local_pf(i,j-1,k) + &
+!                                           0.50 * local_pf(i,j,k)   + &
+!                                           0.25 * local_pf(i,j+1,k)
+!                      ENDDO
+!                   ENDDO
+!                ENDDO
 
              CASE ( 'p', 'p_xy', 'p_xz', 'p_yz' )
@@ -475 +477 @@
              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 )
+!             CALL DVRP_SLICER( m-1, section_mode, slicer_position )
+             CALL DVRP_SLICER( m-1, 2, 1.0 )
+             WRITE (9,*) 'nx_dvrp=', nx_dvrp
+             WRITE (9,*) 'ny_dvrp=', ny_dvrp
+             WRITE (9,*) 'nz_dvrp=', nz_dvrp
+             WRITE (9,*) 'section_mode=', section_mode
+             WRITE (9,*) 'slicer_position=', slicer_position
+             CALL local_flush( 9 )
+
              CALL DVRP_VISUALIZE( m-1, 2, dvrp_filecount )