source: palm/trunk/SOURCE/advec_u_up.f90 @ 757

Last change on this file since 757 was 484, checked in by raasch, 15 years ago

typo in file headers removed

  • Property svn:keywords set to Id
File size: 4.7 KB
Line 
1 MODULE advec_u_up_mod
2
3!------------------------------------------------------------------------------!
4! Current revisions:
5! -----------------
6!
7!
8! Former revisions:
9! -----------------
10! $Id: advec_u_up.f90 484 2010-02-05 07:36:54Z helmke $
11!
12! 106 2007-08-16 14:30:26Z raasch
13! i loop is starting from nxlu (needed for non-cyclic boundary conditions)
14!
15! 75 2007-03-22 09:54:05Z raasch
16! uxrp eliminated
17!
18! RCS Log replace by Id keyword, revision history cleaned up
19!
20! Revision 1.12  2006/02/23 09:45:04  raasch
21! nzb_2d replaced by nzb_u_inner
22!
23! Revision 1.1  1997/08/29 08:55:25  raasch
24! Initial revision
25!
26!
27! Description:
28! ------------
29! Advection term for the u velocity-component using upstream scheme.
30! NOTE: vertical advection at k=1 still has wrong grid spacing for w>0!
31!       The same problem occurs for all topography boundaries!
32!------------------------------------------------------------------------------!
33
34    PRIVATE
35    PUBLIC advec_u_up
36
37    INTERFACE advec_u_up
38       MODULE PROCEDURE advec_u_up
39       MODULE PROCEDURE advec_u_up_ij
40    END INTERFACE advec_u_up
41
42 CONTAINS
43
44
45!------------------------------------------------------------------------------!
46! Call for all grid points
47!------------------------------------------------------------------------------!
48    SUBROUTINE advec_u_up
49
50       USE arrays_3d
51       USE control_parameters
52       USE grid_variables
53       USE indices
54
55       IMPLICIT NONE
56
57       INTEGER ::  i, j, k
58
59       REAL ::  ukomp, vkomp, wkomp
60
61
62       DO  i = nxlu, nxr
63          DO  j = nys, nyn
64             DO  k = nzb_u_inner(j,i)+1, nzt
65!
66!--             x-direction
67                ukomp = u(k,j,i) - u_gtrans
68                IF ( ukomp > 0.0 )  THEN
69                   tend(k,j,i) = tend(k,j,i) - ukomp * &
70                                         ( u(k,j,i) - u(k,j,i-1) ) * ddx
71                ELSE
72                   tend(k,j,i) = tend(k,j,i) - ukomp * &
73                                          ( u(k,j,i+1) - u(k,j,i) ) * ddx
74                ENDIF
75!
76!--             y-direction
77                vkomp = 0.25 * ( v(k,j,i)   + v(k,j+1,i) + &
78                                 v(k,j,i-1) + v(k,j+1,i-1) ) - v_gtrans
79                IF ( vkomp > 0.0 )  THEN
80                   tend(k,j,i) = tend(k,j,i) - vkomp * &
81                                         ( u(k,j,i) - u(k,j-1,i) ) * ddy
82                ELSE
83                   tend(k,j,i) = tend(k,j,i) - vkomp * &
84                                         ( u(k,j+1,i) - u(k,j,i) ) * ddy
85                ENDIF
86!
87!--             z-direction
88                wkomp = 0.25 * ( w(k,j,i)   + w(k-1,j,i) + &
89                                 w(k,j,i-1) + w(k-1,j,i-1) )
90                IF ( wkomp > 0.0 )  THEN
91                   tend(k,j,i) = tend(k,j,i) - wkomp * &
92                                         ( u(k,j,i) - u(k-1,j,i) ) * ddzu(k)
93                ELSE
94                   tend(k,j,i) = tend(k,j,i) - wkomp * &
95                                         ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1)
96                ENDIF
97
98             ENDDO
99          ENDDO
100       ENDDO
101
102    END SUBROUTINE advec_u_up
103
104
105!------------------------------------------------------------------------------!
106! Call for grid point i,j
107!------------------------------------------------------------------------------!
108    SUBROUTINE advec_u_up_ij( i, j )
109
110       USE arrays_3d
111       USE control_parameters
112       USE grid_variables
113       USE indices
114
115       IMPLICIT NONE
116
117       INTEGER ::  i, j, k
118
119       REAL ::  ukomp, vkomp, wkomp
120
121
122       DO  k = nzb_u_inner(j,i)+1, nzt
123!
124!--       x-direction
125          ukomp = u(k,j,i) - u_gtrans
126          IF ( ukomp > 0.0 )  THEN
127             tend(k,j,i) = tend(k,j,i) - ukomp * &
128                                         ( u(k,j,i) - u(k,j,i-1) ) * ddx
129          ELSE
130             tend(k,j,i) = tend(k,j,i) - ukomp * &
131                                         ( u(k,j,i+1) - u(k,j,i) ) * ddx
132          ENDIF
133!
134!--       y-direction
135          vkomp = 0.25 * ( v(k,j,i) + v(k,j+1,i) + v(k,j,i-1) + v(k,j+1,i-1) &
136                         ) - v_gtrans
137          IF ( vkomp > 0.0 )  THEN
138             tend(k,j,i) = tend(k,j,i) - vkomp * &
139                                         ( u(k,j,i) - u(k,j-1,i) ) * ddy
140          ELSE
141             tend(k,j,i) = tend(k,j,i) - vkomp * &
142                                         ( u(k,j+1,i) - u(k,j,i) ) * ddy
143          ENDIF
144!
145!--       z-direction
146          wkomp = 0.25 * ( w(k,j,i) + w(k-1,j,i) + w(k,j,i-1) + w(k-1,j,i-1) )
147          IF ( wkomp > 0.0 )  THEN
148             tend(k,j,i) = tend(k,j,i) - wkomp * &
149                                         ( u(k,j,i) - u(k-1,j,i) ) * ddzu(k)
150          ELSE
151             tend(k,j,i) = tend(k,j,i) - wkomp * &
152                                         ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1)
153          ENDIF
154
155       ENDDO
156
157    END SUBROUTINE advec_u_up_ij
158
159 END MODULE advec_u_up_mod
Note: See TracBrowser for help on using the repository browser.