Ignore:
Timestamp:
Mar 8, 2009 8:56:27 AM (14 years ago)
Author:
letzel
Message:
  • topography_grid_convention moved from userpar to inipar
  • documentation and examples updated
File:
1 edited

Legend:

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

    r241 r256  
    66! Bugfix: pressure included for profile output
    77! Check pressure gradient conditions
     8! topography_grid_convention moved from user_check_parameters
     9! 'single_street_canyon'
    810!
    911! Former revisions:
     
    324326          CALL message( 'check_parameters', 'PA0014', 1, 2, 0, 6, 0 )
    325327       ENDIF
     328!
     329!--    In case of non-flat topography, check whether the convention how to
     330!--    define the topography grid has been set correctly, or whether the default
     331!--    is applicable. If this is not possible, abort.
     332       IF ( TRIM( topography_grid_convention ) == ' ' )  THEN
     333          IF ( TRIM( topography ) /= 'single_building' .AND.  &
     334               TRIM( topography ) /= 'single_street_canyon' .AND.  &
     335               TRIM( topography ) /= 'read_from_file' )  THEN
     336!--          The default value is not applicable here, because it is only valid
     337!--          for the two standard cases 'single_building' and 'read_from_file'
     338!--          defined in init_grid.
     339             WRITE( message_string, * )  &
     340                  'The value for "topography_grid_convention" ',  &
     341                  'is not set. Its default value is & only valid for ',  &
     342                  '"topography" = ''single_building'', ',  &
     343                  '''single_street_canyon'' & or ''read_from_file''.',  &
     344                  ' & Choose ''cell_edge'' or ''cell_center''.'
     345             CALL message( 'user_check_parameters', 'PA0239', 1, 2, 0, 6, 0 )
     346          ELSE
     347!--          The default value is applicable here.
     348!--          Set convention according to topography.
     349             IF ( TRIM( topography ) == 'single_building' .OR.  &
     350                  TRIM( topography ) == 'single_street_canyon' )  THEN
     351                topography_grid_convention = 'cell_edge'
     352             ELSEIF ( TRIM( topography ) == 'read_from_file' )  THEN
     353                topography_grid_convention = 'cell_center'
     354             ENDIF
     355          ENDIF
     356       ELSEIF ( TRIM( topography_grid_convention ) /= 'cell_edge' .AND.  &
     357                TRIM( topography_grid_convention ) /= 'cell_center' )  THEN
     358          WRITE( message_string, * )  &
     359               'The value for "topography_grid_convention" is ', &
     360               'not recognized. & Choose ''cell_edge'' or ''cell_center''.'
     361          CALL message( 'user_check_parameters', 'PA0240', 1, 2, 0, 6, 0 )
     362       ENDIF
     363
    326364    ENDIF
    327365
Note: See TracChangeset for help on using the changeset viewer.