[1682] | 1 | !> @file user_init_3d_model.f90 |
---|
[4498] | 2 | !--------------------------------------------------------------------------------------------------! |
---|
[2696] | 3 | ! This file is part of the PALM model system. |
---|
[1036] | 4 | ! |
---|
[4498] | 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. |
---|
[1036] | 8 | ! |
---|
[4498] | 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. |
---|
[1036] | 12 | ! |
---|
[4498] | 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/>. |
---|
[1036] | 15 | ! |
---|
[4360] | 16 | ! Copyright 1997-2020 Leibniz Universitaet Hannover |
---|
[4498] | 17 | !--------------------------------------------------------------------------------------------------! |
---|
[1036] | 18 | ! |
---|
[4498] | 19 | ! |
---|
[484] | 20 | ! Current revisions: |
---|
[211] | 21 | ! ----------------- |
---|
[1683] | 22 | ! |
---|
[2233] | 23 | ! |
---|
[211] | 24 | ! Former revisions: |
---|
| 25 | ! ----------------- |
---|
| 26 | ! $Id: user_init_3d_model.f90 4498 2020-04-15 14:26:31Z maronga $ |
---|
[4498] | 27 | ! file re-formatted to follow the PALM coding standard |
---|
| 28 | ! |
---|
| 29 | ! |
---|
| 30 | ! 4360 2020-01-07 11:25:50Z suehring |
---|
[4182] | 31 | ! Corrected "Former revisions" section |
---|
[4498] | 32 | ! |
---|
[4182] | 33 | ! 3768 2019-02-27 14:35:58Z raasch |
---|
[3768] | 34 | ! variables commented out to avoid compiler warnings about unused variables |
---|
[4498] | 35 | ! |
---|
[3768] | 36 | ! 3655 2019-01-07 16:51:22Z knoop |
---|
[2716] | 37 | ! Corrected "Former revisions" section |
---|
| 38 | ! |
---|
[4182] | 39 | ! 211 2008-11-11 04:46:24Z raasch |
---|
| 40 | ! Former file user_interface.f90 split into one file per subroutine |
---|
| 41 | ! |
---|
[211] | 42 | ! Description: |
---|
| 43 | ! ------------ |
---|
[1682] | 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! |
---|
[4498] | 48 | !--------------------------------------------------------------------------------------------------! |
---|
[1682] | 49 | SUBROUTINE user_init_3d_model |
---|
[211] | 50 | |
---|
[4498] | 51 | |
---|
[211] | 52 | USE arrays_3d |
---|
[4498] | 53 | |
---|
[211] | 54 | USE control_parameters |
---|
[4498] | 55 | |
---|
[211] | 56 | USE indices |
---|
[4498] | 57 | |
---|
[1320] | 58 | USE kinds |
---|
[2232] | 59 | |
---|
| 60 | USE surface_mod |
---|
[4498] | 61 | |
---|
[211] | 62 | USE user |
---|
| 63 | |
---|
| 64 | IMPLICIT NONE |
---|
| 65 | |
---|
[4498] | 66 | ! INTEGER(iwp) :: l !< running index surface orientation |
---|
| 67 | ! INTEGER(iwp) :: m !< running index surface elements |
---|
[211] | 68 | |
---|
[2618] | 69 | ! |
---|
| 70 | !-- Initialization of surface-related quantities. |
---|
[4498] | 71 | !-- The following example shows required initialization of surface quantitites at default-type |
---|
| 72 | !-- upward-facing surfaces. |
---|
[2618] | 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 |
---|
[4498] | 93 | ! ENDDO |
---|
[2618] | 94 | ! |
---|
| 95 | !-- Same for natural and urban type surfaces |
---|
| 96 | ! DO m = 1, surf_lsm_h%ns |
---|
| 97 | ! ... |
---|
[4498] | 98 | ! ENDDO |
---|
[2618] | 99 | ! DO m = 1, surf_usm_h%ns |
---|
| 100 | ! ... |
---|
| 101 | ! ENDDO |
---|
| 102 | ! |
---|
[4498] | 103 | !-- Also care for vertically aligned surfaces (default-, natural-, and |
---|
[2618] | 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 | |
---|
[211] | 120 | END SUBROUTINE user_init_3d_model |
---|
| 121 | |
---|