Changeset 671 for palm

Timestamp:

Jan 11, 2011 12:04:00 PM (13 years ago)

Author:

heinze

Message:

bugfix: access to ddzu(nzt+2) which is not defined

File:

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

palm/trunk/SOURCE/subsidence.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! bugfix: access to ddzu(nzt+2) which is not defined
 !
 ! Former revisions:
@@ -167 +168 @@
          IF ( w_subs(k) >= 0.0 ) THEN  ! large-scale ascent
             var_mod(k) = var_init(k) - dt_3d * w_subs(k) *  &
-                               ( var_init(k) - var_init(k-1) ) * ddzu(k+1) 
+                               ( var_init(k) - var_init(k-1) ) * ddzu(k) 
          ENDIF
       ENDDO
@@ -241 +242 @@
              IF ( w_subs(k) >= 0.0 ) THEN  ! large-scale ascent
                 var_mod(k) = var_init(k) - dt_3d * w_subs(k) *  &
-                                   ( var_init(k) - var_init(k-1) ) * ddzu(k+1)
+                                   ( var_init(k) - var_init(k-1) ) * ddzu(k)
              ENDIF
           ENDDO