Changeset 20 for palm/trunk/SOURCE/diffusion_s.f90
- Timestamp:
- Feb 26, 2007 12:12:32 AM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diffusion_s.f90
r19 r20 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Bugfix: ddzw dimensioned 1:nzt"+1" 6 7 ! 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, 8 9 ! 9 10 ! Former revisions: … … 50 51 INTEGER :: i, j, k 51 52 REAL :: vertical_gridspace 52 REAL :: ddzu(1:nzt+1), ddzw(1:nzt )53 REAL :: ddzu(1:nzt+1), ddzw(1:nzt+1) 53 54 REAL :: tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) 54 55 REAL, DIMENSION(:,:), POINTER :: s_flux_b, s_flux_t … … 134 135 tend(k,j,i) = tend(k,j,i) & 135 136 + ( - 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) & 139 140 ) * ddzw(k) 140 141 … … 161 162 INTEGER :: i, j, k 162 163 REAL :: vertical_gridspace 163 REAL :: ddzu(1:nzt+1), ddzw(1:nzt )164 REAL :: ddzu(1:nzt+1), ddzw(1:nzt+1) 164 165 REAL :: tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) 165 166 REAL, DIMENSION(:,:), POINTER :: s_flux_b, s_flux_t
Note: See TracChangeset
for help on using the changeset viewer.