OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
plot_curve.F File Reference
#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "scr17_c.inc"
#include "scr03_c.inc"
#include "units_c.inc"
#include "warn_c.inc"
#include "param_c.inc"
#include "remesh_c.inc"
#include "sphcom.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine plot_curve (x, y, nb_points, input_size_x, input_size_y, input_x_minvalue, input_y_minvalue, input_x_maxvalue, input_y_maxvalue, input_symbol, input_curve_type, input_txt_x, input_txt_y)

Function/Subroutine Documentation

◆ plot_curve()

subroutine plot_curve ( x,
y,
integer nb_points,
integer, optional input_size_x,
integer, optional input_size_y,
optional input_x_minvalue,
optional input_y_minvalue,
optional input_x_maxvalue,
optional input_y_maxvalue,
character*1, optional input_symbol,
integer, optional input_curve_type,
character(*), optional input_txt_x,
character(*), optional input_txt_y )

Definition at line 32 of file plot_curve.F.

36c
37C-----------------------------------------------
38 USE message_mod
39 USE names_and_titles_mod, ONLY: ncharline
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C A n a l y s e M o d u l e
46C-----------------------------------------------
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "com01_c.inc"
51#include "com04_c.inc"
52#include "scr17_c.inc"
53#include "scr03_c.inc"
54#include "units_c.inc"
55#include "warn_c.inc"
56#include "param_c.inc"
57#include "remesh_c.inc"
58#include "sphcom.inc"
59C-----------------------------------------------
60C D u m m y A r g u m e n t s
61C-----------------------------------------------
62 my_real x(nb_points),y(nb_points)
63 INTEGER NB_POINTS,INPUT_SIZE_X,INPUT_SIZE_Y, INPUT_CURVE_TYPE
64 my_real input_x_minvalue, input_x_maxvalue, input_y_minvalue, input_y_maxvalue
65 CHARACTER*1 INPUT_SYMBOL
66 CHARACTER(*) INPUT_TXT_X, INPUT_TXT_Y
67C-----------------------------------------------
68 OPTIONAL :: input_size_x, input_size_y, input_x_minvalue, input_y_minvalue,
69 . input_x_maxvalue, input_y_maxvalue, input_symbol,
70 . input_curve_type, input_txt_x, input_txt_y
71C-----------------------------------------------
72C L o c a l V a r i a b l e s
73C-----------------------------------------------
74 CHARACTER(len=60) :: FMTA
75 INTEGER I,J,II,CURVE(NB_POINTS),SIZE_X,SIZE_Y,CURVE_TYPE,COORDX,COORDY,LENA
76 my_real incr_x,incr_y,x_minvalue, x_maxvalue, y_minvalue, y_maxvalue
77 CHARACTER(LEN=NCHARLINE) :: CHAR(100)
78 CHARACTER(LEN=NCHARLINE) :: CHAR1(100)
79 CHARACTER(LEN=NCHARLINE) :: CHAR2
80 CHARACTER(LEN=NCHARLINE) :: TXT_X
81 CHARACTER(LEN=NCHARLINE) :: TXT_Y
82 CHARACTER*1 SYMBOL
83 INTEGER, DIMENSION(:,:), ALLOCATABLE :: GRID
84C--------------------------------------------------------------
85C GENERIC ROUTINE TO PLOT CURVES OR HISTOGRAMS IN .out files
86C--------------------------------------------------------------
87C MANDATORY
88C NB_POINTS : NUMBER OF POINTS OF THE CURVE TO PLOT
89C X : X COORDS
90C Y : Y COORDS
91C OPTIONAL :
92C INPUT_SIZE_X : X SIZE (NUMBER OF CHAR) OF THE PLOT, DEFAULT = 100
93C INPUT_SIZE_Y : Y SIZE (NUMBER OF CHAR) OF THE PLOT, DEFAULT = 30
94C INPUT_X_MINVALUE : MIN VALUE OF THE ABSCISSA , DEFAULT : VALUE IS SET AUTOMATICALLY
95C INPUT_Y_MINVALUE : MIN VALUE OF THE ORDINATE , DEFAULT : VALUE IS SET AUTOMATICALLY
96C INPUT_X_MAXVALUE : MAX VALUE OF THE ABSCISSA , DEFAULT : VALUE IS SET AUTOMATICALLY
97C INPUT_Y_MAXVALUE : MAX VALUE OF THE ORDINATE , DEFAULT : VALUE IS SET AUTOMATICALLY
98C INPUT_SYMBOL : CHARACTER TO PLOT THE CURVE/HISTOGRAM, DEFAULT = '#'
99C INPUT_CURVE_TYPE , DEFAULT = 0'
100C 0 : PLOT CURVE
101C 1 : PLOT HISTOGRAM
102C INPUT_TXT_X : ABSCISSA LEGEND , DEFAULT : NONE
103C INPUT_TXT_Y : ORDINATE LEGEND , DEFAULT : NONE
104C=======================================================================
105 IF ( .NOT. PRESENT(input_size_x)) THEN
106 size_x = 100
107 ELSE
108 size_x = input_size_x
109 ENDIF
110c
111 IF ( .NOT. PRESENT(input_size_y)) THEN
112 size_y = 30
113 ELSE
114 size_y = input_size_y
115 ENDIF
116c
117 IF ( .NOT. PRESENT(input_x_minvalue)) THEN
118 x_minvalue = ep20
119 DO i=1,nb_points
120 x_minvalue = min(x(i),x_minvalue)
121 ENDDO
122 ELSE
123 x_minvalue = input_x_minvalue
124 ENDIF
125c
126 IF ( .NOT. PRESENT(input_y_minvalue)) THEN
127 y_minvalue = ep20
128 DO i=1,nb_points
129 y_minvalue = min(y(i),y_minvalue)
130 ENDDO
131 ELSE
132 y_minvalue = input_y_minvalue
133 ENDIF
134c
135 IF ( .NOT. PRESENT(input_x_maxvalue)) THEN
136 x_maxvalue = -ep20
137 DO i=1,nb_points
138 x_maxvalue = max(x(i),x_maxvalue)
139 ENDDO
140 ELSE
141 x_maxvalue = input_x_maxvalue
142 ENDIF
143c
144 IF ( .NOT. PRESENT(input_y_maxvalue)) THEN
145 y_maxvalue = -ep20
146 DO i=1,nb_points
147 y_maxvalue = max(y(i),y_maxvalue)
148 ENDDO
149 ELSE
150 y_maxvalue = input_y_maxvalue
151 ENDIF
152c
153 IF ( .NOT. PRESENT(input_symbol)) THEN
154 symbol(1:1) = '#'
155 ELSE
156 symbol(1:1) = input_symbol(1:1)
157 ENDIF
158c
159 IF ( .NOT. PRESENT(input_curve_type)) THEN
160 curve_type = 0
161 ELSE
162 curve_type = input_curve_type
163 ENDIF
164c
165 IF ( .NOT. PRESENT(input_txt_x)) THEN
166 txt_x = ''
167 ELSE
168 txt_x = ''
169 txt_x(1:len_trim(input_txt_x)) = input_txt_x(1:len_trim(input_txt_x))
170 ENDIF
171c
172 IF ( .NOT. PRESENT(input_txt_y)) THEN
173 txt_y = ''
174 ELSE
175 txt_y = ''
176 txt_y(1:len_trim(input_txt_y)) = input_txt_y(1:len_trim(input_txt_y))
177 ENDIF
178C=======================================================================
179 char=''
180 char1=''
181 char2=''
182 curve(1:nb_points) = 0
183 incr_x = (x_maxvalue-x_minvalue)/(size_x-1)
184 incr_y = (y_maxvalue-y_minvalue)/(size_y-1)
185 IF(incr_x == 0) incr_x = 1.0e-6
186 IF(incr_y == 0) incr_y = 1.0e-6
187c
188 ALLOCATE(grid(size_x,size_y))
189 DO j = 1, size_y
190 DO i = 1, size_x
191 grid(i,j) = 0
192 ENDDO
193 ENDDO
194c
195 IF(curve_type == 0)THEN
196 DO i=1,nb_points
197 coordx = min(size_x,max(1,int((x(i) - x_minvalue)/incr_x)))
198 coordy = min(size_y,max(1,int((y(i) - y_minvalue)/incr_y)))
199 grid(coordx,coordy) = 1
200 ENDDO
201 ELSEIF(curve_type == 1)THEN
202 DO i=1,nb_points
203 coordx = min(size_x,max(1,int((x(i) - x_minvalue)/incr_x)))
204 coordy = min(size_y,max(1,int((y(i) - y_minvalue)/incr_y)))
205 DO j=1,coordy
206 grid(coordx,j) = 1
207 ENDDO
208 ENDDO
209 ENDIF
210c
211 DO j = 1, size_y
212 char(j)= " "
213 ENDDO
214 DO j = 1, size_y
215 DO i = 1, size_x
216 IF( grid(i,size_y+1-j) == 1) char(j)(i:i)= symbol(1:1)
217 ENDDO
218 ENDDO
219
220 DEALLOCATE(grid)
221
222 IF(PRESENT(input_txt_y)) THEN
223 WRITE(iout,'(12X,A)') txt_y(1:len_trim(txt_y))
224 ENDIF
225C
226 WRITE(iout,'(19X,A)') '^'
227 WRITE(iout,'(19X,A,A)') '|',char(1)(1:size_x)
228 WRITE(iout,'(9X,F7.3,1X,A,A,A)') y_maxvalue,'--','|',char(2)(1:size_x)
229 DO i = 3,size_y-1
230 WRITE(iout,'(19X,A,A)') '|',char(i)(1:size_x)
231 ENDDO
232 WRITE(iout,'(9X,F7.3,1X,A,A,A)') y_minvalue,'--','|',char(size_y)(1:size_x)
233C
234 DO ii = 1,size_x+2
235 char2(ii:ii) = '-'
236 ENDDO
237 char2(size_x+2:size_x+3) = '>'
238 IF(PRESENT(input_txt_x) ) THEN
239 char2(size_x+4:size_x+3+len_trim(txt_x)) = txt_x(1:len_trim(txt_x))
240 ENDIF
241C
242 lena = len_trim(char2)
243 WRITE(iout,'(20X,A)') char2(1:lena)
244c
245 char2 = ' '
246 char2(1:1)= '|'
247 char2(size_x:size_x)= '|'
248 lena = len_trim(char2)
249 WRITE(iout,'(20X,A)') char2(1:lena)
250C
251 IF (size_x <= 10) THEN
252 fmta='(14X,1PG11.3,5X,1PG11.3)'
253 ELSEIF (size_x <= 20) THEN
254 fmta='(14X,1PG11.3,8X,1PG11.3)'
255 ELSEIF (size_x <= 30) THEN
256 fmta='(14X,1PG11.3,18X,1PG11.3)'
257 ELSEIF (size_x <= 40) THEN
258 fmta='(14X,1PG11.3,28X,1PG11.3)'
259 ELSEIF (size_x <= 50) THEN
260 fmta='(14X,1PG11.3,38X,1PG11.3)'
261 ELSEIF (size_x <= 60) THEN
262 fmta='(14X,1PG11.3,48X,1PG11.3)'
263 ELSEIF (size_x <= 70) THEN
264 fmta='(14X,1PG11.3,58X,1PG11.3)'
265 ELSEIF (size_x <= 80) THEN
266 fmta='(14X,1PG11.3,68X,1PG11.3)'
267 ELSEIF (size_x <= 90) THEN
268 fmta='(14X,1PG11.3,78X,1PG11.3)'
269 ELSE
270 fmta='(14 X,1PG11.3,85X,1PG11.3)'
271 ENDIF
272c
273 WRITE(iout,fmt=fmta) x_minvalue,x_maxvalue
274c
275 WRITE(iout,*) ' '
276 WRITE(iout,*) ' '
277C------------------------------
278 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::ncharline
integer, parameter ncharline
Definition names_and_titles_mod.F:1006