source: palm/trunk/SOURCE/user_data_output_2d.f90 @ 375

Last change on this file since 375 was 343, checked in by maronga, 16 years ago

adjustments for lcxt4 and ibmy, allow user 2d xy cross section output at z=nzb+1

  • Property svn:keywords set to Id
File size: 1.9 KB
Line 
1 SUBROUTINE user_data_output_2d( av, variable, found, grid, local_pf, two_d )
2
3!------------------------------------------------------------------------------!
4! Actual revisions:
5! -----------------
6!
7!
8! Former revisions:
9! -----------------
10! $Id: user_data_output_2d.f90 343 2009-06-24 12:59:09Z heinze $
11!
12! 211 2008-11-11 04:46:24Z raasch
13! Former file user_interface.f90 split into one file per subroutine
14!
15! Description:
16! ------------
17! Resorts the user-defined output quantity with indices (k,j,i) to a
18! temporary array with indices (i,j,k) and sets the grid on which it is defined.
19! Allowed values for grid are "zu" and "zw".
20!------------------------------------------------------------------------------!
21
22    USE indices
23    USE user
24
25    IMPLICIT NONE
26
27    CHARACTER (LEN=*) ::  grid, variable
28
29    INTEGER ::  av, i, j, k
30
31    LOGICAL ::  found, two_d
32
33    REAL, DIMENSION(nxl-1:nxr+1,nys-1:nyn+1,nzb:nzt+1) ::  local_pf
34
35
36    found = .TRUE.
37
38    SELECT CASE ( TRIM( variable ) )
39
40!
41!--    Uncomment and extend the following lines, if necessary.
42!--    The arrays for storing the user defined quantities (here u2 and u2_av)
43!--    have to be declared and defined by the user!
44!--    Sample for user-defined output:
45!       CASE ( 'u2_xy', 'u2_xz', 'u2_yz' )
46!          IF ( av == 0 )  THEN
47!             DO  i = nxl-1, nxr+1
48!                DO  j = nys-1, nyn+1
49!                   DO  k = nzb, nzt+1
50!                      local_pf(i,j,k) = u2(k,j,i)
51!                   ENDDO
52!                ENDDO
53!             ENDDO
54!          ELSE
55!             DO  i = nxl-1, nxr+1
56!                DO  j = nys-1, nyn+1
57!                   DO  k = nzb, nzt+1
58!                      local_pf(i,j,k) = u2_av(k,j,i)
59!                   ENDDO
60!                ENDDO
61!             ENDDO
62!          ENDIF
63!
64!          grid = 'zu'
65
66
67       CASE DEFAULT
68          found = .FALSE.
69          grid  = 'none'
70
71    END SELECT
72
73
74 END SUBROUTINE user_data_output_2d
75
Note: See TracBrowser for help on using the repository browser.