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_v.f90

    r53 r56  
    3232!------------------------------------------------------------------------------!
    3333
     34    USE wall_fluxes_mod
     35
    3436    PRIVATE
    3537    PUBLIC diffusion_v
     
    5961       REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
    6062       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
    61        REAL, DIMENSION(nzb:nzt+1)      ::  vsus
    6263       REAL, DIMENSION(:,:),   POINTER ::  vsws
    6364       REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
     65       REAL, DIMENSION(nzb:nzt+1,nys:nyn+vynp,nxl:nxr) ::  vsus
     66
     67!
     68!--    First calculate horizontal momentum flux v'u' at vertical walls,
     69!--    if neccessary
     70       IF ( topography /= 'flat' )  THEN
     71          CALL wall_fluxes( vsus, 0.0, 1.0, 0.0, 0.0, 0, vynp, nzb_v_inner, &
     72                            nzb_v_outer, wall_v )
     73       ENDIF
    6474
    6575       DO  i = nxl, nxr
     
    101111!--          Wall functions at the left and right walls, respectively
    102112             IF ( wall_v(j,i) /= 0.0 )  THEN
    103 
    104 !
    105 !--             Calculate the horizontal momentum flux v'u'
    106                 CALL wall_fluxes( i, j, nzb_v_inner(j,i)+1, nzb_v_outer(j,i), &
    107                                   vsus, 0.0, 1.0, 0.0, 0.0 )
    108113
    109114                DO  k = nzb_v_inner(j,i)+1, nzb_v_outer(j,i)
     
    137142                                + kmxm_y * ( u(k,j,i) - u(k,j-1,i) ) * ddy     &
    138143                                                  )                            &
    139                                      + wall_v(j,i) * vsus(k)                   &
     144                                     + wall_v(j,i) * vsus(k,j,i)               &
    140145                                   ) * ddx
    141146                ENDDO
Note: See TracChangeset for help on using the changeset viewer.