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

Last change on this file since 1320 was 1320, checked in by raasch, 7 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: 5.4 KB
Line 
1 MODULE advec_u_pw_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_pw.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/11 06:09:21  raasch
37! Initial revision
38!
39!
40! Description:
41! ------------
42! Advection term for u velocity-component using Piacsek and Williams.
43! Vertical advection at the first grid point above the surface is done with
44! normal centred differences, because otherwise no information from the surface
45! would be communicated upwards due to w=0 at K=nzb.
46!------------------------------------------------------------------------------!
47
48    PRIVATE
49    PUBLIC advec_u_pw
50
51    INTERFACE advec_u_pw
52       MODULE PROCEDURE advec_u_pw
53       MODULE PROCEDURE advec_u_pw_ij
54    END INTERFACE advec_u_pw
55 
56 CONTAINS
57
58
59!------------------------------------------------------------------------------!
60! Call for all grid points
61!------------------------------------------------------------------------------!
62    SUBROUTINE advec_u_pw
63
64       USE arrays_3d,                                                          &
65           ONLY:  ddzw, tend, u, v, w
66
67       USE control_parameters,                                                 &
68           ONLY:  u_gtrans, v_gtrans
69
70       USE grid_variables,                                                     &
71           ONLY:  ddx, ddy
72
73       USE indices,                                                            &
74           ONLY:  nxlu, nxr, nyn, nys, nzb_u_inner, nzt
75
76       USE kinds
77
78
79       IMPLICIT NONE
80
81       INTEGER(iwp) ::  i !:
82       INTEGER(iwp) ::  j !:
83       INTEGER(iwp) ::  k !:
84       
85       REAL(wp)    ::  gu !:
86       REAL(wp)    ::  gv !:
87 
88       gu = 2.0 * u_gtrans
89       gv = 2.0 * v_gtrans
90       DO  i = nxlu, nxr
91          DO  j = nys, nyn
92             DO  k = nzb_u_inner(j,i)+1, nzt
93                tend(k,j,i) = tend(k,j,i) - 0.25 * (                           &
94                         ( u(k,j,i+1) * ( u(k,j,i+1) + u(k,j,i) - gu )         &
95                         - u(k,j,i-1) * ( u(k,j,i) + u(k,j,i-1) - gu ) ) * ddx &
96                       + ( u(k,j+1,i) * ( v(k,j+1,i) + v(k,j+1,i-1) - gv )     &
97                         - u(k,j-1,i) * ( v(k,j,i) + v(k,j,i-1) - gv ) ) * ddy &
98                       + ( u(k+1,j,i) * ( w(k,j,i) + w(k,j,i-1) )              &
99                         - u(k-1,j,i) * ( w(k-1,j,i) + w(k-1,j,i-1) ) )        &
100                                                                  * ddzw(k)    &
101                                                   )
102             ENDDO
103          ENDDO
104       ENDDO
105
106    END SUBROUTINE advec_u_pw
107
108
109!------------------------------------------------------------------------------!
110! Call for grid point i,j
111!------------------------------------------------------------------------------!
112    SUBROUTINE advec_u_pw_ij( i, j )
113
114       USE arrays_3d,                                                          &
115           ONLY:  ddzw, tend, u, v, w
116
117       USE control_parameters,                                                 &
118           ONLY:  u_gtrans, v_gtrans
119
120       USE grid_variables,                                                     &
121           ONLY:  ddx, ddy
122
123       USE indices,                                                            &
124           ONLY:  nzb_u_inner, nzt
125
126       USE kinds
127
128
129       IMPLICIT NONE
130
131       INTEGER(iwp) ::  i !:
132       INTEGER(iwp) ::  j !:
133       INTEGER(iwp) ::  k !:
134       
135       REAL(wp)    ::  gu !:
136       REAL(wp)    ::  gv !:
137
138       gu = 2.0 * u_gtrans
139       gv = 2.0 * v_gtrans
140       DO  k = nzb_u_inner(j,i)+1, nzt
141          tend(k,j,i) = tend(k,j,i) - 0.25 * (                                 &
142                         ( u(k,j,i+1) * ( u(k,j,i+1) + u(k,j,i) - gu )         &
143                         - u(k,j,i-1) * ( u(k,j,i) + u(k,j,i-1) - gu ) ) * ddx &
144                       + ( u(k,j+1,i) * ( v(k,j+1,i) + v(k,j+1,i-1) - gv )     &
145                         - u(k,j-1,i) * ( v(k,j,i) + v(k,j,i-1) - gv ) ) * ddy &
146                       + ( u(k+1,j,i) * ( w(k,j,i) + w(k,j,i-1) )              &
147                         - u(k-1,j,i) * ( w(k-1,j,i) + w(k-1,j,i-1) ) )        &
148                                                                  * ddzw(k)    &
149                                             )
150       ENDDO
151
152    END SUBROUTINE advec_u_pw_ij
153
154 END MODULE advec_u_pw_mod
Note: See TracBrowser for help on using the repository browser.