Ignore:
Timestamp:
Mar 8, 2007 1:57:07 PM (15 years ago)
Author:
raasch
Message:

further checkin of preliminary changes

File:
1 edited

Legend:

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

    r53 r56  
    3434!------------------------------------------------------------------------------!
    3535
     36    USE wall_fluxes_mod
     37
    3638    PRIVATE
    3739    PUBLIC diffusion_u
     
    6163       REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
    6264       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
    63        REAL, DIMENSION(nzb:nzt+1)      ::  usvs
    6465       REAL, DIMENSION(:,:),   POINTER ::  usws
    6566       REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
     67       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr+uxrp) ::  usvs
     68
     69!
     70!--    First calculate horizontal momentum flux u'v' at vertical walls,
     71!--    if neccessary
     72       IF ( topography /= 'flat' )  THEN
     73          CALL wall_fluxes( usvs, 1.0, 0.0, 0.0, 0.0, uxrp, 0, nzb_u_inner, &
     74                            nzb_u_outer, wall_u )
     75       ENDIF
    6676
    6777       DO  i = nxl, nxr+uxrp
     
    103113!--          Wall functions at the north and south walls, respectively
    104114             IF ( wall_u(j,i) /= 0.0 )  THEN
    105 
    106 !
    107 !--             Calculate the horizontal momentum flux u'v'
    108                 CALL wall_fluxes( i, j, nzb_u_inner(j,i)+1, nzb_u_outer(j,i),  &
    109                                   usvs, 1.0, 0.0, 0.0, 0.0 )
    110115
    111116                DO  k = nzb_u_inner(j,i)+1, nzb_u_outer(j,i)
     
    139144                                + kmym_x * ( v(k,j,i) - v(k,j,i-1) ) * ddx     &
    140145                                                  )                            &
    141                                      + wall_u(j,i) * usvs(k)                   &
     146                                     + wall_u(j,i) * usvs(k,j,i)               &
    142147                                   ) * ddy
    143148                ENDDO
Note: See TracChangeset for help on using the changeset viewer.