Changeset 1382 for palm/trunk/SOURCE/subsidence.f90

Timestamp:

Apr 30, 2014 12:15:41 PM (10 years ago)

Author:

boeske

Message:

minor changes in profile data output of lsf tendencies, variables renamed

File:

• palm/trunk/SOURCE/subsidence.f90 (modified) (6 diffs)

palm/trunk/SOURCE/subsidence.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Changed the weighting factor that is used in the summation of subsidence 
+! tendencies for profile data output from weight_pres to weight_substep
+! added Neumann boundary conditions for profile data output of subsidence terms
+! at nzt+1
 ! 
 ! Former revisions:
@@ -163 +166 @@
 
        USE statistics,                                                         &
-           ONLY:  sums_ls_l, weight_pres
+           ONLY:  sums_ls_l, weight_substep
 
        IMPLICIT NONE
@@ -186 +189 @@
        DO  i = nxl, nxr
           DO  j = nys, nyn
+
              DO  k = nzb_s_inner(j,i)+1, nzt 
                 IF ( w_subs(k) < 0.0_wp )  THEN    ! large-scale subsidence
@@ -199 +203 @@
                 IF ( large_scale_forcing )  THEN
                    sums_ls_l(k,ls_index) = sums_ls_l(k,ls_index) + tmp_tend    &
-                                      * weight_pres(intermediate_timestep_count)
+                                   * weight_substep(intermediate_timestep_count)
                 ENDIF
              ENDDO
+
+             sums_ls_l(nzt+1,ls_index) = sums_ls_l(nzt,ls_index)
+
           ENDDO
        ENDDO
@@ -259 +266 @@
 
        USE statistics,                                                         &
-           ONLY:  sums_ls_l, weight_pres
+           ONLY:  sums_ls_l, weight_substep
 
        IMPLICIT NONE
@@ -291 +298 @@
           IF ( large_scale_forcing )  THEN
              sums_ls_l(k,ls_index) = sums_ls_l(k,ls_index) + tmp_tend          &
-                                     * weight_pres(intermediate_timestep_count)
+                                   * weight_substep(intermediate_timestep_count)
           ENDIF
        ENDDO
+
+       sums_ls_l(nzt+1,ls_index) = sums_ls_l(nzt,ls_index)