Changeset 388 for palm/trunk/SOURCE/flow_statistics.f90
- Timestamp:
- Sep 23, 2009 9:40:33 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/flow_statistics.f90
r343 r388 4 4 ! Current revisions: 5 5 ! ----------------- 6 ! Vertical profiles of potential density and hydrostatic pressure are 7 ! calculated. 6 8 ! Added missing timeseries calculation of w"q"(0), moved timeseries q* to the 7 9 ! end. … … 531 533 IF ( humidity ) THEN 532 534 sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + & 533 qswst(j,i) * rmask(j,i,sr) 535 qswst(j,i) * rmask(j,i,sr) ! w"q" (w"qv") 534 536 IF ( cloud_physics ) THEN 535 537 sums_l(nzt,45,tn) = sums_l(nzt,45,tn) + ( & … … 546 548 IF ( passive_scalar ) THEN 547 549 sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + & 548 qswst(j,i) * rmask(j,i,sr) 550 qswst(j,i) * rmask(j,i,sr) ! w"q" (w"qv") 549 551 ENDIF 550 552 ENDIF … … 597 599 rmask(j,i,sr) 598 600 sums_l(k,64,tn) = sums_l(k,64,tn) + rho(k,j,i) * & 601 rmask(j,i,sr) 602 sums_l(k,71,tn) = sums_l(k,71,tn) + prho(k,j,i) * & 599 603 rmask(j,i,sr) 600 604 ENDIF … … 644 648 ! 645 649 !-- Density at top follows Neumann condition 646 IF ( ocean ) sums_l(nzt+1,64,tn) = sums_l(nzt,64,tn) 650 IF ( ocean ) THEN 651 sums_l(nzt+1,64,tn) = sums_l(nzt,64,tn) 652 sums_l(nzt+1,71,tn) = sums_l(nzt,71,tn) 653 ENDIF 647 654 648 655 ! … … 897 904 hom(:,1,69,sr) = sums(:,69) ! w"e + w"p"/rho 898 905 hom(:,1,70,sr) = sums(:,70) ! q*2 906 hom(:,1,71,sr) = sums(:,71) ! prho 907 hom(:,1,72,sr) = hyp * 1E-4 ! hyp in kPa 899 908 900 909 hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1) … … 916 925 !-- is less than 1.5 times the height where the heat flux becomes negative 917 926 !-- (positive) for the first time. 918 !-- NOTE: This criterion is still capable of improving!919 927 z_i(1) = 0.0 920 928 first = .TRUE.
Note: See TracChangeset
for help on using the changeset viewer.