Changeset 119 for palm/trunk/SOURCE
 Timestamp:
 Oct 17, 2007 10:27:13 AM (17 years ago)
 Location:
 palm/trunk/SOURCE
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r117 r119 22 22 23 23 Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail 24 numbers along y (advec_particles)24 numbers along y, small bugfixes in the SGS part (advec_particles) 25 25 26 26 Bugfix: model_string needed a default value (combine_plot_fields) 
palm/trunk/SOURCE/advec_particles.f90
r116 r119 8 8 ! Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail 9 9 ! numbers along y 10 ! Small bugfixes in the SGS part 10 11 ! TEST: PRINT statements on unit 9 (commented out) 11 12 ! … … 1056 1057 gp_outside_of_building(1) = 1 1057 1058 location(num_gp,1) = i * dx 1058 location(num_gp,2) = j * d x1059 location(num_gp,2) = j * dy 1059 1060 location(num_gp,3) = k * dz  0.5 * dz 1060 1061 ei(num_gp) = e(k,j,i) … … 1070 1071 gp_outside_of_building(2) = 1 1071 1072 location(num_gp,1) = i * dx 1072 location(num_gp,2) = (j+1) * d x1073 location(num_gp,2) = (j+1) * dy 1073 1074 location(num_gp,3) = k * dz  0.5 * dz 1074 1075 ei(num_gp) = e(k,j+1,i) … … 1513 1514 ! Neumann boundary condition has to be applied 1514 1515 ! (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) THEN1516 IF ( gp_outside_of_building(1) == 0 .AND. & 1517 gp_outside_of_building(3) == 1 ) THEN 1517 1518 num_gp = num_gp + 1 1518 1519 location(num_gp,1) = i * dx … … 1530 1531 ! Neumann boundary condition has to be applied 1531 1532 ! (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) THEN1533 IF ( gp_outside_of_building(5) == 0 .AND. & 1534 gp_outside_of_building(7) == 1 ) THEN 1534 1535 num_gp = num_gp + 1 1535 1536 location(num_gp,1) = (i+1) * dx … … 1547 1548 ! Neumann boundary condition has to be applied 1548 1549 ! (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) THEN1550 IF ( gp_outside_of_building(2) == 0 .AND. & 1551 gp_outside_of_building(4) == 1 ) THEN 1551 1552 num_gp = num_gp + 1 1552 1553 location(num_gp,1) = i * dx … … 1564 1565 ! Neumann boundary condition has to be applied 1565 1566 ! (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) THEN1567 IF ( gp_outside_of_building(6) == 0 .AND. & 1568 gp_outside_of_building(8) == 1 ) THEN 1568 1569 num_gp = num_gp + 1 1569 1570 location(num_gp,1) = (i+1) * dx … … 1584 1585 ! building, it follows that the values at the particle 1585 1586 ! 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) 1588 1592 ELSE IF ( num_gp > 1 ) THEN 1589 1593
Note: See TracChangeset
for help on using the changeset viewer.