Changeset 119 for palm/trunk/SOURCE


Ignore:
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:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r117 r119  
    2222
    2323Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
    24 numbers along y (advec_particles)
     24numbers along y, small bugfixes in the SGS part (advec_particles)
    2525
    2626Bugfix: model_string needed a default value (combine_plot_fields)
  • palm/trunk/SOURCE/advec_particles.f90

    r116 r119  
    88! Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
    99! numbers along y
     10! Small bugfixes in the SGS part
    1011! TEST: PRINT statements on unit 9 (commented out)
    1112!
     
    10561057                   gp_outside_of_building(1) = 1
    10571058                   location(num_gp,1) = i * dx
    1058                    location(num_gp,2) = j * dx
     1059                   location(num_gp,2) = j * dy
    10591060                   location(num_gp,3) = k * dz - 0.5 * dz
    10601061                   ei(num_gp)     = e(k,j,i)
     
    10701071                   gp_outside_of_building(2) = 1
    10711072                   location(num_gp,1) = i * dx
    1072                    location(num_gp,2) = (j+1) * dx
     1073                   location(num_gp,2) = (j+1) * dy
    10731074                   location(num_gp,3) = k * dz - 0.5 * dz
    10741075                   ei(num_gp)     = e(k,j+1,i)
     
    15131514!--                Neumann boundary condition has to be applied
    15141515!--                (only one case as only building beneath is possible)
    1515                    IF ( gp_outside_of_building(1) == 1  .AND. &
    1516                         gp_outside_of_building(3) == 0 )  THEN
     1516                   IF ( gp_outside_of_building(1) == 0  .AND. &
     1517                        gp_outside_of_building(3) == 1 )  THEN
    15171518                      num_gp = num_gp + 1
    15181519                      location(num_gp,1) = i * dx
     
    15301531!--                Neumann boundary condition has to be applied
    15311532!--                (only one case as only building beneath is possible)
    1532                    IF ( gp_outside_of_building(5) == 1  .AND. &
    1533                         gp_outside_of_building(7) == 0 )  THEN
     1533                   IF ( gp_outside_of_building(5) == 0  .AND. &
     1534                        gp_outside_of_building(7) == 1 )  THEN
    15341535                      num_gp = num_gp + 1
    15351536                      location(num_gp,1) = (i+1) * dx
     
    15471548!--                Neumann boundary condition has to be applied
    15481549!--                (only one case as only building beneath is possible)
    1549                    IF ( gp_outside_of_building(2) == 1  .AND. &
    1550                         gp_outside_of_building(4) == 0 )  THEN
     1550                   IF ( gp_outside_of_building(2) == 0  .AND. &
     1551                        gp_outside_of_building(4) == 1 )  THEN
    15511552                      num_gp = num_gp + 1
    15521553                      location(num_gp,1) = i * dx
     
    15641565!--                Neumann boundary condition has to be applied
    15651566!--                (only one case as only building beneath is possible)
    1566                    IF ( gp_outside_of_building(6) == 1  .AND. &
    1567                         gp_outside_of_building(8) == 0 )  THEN
     1567                   IF ( gp_outside_of_building(6) == 0  .AND. &
     1568                        gp_outside_of_building(8) == 1 )  THEN
    15681569                      num_gp = num_gp + 1
    15691570                      location(num_gp,1) = (i+1) * dx
     
    15841585!--                   building, it follows that the values at the particle
    15851586!--                   location are the same as the gridpoint values
    1586                       e_int = ei(num_gp)
    1587 
     1587                      e_int     = ei(num_gp)
     1588                      diss_int  = dissi(num_gp)
     1589                      de_dx_int = de_dxi(num_gp)
     1590                      de_dy_int = de_dyi(num_gp)
     1591                      de_dz_int = de_dzi(num_gp)
    15881592                   ELSE IF ( num_gp > 1 )  THEN
    15891593
Note: See TracChangeset for help on using the changeset viewer.