Changeset 1322 for palm/trunk/SOURCE/init_grid.f90
- Timestamp:
- Mar 20, 2014 4:38:49 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_grid.f90
r1321 r1322 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants defined as wp-kind 23 23 ! 24 24 ! Former revisions: … … 377 377 !-- Compute the grid-dependent mixing length. 378 378 DO k = 1, nzt 379 l_grid(k) = ( dx * dy * dzw(k) )**0.33333333333333 379 l_grid(k) = ( dx * dy * dzw(k) )**0.33333333333333_wp 380 380 ENDDO 381 381 … … 516 516 bh = NINT( building_height / dz ) 517 517 518 IF ( building_wall_left == 9999999.9 ) THEN518 IF ( building_wall_left == 9999999.9_wp ) THEN 519 519 building_wall_left = ( nx + 1 - blx ) / 2 * dx 520 520 ENDIF … … 522 522 bxr = bxl + blx 523 523 524 IF ( building_wall_south == 9999999.9 ) THEN524 IF ( building_wall_south == 9999999.9_wp ) THEN 525 525 building_wall_south = ( ny + 1 - bly ) / 2 * dy 526 526 ENDIF … … 547 547 !-- Single quasi-2D street canyon of infinite length in x or y direction. 548 548 !-- The canyon is centered in the other direction by default. 549 IF ( canyon_width_x /= 9999999.9 ) THEN549 IF ( canyon_width_x /= 9999999.9_wp ) THEN 550 550 ! 551 551 !-- Street canyon in y direction 552 552 cwx = NINT( canyon_width_x / dx ) 553 IF ( canyon_wall_left == 9999999.9 ) THEN553 IF ( canyon_wall_left == 9999999.9_wp ) THEN 554 554 canyon_wall_left = ( nx + 1 - cwx ) / 2 * dx 555 555 ENDIF … … 557 557 cxr = cxl + cwx 558 558 559 ELSEIF ( canyon_width_y /= 9999999.9 ) THEN559 ELSEIF ( canyon_width_y /= 9999999.9_wp ) THEN 560 560 ! 561 561 !-- Street canyon in x direction 562 562 cwy = NINT( canyon_width_y / dy ) 563 IF ( canyon_wall_south == 9999999.9 ) THEN563 IF ( canyon_wall_south == 9999999.9_wp ) THEN 564 564 canyon_wall_south = ( ny + 1 - cwy ) / 2 * dy 565 565 ENDIF … … 578 578 ! 579 579 !-- Street canyon size has to meet some requirements 580 IF ( canyon_width_x /= 9999999.9 ) THEN580 IF ( canyon_width_x /= 9999999.9_wp ) THEN 581 581 IF ( ( cxl < 1 ) .OR. ( cxr > nx-1 ) .OR. ( cwx < 3 ) .OR. & 582 582 ( ch < 3 ) ) THEN … … 587 587 CALL message( 'init_grid', 'PA0205', 1, 2, 0, 6, 0 ) 588 588 ENDIF 589 ELSEIF ( canyon_width_y /= 9999999.9 ) THEN589 ELSEIF ( canyon_width_y /= 9999999.9_wp ) THEN 590 590 IF ( ( cys < 1 ) .OR. ( cyn > ny-1 ) .OR. ( cwy < 3 ) .OR. & 591 591 ( ch < 3 ) ) THEN … … 597 597 ENDIF 598 598 ENDIF 599 IF ( canyon_width_x /= 9999999.9 .AND. canyon_width_y /= 9999999.9) &599 IF ( canyon_width_x /= 9999999.9_wp .AND. canyon_width_y /= 9999999.9_wp ) & 600 600 THEN 601 601 message_string = 'inconsistent canyon parameters:' // & … … 606 606 607 607 nzb_local = ch 608 IF ( canyon_width_x /= 9999999.9 ) THEN608 IF ( canyon_width_x /= 9999999.9_wp ) THEN 609 609 nzb_local(:,cxl+1:cxr-1) = 0 610 ELSEIF ( canyon_width_y /= 9999999.9 ) THEN610 ELSEIF ( canyon_width_y /= 9999999.9_wp ) THEN 611 611 nzb_local(cys+1:cyn-1,:) = 0 612 612 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.