Ignore:
Timestamp:
Mar 14, 2018 4:20:52 PM (3 years ago)
Author:
suehring
Message:

Revise informative message concerning filtered topography.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_grid.f90

    r2892 r2893  
    2525! -----------------
    2626! $Id$
     27! Revise informative message concerning filtered topography (1 grid-point
     28! holes).
     29!
     30! 2892 2018-03-14 15:06:29Z suehring
    2731! Bugfix, uninitialized array in case of single_building.
    2832!
     
    11601164    IMPLICIT NONE
    11611165
     1166    LOGICAL      ::  filled = .FALSE. !< flag indicating if holes were filled
     1167
    11621168    INTEGER(iwp) ::  i          !< running index along x-direction
    1163     INTEGER(iwp) ::  it         !< counter for number of iterations
    11641169    INTEGER(iwp) ::  j          !< running index along y-direction
    11651170    INTEGER(iwp) ::  k          !< running index along z-direction
     
    11791184
    11801185    num_hole = 99999
    1181     it = 0
    11821186    DO WHILE ( num_hole > 0 )       
    11831187
     
    12001204       ENDIF
    12011205
    1202        it = it + 1
    12031206       num_hole_l = 0
    12041207       DO i = nxl, nxr
     
    13011304       num_hole = num_hole_l
    13021305#endif   
    1303 
    1304 !
    1305 !--    Create an informative message if any holes were filled.
    1306        IF ( num_hole > 0 )  THEN
    1307           WRITE( message_string, * ) num_hole, 'hole(s) resolved by only '//   &
    1308                                      'one grid point were filled at ', it,     &
    1309                                      'th iteration '
    1310           CALL message( 'init_grid', 'PA0430', 0, 0, 0, 6, 0 )
    1311        ENDIF
    1312 
     1306       IF ( num_hole > 0  .AND.  .NOT. filled )  filled = .TRUE.
    13131307
    13141308    ENDDO
     1309!
     1310!-- Create an informative message if any holes were filled.
     1311    IF ( filled )  THEN
     1312       WRITE( message_string, * ) 'Topography was filtered, i.e. holes ' //    &
     1313                                  'resolved by only one grid point '     //    &
     1314                                  'were filled during initialization.'
     1315       CALL message( 'init_grid', 'PA0430', 0, 0, 0, 6, 0 )
     1316    ENDIF
    13151317
    13161318    DEALLOCATE( topo_tmp )
     
    15341536    INTEGER(iwp) ::  cys           !< index for south canyon wall
    15351537    INTEGER(iwp) ::  i             !< index variable along x
    1536     INTEGER(iwp) ::  it            !< Number of iterations required to fill all problematic holes
    15371538    INTEGER(iwp) ::  j             !< index variable along y
    15381539    INTEGER(iwp) ::  k             !< index variable along z
Note: See TracChangeset for help on using the changeset viewer.