Ignore:
Timestamp:
Feb 26, 2007 12:12:32 AM (15 years ago)
Author:
raasch
Message:

new parameter use_top_fluxes, Bugfix: ddzw dimensioned 1:nzt+1

File:
1 edited

Legend:

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

    r19 r20  
    44! Actual revisions:
    55! -----------------
     6! Bugfix: ddzw dimensioned 1:nzt"+1"
    67! Calculation extended for gridpoint nzt, fluxes can be given at top,
    7 ! +s_flux_t in parameter list, s_flux renamed s_flux_b
     8! +s_flux_t in parameter list, s_flux renamed s_flux_b,
    89!
    910! Former revisions:
     
    5051       INTEGER ::  i, j, k
    5152       REAL    ::  vertical_gridspace
    52        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt)
     53       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1)
    5354       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
    5455       REAL, DIMENSION(:,:),   POINTER ::  s_flux_b, s_flux_t
     
    134135                tend(k,j,i) = tend(k,j,i)                                     &
    135136                                       + ( - s_flux_t(j,i)                    &
    136                                          - 0.5 * ( kh(k-1,j,i)+kh(k,j,i) )    &
    137                                                * ( s(k,j,i)-s(k-1,j,i) )      &
    138                                                * ddzu(k)                      &
     137                                           - 0.5 * ( kh(k-1,j,i)+kh(k,j,i) )  &
     138                                                 * ( s(k,j,i)-s(k-1,j,i) )    &
     139                                                 * ddzu(k)                    &
    139140                                         ) * ddzw(k)
    140141
     
    161162       INTEGER ::  i, j, k
    162163       REAL    ::  vertical_gridspace
    163        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt)
     164       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1)
    164165       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
    165166       REAL, DIMENSION(:,:),   POINTER ::  s_flux_b, s_flux_t
Note: See TracChangeset for help on using the changeset viewer.