Changeset 134 for palm/trunk/SOURCE

Timestamp:

Nov 21, 2007 7:28:38 AM (16 years ago)

Author:

letzel

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/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (2 diffs)
• init_grid.f90 (modified) (3 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -21 +21 @@
 -------
 
+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
@@ -40 +45 @@
 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

@@ -4 +4 @@
 ! 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
 !
@@ -365 +369 @@
 
 !
-!--       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' )
@@ -468 +468 @@
           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,:)
 
 !