Ignore:
Timestamp:
Aug 16, 2007 2:30:26 PM (14 years ago)
Author:
raasch
Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/data_output_dvrp.f90

    r86 r106  
    151151!
    152152!--       Definition of characteristics of particle material
    153           CALL DVRP_MATERIAL_RGB( m-1, 1, 0.1, 0.7, 0.1, 0.0 )
     153!          CALL DVRP_MATERIAL_RGB( m-1, 1, 0.1, 0.7, 0.1, 0.0 )
     154          CALL DVRP_MATERIAL_RGB( m-1, 1, 0.0, 0.0, 0.0, 0.0 )
    154155
    155156!
     
    318319                   ENDDO
    319320                ENDDO
    320                 DO  k = nzb, nz_do3d
    321                    DO  j = nys+1, nyn
    322                       DO  i = nxl, nxr+1
    323                          local_pf(i,j,k) = 0.25 * local_pf(i,j-1,k) + &
    324                                            0.50 * local_pf(i,j,k)   + &
    325                                            0.25 * local_pf(i,j+1,k)
    326                       ENDDO
    327                    ENDDO
    328                 ENDDO
     321! Averaging for Langmuir circulation
     322!                DO  k = nzb, nz_do3d
     323!                   DO  j = nys+1, nyn
     324!                      DO  i = nxl, nxr+1
     325!                         local_pf(i,j,k) = 0.25 * local_pf(i,j-1,k) + &
     326!                                           0.50 * local_pf(i,j,k)   + &
     327!                                           0.25 * local_pf(i,j+1,k)
     328!                      ENDDO
     329!                   ENDDO
     330!                ENDDO
    329331
    330332             CASE ( 'p', 'p_xy', 'p_xz', 'p_yz' )
     
    475477             CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, &
    476478                             cyclic_dvrp, cyclic_dvrp, cyclic_dvrp )
    477              CALL DVRP_SLICER( m-1, section_mode, slicer_position )
     479!             CALL DVRP_SLICER( m-1, section_mode, slicer_position )
     480             CALL DVRP_SLICER( m-1, 2, 1.0 )
     481             WRITE (9,*) 'nx_dvrp=', nx_dvrp
     482             WRITE (9,*) 'ny_dvrp=', ny_dvrp
     483             WRITE (9,*) 'nz_dvrp=', nz_dvrp
     484             WRITE (9,*) 'section_mode=', section_mode
     485             WRITE (9,*) 'slicer_position=', slicer_position
     486             CALL local_flush( 9 )
     487
    478488             CALL DVRP_VISUALIZE( m-1, 2, dvrp_filecount )
    479489
Note: See TracChangeset for help on using the changeset viewer.