Changeset 56 for palm/trunk/SOURCE/diffusion_v.f90
- Timestamp:
- Mar 8, 2007 1:57:07 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diffusion_v.f90
r53 r56 32 32 !------------------------------------------------------------------------------! 33 33 34 USE wall_fluxes_mod 35 34 36 PRIVATE 35 37 PUBLIC diffusion_v … … 59 61 REAL :: z0(nys-1:nyn+1,nxl-1:nxr+1) 60 62 REAL :: tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) 61 REAL, DIMENSION(nzb:nzt+1) :: vsus62 63 REAL, DIMENSION(:,:), POINTER :: vsws 63 64 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 64 74 65 75 DO i = nxl, nxr … … 101 111 !-- Wall functions at the left and right walls, respectively 102 112 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 )108 113 109 114 DO k = nzb_v_inner(j,i)+1, nzb_v_outer(j,i) … … 137 142 + kmxm_y * ( u(k,j,i) - u(k,j-1,i) ) * ddy & 138 143 ) & 139 + wall_v(j,i) * vsus(k )&144 + wall_v(j,i) * vsus(k,j,i) & 140 145 ) * ddx 141 146 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.