Changeset 1931 for palm/trunk/SOURCE/init_grid.f90
 Timestamp:
 Jun 10, 2016 12:06:59 PM (8 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/init_grid.f90
r1911 r1931 19 19 ! Current revisions: 20 20 !  21 ! 21 ! Rename multigrid into multigrid_noopt and multigrid_fast into multigrid 22 22 ! 23 23 ! Former revisions: … … 225 225 INTEGER(iwp) :: cyn !< index for north canyon wall 226 226 INTEGER(iwp) :: cys !< index for south canyon wall 227 INTEGER(iwp) :: gls !< number of lateral ghost points at total model 228 !< domain boundaries required for multigrid solver 227 INTEGER(iwp) :: gls !< number of lateral ghost points at total model domain boundaries required for multigrid solver 229 228 INTEGER(iwp) :: i !< index variable along x 230 229 INTEGER(iwp) :: ii !< loop variable for reading topography file … … 242 241 243 242 INTEGER(iwp), DIMENSION(:), ALLOCATABLE :: & 244 vertical_influence !< number of vertical grid points above 245 !< obstacle where adjustment of near 246 !< wall mixing length is required 243 vertical_influence !< number of vertical grid points above obstacle where adjustment of nearwall mixing length is required 247 244 248 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: corner_nl !< index of 249 !< northleft corner location to limit 250 !< nearwall mixing length 245 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: corner_nl !< index of northleft corner location to limit nearwall mixing length 251 246 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: corner_nr !< northright 252 247 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: corner_sl !< southleft 253 248 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: corner_sr !< southright 254 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: wall_l !< distance to 255 !< adjacent leftfacing 249 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: wall_l !< distance to adjacent leftfacing 256 250 !< wall 257 251 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: wall_n !< northfacing … … 260 254 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_local !< index for topography 261 255 !< top at cellcenter 262 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_tmp !< dummy to calculate 263 !< topography indices 264 !< on u and vgrid 256 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_tmp !< dummy to calculate topography indices on u and vgrid 265 257 266 258 LOGICAL :: flag_set = .FALSE. !< steering variable for advection flags … … 270 262 REAL(wp) :: dz_stretched !< stretched vertical grid spacing 271 263 272 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: topo_height !< input variable for 273 !< topography height 264 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: topo_height !< input variable for topography height 274 265 275 266 … … 1138 1129 1139 1130 ! 1140 ! Calculate wall flag arrays for the multigrid method 1141 IF ( psolver(1:9) == 'multigrid' ) THEN 1131 ! Calculate wall flag arrays for the multigrid method. 1132 ! Please note, wall flags are only applied in the not cacheoptimized 1133 ! version. 1134 IF ( psolver == 'multigrid_noopt' ) THEN 1142 1135 ! 1143 1136 ! Gridpoint increment of the current level … … 1193 1186 ! In case of masking method, flags are not set and multigrid method 1194 1187 ! works like FFTsolver 1195 IF ( psolver == 'multigrid' .AND..NOT. masking_method ) THEN1188 IF ( .NOT. masking_method ) THEN 1196 1189 1197 1190 DO i = nxl_l1, nxr_l+1
Note: See TracChangeset
for help on using the changeset viewer.