Ignore:
Timestamp:
Mar 6, 2017 6:16:35 PM (8 years ago)
Author:
suehring
Message:

Bugfix, move setting of topography_grid_convention to init_grid, else, in case of restarts generic topography may not set properly if no value is prescribed

File:
1 edited

Legend:

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

    r2143 r2169  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Bugfix, move setting for topography grid convention to init_grid
    2323!
    2424! Former revisions:
     
    826826                           TRIM( action )
    827827          CALL message( 'check_parameters', 'PA0014', 1, 2, 0, 6, 0 )
    828        ENDIF
    829 !
    830 !--    In case of non-flat topography, check whether the convention how to
    831 !--    define the topography grid has been set correctly, or whether the default
    832 !--    is applicable. If this is not possible, abort.
    833        IF ( TRIM( topography_grid_convention ) == ' ' )  THEN
    834           IF ( TRIM( topography ) /= 'single_building' .AND.                   &
    835                TRIM( topography ) /= 'single_street_canyon' .AND.              &
    836                TRIM( topography ) /= 'read_from_file' )  THEN
    837 !--          The default value is not applicable here, because it is only valid
    838 !--          for the two standard cases 'single_building' and 'read_from_file'
    839 !--          defined in init_grid.
    840              WRITE( message_string, * )                                        &
    841                   'The value for "topography_grid_convention" ',               &
    842                   'is not set. Its default value is & only valid for ',        &
    843                   '"topography" = ''single_building'', ',                      &
    844                   '''single_street_canyon'' & or ''read_from_file''.',         &
    845                   ' & Choose ''cell_edge'' or ''cell_center''.'
    846              CALL message( 'check_parameters', 'PA0239', 1, 2, 0, 6, 0 )
    847           ELSE
    848 !--          The default value is applicable here.
    849 !--          Set convention according to topography.
    850              IF ( TRIM( topography ) == 'single_building' .OR.                 &
    851                   TRIM( topography ) == 'single_street_canyon' )  THEN
    852                 topography_grid_convention = 'cell_edge'
    853              ELSEIF ( TRIM( topography ) == 'read_from_file' )  THEN
    854                 topography_grid_convention = 'cell_center'
    855              ENDIF
    856           ENDIF
    857        ELSEIF ( TRIM( topography_grid_convention ) /= 'cell_edge' .AND.        &
    858                 TRIM( topography_grid_convention ) /= 'cell_center' )  THEN
    859           WRITE( message_string, * )                                           &
    860                'The value for "topography_grid_convention" is ',               &
    861                'not recognized. & Choose ''cell_edge'' or ''cell_center''.'
    862           CALL message( 'check_parameters', 'PA0240', 1, 2, 0, 6, 0 )
    863828       ENDIF
    864829
Note: See TracChangeset for help on using the changeset viewer.