Home

Context Navigation

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

Last change on this file since 1682 was 1682, checked in by knoop, 9 years ago
Code annotations made doxygen readable
Property svn:keywords set to `Id`
File size: 7.2 KB

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