Changeset 119 for palm/trunk/SOURCE

Timestamp:

Oct 17, 2007 10:27:13 AM (17 years ago)

Author:

raasch

Message:

small bugfixes in SGS part of adved_particles

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (1 diff)
• advec_particles.f90 (modified) (8 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -22 +22 @@
 
 Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
-numbers along y (advec_particles)
+numbers along y, small bugfixes in the SGS part (advec_particles)
 
 Bugfix: model_string needed a default value (combine_plot_fields)

palm/trunk/SOURCE/advec_particles.f90

@@ -8 +8 @@
 ! Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
 ! numbers along y
+! Small bugfixes in the SGS part
 ! TEST: PRINT statements on unit 9 (commented out)
 !
@@ -1056 +1057 @@
                    gp_outside_of_building(1) = 1
                    location(num_gp,1) = i * dx
-                   location(num_gp,2) = j * dx
+                   location(num_gp,2) = j * dy
                    location(num_gp,3) = k * dz - 0.5 * dz
                    ei(num_gp)     = e(k,j,i)
@@ -1070 +1071 @@
                    gp_outside_of_building(2) = 1
                    location(num_gp,1) = i * dx
-                   location(num_gp,2) = (j+1) * dx
+                   location(num_gp,2) = (j+1) * dy
                    location(num_gp,3) = k * dz - 0.5 * dz
                    ei(num_gp)     = e(k,j+1,i)
@@ -1513 +1514 @@
 !--                Neumann boundary condition has to be applied
 !--                (only one case as only building beneath is possible)
-                   IF ( gp_outside_of_building(1) == 1  .AND. &
-                        gp_outside_of_building(3) == 0 )  THEN
+                   IF ( gp_outside_of_building(1) == 0  .AND. &
+                        gp_outside_of_building(3) == 1 )  THEN
                       num_gp = num_gp + 1
                       location(num_gp,1) = i * dx
@@ -1530 +1531 @@
 !--                Neumann boundary condition has to be applied
 !--                (only one case as only building beneath is possible)
-                   IF ( gp_outside_of_building(5) == 1  .AND. &
-                        gp_outside_of_building(7) == 0 )  THEN
+                   IF ( gp_outside_of_building(5) == 0  .AND. &
+                        gp_outside_of_building(7) == 1 )  THEN
                       num_gp = num_gp + 1
                       location(num_gp,1) = (i+1) * dx
@@ -1547 +1548 @@
 !--                Neumann boundary condition has to be applied
 !--                (only one case as only building beneath is possible)
-                   IF ( gp_outside_of_building(2) == 1  .AND. &
-                        gp_outside_of_building(4) == 0 )  THEN
+                   IF ( gp_outside_of_building(2) == 0  .AND. &
+                        gp_outside_of_building(4) == 1 )  THEN
                       num_gp = num_gp + 1
                       location(num_gp,1) = i * dx
@@ -1564 +1565 @@
 !--                Neumann boundary condition has to be applied
 !--                (only one case as only building beneath is possible)
-                   IF ( gp_outside_of_building(6) == 1  .AND. &
-                        gp_outside_of_building(8) == 0 )  THEN
+                   IF ( gp_outside_of_building(6) == 0  .AND. &
+                        gp_outside_of_building(8) == 1 )  THEN
                       num_gp = num_gp + 1
                       location(num_gp,1) = (i+1) * dx
@@ -1584 +1585 @@
 !--                   building, it follows that the values at the particle 
 !--                   location are the same as the gridpoint values
-                      e_int = ei(num_gp)
-
+                      e_int     = ei(num_gp)
+                      diss_int  = dissi(num_gp)
+                      de_dx_int = de_dxi(num_gp)
+                      de_dy_int = de_dyi(num_gp)
+                      de_dz_int = de_dzi(num_gp)
                    ELSE IF ( num_gp > 1 )  THEN