source: palm/trunk/SOURCE/user_init_grid.f90 @ 357

Last change on this file since 357 was 258, checked in by heinze, 16 years ago

Output of messages replaced by message handling routine.

  • Property svn:keywords set to Id
File size: 2.0 KB
RevLine 
[211]1 SUBROUTINE user_init_grid( gls, nzb_local )
2
3!------------------------------------------------------------------------------!
[258]4! Current revisions:
[211]5! -----------------
[258]6! Output of messages replaced by message handling routine.
[240]7! add 'single_street_canyon' as standard topography case
[256]8! -topography_grid_convention
[211]9!
10! Former revisions:
11! -----------------
12! $Id: user_init_grid.f90 258 2009-03-13 12:36:03Z heinze $
13!
[226]14! 217 2008-12-09 18:00:48Z letzel
15! +topography_grid_convention
16! Former file user_interface.f90 split into one file per subroutine
17!
[211]18! Description:
19! ------------
20! Execution of user-defined grid initializing actions
21! First argument gls contains the number of ghost layers, which is > 1 if the
22! multigrid method for the pressure solver is used
23!------------------------------------------------------------------------------!
24
25    USE control_parameters
26    USE indices
27    USE user
28
29    IMPLICIT NONE
30
31    INTEGER ::  gls
32
33    INTEGER, DIMENSION(-gls:ny+gls,-gls:nx+gls) ::  nzb_local
34
35!
36!-- Here the user-defined grid initializing actions follow:
37
38!
39!-- Set the index array nzb_local for non-flat topography.
40!-- Here consistency checks concerning domain size and periodicity are necessary
41    SELECT CASE ( TRIM( topography ) )
42
[240]43       CASE ( 'flat', 'single_building', 'single_street_canyon' )
[211]44!
45!--       Not allowed here since these are the standard cases used in init_grid.
46
47       CASE ( 'user_defined_topography_1' )
48!
[256]49!--       Here the user can define his own topography.
[217]50!--       After definition, please remove the following three lines!
[258]51          message_string = 'topography "' // topography // '" not available yet'
52          CALL message( 'user_init_grid', 'UI0005', 1, 2, 0, 6, 0 )
[211]53
54       CASE DEFAULT
55!
56!--       The DEFAULT case is reached if the parameter topography contains a
57!--       wrong character string that is neither recognized in init_grid nor
58!--       here in user_init_grid.
[258]59          message_string = 'unknown topography "' // topography // '"'
60          CALL message( 'user_init_grid', 'UI0006', 1, 2, 0, 6, 0 )
[211]61
62    END SELECT
63
64
65 END SUBROUTINE user_init_grid
66
Note: See TracBrowser for help on using the repository browser.