ETS  \$Id: Doxyfile 2162 2020-02-26 14:16:09Z g2dpc $
 All Classes Files Functions Variables Pages
plot_data.f90
Go to the documentation of this file.
1 module plot_data
2 
3  use itm_types
4 
5  implicit none
6 
7  character(len = 132) :: path
8 
9  contains
10 
11  subroutine plot_data_1d(type_plot, x, y, np, printname, z)
12 !----------------------------------------------------------|
13 ! writing of plot data for xmgr or idl |
14 ! |
15 ! input: |
16 ! type_plot : type of plot |
17 ! line: n sets of lenght np of the ordinate y vs. the |
18 ! same abscissa x |
19 ! 2d: n sets of lenght np of the ordinate y vs. x |
20 ! 2d_t: n sets of lenght np of the ordinate y vs. |
21 ! x (transposed) |
22 ! grid: n sets of lenght np of the ordinate y vs. x |
23 ! and x (transposed) to generate a xy-grid |
24 ! 3d: writes data (x, y, z) for 3d plots |
25 ! x, y, z (optional, only if type_plot == "3d") : |
26 ! arrays of reals |
27 ! np: length of a part of x, y and z |
28 ! printname: name of output file |
29 ! |
30 ! REG 12.04.2007 |
31 !----------------------------------------------------------|
32 
33  character (len=*), intent(in) :: type_plot
34  real(r8), dimension(:), intent(in) :: x, y
35  integer(itm_i4), intent(in) :: np
36  character (len=*), intent(in) :: printname
37  real(r8), dimension(:), intent(in), optional :: z
38 
39  character (len=20) :: type_pl
40  character (len=100) :: filename
41  integer(itm_i4) :: i, j, n, n_y_dim, i_error
42  integer(itm_i4), parameter :: out_ppc = 95
43 
44  n_y_dim = size(y)
45  n = n_y_dim / np
46 
47  type_pl = adjustl(trim(type_plot))
48 
49  filename = trim(printname) // '.' // adjustl('dat')
50  open (unit = out_ppc, file = trim(adjustl(path)) // 'plots/' &
51  // filename, status = 'replace', form = 'formatted', action = 'write', &
52  iostat = i_error)
53 
54  if (type_pl == "line") then
55  if (n == 1) then
56  do i = 1, np
57  write(out_ppc, 1) x(i), y(i)
58  end do
59  else
60  do i = 1, n
61  do j = 1, np
62  write(out_ppc, 1) x(j), y((j - 1) * n + i)
63  end do
64  write(out_ppc, 2)
65  end do
66  end if
67  write(out_ppc, 3)
68  else
69  if (type_pl == "2d" .or. type_pl == "grid") then
70  do i = 1, n
71  do j = 1, np
72  write(out_ppc, 1) x((i - 1) * np + j), y((i - 1) * np + j)
73  end do
74  write(out_ppc, 1) x((i - 1) * np + 1), y((i - 1) * np + 1)
75  write(out_ppc, 2)
76  end do
77  end if
78  if (type_pl == "2d_t" .or. type_pl == "grid") then
79  do j = 1, np
80  do i = 1, n
81  write(out_ppc, 1) x((i - 1) * np + j), y((i - 1) * np + j)
82  end do
83  write(out_ppc, 2)
84  end do
85  end if
86  if (type_pl == "3d") then
87  write(out_ppc, 4) n, np
88  do i = 1, n * np
89  write(out_ppc, 5) x(i), y(i), z(i)
90  end do
91  else
92  write(out_ppc, 3)
93  end if
94  end if
95 
96  close(out_ppc)
97 
98  1 format(1x, e16.8, 3x, e16.8)
99  2 format('&')
100  3 format('&&')
101  4 format(1x, i4.4, 1x, i4.4)
102  5 format(1x, e16.8, 3x, e16.8, 3x, e16.8)
103 
104  return
105 
106  end subroutine plot_data_1d
107 
108 
109  subroutine plot_data_2d(type_plot, x, y, printname, z)
110 !----------------------------------------------------------
111 ! writing of plot data for xmgr or idl
112 !
113 ! input:
114 ! type_plot : type of plot
115 ! line: n_chi sets of lenght n_r of the ordinate y vs. the
116 ! same abscissa x
117 ! 3d: writes data (x, y, z) for 3d plots
118 ! x, y, z (optional, only if type_plot == "3d") :
119 ! arrays of reals
120 ! printname: name of output file
121 !
122 ! REG 19.02.2008
123 !----------------------------------------------------------
124 
125  character (len = *), intent(in) :: type_plot
126  real(r8), dimension(:, :), intent(in) :: x, y
127  character (len = *), intent(in) :: printname
128  real(r8), dimension(:, :), intent(in), optional :: z
129 
130  character (len = 20) :: type_pl
131  character (len = 100) :: filename
132  integer(itm_i4) :: i, j, i_error
133  integer(itm_i4) :: n_chi, n_r
134  integer(itm_i4), parameter :: out_ppc = 95
135 
136  n_r = size(y, 1)
137  n_chi = size(y, 2)
138 
139  type_pl = adjustl(trim(type_plot))
140 
141  filename = trim(printname) // '.' // adjustl('dat')
142  open (unit = out_ppc, file = trim(adjustl(path)) // 'plots/' &
143  // filename, status = 'replace', form = 'formatted', action = 'write', &
144  iostat = i_error)
145 
146  if (type_pl == "line") then
147  do j = 1, n_chi
148  do i = 1, n_r
149  write(out_ppc, 1) x(i, 1), y(i, j)
150  end do
151  write(out_ppc, 2)
152  end do
153  write(out_ppc, 3)
154  else
155  if (type_pl == "2d" .or. type_pl == "grid") then
156  do i = 1, n_r
157  do j = 1, n_chi
158  write(out_ppc, 1) x(i, j), y(i, j)
159  end do
160  write(out_ppc, 1) x(i, 1), y(i, 1)
161  write(out_ppc, 2)
162  end do
163  end if
164  if (type_pl == "2d_t" .or. type_pl == "grid") then
165  do j = 1, n_chi
166  do i = 1, n_r
167  write(out_ppc, 1) x(i, j), y(i, j)
168  end do
169  write(out_ppc, 2)
170  end do
171  end if
172  if (type_pl == "3d") then
173  write(out_ppc, 4) n_chi, n_r
174  do i = 1, n_r
175  do j = 1, n_chi
176  write(out_ppc, 5) x(i, j), y(i, j), z(i, j)
177  end do
178  end do
179  else
180  write(out_ppc, 3)
181  end if
182  end if
183 
184  close(out_ppc)
185 
186  1 format(1x, e16.8, 3x, e16.8)
187  2 format('&')
188  3 format('&&')
189  4 format(1x, i4.4, 1x, i4.4)
190  5 format(1x, e16.8, 3x, e16.8, 3x, e16.8)
191 
192  return
193 
194  end subroutine plot_data_2d
195 
196 
197 end module plot_data
plot_data::plot_data_2d
subroutine plot_data_2d(type_plot, x, y, printname, z)
Definition: plot_data.f90:109
plot_data
Definition: plot_data.f90:1
plot_data::plot_data_1d
subroutine plot_data_1d(type_plot, x, y, np, printname, z)
Definition: plot_data.f90:11