Changeset 1660 for palm/trunk/SOURCE/init_grid.f90
- Timestamp:
- Sep 21, 2015 8:15:16 AM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_grid.f90
r1581 r1660 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! Bugfix: Definition of topography grid levels if vertical grid stretching 23 ! starts below the maximum topography height. 22 24 ! 23 25 ! … … 541 543 blx = NINT( building_length_x / dx ) 542 544 bly = NINT( building_length_y / dy ) 543 bh = NINT( building_height / dz ) 545 bh = 0 546 DO k = 0, nzt 547 IF ( building_height .GT. zw(k) ) THEN 548 bh = k+1 549 ENDIF 550 ENDDO 544 551 545 552 IF ( building_wall_left == 9999999.9_wp ) THEN … … 601 608 ENDIF 602 609 603 ch = NINT( canyon_height / dz ) 610 ch = 0 611 DO k = 0, nzt 612 IF ( canyon_height .GT. zw(k) ) THEN 613 ch = k+1 614 ENDIF 615 ENDDO 604 616 dp_level_ind_b = ch 605 617 ! … … 675 687 DO i = 0, nx 676 688 DO j = 0, ny 677 nzb_local(j,i) = NINT( topo_height(j,i) / dz ) 689 nzb_local(j,i) = 0 690 DO k = 0, nzt 691 IF ( topo_height(j,i) .GT. zw(k) ) THEN 692 nzb_local(j,i) = k+1 693 ENDIF 694 ENDDO 678 695 ENDDO 679 696 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.