Changeset 187 for palm/trunk/SOURCE/production_e.f90
- Timestamp:
- Aug 6, 2008 4:25:09 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/production_e.f90
r139 r187 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! Change: add 'minus' sign to fluxes obtained from subroutine wall_fluxes_e for 7 ! consistency with subroutine wall_fluxes 7 8 ! 8 9 ! Former revisions: … … 164 165 CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, & 165 166 usvs, 1.0, 0.0, 0.0, 0.0 ) 166 dudy = wall_e_y(j,i) * usvs(k) / km(k,j,i)167 ! dudy = wall_e_y(j,i) * usvs(k,j,i) / km(k,j,i)167 dudy = - wall_e_y(j,i) * usvs(k) / km(k,j,i) 168 ! dudy = - wall_e_y(j,i) * usvs(k,j,i) / km(k,j,i) 168 169 CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, & 169 170 wsvs, 0.0, 0.0, 1.0, 0.0 ) 170 dwdy = wall_e_y(j,i) * wsvs(k) / km(k,j,i)171 ! dwdy = wall_e_y(j,i) * wsvs(k,j,i) / km(k,j,i)171 dwdy = - wall_e_y(j,i) * wsvs(k) / km(k,j,i) 172 ! dwdy = - wall_e_y(j,i) * wsvs(k,j,i) / km(k,j,i) 172 173 ELSE 173 174 dudy = 0.25 * ( u(k,j+1,i) + u(k,j+1,i+1) - & … … 180 181 CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, & 181 182 vsus, 0.0, 1.0, 0.0, 0.0 ) 182 dvdx = wall_e_x(j,i) * vsus(k) / km(k,j,i)183 ! dvdx = wall_e_x(j,i) * vsus(k,j,i) / km(k,j,i)183 dvdx = - wall_e_x(j,i) * vsus(k) / km(k,j,i) 184 ! dvdx = - wall_e_x(j,i) * vsus(k,j,i) / km(k,j,i) 184 185 CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, & 185 186 wsus, 0.0, 0.0, 0.0, 1.0 ) 186 dwdx = wall_e_x(j,i) * wsus(k) / km(k,j,i)187 ! dwdx = wall_e_x(j,i) * wsus(k,j,i) / km(k,j,i)187 dwdx = - wall_e_x(j,i) * wsus(k) / km(k,j,i) 188 ! dwdx = - wall_e_x(j,i) * wsus(k,j,i) / km(k,j,i) 188 189 ELSE 189 190 dvdx = 0.25 * ( v(k,j,i+1) + v(k,j+1,i+1) - & … … 220 221 221 222 IF ( wall_e_y(j,i) /= 0.0 ) THEN 222 dudy = wall_e_y(j,i) * usvs(k) / km(k,j,i)223 ! dudy = wall_e_y(j,i) * usvs(k,j,i) / km(k,j,i)224 dwdy = wall_e_y(j,i) * wsvs(k) / km(k,j,i)225 ! dwdy = wall_e_y(j,i) * wsvs(k,j,i) / km(k,j,i)223 dudy = - wall_e_y(j,i) * usvs(k) / km(k,j,i) 224 ! dudy = - wall_e_y(j,i) * usvs(k,j,i) / km(k,j,i) 225 dwdy = - wall_e_y(j,i) * wsvs(k) / km(k,j,i) 226 ! dwdy = - wall_e_y(j,i) * wsvs(k,j,i) / km(k,j,i) 226 227 ELSE 227 228 dudy = 0.25 * ( u(k,j+1,i) + u(k,j+1,i+1) - & … … 232 233 233 234 IF ( wall_e_x(j,i) /= 0.0 ) THEN 234 dvdx = wall_e_x(j,i) * vsus(k) / km(k,j,i)235 ! dvdx = wall_e_x(j,i) * vsus(k,j,i) / km(k,j,i)236 dwdx = wall_e_x(j,i) * wsus(k) / km(k,j,i)237 ! dwdx = wall_e_x(j,i) * wsus(k,j,i) / km(k,j,i)235 dvdx = - wall_e_x(j,i) * vsus(k) / km(k,j,i) 236 ! dvdx = - wall_e_x(j,i) * vsus(k,j,i) / km(k,j,i) 237 dwdx = - wall_e_x(j,i) * wsus(k) / km(k,j,i) 238 ! dwdx = - wall_e_x(j,i) * wsus(k,j,i) / km(k,j,i) 238 239 ELSE 239 240 dvdx = 0.25 * ( v(k,j,i+1) + v(k,j+1,i+1) - & … … 629 630 CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, & 630 631 usvs, 1.0, 0.0, 0.0, 0.0 ) 631 dudy = wall_e_y(j,i) * usvs(k) / km(k,j,i)632 dudy = - wall_e_y(j,i) * usvs(k) / km(k,j,i) 632 633 CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, & 633 634 wsvs, 0.0, 0.0, 1.0, 0.0 ) 634 dwdy = wall_e_y(j,i) * wsvs(k) / km(k,j,i)635 dwdy = - wall_e_y(j,i) * wsvs(k) / km(k,j,i) 635 636 ELSE 636 637 dudy = 0.25 * ( u(k,j+1,i) + u(k,j+1,i+1) - & … … 643 644 CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, & 644 645 vsus, 0.0, 1.0, 0.0, 0.0 ) 645 dvdx = wall_e_x(j,i) * vsus(k) / km(k,j,i)646 dvdx = - wall_e_x(j,i) * vsus(k) / km(k,j,i) 646 647 CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, & 647 648 wsus, 0.0, 0.0, 0.0, 1.0 ) 648 dwdx = wall_e_x(j,i) * wsus(k) / km(k,j,i)649 dwdx = - wall_e_x(j,i) * wsus(k) / km(k,j,i) 649 650 ELSE 650 651 dvdx = 0.25 * ( v(k,j,i+1) + v(k,j+1,i+1) - & … … 679 680 680 681 IF ( wall_e_y(j,i) /= 0.0 ) THEN 681 dudy = wall_e_y(j,i) * usvs(k) / km(k,j,i)682 dwdy = wall_e_y(j,i) * wsvs(k) / km(k,j,i)682 dudy = - wall_e_y(j,i) * usvs(k) / km(k,j,i) 683 dwdy = - wall_e_y(j,i) * wsvs(k) / km(k,j,i) 683 684 ELSE 684 685 dudy = 0.25 * ( u(k,j+1,i) + u(k,j+1,i+1) - & … … 689 690 690 691 IF ( wall_e_x(j,i) /= 0.0 ) THEN 691 dvdx = wall_e_x(j,i) * vsus(k) / km(k,j,i)692 dwdx = wall_e_x(j,i) * wsus(k) / km(k,j,i)692 dvdx = - wall_e_x(j,i) * vsus(k) / km(k,j,i) 693 dwdx = - wall_e_x(j,i) * wsus(k) / km(k,j,i) 693 694 ELSE 694 695 dvdx = 0.25 * ( v(k,j,i+1) + v(k,j+1,i+1) - &
Note: See TracChangeset
for help on using the changeset viewer.