Changeset 671 for palm


Ignore:
Timestamp:
Jan 11, 2011 12:04:00 PM (14 years ago)
Author:
heinze
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/subsidence.f90

    r668 r671  
    44! Current revisions:
    55! -----------------
     6! bugfix: access to ddzu(nzt+2) which is not defined
    67!
    78! Former revisions:
     
    167168         IF ( w_subs(k) >= 0.0 ) THEN  ! large-scale ascent
    168169            var_mod(k) = var_init(k) - dt_3d * w_subs(k) *  &
    169                                ( var_init(k) - var_init(k-1) ) * ddzu(k+1)
     170                               ( var_init(k) - var_init(k-1) ) * ddzu(k)
    170171         ENDIF
    171172      ENDDO
     
    241242             IF ( w_subs(k) >= 0.0 ) THEN  ! large-scale ascent
    242243                var_mod(k) = var_init(k) - dt_3d * w_subs(k) *  &
    243                                    ( var_init(k) - var_init(k-1) ) * ddzu(k+1)
     244                                   ( var_init(k) - var_init(k-1) ) * ddzu(k)
    244245             ENDIF
    245246          ENDDO
Note: See TracChangeset for help on using the changeset viewer.