Ignore:
Timestamp:
Aug 6, 2008 4:25:09 PM (13 years ago)
Author:
letzel
Message:
  • new: descriptions of plant canopy model and prandtl layer (trunk/DOC/misc)
  • changed: more consistent flux definitions; modification of the integrated version of the profile function for momentum for unstable stratification (wall_fluxes, production_e)
  • bugfix: change definition of us_wall from 1D to 2D (prandtl_fluxes, wall_fluxes)
File:
1 edited

Legend:

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

    r139 r187  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Change: add 'minus' sign to fluxes obtained from subroutine wall_fluxes_e for
     7!         consistency with subroutine wall_fluxes
    78!
    89! Former revisions:
     
    164165                      CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, &
    165166                                          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)
    168169                      CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, &
    169170                                          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)
    172173                   ELSE
    173174                      dudy = 0.25 * ( u(k,j+1,i) + u(k,j+1,i+1) - &
     
    180181                      CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, &
    181182                                          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)
    184185                      CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, &
    185186                                          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)
    188189                   ELSE
    189190                      dvdx = 0.25 * ( v(k,j,i+1) + v(k,j+1,i+1) - &
     
    220221
    221222                      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)
    226227                      ELSE
    227228                         dudy = 0.25 * ( u(k,j+1,i) + u(k,j+1,i+1) - &
     
    232233
    233234                      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)
    238239                      ELSE
    239240                         dvdx = 0.25 * ( v(k,j,i+1) + v(k,j+1,i+1) - &
     
    629630                CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, &
    630631                                    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)
    632633                CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, &
    633634                                    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)
    635636             ELSE
    636637                dudy = 0.25 * ( u(k,j+1,i) + u(k,j+1,i+1) - &
     
    643644                CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, &
    644645                                    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)
    646647                CALL wall_fluxes_e( i, j, k, nzb_diff_s_outer(j,i)-2, &
    647648                                    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)
    649650             ELSE
    650651                dvdx = 0.25 * ( v(k,j,i+1) + v(k,j+1,i+1) - &
     
    679680
    680681                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)
    683684                ELSE
    684685                   dudy = 0.25 * ( u(k,j+1,i) + u(k,j+1,i+1) - &
     
    689690
    690691                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)
    693694                ELSE
    694695                   dvdx = 0.25 * ( v(k,j,i+1) + v(k,j+1,i+1) - &
Note: See TracChangeset for help on using the changeset viewer.