Home

Context Navigation

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

Last change on this file since 1320 was 1320, checked in by raasch, 10 years ago
ONLY-attribute added to USE-statements, kind-parameters added to all INTEGER and REAL declaration statements, kinds are defined in new module kinds, old module precision_kind is removed, revision history before 2012 removed, comment fields (!:) to be used for variable explanations added to all variable declaration statements
Property svn:keywords set to `Id`
File size: 6.9 KB

Line
1	MODULE advec_u_up_mod
2
3	!--------------------------------------------------------------------------------!
4	! This file is part of PALM.
5	!
6	! PALM is free software: you can redistribute it and/or modify it under the terms
7	! of the GNU General Public License as published by the Free Software Foundation,
8	! either version 3 of the License, or (at your option) any later version.
9	!
10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
13	!
14	! You should have received a copy of the GNU General Public License along with
15	! PALM. If not, see <http://www.gnu.org/licenses/>.
16	!
17	! Copyright 1997-2014 Leibniz Universitaet Hannover
18	!--------------------------------------------------------------------------------!
19	!
20	! Current revisions:
21	! -----------------
22	! ONLY-attribute added to USE-statements,
23	! kind-parameters added to all INTEGER and REAL declaration statements,
24	! kinds are defined in new module kinds,
25	! revision history before 2012 removed,
26	! comment fields (!:) to be used for variable explanations added to
27	! all variable declaration statements
28	!
29	! Former revisions:
30	! -----------------
31	! $Id: advec_u_up.f90 1320 2014-03-20 08:40:49Z raasch $
32	!
33	! 1036 2012-10-22 13:43:42Z raasch
34	! code put under GPL (PALM 3.9)
35	!
36	! Revision 1.1 1997/08/29 08:55:25 raasch
37	! Initial revision
38	!
39	!
40	! Description:
41	! ------------
42	! Advection term for the u velocity-component using upstream scheme.
43	! NOTE: vertical advection at k=1 still has wrong grid spacing for w>0!
44	! The same problem occurs for all topography boundaries!
45	!------------------------------------------------------------------------------!
46
47	PRIVATE
48	PUBLIC advec_u_up
49
50	INTERFACE advec_u_up
51	MODULE PROCEDURE advec_u_up
52	MODULE PROCEDURE advec_u_up_ij
53	END INTERFACE advec_u_up
54
55	CONTAINS
56
57
58	!------------------------------------------------------------------------------!
59	! Call for all grid points
60	!------------------------------------------------------------------------------!
61	SUBROUTINE advec_u_up
62
63	USE arrays_3d, &
64	ONLY: ddzu, tend, u, v, w
65
66	USE control_parameters, &
67	ONLY: u_gtrans, v_gtrans
68
69	USE grid_variables, &
70	ONLY: ddx, ddy
71
72	USE indices, &
73	ONLY: nxlu, nxr, nyn, nys, nzb_u_inner, nzt
74
75	USE kinds
76
77
78	IMPLICIT NONE
79
80	INTEGER(iwp) :: i !:
81	INTEGER(iwp) :: j !:
82	INTEGER(iwp) :: k !:
83
84	REAL(wp) :: ukomp !:
85	REAL(wp) :: vkomp !:
86	REAL(wp) :: wkomp !:
87
88
89	DO i = nxlu, nxr
90	DO j = nys, nyn
91	DO k = nzb_u_inner(j,i)+1, nzt
92	!
93	!-- x-direction
94	ukomp = u(k,j,i) - u_gtrans
95	IF ( ukomp > 0.0 ) THEN
96	tend(k,j,i) = tend(k,j,i) - ukomp * &
97	( u(k,j,i) - u(k,j,i-1) ) * ddx
98	ELSE
99	tend(k,j,i) = tend(k,j,i) - ukomp * &
100	( u(k,j,i+1) - u(k,j,i) ) * ddx
101	ENDIF
102	!
103	!-- y-direction
104	vkomp = 0.25 * ( v(k,j,i) + v(k,j+1,i) + &
105	v(k,j,i-1) + v(k,j+1,i-1) ) - v_gtrans
106	IF ( vkomp > 0.0 ) THEN
107	tend(k,j,i) = tend(k,j,i) - vkomp * &
108	( u(k,j,i) - u(k,j-1,i) ) * ddy
109	ELSE
110	tend(k,j,i) = tend(k,j,i) - vkomp * &
111	( u(k,j+1,i) - u(k,j,i) ) * ddy
112	ENDIF
113	!
114	!-- z-direction
115	wkomp = 0.25 * ( w(k,j,i) + w(k-1,j,i) + &
116	w(k,j,i-1) + w(k-1,j,i-1) )
117	IF ( wkomp > 0.0 ) THEN
118	tend(k,j,i) = tend(k,j,i) - wkomp * &
119	( u(k,j,i) - u(k-1,j,i) ) * ddzu(k)
120	ELSE
121	tend(k,j,i) = tend(k,j,i) - wkomp * &
122	( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1)
123	ENDIF
124
125	ENDDO
126	ENDDO
127	ENDDO
128
129	END SUBROUTINE advec_u_up
130
131
132	!------------------------------------------------------------------------------!
133	! Call for grid point i,j
134	!------------------------------------------------------------------------------!
135	SUBROUTINE advec_u_up_ij( i, j )
136
137	USE arrays_3d, &
138	ONLY: ddzu, tend, u, v, w
139
140	USE control_parameters, &
141	ONLY: u_gtrans, v_gtrans
142
143	USE grid_variables, &
144	ONLY: ddx, ddy
145
146	USE indices, &
147	ONLY: nzb_u_inner, nzt
148
149	USE kinds
150
151
152	IMPLICIT NONE
153
154	INTEGER(iwp) :: i !:
155	INTEGER(iwp) :: j !:
156	INTEGER(iwp) :: k !:
157
158	REAL(wp) :: ukomp !:
159	REAL(wp) :: vkomp !:
160	REAL(wp) :: wkomp !:
161
162
163	DO k = nzb_u_inner(j,i)+1, nzt
164	!
165	!-- x-direction
166	ukomp = u(k,j,i) - u_gtrans
167	IF ( ukomp > 0.0 ) THEN
168	tend(k,j,i) = tend(k,j,i) - ukomp * &
169	( u(k,j,i) - u(k,j,i-1) ) * ddx
170	ELSE
171	tend(k,j,i) = tend(k,j,i) - ukomp * &
172	( u(k,j,i+1) - u(k,j,i) ) * ddx
173	ENDIF
174	!
175	!-- y-direction
176	vkomp = 0.25 * ( v(k,j,i) + v(k,j+1,i) + v(k,j,i-1) + v(k,j+1,i-1) &
177	) - v_gtrans
178	IF ( vkomp > 0.0 ) THEN
179	tend(k,j,i) = tend(k,j,i) - vkomp * &
180	( u(k,j,i) - u(k,j-1,i) ) * ddy
181	ELSE
182	tend(k,j,i) = tend(k,j,i) - vkomp * &
183	( u(k,j+1,i) - u(k,j,i) ) * ddy
184	ENDIF
185	!
186	!-- z-direction
187	wkomp = 0.25 * ( w(k,j,i) + w(k-1,j,i) + w(k,j,i-1) + w(k-1,j,i-1) )
188	IF ( wkomp > 0.0 ) THEN
189	tend(k,j,i) = tend(k,j,i) - wkomp * &
190	( u(k,j,i) - u(k-1,j,i) ) * ddzu(k)
191	ELSE
192	tend(k,j,i) = tend(k,j,i) - wkomp * &
193	( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1)
194	ENDIF
195
196	ENDDO
197
198	END SUBROUTINE advec_u_up_ij
199
200	END MODULE advec_u_up_mod

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |