1 | !> @file user_init_3d_model.f90 |
---|
2 | !--------------------------------------------------------------------------------------------------! |
---|
3 | ! This file is part of the PALM model system. |
---|
4 | ! |
---|
5 | ! PALM is free software: you can redistribute it and/or modify it under the terms of the GNU General |
---|
6 | ! Public License as published by the Free Software Foundation, either version 3 of the License, or |
---|
7 | ! (at your option) any later version. |
---|
8 | ! |
---|
9 | ! PALM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the |
---|
10 | ! implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General |
---|
11 | ! Public License for more details. |
---|
12 | ! |
---|
13 | ! You should have received a copy of the GNU General Public License along with PALM. If not, see |
---|
14 | ! <http://www.gnu.org/licenses/>. |
---|
15 | ! |
---|
16 | ! Copyright 1997-2021 Leibniz Universitaet Hannover |
---|
17 | !--------------------------------------------------------------------------------------------------! |
---|
18 | ! |
---|
19 | ! |
---|
20 | ! Current revisions: |
---|
21 | ! ----------------- |
---|
22 | ! |
---|
23 | ! |
---|
24 | ! Former revisions: |
---|
25 | ! ----------------- |
---|
26 | ! $Id: user_init_3d_model.f90 4828 2021-01-05 11:21:41Z moh.hefny $ |
---|
27 | ! file re-formatted to follow the PALM coding standard |
---|
28 | ! |
---|
29 | ! |
---|
30 | ! 4360 2020-01-07 11:25:50Z suehring |
---|
31 | ! Corrected "Former revisions" section |
---|
32 | ! |
---|
33 | ! 3768 2019-02-27 14:35:58Z raasch |
---|
34 | ! variables commented out to avoid compiler warnings about unused variables |
---|
35 | ! |
---|
36 | ! 3655 2019-01-07 16:51:22Z knoop |
---|
37 | ! Corrected "Former revisions" section |
---|
38 | ! |
---|
39 | ! 211 2008-11-11 04:46:24Z raasch |
---|
40 | ! Former file user_interface.f90 split into one file per subroutine |
---|
41 | ! |
---|
42 | ! Description: |
---|
43 | ! ------------ |
---|
44 | !> Allows the complete initialization of the 3d model. |
---|
45 | !> |
---|
46 | !> @attention The user is responsible to set at least all those quantities which |
---|
47 | !> are normally set within init_3d_model! |
---|
48 | !--------------------------------------------------------------------------------------------------! |
---|
49 | SUBROUTINE user_init_3d_model |
---|
50 | |
---|
51 | |
---|
52 | USE arrays_3d |
---|
53 | |
---|
54 | USE control_parameters |
---|
55 | |
---|
56 | USE indices |
---|
57 | |
---|
58 | USE kinds |
---|
59 | |
---|
60 | USE surface_mod |
---|
61 | |
---|
62 | USE user |
---|
63 | |
---|
64 | IMPLICIT NONE |
---|
65 | |
---|
66 | ! INTEGER(iwp) :: l !< running index surface orientation |
---|
67 | ! INTEGER(iwp) :: m !< running index surface elements |
---|
68 | |
---|
69 | ! |
---|
70 | !-- Initialization of surface-related quantities. |
---|
71 | !-- The following example shows required initialization of surface quantitites at default-type |
---|
72 | !-- upward-facing surfaces. |
---|
73 | ! DO m = 1, surf_def_h(0)%ns |
---|
74 | ! surf_def_h(0)%ol(m) = ... ! Obukhov length |
---|
75 | ! surf_def_h(0)%us(m ) = ... ! friction velocity |
---|
76 | ! surf_def_h(0)%usws(m) = ... ! vertical momentum flux, u-component |
---|
77 | ! surf_def_h(0)%vsws(m) = ... ! vertical momentum flux, v-component |
---|
78 | ! surf_def_h(0)%z0(m) = ... ! roughness length for momentum |
---|
79 | ! IF ( .NOT. neutral ) THEN |
---|
80 | ! surf_def_h(0)%ts(m) = ... ! scaling parameter |
---|
81 | ! surf_def_h(0)%shf(m) = ... ! surface sensible heat flux |
---|
82 | ! surf_def_h(0)%z0h(m) = ... ! roughness length for heat |
---|
83 | ! ENDIF |
---|
84 | ! IF ( humditiy ) THEN |
---|
85 | ! surf_def_h(0)%qs(m) = ... ! scaling parameter |
---|
86 | ! surf_def_h(0)%qsws(m) = ... ! surface latent heat flux |
---|
87 | ! surf_def_h(0)%z0q(m) = ... ! roughness length for moisture |
---|
88 | ! ENDIF |
---|
89 | ! IF ( passive_scalar ) THEN |
---|
90 | ! surf_def_h(0)%ss(m) = ... ! scaling parameter |
---|
91 | ! surf_def_h(0)%ssws(m) = ... ! surface latent heat flux |
---|
92 | ! ENDIF |
---|
93 | ! ENDDO |
---|
94 | ! |
---|
95 | !-- Same for natural and urban type surfaces |
---|
96 | ! DO m = 1, surf_lsm_h%ns |
---|
97 | ! ... |
---|
98 | ! ENDDO |
---|
99 | ! DO m = 1, surf_usm_h%ns |
---|
100 | ! ... |
---|
101 | ! ENDDO |
---|
102 | ! |
---|
103 | !-- Also care for vertically aligned surfaces (default-, natural-, and |
---|
104 | !-- urban-type). |
---|
105 | ! DO l = 0, 3 |
---|
106 | ! DO m = 1, surf_def_v(l)%ns |
---|
107 | ! ... |
---|
108 | ! ENDDO |
---|
109 | ! DO m = 1, surf_lsm_v(l)%ns |
---|
110 | ! ... |
---|
111 | ! ENDDO |
---|
112 | ! DO m = 1, surf_usm_v(l)%ns |
---|
113 | ! ... |
---|
114 | ! ENDDO |
---|
115 | ! ENDDO |
---|
116 | ! |
---|
117 | ! |
---|
118 | !-- In the following, initialize 3D quantities, e.g. u, v, w, pt, etc.. |
---|
119 | |
---|
120 | END SUBROUTINE user_init_3d_model |
---|
121 | |
---|