# Changeset 134

Ignore:
Timestamp:
Nov 21, 2007 7:28:38 AM (14 years ago)
Message:

Redefine initial nzb_local as the actual total size of topography (later the
extent of topography in nzb_local is reduced by 1dx at the E topography walls
and by 1dy at the N topography walls to form the basis for nzb_s_inner);
for consistency redefine 'single_building' case.

Location:
palm/trunk
Files:
4 edited

Unmodified
Removed
• ## palm/trunk/DOC/app/chapter_3.4.html

 r108 PALM chapter 3.4PALM chapter 3.4                        S

These data must exactly match the horizontal grid.
S

Example for a 50m building surrounded by a 12.5m podium on flat ground:
0    0    0    0    0    0    0    0    0    0    0    0    0
0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5    0    0
0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5    0    0
0 12.5 12.5   50   50   50   50   50   50 12.5 12.5    0    0
0 12.5 12.5   50   50   50   50   50   50 12.5 12.5    0    0
0 12.5 12.5   50   50   50   50   50   50 12.5 12.5    0    0
0 12.5 12.5   50   50   50   50   50   50 12.5 12.5    0    0
0 12.5 12.5   50   50   50   50   50   50 12.5 12.5    0    0
0 12.5 12.5   50   50   50   50   50   50 12.5 12.5    0    0
0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5    0    0
0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5    0    0
0    0    0    0    0    0    0    0    0    0    0    0    0

These data must exactly match the horizontal grid. Due to the staggered grid the topography will be displaced by -0.5 dx in x-direction and -0.5 dy in y-direction.

• ## palm/trunk/DOC/app/chapter_4.1.html

 r132 = ( ( nx + 1 ) * dxbuilding_length_x ) / 2 centers the building in x-direction. centers the building in x-direction. Due to the staggered grid the building will be displaced by -0.5 dx in x-direction and -0.5 dy in y-direction. = ( ( ny + 1 ) * dybuilding_length_y ) / 2 centers the building in y-direction. centers the building in y-direction. Due to the staggered grid the building will be displaced by -0.5 dx in x-direction and -0.5 dy in y-direction. The building size and location can be specified with the parameters building_height, building_length_x, building_length_y, building_wall_left and building_wall_south. The building size and location can be specified with the parameters building_height, building_length_x, building_length_y, building_wall_left and building_wall_south.
Due to the staggered grid the building will be displaced by -0.5 dx in x-direction and -0.5 dy in y-direction. information in m. These data must exactly match the horizontal grid. exactly match the horizontal grid.
Due to the staggered grid the topography will be displaced by -0.5 dx in x-direction and -0.5 dy in y-direction. Alternatively, the user may add code to the user interface subroutine user_init_grid to allow further topography modes.
to allow further topography modes. Again,
due to the staggered grid the topography will be displaced by -0.5 dx in x-direction and -0.5 dy in y-direction.
• ## palm/trunk/SOURCE/CURRENT_MODIFICATIONS

 r133 ------- Redefine initial nzb_local as the actual total size of topography (later the extent of topography in nzb_local is reduced by 1dx at the E topography walls and by 1dy at the N topography walls to form the basis for nzb_s_inner); for consistency redefine 'single_building' case. Vertical profiles now based on nzb_s_inner; they are divided by ngp_2dh_s_inner (scalars, procucts of scalars) and ngp_2dh (staggered velocity humidity/scalar/salinity in a future release. buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, modules, prognostic_equations buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, init_grid, modules, prognostic_equations
• ## palm/trunk/SOURCE/init_grid.f90

 r116 ! Actual revisions: ! ----------------- ! Redefine initial nzb_local as the actual total size of topography (later the ! extent of topography in nzb_local is reduced by 1dx at the E topography walls ! and by 1dy at the N topography walls to form the basis for nzb_s_inner); ! for consistency redefine 'single_building' case. ! Calculation of wall flag arrays ! ! !--       Set the individual index arrays for all velocity components and !--       scalars, taking into account the staggered grid. The horizontal !--       wind component normal to a wall defines the position of the wall, and !--       in the respective direction the building is as long as specified in !--       building_length_?, but in the other horizontal direction (for w and s !--       in both horizontal directions) the building appears shortened by one !--       grid length due to the staggered grid. !--       Set the actual total size of the building. Due to the staggered grid, !--       the building will be displaced by -0.5dx in x-direction and by -0.5dy !--       in y-direction compared to the scalar grid. nzb_local = 0 nzb_local(bys:byn-1,bxl:bxr-1) = bh nzb_local(bys:byn,bxl:bxr) = bh CASE ( 'read_from_file' ) ENDIF ENDIF ! !--    The array nzb_local as defined above describes the actual total size of !--    topography which is defined by u=0 on the topography walls in x-direction !--    and by v=0 on the topography walls in y-direction. However, PALM uses !--    individual arrays nzb_u|v|w|s_inner|outer that are based on nzb_s_inner. !--    Therefore, the extent of topography in nzb_local is now reduced by 1dx !--    at the E topography walls and by 1dy at the N topography walls to form !--    the basis for nzb_s_inner. DO  j = -gls, ny + gls DO  i = -gls, nx nzb_local(j,i) = MIN( nzb_local(j,i), nzb_local(j,i+1) ) ENDDO ENDDO !--    apply cyclic boundary conditions in x-direction nzb_local(:,nx+1:nx+gls) = nzb_local(:,0:gls-1) DO  i = -gls, nx + gls DO  j = -gls, ny nzb_local(j,i) = MIN( nzb_local(j,i), nzb_local(j+1,i) ) ENDDO ENDDO !--    apply cyclic boundary conditions in y-direction nzb_local(ny+1:ny+gls,:) = nzb_local(0:gls-1,:) !
Note: See TracChangeset for help on using the changeset viewer.