Changeset 699
 Timestamp:
 Mar 22, 2011 5:52:22 PM (13 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/flow_statistics.f90
r697 r699 4 4 ! Current revisions: 5 5 !  6 ! 6 ! Bugfix in calculation of vertical velocity skewness. The added absolute value 7 ! avoid negative values in the root. Negative values of w'w' can occur at the 8 ! top or bottom of the model domain due to degrading the order of advection 9 ! scheme. Furthermore the calculation will be the same for all advection 10 ! schemes. 7 11 ! 8 12 ! Former revisions: … … 437 441 ( q(k,j,i)hom(k,1,41,sr) )**2 * rmask(j,i,sr) 438 442 ENDIF 443 444 ! 445 ! Higher moments 446 ! (Computation of the skewness of w further below) 447 sums_l(k,38,tn) = sums_l(k,38,tn) + w(k,j,i)**3 * rmask(j,i,sr) 439 448 440 449 sums_l_etot = sums_l_etot + & … … 484 493 sums_l(k,32,tn) = sums_l(k,32,tn) + w2 * rmask(j,i,sr) 485 494 ! 486 ! Higher moments487 ! (Computation of the skewness of w further below)488 sums_l(k,38,tn) = sums_l(k,38,tn) + w(k,j,i) * w2 * &489 rmask(j,i,sr)490 !491 495 ! Perturbation energy 492 496 … … 501 505 ENDDO 502 506 ENDDO 503 ELSE 504 !$OMP DO 505 DO i = nxl, nxr 506 DO j = nys, nyn 507 DO k = nzb_s_inner(j,i), nzt + 1 508 w2 = w(k,j,i)**2 509 ! 510 ! Higher moments 511 ! (Computation of the skewness of w further below) 512 sums_l(k,38,tn) = sums_l(k,38,tn) + w(k,j,i) * w2 * & 513 rmask(j,i,sr) 514 ENDDO 515 ENDDO 516 ENDDO 517 ENDIF 518 507 ENDIF 519 508 ! 520 509 ! Horizontally averaged profiles of the vertical fluxes … … 1031 1020 hom(:,1,37,sr) = sums(:,37) ! w*e* 1032 1021 hom(:,1,38,sr) = sums(:,38) ! w*3 1033 hom(:,1,39,sr) = sums(:,38) / ( sums(:,32) + 1E20 )**1.5! Sw1022 hom(:,1,39,sr) = sums(:,38) / ( abs( sums(:,32) ) + 1E20 )**1.5 ! Sw 1034 1023 hom(:,1,40,sr) = sums(:,40) ! p 1035 1024 hom(:,1,45,sr) = sums(:,45) ! w"vpt"
Note: See TracChangeset
for help on using the changeset viewer.