Changeset 1007 for palm/trunk/SOURCE/data_output_2d.f90
- Timestamp:
- Sep 19, 2012 2:30:36 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_2d.f90
r979 r1007 4 4 ! Current revisions: 5 5 ! ----------------- 6 ! 6 ! Bugfix: missing calculation of ql_vp added 7 7 ! 8 8 ! Former revisions: … … 428 428 CASE ( 'ql_vp_xy', 'ql_vp_xz', 'ql_vp_yz' ) 429 429 IF ( av == 0 ) THEN 430 IF ( simulated_time >= particle_advection_start ) THEN 431 DO i = nxl, nxr 432 DO j = nys, nyn 433 DO k = nzb, nzt+1 434 psi = prt_start_index(k,j,i) 435 DO n = psi, psi+prt_count(k,j,i)-1 436 tend(k,j,i) = tend(k,j,i) + & 437 particles(n)%weight_factor / & 438 prt_count(k,j,i) 439 ENDDO 440 ENDDO 441 ENDDO 442 ENDDO 443 CALL exchange_horiz( tend, nbgp ) 444 ELSE 445 tend = 0.0 446 END IF 447 DO i = nxlg, nxrg 448 DO j = nysg, nyng 449 DO k = nzb, nzt+1 450 local_pf(i,j,k) = tend(k,j,i) 451 ENDDO 452 ENDDO 453 ENDDO 454 resorted = .TRUE. 455 ELSE 456 CALL exchange_horiz( ql_vp_av, nbgp ) 430 457 to_be_resorted => ql_vp 431 ELSE432 to_be_resorted => ql_vp_av433 458 ENDIF 434 459 IF ( mode == 'xy' ) level_z = zu
Note: See TracChangeset
for help on using the changeset viewer.