source: palm/trunk/SOURCE/user_check_data_output_pr.f90 @ 1836

Last change on this file since 1836 was 1818, checked in by maronga, 9 years ago

last commit documented / copyright update

  • Property svn:keywords set to Id
File size: 3.7 KB
RevLine 
[1682]1!> @file user_check_data_output_pr.f90
[1036]2!--------------------------------------------------------------------------------!
3! This file is part of PALM.
4!
5! PALM is free software: you can redistribute it and/or modify it under the terms
6! of the GNU General Public License as published by the Free Software Foundation,
7! either version 3 of the License, or (at your option) any later version.
8!
9! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
10! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
11! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
12!
13! You should have received a copy of the GNU General Public License along with
14! PALM. If not, see <http://www.gnu.org/licenses/>.
15!
[1818]16! Copyright 1997-2016 Leibniz Universitaet Hannover
[1036]17!--------------------------------------------------------------------------------!
18!
[484]19! Current revisions:
[211]20! -----------------
[1321]21!
[1784]22!
[1321]23! Former revisions:
24! -----------------
25! $Id: user_check_data_output_pr.f90 1818 2016-04-06 15:53:27Z raasch $
26!
[1784]27! 1783 2016-03-06 18:36:17Z raasch
28! netcdf module name changed + relatec changes
29!
[1683]30! 1682 2015-10-07 23:56:08Z knoop
31! Code annotations made doxygen readable
32!
[1321]33! 1320 2014-03-20 08:40:49Z raasch
[1320]34! kind-parameters added to all INTEGER and REAL declaration statements,
35! kinds are defined in new module kinds,
36! revision history before 2012 removed,
37! comment fields (!:) to be used for variable explanations added to
38! all variable declaration statements
[211]39!
[1037]40! 1036 2012-10-22 13:43:42Z raasch
41! code put under GPL (PALM 3.9)
42!
[226]43! 211 2008-11-11 04:46:24Z raasch
44! Former file user_interface.f90 split into one file per subroutine
45!
[211]46! Description:
47! ------------
[1682]48!> Set the unit of user defined profile output quantities. For those variables
49!> not recognized by the user, the parameter unit is set to "illegal", which
50!> tells the calling routine that the output variable is not defined and leads
51!> to a program abort.
[211]52!------------------------------------------------------------------------------!
[1682]53 SUBROUTINE user_check_data_output_pr( variable, var_count, unit )
54 
[211]55
56    USE arrays_3d
[1320]57
[211]58    USE indices
[1320]59
60    USE kinds
61
[1783]62    USE netcdf_interface,                                                      &
63        ONLY:  dopr_unit
[1320]64
[211]65    USE profil_parameter
[1320]66
[211]67    USE statistics
[1320]68
[211]69    USE user
70
71    IMPLICIT NONE
72
[1682]73    CHARACTER (LEN=*) ::  unit     !<
74    CHARACTER (LEN=*) ::  variable !<
[211]75
[1682]76    INTEGER(iwp) ::  user_pr_index !<
77    INTEGER(iwp) ::  var_count     !<
[211]78
79    SELECT CASE ( TRIM( variable ) )
80
81!
82!--    Uncomment and extend the following lines, if necessary.
83!--    Add additional CASE statements depending on the number of quantities
84!--    for which profiles are to be calculated. The respective calculations
85!--    to be performed have to be added in routine user_statistics.
86!--    The quantities are (internally) identified by a user-profile-number
[1320]87!--    (see variable "user_pr_index" below). The first user-profile must be assigned
[211]88!--    the number "pr_palm+1", the second one "pr_palm+2", etc. The respective
89!--    user-profile-numbers have also to be used in routine user_statistics!
90!       CASE ( 'u*v*' )                      ! quantity string as given in
91!                                            ! data_output_pr_user
[1320]92!          user_pr_index = pr_palm + 1
93!          dopr_index(var_count)  = user_pr_index    ! quantities' user-profile-number
[211]94!          dopr_unit(var_count)   = 'm2/s2'  ! quantity unit
[1320]95!          hom(:,2,user_pr_index,:)       = SPREAD( zu, 2, statistic_regions+1 )
[211]96!                                            ! grid on which the quantity is
97!                                            ! defined (use zu or zw)
98
99       CASE DEFAULT
100          unit = 'illegal'
101
102    END SELECT
103
104
105 END SUBROUTINE user_check_data_output_pr
106
Note: See TracBrowser for help on using the repository browser.