Changeset 256 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Mar 8, 2009 8:56:27 AM (15 years ago)

Author:

letzel

Message:

• topography_grid_convention moved from userpar to inipar
• documentation and examples updated

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (2 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -6 +6 @@
 ! Bugfix: pressure included for profile output
 ! Check pressure gradient conditions
+! topography_grid_convention moved from user_check_parameters
+! 'single_street_canyon'
 !
 ! Former revisions:
@@ -324 +326 @@
           CALL message( 'check_parameters', 'PA0014', 1, 2, 0, 6, 0 )
        ENDIF
+!
+!--    In case of non-flat topography, check whether the convention how to 
+!--    define the topography grid has been set correctly, or whether the default
+!--    is applicable. If this is not possible, abort.
+       IF ( TRIM( topography_grid_convention ) == ' ' )  THEN
+          IF ( TRIM( topography ) /= 'single_building' .AND.  &
+               TRIM( topography ) /= 'single_street_canyon' .AND.  &
+               TRIM( topography ) /= 'read_from_file' )  THEN
+!--          The default value is not applicable here, because it is only valid
+!--          for the two standard cases 'single_building' and 'read_from_file'
+!--          defined in init_grid.
+             WRITE( message_string, * )  &
+                  'The value for "topography_grid_convention" ',  &
+                  'is not set. Its default value is & only valid for ',  &
+                  '"topography" = ''single_building'', ',  &
+                  '''single_street_canyon'' & or ''read_from_file''.',  &
+                  ' & Choose ''cell_edge'' or ''cell_center''.'
+             CALL message( 'user_check_parameters', 'PA0239', 1, 2, 0, 6, 0 )
+          ELSE
+!--          The default value is applicable here.
+!--          Set convention according to topography.
+             IF ( TRIM( topography ) == 'single_building' .OR.  &
+                  TRIM( topography ) == 'single_street_canyon' )  THEN
+                topography_grid_convention = 'cell_edge'
+             ELSEIF ( TRIM( topography ) == 'read_from_file' )  THEN
+                topography_grid_convention = 'cell_center'
+             ENDIF
+          ENDIF
+       ELSEIF ( TRIM( topography_grid_convention ) /= 'cell_edge' .AND.  &
+                TRIM( topography_grid_convention ) /= 'cell_center' )  THEN
+          WRITE( message_string, * )  &
+               'The value for "topography_grid_convention" is ', &
+               'not recognized. & Choose ''cell_edge'' or ''cell_center''.'
+          CALL message( 'user_check_parameters', 'PA0240', 1, 2, 0, 6, 0 )
+       ENDIF
+
     ENDIF