Changeset 75 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Mar 22, 2007 9:54:05 AM (17 years ago)

Author:

raasch

Message:

preliminary update for changes concerning non-cyclic boundary conditions

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (10 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -5 +5 @@
 ! -----------------
 ! Collection of time series quantities moved from routine flow_statistics to
-! here, routine user_statistics is called for each statistic region
+! here, routine user_statistics is called for each statistic region,
+! moisture renamed humidity
 !
 ! Former revisions:
@@ -115 +116 @@
 !--    total water content, specific humidity and liquid water potential
 !--    temperature
-       IF ( moisture )  THEN
+       IF ( humidity )  THEN
           !$OMP DO
           DO  i = nxl, nxr
@@ -164 +165 @@
              sums_l(:,2,0) = sums_l(:,2,0) + sums_l(:,2,i)
              sums_l(:,4,0) = sums_l(:,4,0) + sums_l(:,4,i)
-             IF ( moisture )  THEN
+             IF ( humidity )  THEN
                 sums_l(:,41,0) = sums_l(:,41,0) + sums_l(:,41,i)
                 sums_l(:,44,0) = sums_l(:,44,0) + sums_l(:,44,i)
@@ -187 +188 @@
        CALL MPI_ALLREDUCE( sums_l(nzb,4,0), sums(nzb,4), nzt+2-nzb, MPI_REAL, &
                            MPI_SUM, comm2d, ierr )
-       IF ( moisture ) THEN
+       IF ( humidity ) THEN
           CALL MPI_ALLREDUCE( sums_l(nzb,44,0), sums(nzb,44), nzt+2-nzb, &
                               MPI_REAL, MPI_SUM, comm2d, ierr )
@@ -208 +209 @@
        sums(:,2) = sums_l(:,2,0)
        sums(:,4) = sums_l(:,4,0)
-       IF ( moisture ) THEN
+       IF ( humidity ) THEN
           sums(:,44) = sums_l(:,44,0)
           sums(:,41) = sums_l(:,41,0)
@@ -231 +232 @@
 !
 !--    Humidity and cloud parameters
-       IF ( moisture ) THEN
+       IF ( humidity ) THEN
           sums(:,44) = sums(:,44) / ngp_2dh_outer(:,sr)
           sums(:,41) = sums(:,41) / ngp_2dh_outer(:,sr)
@@ -365 +366 @@
 !
 !--             Buoyancy flux, water flux (humidity flux) w"q"
-                IF ( moisture ) THEN
+                IF ( humidity ) THEN
                    sums_l(k,45,tn) = sums_l(k,45,tn)                           &
                                          - 0.5 * ( kh(k,j,i) + kh(k+1,j,i) )   &
@@ -407 +408 @@
                 sums_l(nzb,61,tn) = sums_l(nzb,61,tn) + &
                                     0.0 * rmask(j,i,sr)           ! v"pt"
-                IF ( moisture )  THEN
+                IF ( humidity )  THEN
                    sums_l(nzb,48,tn) = sums_l(nzb,48,tn) + &
                                        qsws(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
@@ -437 +438 @@
                 sums_l(nzt,61,tn) = sums_l(nzt,61,tn) + &
                                     0.0 * rmask(j,i,sr)           ! v"pt"
-                IF ( moisture )  THEN
+                IF ( humidity )  THEN
                    sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
                                        qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
@@ -499 +500 @@
 !--             Buoyancy flux, water flux, humidity flux and liquid water
 !--             content
-                IF ( moisture )  THEN
+                IF ( humidity )  THEN
                    pts = 0.5 * ( vpt(k,j,i)   - hom(k,1,44,sr) + &
                                  vpt(k+1,j,i) - hom(k+1,1,44,sr) )