Ignore:
Timestamp:
Nov 18, 2008 9:54:31 AM (13 years ago)
Author:
raasch
Message:

precompilation mechanism completely revised: now one depository per configuration block, further change of output messages

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE

    • Property svn:ignore set to
      palm_sources.tar
  • palm/trunk/SOURCE/data_output_2d.f90

    r108 r215  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: zero output of particle concentration and radius unless particles
     7! have been started
    78!
    89! Former revisions:
     
    215216             CASE ( 'pc_xy', 'pc_xz', 'pc_yz' )  ! particle concentration
    216217                IF ( av == 0 )  THEN
    217                    tend = prt_count
    218                    CALL exchange_horiz( tend )
     218                   IF ( simulated_time >= particle_advection_start )  THEN
     219                      tend = prt_count
     220                      CALL exchange_horiz( tend )
     221                   ELSE
     222                      tend = 0.0
     223                   ENDIF
    219224                   DO  i = nxl-1, nxr+1
    220225                      DO  j = nys-1, nyn+1
     
    232237             CASE ( 'pr_xy', 'pr_xz', 'pr_yz' )  ! mean particle radius
    233238                IF ( av == 0 )  THEN
    234                    DO  i = nxl, nxr
    235                       DO  j = nys, nyn
    236                          DO  k = nzb, nzt+1
    237                             psi = prt_start_index(k,j,i)
    238                             s_r3 = 0.0
    239                             s_r4 = 0.0
    240                             DO  n = psi, psi+prt_count(k,j,i)-1
    241                                s_r3 = s_r3 + particles(n)%radius**3
    242                                s_r4 = s_r4 + particles(n)%radius**4
     239                   IF ( simulated_time >= particle_advection_start )  THEN
     240                      DO  i = nxl, nxr
     241                         DO  j = nys, nyn
     242                            DO  k = nzb, nzt+1
     243                               psi = prt_start_index(k,j,i)
     244                               s_r3 = 0.0
     245                               s_r4 = 0.0
     246                               DO  n = psi, psi+prt_count(k,j,i)-1
     247                                  s_r3 = s_r3 + particles(n)%radius**3
     248                                  s_r4 = s_r4 + particles(n)%radius**4
     249                               ENDDO
     250                               IF ( s_r3 /= 0.0 )  THEN
     251                                  mean_r = s_r4 / s_r3
     252                               ELSE
     253                                  mean_r = 0.0
     254                               ENDIF
     255                               tend(k,j,i) = mean_r
    243256                            ENDDO
    244                             IF ( s_r3 /= 0.0 )  THEN
    245                                mean_r = s_r4 / s_r3
    246                             ELSE
    247                                mean_r = 0.0
    248                             ENDIF
    249                             tend(k,j,i) = mean_r
    250257                         ENDDO
    251258                      ENDDO
    252                    ENDDO
    253                    CALL exchange_horiz( tend )
     259                      CALL exchange_horiz( tend )
     260                   ELSE
     261                      tend = 0.0
     262                   ENDIF
    254263                   DO  i = nxl-1, nxr+1
    255264                      DO  j = nys-1, nyn+1
Note: See TracChangeset for help on using the changeset viewer.