SFEMaNS  version 5.3
Reference documentation for SFEMaNS
sub_plot.f90
Go to the documentation of this file.
1 !
2 !Authors: Jean-Luc Guermond, Lugi Quartapelle, Copyright 1994
3 !
4 MODULE sub_plot
5 
6 CONTAINS
7 
8  SUBROUTINE plot_vit_2d(jj, rr, uu)
9  !---FORMAT PLTMTV
10 
11  IMPLICIT NONE
12 
13  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
14  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
15  REAL(KIND=8), DIMENSION(:), INTENT(INOUT) :: uu ! rigid translation
16 
17  INTEGER :: m, n1, n2, n3
18 
19  OPEN(unit = 20, file ='vit_u', form ='formatted', status ='unknown')
20 
21  WRITE (20, *) '$ DATA = CONTCURVE'
22  WRITE (20, *) '% contstyle = 2'
23  WRITE (20,*) '% meshplot = False'
24 
25  DO m = 1, SIZE(jj, 2)
26  ! n1 = jj(1, m)
27  ! n2 = jj(2, m)
28  ! n3 = jj(3, m)
29  ! WRITE (20, *) rr(1,n1), rr(2,n1), uu(n1), n1
30  ! WRITE (20, *) rr(1,n2), rr(2,n2), uu(n2), n2
31  ! WRITE (20, *) rr(1,n3), rr(2,n3), uu(n3), n3
32  ! WRITE (20, *)
33 
34  n1 = jj(1, m)
35  n2 = jj(6, m)
36  n3 = jj(5, m)
37  WRITE (20, *) rr(1,n1), rr(2,n1), uu(n1), n1
38  WRITE (20, *) rr(1,n2), rr(2,n2), uu(n2), n2
39  WRITE (20, *) rr(1,n3), rr(2,n3), uu(n3), n3
40  WRITE (20, *)
41 
42 
43  n1 = jj(2, m)
44  n2 = jj(4, m)
45  n3 = jj(6, m)
46  WRITE (20, *) rr(1,n1), rr(2,n1), uu(n1), n1
47  WRITE (20, *) rr(1,n2), rr(2,n2), uu(n2), n2
48  WRITE (20, *) rr(1,n3), rr(2,n3), uu(n3), n3
49  WRITE (20, *)
50 
51  n1 = jj(4, m)
52  n2 = jj(3, m)
53  n3 = jj(5, m)
54  WRITE (20, *) rr(1,n1), rr(2,n1), uu(n1), n1
55  WRITE (20, *) rr(1,n2), rr(2,n2), uu(n2), n2
56  WRITE (20, *) rr(1,n3), rr(2,n3), uu(n3), n3
57  WRITE (20, *)
58 
59  n1 = jj(5, m)
60  n2 = jj(6, m)
61  n3 = jj(4, m)
62  WRITE (20, *) rr(1,n1), rr(2,n1), uu(n1), n1
63  WRITE (20, *) rr(1,n2), rr(2,n2), uu(n2), n2
64  WRITE (20, *) rr(1,n3), rr(2,n3), uu(n3), n3
65  WRITE (20, *)
66  ENDDO
67 
68  END SUBROUTINE plot_vit_2d
69 
70  !----------------------------------------------------------------
71 
72  SUBROUTINE plot_arrow(jj, rr, vv)
73  !---FORMAT PLTMTV
74 
75  IMPLICIT NONE
76 
77  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
78  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr, vv
79 
80  INTEGER :: n
81  REAL(KIND=8) :: z_d, w_d
82 
83  OPEN(unit = 20,file ='velocity',form ='formatted',status ='unknown')
84 
85  WRITE (20, 9000) SIZE(rr,2), SIZE(jj,1), SIZE(jj,2)
86 9000 FORMAT('#',3i10)
87  WRITE (20, *) '$ DATA = VECTOR'
88  WRITE (20, *) '% axisscale = FALSE'
89  WRITE (20, *) '% vscale = 1.'
90 
91  ! WRITE (20, *) '% contstyle = 2'
92  ! WRITE (20, *) '% meshplot = True'
93 
94  SELECT CASE(SIZE(rr,1))
95 
96  CASE(2)
97  z_d = 0
98  w_d = 0
99  DO n = 1, SIZE(rr, 2)
100  WRITE (20,*) rr(1,n), rr(2,n), z_d, vv(1,n), vv(2,n), w_d
101  ENDDO
102 
103  CASE(3)
104  DO n = 1, SIZE(rr, 2)
105  WRITE (20,*) rr(:,n), vv(:,n)
106  ENDDO
107 
108  END SELECT
109 
110  CLOSE(20)
111 
112  END SUBROUTINE plot_arrow
113 
114  !------------------------------------------------------------------------------
115 
116  SUBROUTINE plot_pressure_2d(jj, rr, uu)
117  !---FORMAT PLTMTV
118 
119  IMPLICIT NONE
120 
121  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
122  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
123  REAL(KIND=8), DIMENSION(:), INTENT(INOUT) :: uu ! rigid translation
124 
125  INTEGER :: na, nb, n, m, n1, n2, n3
126  REAL(KIND=8) :: xa, xb, xm, r, p_midpoint
127 
128  OPEN (unit = 20, file = 'pressure_2D.plt', form = 'formatted', status = 'unknown')
129 
130  ! set pressure field so as to have a zero pressure at
131  ! the mid-point of the bottom unit side of the cavity
132 
133  ! find the mid-point
134 
135  xa = 0; xb = 1; xm = 0.5
136  na = 1; nb = 1
137  DO n = 1, SIZE(uu)
138 
139  IF (rr(2,n) < 1.0d-4) THEN
140 
141  r = rr(1,n)
142  IF (r < xm .AND. r > xa) then; xa = r; na = n
143  ENDIF
144  IF (r > xm .AND. r < xb) then; xb = r; nb = n
145  ENDIF
146 
147  ENDIF
148 
149  ENDDO
150 
151  p_midpoint = uu(na) + (uu(nb) - uu(na)) * (xm-xa)/(xb-xa)
152 
153  ! reset pressure field
154 
155  uu = uu - p_midpoint
156 
157  WRITE (20, *) '$ DATA = CONTCURVE'
158  WRITE (20, *) '% contstyle = 1'
159  WRITE (20, *) '% contours = "( -.12 -.06 -.03 )"'
160  WRITE (20, *) '% contours = "( -.015 -.0075 -.00375 )"'
161  WRITE (20, *) '% contours = "( 0.0 )"'
162  WRITE (20, *) '% contours = "( .12 .06 .03 )"'
163  WRITE (20, *) '% contours = "( .015 .0075 .00375 )"'
164 
165  WRITE (20,*) '% meshplot = False'
166 
167  DO m = 1, SIZE(jj, 2)
168  n1 = jj(1, m)
169  n2 = jj(2, m)
170  n3 = jj(3, m)
171  WRITE (20, *) rr(1,n1), rr(2,n1), uu(n1), n1
172  WRITE (20, *) rr(1,n2), rr(2,n2), uu(n2), n2
173  WRITE (20, *) rr(1,n3), rr(2,n3), uu(n3), n3
174  WRITE (20, *)
175  ENDDO
176  CLOSE(20)
177 
178  END SUBROUTINE plot_pressure_2d
179 
180  !------------------------------------------------------------------------------
181 
182  SUBROUTINE plot_loc_rel_var(jj, rr, vo, vv)
183  !---FORMAT PLTMTV
184 
185  ! plot local relative variation of speed
186 
187 
188  IMPLICIT NONE
189 
190  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
191  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr, vo, vv
192 
193  REAL(KIND=8), DIMENSION(SIZE(vv,2)) :: uu
194  REAL(KIND=8) :: vo_n_mod
195  INTEGER :: n, m, n1, n2, n3
196 
197  OPEN (unit=20, file ='speed_var',form ='formatted',status ='unknown')
198 
199 
200  DO n = 1, SIZE(vv,2)
201 
202  vo_n_mod = sqrt(sum(vo(:,n)**2))
203 
204  IF (vo_n_mod < 1.0d-15) THEN
205 
206  uu(n) = 0
207 
208  ELSE
209 
210  uu(n) = sqrt(sum( (vo(:,n) - vv(:,n))**2 ))/vo_n_mod
211 
212  ENDIF
213 
214  ENDDO
215 
216 
217  WRITE (20, *) '$ DATA = CONTCURVE'
218  WRITE (20, *) '% contstyle = 2'
219  WRITE (20, *) '% meshplot = True'
220 
221  DO m = 1, SIZE(jj, 2)
222  n1 = jj(1, m)
223  n2 = jj(2, m)
224  n3 = jj(3, m)
225  WRITE (20, *) rr(1,n1), rr(2,n1), uu(n1), n1
226  WRITE (20, *) rr(1,n2), rr(2,n2), uu(n2), n2
227  WRITE (20, *) rr(1,n3), rr(2,n3), uu(n3), n3
228  WRITE (20, *)
229  ENDDO
230 
231  END SUBROUTINE plot_loc_rel_var
232 
233  !------------------------------------------------------------------------------
234 
235  SUBROUTINE plot_pressure(jj, rr, uu)
236  !---FORMAT PLTMTV
237 
238  IMPLICIT NONE
239 
240  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
241  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
242  REAL(KIND=8), DIMENSION(:), INTENT(INOUT) :: uu ! rigid translation
243  INTEGER :: m, n1, n2, n3
244 
245  OPEN (unit=20, file='pressure', form='formatted', status='unknown')
246 
247  WRITE (20, *) '$ DATA = CONTCURVE'
248  WRITE (20, *) '% contstyle = 2'
249  WRITE (20, *) '% cstep = 20'
250  WRITE (20, *) '% meshplot = False'
251 
252  DO m = 1, SIZE(jj, 2)
253  n1 = jj(1, m)
254  n2 = jj(2, m)
255  n3 = jj(3, m)
256  WRITE (20, *) rr(1,n1), rr(2,n1), uu(n1), n1
257  WRITE (20, *) rr(1,n2), rr(2,n2), uu(n2), n2
258  WRITE (20, *) rr(1,n3), rr(2,n3), uu(n3), n3
259  WRITE (20, *)
260  ENDDO
261 
262  END SUBROUTINE plot_pressure
263 
264  SUBROUTINE plot_const_p1_label(jj, rr, uu, file_name )
265  !---FORMAT PLTMTV
266 
267  IMPLICIT NONE
268 
269  CHARACTER(*) :: file_name
270  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
271  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
272  REAL(KIND=8), DIMENSION(:), INTENT(IN) :: uu
273  INTEGER :: m, n1, n2, n3
274 
275  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
276 
277  WRITE (20, *) '$ DATA = CONTCURVE'
278  WRITE (20, *) '% contstyle = 2'
279  WRITE (20, *) '% nsteps = 50'
280  WRITE (20, *) '% meshplot = true'
281  WRITE (20, *)
282 
283  DO m = 1, SIZE(jj, 2)
284  n1 = jj(1, m)
285  n2 = jj(2, m)
286  n3 = jj(3, m)
287  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(m), n1
288  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(m), n2
289  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(m), n3
290  WRITE (20, *)
291  ENDDO
292 100 FORMAT(3(e12.5,3x),i5)
293  CLOSE(20)
294 
295  END SUBROUTINE plot_const_p1_label
296 
297  SUBROUTINE plot_pressure_p1_label(jj, rr, uu, file_name )
298  !---FORMAT PLTMTV
299 
300  IMPLICIT NONE
301 
302  CHARACTER(*) :: file_name
303  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
304  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
305  REAL(KIND=8), DIMENSION(:), INTENT(INOUT) :: uu ! rigid translation
306  INTEGER :: m, n1, n2, n3
307 
308  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
309 
310  WRITE (20, *) '$ DATA = CONTCURVE'
311  WRITE (20, *) '% contstyle = 2'
312  WRITE (20, *) '% nsteps = 50'
313  WRITE (20, *) '% meshplot = false'
314  WRITE (20, *)
315 
316  DO m = 1, SIZE(jj, 2)
317 
318  n1 = jj(1, m)
319  n2 = jj(2, m)
320  n3 = jj(3, m)
321  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
322  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
323  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
324  WRITE (20, *)
325  ENDDO
326 100 FORMAT(3(e12.5,3x),i5)
327  CLOSE(20)
328 
329  END SUBROUTINE plot_pressure_p1_label
330 
331  SUBROUTINE plot_p1_cont_label(jj, jjs, sides, list, rr, uu, file_name )
332  !---FORMAT PLOTMTV
333 
334  IMPLICIT NONE
335 
336  CHARACTER(*) :: file_name
337  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj, jjs
338  INTEGER, DIMENSION(:), INTENT(IN) :: sides, list
339  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
340  REAL(KIND=8), DIMENSION(:), INTENT(INOUT) :: uu ! rigid translation
341 
342  INTEGER :: m, n1, n2, n3, ms
343 
344  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
345 
346  WRITE (20, *) '$ DATA = CONTCURVE'
347  WRITE (20, *) '% contstyle = 2'
348  WRITE (20, *) '% nsteps = 50'
349  WRITE (20, *) '% meshplot = f'
350  WRITE (20, *)
351 
352  DO ms = 1, SIZE(sides)
353  IF(minval(abs(list-sides(ms))) /= 0 ) cycle
354  n1 = jjs(1,ms)
355  n2 = jjs(2,ms)
356  WRITE(20,110) '@line x1=',rr(1,n1), 'y1=',rr(2,n1), 'z1=',0., &
357  'x2=',rr(1,n2), 'y2=',rr(2,n2), 'z2=',0.
358  END DO
359 110 FORMAT(6(a,x,e12.5,3x))
360 
361  WRITE (20, *)
362  DO m = 1, SIZE(jj, 2)
363  n1 = jj(1, m)
364  n2 = jj(2, m)
365  n3 = jj(3, m)
366  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
367  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
368  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
369  WRITE (20, *)
370  ENDDO
371 100 FORMAT(3(e12.5,3x),i5)
372  CLOSE(20)
373 
374  END SUBROUTINE plot_p1_cont_label
375 
376  SUBROUTINE plot_p1_matiere_label(jj, neigh, i_d, rr, uu, file_name )
377  !---FORMAT PLOTMTV
378 
379  IMPLICIT NONE
380 
381  CHARACTER(*) :: file_name
382  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj, neigh
383  INTEGER, DIMENSION(:), INTENT(IN) :: i_d
384  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
385  REAL(KIND=8), DIMENSION(:), INTENT(INOUT) :: uu ! rigid translation
386 
387  INTEGER :: n, m, n1, n2, n3, nghm
388 
389  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
390 
391  WRITE (20, *) '$ DATA = CONTCURVE'
392  WRITE (20, *) '% contstyle = 2'
393  WRITE (20, *) '% nsteps = 50'
394  WRITE (20, *) '% meshplot = f'
395  WRITE (20, *)
396 
397  DO m = 1, SIZE(neigh,2)
398  DO n = 1, 3
399  nghm = neigh(n,m)
400 
401  IF (nghm == 0) THEN
402  CONTINUE
403  ELSE IF (i_d(nghm) /= i_d(m)) THEN
404  CONTINUE
405  ELSE
406  cycle
407  END IF
408 
409  n1 = jj(modulo(n,3) + 1,m)
410  n2 = jj(modulo(n+1,3) + 1,m)
411  WRITE(20,110) '@line x1=',rr(1,n1), 'y1=',rr(2,n1), 'z1=',0., &
412  'x2=',rr(1,n2), 'y2=',rr(2,n2), 'z2=',0.
413  END DO
414  END DO
415 110 FORMAT(6(a,x,e12.5,3x))
416 
417  WRITE (20, *)
418  DO m = 1, SIZE(jj, 2)
419  n1 = jj(1, m)
420  n2 = jj(2, m)
421  n3 = jj(3, m)
422  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
423  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
424  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
425  WRITE (20, *)
426  ENDDO
427 100 FORMAT(3(e12.5,3x),i5)
428  CLOSE(20)
429 
430  END SUBROUTINE plot_p1_matiere_label
431 
432  SUBROUTINE plot_pressure_label(jj, rr, uu, file_name )
433  !---FORMAT PLTMTV
434 
435  IMPLICIT NONE
436 
437  CHARACTER(*) :: file_name
438  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
439  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
440  REAL(KIND=8), DIMENSION(:), INTENT(INOUT) :: uu ! rigid translation
441 
442  INTEGER :: m, n1, n2, n3
443 
444  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
445 
446  WRITE (20, *) '$ DATA = CONTCURVE'
447  WRITE (20, *) '% contstyle = 2'
448  WRITE (20, *) '% nsteps = 50'
449  WRITE (20, *) '% meshplot = False'
450  WRITE (20, *)
451 
452  DO m = 1, SIZE(jj, 2)
453  n1 = jj(1, m)
454  n2 = jj(2, m)
455  n3 = jj(3, m)
456  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
457  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
458  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
459  WRITE (20, *)
460  ENDDO
461 100 FORMAT(3(e12.5,3x),i5)
462 
463  CLOSE(20)
464 
465  END SUBROUTINE plot_pressure_label
466 
467  SUBROUTINE plot_arrow_label(jj, rr, vv, file_name)
468  !---FORMAT PLTMTV
469 
470  IMPLICIT NONE
471 
472  CHARACTER(*) :: file_name
473  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
474  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr, vv
475 
476  INTEGER :: n
477  REAL(KIND=8) :: z_d, w_d
478 
479  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
480 
481  WRITE (20, 9000) SIZE(rr,2), SIZE(jj,1), SIZE(jj,2)
482 9000 FORMAT('#',3i10)
483  WRITE (20, *) '$ DATA = VECTOR'
484  WRITE (20, *) '% axisscale = FALSE'
485  WRITE (20, *) '% vscale = 1.'
486 
487  ! WRITE (20, *) '% contstyle = 2'
488  ! WRITE (20, *) '% meshplot = True'
489 
490  SELECT CASE(SIZE(rr,1))
491 
492  CASE(2)
493  z_d = 0
494  w_d = 0
495  IF (SIZE(vv,1)/=2) THEN
496  DO n = 1, SIZE(rr, 2)
497  WRITE (20,100) rr(1,n), rr(2,n), z_d, vv(n,1), vv(n,2), w_d
498  ENDDO
499  ELSE
500  DO n = 1, SIZE(rr, 2)
501  WRITE (20,100) rr(1,n), rr(2,n), z_d, vv(1,n), vv(2,n), w_d
502  ENDDO
503  END IF
504  CASE(3)
505  IF (SIZE(vv,1)/=3) THEN
506  DO n = 1, SIZE(rr, 2)
507  WRITE (20,100) rr(:,n), vv(:,n)
508  ENDDO
509  ELSE
510  DO n = 1, SIZE(rr, 2)
511  WRITE (20,100) rr(:,n), vv(:,n)
512  ENDDO
513  END IF
514  END SELECT
515 
516  CLOSE(20)
517 100 FORMAT(6(e12.5,3x))
518 
519  END SUBROUTINE plot_arrow_label
520 
521  SUBROUTINE plot_pressure_p2_label(jj, rr, uu, file_name)
522  !---FORMAT PLTMTV
523 
524  IMPLICIT NONE
525 
526  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
527  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
528  REAL(KIND=8), DIMENSION(:), INTENT(INOUT) :: uu ! rigid translation
529  CHARACTER(*) :: file_name
530 
531  INTEGER :: m, n1, n2, n3
532 
533  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
534 
535  WRITE (20, *) '$ DATA = CONTCURVE'
536  WRITE (20, *) '% contstyle = 2'
537  WRITE (20, *) '% nsteps = 50'
538  WRITE (20,*) '% meshplot = false'
539 
540  DO m = 1, SIZE(jj, 2)
541 
542  n1 = jj(1, m)
543  n2 = jj(6, m)
544  n3 = jj(5, m)
545  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
546  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
547  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
548  WRITE (20, 100)
549 
550 
551  n1 = jj(2, m)
552  n2 = jj(4, m)
553  n3 = jj(6, m)
554  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
555  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
556  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
557  WRITE (20, 100)
558 
559  n1 = jj(4, m)
560  n2 = jj(3, m)
561  n3 = jj(5, m)
562  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
563  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
564  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
565  WRITE (20, 100)
566 
567  n1 = jj(5, m)
568  n2 = jj(6, m)
569  n3 = jj(4, m)
570  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
571  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
572  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
573  WRITE (20, 100)
574  ENDDO
575 100 FORMAT(3(e12.5,3x),i5)
576 
577  CLOSE(20)
578 
579  END SUBROUTINE plot_pressure_p2_label
580 
581  !----------------------------------------------------------------
582 
583  SUBROUTINE plot_ensight_vecteur(u8, vit)
584  !--FORMAT ENSIGHT
585 
586  IMPLICIT NONE
587 
588  CHARACTER(LEN=*), INTENT(IN) :: vit
589  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: u8
590  CHARACTER(LEN=80) :: sketuve
591  INTEGER :: i, j
592 
593  OPEN(unit=33,file=vit,status='unknown',form='UNFORMATTED')
594 
595  sketuve = 'rien'
596  WRITE(33) sketuve
597 
598  WRITE(33) ((REAL(u8(i,j),KIND=4),i=1,size(u8,1)),j=1,SIZE(u8,2))
599 
600  CLOSE(33)
601 
602  END SUBROUTINE plot_ensight_vecteur
603 
604  !----------------------------------------------------------------
605 
606  SUBROUTINE plot_ensight_scalaire(p8, pres)
607  !--FORMAT ENSIGHT
608 
609  IMPLICIT NONE
610 
611  CHARACTER(LEN=*), INTENT(IN) :: pres
612  REAL(KIND=8), DIMENSION(:), INTENT(IN) :: p8
613  CHARACTER(LEN=80) :: sketuve
614  INTEGER :: i
615 
616  OPEN(unit=44,file=pres,status='unknown',form='UNFORMATTED')
617 
618  sketuve = 'rien'
619  WRITE(44) sketuve
620 
621  WRITE(44) (REAL(p8(i),KIND =4),i=1,size(p8))
622 
623  CLOSE(44)
624 
625  END SUBROUTINE plot_ensight_scalaire
626 
627  SUBROUTINE plot_vorticity (jj, rr, zz, i)
628 
629  IMPLICIT NONE
630 
631  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
632  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
633  REAL(KIND=8), DIMENSION(:), INTENT(IN) :: zz
634  INTEGER, INTENT(IN) :: i
635 
636  INTEGER :: m, n1, n2, n3
637 
638  ! OPEN (UNIT = 20, FILE = 'vorticity', FORM = 'formatted', STATUS = 'unknown')
639 
640  WRITE (i, *) '$ DATA = CONTCURVE'
641  WRITE (i, *) '% contstyle = 1'
642  WRITE (i, *) '% nsteps = 40'
643 
644  ! WRITE (i, *) '% contours = "(-6.0 -5.0 -4.0 -3.0 -2.0 -1.0)"'
645  ! WRITE (i, *) '% contours = "(-0.5 0.0 0.5)"'
646  ! WRITE (i, *) '% contours = "( 1.0 2.0 3.0 4.0 5.0 6.0)"'
647 
648  WRITE (i, *) '% meshplot = False'
649 
650  DO m = 1, SIZE(jj, 2)
651  n1 = jj(1, m)
652  n2 = jj(2, m)
653  n3 = jj(3, m)
654  WRITE (i, *) rr(1,n1), rr(2,n1), zz(n1), n1
655  WRITE (i, *) rr(1,n2), rr(2,n2), zz(n2), n2
656  WRITE (i, *) rr(1,n3), rr(2,n3), zz(n3), n3
657  WRITE (i, *)
658  ENDDO
659 
660  CLOSE (unit = i)
661 
662  END SUBROUTINE plot_vorticity
663 
664  SUBROUTINE plot_stream (jj, rr, pp, i)
665 
666  IMPLICIT NONE
667 
668  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
669  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
670  REAL(KIND=8), DIMENSION(:), INTENT(IN) :: pp
671  INTEGER, INTENT(IN) :: i
672 
673  INTEGER :: m, n1, n2, n3
674 
675  ! OPEN (UNIT = 20, FILE = 'stream', FORM = 'formatted', STATUS = 'unknown')
676 
677  WRITE (i, *) '$ DATA = CONTCURVE'
678  WRITE (i, *) '% contstyle = 1'
679  WRITE (i, *) '% nsteps = 30'
680 
681  ! WRITE (i, *) '% contours = "(-1.e-10 -1.e-7 -1.e-5 -1.e-4 )"'
682  ! WRITE (i, *) '% contours = "(-0.01 -0.03 -0.05 -0.07 -0.09 )"'
683  ! WRITE (i, *) '% contours = "(-0.100 -0.110 -0.115 -0.1175 )"'
684  ! WRITE (i, *) '% contours = "(+1.0e-8 +1.0e-7 +1.0e-6 +1.0e-5 )"'
685  ! WRITE (i, *) '% contours = "(+5.0e-5 +1.0e-4 +2.5e-4 +5.0e-4 )"'
686  ! WRITE (i, *) '% contours = "(+1.0e-3 +1.5e-3 +3.0e-3 )"'
687  WRITE (i, *) '% meshplot = False'
688 
689  DO m = 1, SIZE(jj, 2)
690  n1 = jj(1, m)
691  n2 = jj(2, m)
692  n3 = jj(3, m)
693  WRITE (i, *) rr(1,n1), rr(2,n1), pp(n1), n1
694  WRITE (i, *) rr(1,n2), rr(2,n2), pp(n2), n2
695  WRITE (i, *) rr(1,n3), rr(2,n3), pp(n3), n3
696  WRITE (i, *)
697  ENDDO
698 
699  CLOSE (unit = i)
700 
701  END SUBROUTINE plot_stream
702 
703  SUBROUTINE plot_scalar_field(jj, rr, uu, file_name)
704  !---FORMAT PLTMTV
705 
706  IMPLICIT NONE
707 
708  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
709  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
710  REAL(KIND=8), DIMENSION(:), INTENT(IN) :: uu ! rigid translation
711  CHARACTER(*) :: file_name
712 
713  INTEGER :: m, n1, n2, n3, unit_w=47
714 
715 100 FORMAT(3(e15.8,3x),i5)
716  OPEN (unit=unit_w, file=file_name, form='formatted', status='unknown')
717 
718  WRITE (unit_w, '(A)') '$ DATA = CONTCURVE'
719  WRITE (unit_w, '(A)') '% contstyle = 2'
720  WRITE (unit_w, '(A)') '% nsteps = 50'
721  WRITE (unit_w, '(A)') '% meshplot = false'
722  WRITE (unit_w, '(A)')
723 
724  IF (SIZE(jj,1)==3) THEN
725  DO m = 1, SIZE(jj, 2)
726  n1 = jj(1, m)
727  n2 = jj(2, m)
728  n3 = jj(3, m)
729  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
730  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
731  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
732  WRITE (unit_w, 100)
733  ENDDO
734  ELSE IF (SIZE(jj,1)==6) THEN
735  DO m = 1, SIZE(jj, 2)
736  n1 = jj(1, m)
737  n2 = jj(6, m)
738  n3 = jj(5, m)
739  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
740  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
741  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
742  WRITE (unit_w, 100)
743 
744  n1 = jj(2, m)
745  n2 = jj(4, m)
746  n3 = jj(6, m)
747  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
748  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
749  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
750  WRITE (unit_w, 100)
751 
752  n1 = jj(4, m)
753  n2 = jj(3, m)
754  n3 = jj(5, m)
755  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
756  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
757  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
758  WRITE (unit_w, 100)
759 
760  n1 = jj(5, m)
761  n2 = jj(6, m)
762  n3 = jj(4, m)
763  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
764  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
765  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
766  WRITE (unit_w, 100)
767  ENDDO
768 
769  ELSE IF (SIZE(jj,1)==10) THEN
770 
771  DO m = 1, SIZE(jj, 2)
772  n1 = jj(1, m)
773  n2 = jj(8, m)
774  n3 = jj(6, m)
775  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
776  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
777  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
778  WRITE (unit_w, 100)
779 
780  n1 = jj(8, m)
781  n2 = jj(10, m)
782  n3 = jj(6, m)
783  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
784  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
785  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
786  WRITE (unit_w, 100)
787 
788  n1 = jj(8, m)
789  n2 = jj(9, m)
790  n3 = jj(10, m)
791  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
792  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
793  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
794  WRITE (unit_w, 100)
795 
796  n1 = jj(9, m)
797  n2 = jj(4, m)
798  n3 = jj(10, m)
799  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
800  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
801  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
802  WRITE (unit_w, 100)
803 
804  n1 = jj(9, m)
805  n2 = jj(2, m)
806  n3 = jj(4, m)
807  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
808  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
809  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
810  WRITE (unit_w, 100)
811 
812  n1 = jj(6, m)
813  n2 = jj(10, m)
814  n3 = jj(7, m)
815  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
816  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
817  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
818  WRITE (unit_w, 100)
819 
820  n1 = jj(10, m)
821  n2 = jj(5, m)
822  n3 = jj(7, m)
823  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
824  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
825  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
826  WRITE (unit_w, 100)
827 
828  n1 = jj(10, m)
829  n2 = jj(4, m)
830  n3 = jj(5, m)
831  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
832  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
833  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
834  WRITE (unit_w, 100)
835 
836  n1 = jj(7, m)
837  n2 = jj(5, m)
838  n3 = jj(3, m)
839  WRITE (unit_w, 100) rr(1,n1), rr(2,n1), uu(n1), n1
840  WRITE (unit_w, 100) rr(1,n2), rr(2,n2), uu(n2), n2
841  WRITE (unit_w, 100) rr(1,n3), rr(2,n3), uu(n3), n3
842  WRITE (unit_w, 100)
843  ENDDO
844  ELSE
845  WRITE(*,*) ' Problem in plot_scalar_field ', SIZE(jj,1)
846  stop
847  END IF
848 
849  CLOSE(unit=unit_w)
850 
851  END SUBROUTINE plot_scalar_field
852 
853  SUBROUTINE plot_two_scalar_field(jj, rr, uu, jj2, rr2, uu2, file_name)
854  !---FORMAT PLTMTV
855 
856  IMPLICIT NONE
857 
858  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj, jj2
859  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr, rr2
860  REAL(KIND=8), DIMENSION(:) :: uu, uu2
861  CHARACTER(*) :: file_name
862 
863  INTEGER :: m, n1, n2, n3
864 
865  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
866 
867  WRITE (20, *) '$ DATA = CONTCURVE'
868  WRITE (20, *) '% contstyle = 2'
869  WRITE (20, *) '% nsteps = 50'
870  WRITE (20, *) '% meshplot = false'
871  WRITE (20, *)
872 
873  IF (SIZE(jj,1)==3) THEN
874  DO m = 1, SIZE(jj, 2)
875 
876  n1 = jj(1, m)
877  n2 = jj(2, m)
878  n3 = jj(3, m)
879  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
880  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
881  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
882  WRITE (20, *)
883  ENDDO
884  DO m = 1, SIZE(jj2, 2)
885 
886  n1 = jj2(1, m)
887  n2 = jj2(2, m)
888  n3 = jj2(3, m)
889  WRITE (20, 100) rr2(1,n1), rr2(2,n1), uu2(n1), n1
890  WRITE (20, 100) rr2(1,n2), rr2(2,n2), uu2(n2), n2
891  WRITE (20, 100) rr2(1,n3), rr2(2,n3), uu2(n3), n3
892  WRITE (20, *)
893  ENDDO
894 
895 
896  ELSE IF (SIZE(jj,1)==6) THEN
897 
898  DO m = 1, SIZE(jj, 2)
899 
900  n1 = jj(1, m)
901  n2 = jj(6, m)
902  n3 = jj(5, m)
903  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
904  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
905  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
906  WRITE (20, 100)
907 
908 
909  n1 = jj(2, m)
910  n2 = jj(4, m)
911  n3 = jj(6, m)
912  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
913  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
914  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
915  WRITE (20, 100)
916 
917  n1 = jj(4, m)
918  n2 = jj(3, m)
919  n3 = jj(5, m)
920  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
921  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
922  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
923  WRITE (20, 100)
924 
925  n1 = jj(5, m)
926  n2 = jj(6, m)
927  n3 = jj(4, m)
928  WRITE (20, 100) rr(1,n1), rr(2,n1), uu(n1), n1
929  WRITE (20, 100) rr(1,n2), rr(2,n2), uu(n2), n2
930  WRITE (20, 100) rr(1,n3), rr(2,n3), uu(n3), n3
931  WRITE (20, 100)
932  ENDDO
933  ELSE IF (SIZE(jj,1)==10) THEN
934  DO m = 1, SIZE(jj2, 2)
935 
936  n1 = jj2(1, m)
937  n2 = jj2(6, m)
938  n3 = jj2(5, m)
939  WRITE (20, 100) rr2(1,n1), rr2(2,n1), uu2(n1), n1
940  WRITE (20, 100) rr2(1,n2), rr2(2,n2), uu2(n2), n2
941  WRITE (20, 100) rr2(1,n3), rr2(2,n3), uu2(n3), n3
942  WRITE (20, 100)
943 
944 
945  n1 = jj2(2, m)
946  n2 = jj2(4, m)
947  n3 = jj2(6, m)
948  WRITE (20, 100) rr2(1,n1), rr2(2,n1), uu2(n1), n1
949  WRITE (20, 100) rr2(1,n2), rr2(2,n2), uu2(n2), n2
950  WRITE (20, 100) rr2(1,n3), rr2(2,n3), uu2(n3), n3
951  WRITE (20, 100)
952 
953  n1 = jj2(4, m)
954  n2 = jj2(3, m)
955  n3 = jj2(5, m)
956  WRITE (20, 100) rr2(1,n1), rr2(2,n1), uu2(n1), n1
957  WRITE (20, 100) rr2(1,n2), rr2(2,n2), uu2(n2), n2
958  WRITE (20, 100) rr2(1,n3), rr2(2,n3), uu2(n3), n3
959  WRITE (20, 100)
960 
961  n1 = jj2(5, m)
962  n2 = jj2(6, m)
963  n3 = jj2(4, m)
964  WRITE (20, 100) rr2(1,n1), rr2(2,n1), uu2(n1), n1
965  WRITE (20, 100) rr2(1,n2), rr2(2,n2), uu2(n2), n2
966  WRITE (20, 100) rr2(1,n3), rr2(2,n3), uu2(n3), n3
967  WRITE (20, 100)
968  ENDDO
969  ELSE
970  WRITE(*,*) ' Problem in plot_two_scalar_field ', SIZE(jj,1)
971  stop
972  END IF
973 
974 100 FORMAT(3(e12.5,3x),i5)
975 
976  CLOSE(20)
977 
978  END SUBROUTINE plot_two_scalar_field
979 
980  SUBROUTINE plot_scalar_field_domain(jj, rr, id, index_dom, file_name)
981  !---FORMAT PLTMTV
982 
983  IMPLICIT NONE
984 
985  INTEGER, DIMENSION(:,:), INTENT(IN) :: jj
986  REAL(KIND=8), DIMENSION(:,:), INTENT(IN) :: rr
987  INTEGER, DIMENSION(:), INTENT(IN) :: id
988  INTEGER, INTENT(IN) :: index_dom
989  CHARACTER(*) :: file_name
990 
991  INTEGER :: m, n1, n2, n3
992 
993  OPEN (unit=20, file=file_name, form='formatted', status='unknown')
994 
995  WRITE (20, *) '$ DATA = CONTCURVE'
996  WRITE (20, *) '% contstyle = 2'
997  WRITE (20, *) '% nsteps = 50'
998  WRITE (20, *) '% meshplot = true'
999  WRITE (20, *)
1000 
1001  IF (SIZE(jj,1)==3) THEN
1002  DO m = 1, SIZE(jj, 2)
1003  IF (id(m)/=index_dom) cycle
1004  n1 = jj(1, m)
1005  n2 = jj(2, m)
1006  n3 = jj(3, m)
1007  WRITE (20, 100) rr(1,n1), rr(2,n1), id(m), n1
1008  WRITE (20, 100) rr(1,n2), rr(2,n2), id(m), n2
1009  WRITE (20, 100) rr(1,n3), rr(2,n3), id(m), n3
1010  WRITE (20, *)
1011  ENDDO
1012 
1013  ELSE IF (SIZE(jj,1)==6) THEN
1014 
1015  DO m = 1, SIZE(jj, 2)
1016  IF (id(m)/=index_dom) cycle
1017  n1 = jj(1, m)
1018  n2 = jj(6, m)
1019  n3 = jj(5, m)
1020  WRITE (20, 100) rr(1,n1), rr(2,n1), id(m), n1
1021  WRITE (20, 100) rr(1,n2), rr(2,n2), id(m), n2
1022  WRITE (20, 100) rr(1,n3), rr(2,n3), id(m), n3
1023  WRITE (20, 100)
1024 
1025  n1 = jj(2, m)
1026  n2 = jj(4, m)
1027  n3 = jj(6, m)
1028  WRITE (20, 100) rr(1,n1), rr(2,n1), id(m), n1
1029  WRITE (20, 100) rr(1,n2), rr(2,n2), id(m), n2
1030  WRITE (20, 100) rr(1,n3), rr(2,n3), id(m), n3
1031  WRITE (20, 100)
1032 
1033  n1 = jj(4, m)
1034  n2 = jj(3, m)
1035  n3 = jj(5, m)
1036  WRITE (20, 100) rr(1,n1), rr(2,n1), id(m), n1
1037  WRITE (20, 100) rr(1,n2), rr(2,n2), id(m), n2
1038  WRITE (20, 100) rr(1,n3), rr(2,n3), id(m), n3
1039  WRITE (20, 100)
1040 
1041  n1 = jj(5, m)
1042  n2 = jj(6, m)
1043  n3 = jj(4, m)
1044  WRITE (20, 100) rr(1,n1), rr(2,n1), id(m), n1
1045  WRITE (20, 100) rr(1,n2), rr(2,n2), id(m), n2
1046  WRITE (20, 100) rr(1,n3), rr(2,n3), id(m), n3
1047  WRITE (20, 100)
1048  ENDDO
1049 
1050  ELSE
1051  WRITE(*,*) ' Problem in plot_scalar_field_domain ', SIZE(jj,1)
1052  stop
1053  END IF
1054 
1055 100 FORMAT(2(e12.5,3x),i5,3x,i5)
1056 
1057  CLOSE(20)
1058 
1059  END SUBROUTINE plot_scalar_field_domain
1060 
1061  SUBROUTINE vtk_p1_2d(mesh, champ, unit_file, file_name)
1062 
1063  USE def_type_mesh
1064 
1065  IMPLICIT NONE
1066 
1067  TYPE(mesh_type), TARGET :: mesh
1068  REAL(KIND=8), DIMENSION(:), INTENT(IN) :: champ
1069  INTEGER, INTENT(IN) :: unit_file
1070  CHARACTER(*), INTENT(IN) :: file_name
1071  INTEGER :: i
1072 
1073  write(*,*) 'ecriture header'
1074  OPEN (unit=unit_file,file=file_name,form = 'formatted',&
1075  status = 'unknown')
1076 
1077  write(unit_file,'(A)') '# vtk DataFile Version 3.0'
1078  write(unit_file,'(A)') 'vtk '//file_name//''
1079  write(unit_file,'(A)')'ASCII'
1080  write(unit_file,'(A)')'DATASET UNSTRUCTURED_GRID'
1081  write(unit_file,'(A,I7,A)')'POINTS ', mesh%np, ' float'
1082  write(*,*) 'points ...'
1083  DO i=1, mesh%np
1084  write(unit_file,'(2(e14.7,2x),A)') mesh%rr(1,i), &
1085  mesh%rr(2,i), ' 0.0 '
1086  ENDDO
1087  write(*,*) 'cells ...'
1088  write(unit_file,'(A,I7,I8)') 'CELLS ', mesh%me, mesh%me*4
1089  DO i=1, mesh%me
1090  write(unit_file,'(A,3(I8,1x))') '3 ', mesh%jj(1,i)-1, &
1091  mesh%jj(2,i)-1, mesh%jj(3,i)-1
1092  ENDDO
1093  write(unit_file,'(A,I7)') 'CELL_TYPES ', mesh%me
1094  DO i=1, mesh%me
1095  write(unit_file,'(A)') '5'
1096  ENDDO
1097 
1098  write(*,*) 'data ...'
1099  write(unit_file,'(A,I7)') 'POINT_DATA ',mesh%np
1100  write(unit_file,'(A)') 'SCALARS scalars float 1'
1101  write(unit_file,'(A)') 'LOOKUP_TABLE default'
1102  DO i=1, mesh%np
1103  write(unit_file,'(e14.7,2x)') champ(i)
1104  ENDDO
1105  CLOSE(unit_file)
1106  END SUBROUTINE vtk_p1_2d
1107 
1108  SUBROUTINE vtk_2d(mesh, champ, unit_file, file_name)
1109 
1110  USE def_type_mesh
1111 
1112  IMPLICIT NONE
1113 
1114  TYPE(mesh_type), TARGET :: mesh
1115  REAL(KIND=8), DIMENSION(:), INTENT(IN) :: champ
1116  INTEGER, INTENT(IN) :: unit_file
1117  CHARACTER(*), INTENT(IN) :: file_name
1118  INTEGER :: i, vtk_cell, n, m
1119 
1120 
1121  write(*,*) 'ecriture header'
1122  OPEN (unit=unit_file,file=file_name,form = 'formatted',&
1123  status = 'unknown')
1124 
1125  write(unit_file,'(A)') '# vtk DataFile Version 3.0'
1126  write(unit_file,'(A)') 'vtk '//file_name//''
1127  write(unit_file,'(A)')'ASCII'
1128  write(unit_file,'(A)')'DATASET UNSTRUCTURED_GRID'
1129  write(unit_file,'(A,I7,A)')'POINTS ', mesh%np, ' float'
1130  write(*,*) 'points ...'
1131  DO i=1, mesh%np
1132  write(unit_file,'(2(e14.7,2x),A)') mesh%rr(1,i), &
1133  mesh%rr(2,i), ' 0.0 '
1134  ENDDO
1135  write(*,*) 'cells ...'
1136 
1137  write(unit_file,'(A,I7,I8)') 'CELLS ', mesh%me, mesh%me*(mesh%gauss%n_w+1)
1138  IF (mesh%gauss%n_w==3) THEN
1139  vtk_cell = 5
1140  DO m=1, mesh%me
1141  write(unit_file,'(I2,6(I8,1x))') mesh%gauss%n_w, (mesh%jj(n,m)-1, n=1,mesh%gauss%n_w)
1142  ENDDO
1143  ELSE IF (mesh%gauss%n_w==6) THEN
1144  vtk_cell = 22
1145  DO m=1, mesh%me
1146  write(unit_file,'(I2,6(I8,1x))') mesh%gauss%n_w, &
1147  mesh%jj(1,m)-1, mesh%jj(2,m)-1, mesh%jj(3,m)-1, mesh%jj(6,m)-1, mesh%jj(4,m)-1, mesh%jj(5,m)-1
1148  ENDDO
1149  ELSE
1150  WRITE(*,*) ' BUG in vtk_2d '
1151  END IF
1152  write(unit_file,'(A,I7)') 'CELL_TYPES ', mesh%me
1153  DO i=1, mesh%me
1154  write(unit_file,'(i2)') vtk_cell
1155  ENDDO
1156 
1157  write(*,*) 'data ...'
1158  write(unit_file,'(A,I7)') 'POINT_DATA ',mesh%np
1159  write(unit_file,'(A)') 'SCALARS scalars float 1'
1160  write(unit_file,'(A)') 'LOOKUP_TABLE default'
1161  DO i=1, mesh%np
1162  write(unit_file,'(e14.7,2x)') champ(i)
1163  ENDDO
1164  CLOSE(unit_file)
1165 
1166 
1167  write(*,*) 'ecriture header'
1168  OPEN (unit=unit_file,file=file_name,form = 'formatted',&
1169  status = 'unknown')
1170 
1171 !!$ write(unit_file,'(A)') '# vtk DataFile Version 3.0'
1172 !!$ write(unit_file,'(A)') 'vtk '//file_name//''
1173 !!$ write(unit_file,'(A)')'ASCII'
1174 !!$ write(unit_file,'(A)')'DATASET UNSTRUCTURED_GRID'
1175 !!$ write(unit_file,'(A,I7,A)')'POINTS ', mesh%np, ' float'
1176 !!$ write(*,*) 'points ...'
1177 !!$ DO i = 1, mesh%np
1178 !!$ write(unit_file,'(2(e14.7,2x),A)') mesh%rr(1,i), &
1179 !!$ mesh%rr(2,i), ' 0.0 '
1180 !!$ ENDDO
1181 !!$ write(*,*) 'cells ...'
1182 !!$ write(unit_file,'(A,I7,I8)') 'CELLS ', mesh%me, mesh%me*(mesh%gauss%n_w+1)
1183 !!$ DO m = 1, mesh%me
1184 !!$ write(unit_file,'(A,3(I8,1x))') mesh%gauss%n_w, (mesh%jj(n,m)-1, n=1,mesh%gauss%n_w)
1185 !!$ ENDDO
1186 !!$ IF (mesh%gauss%n_w==3) THEN
1187 !!$ vtK_cell = 5
1188 !!$ ELSE IF (mesh%gauss%n_w==6) THEN
1189 !!$ vtK_cell = 22
1190 !!$ ELSE
1191 !!$ WRITE(*,*) ' BUG in vtk_2d '
1192 !!$ END IF
1193 !!$ write(unit_file,'(A,I7)') 'CELL_TYPES ', mesh%me
1194 !!$ DO m = 1, mesh%me
1195 !!$ write(unit_file,'(I2)') vtk_cell
1196 !!$ ENDDO
1197 !!$
1198 !!$ write(*,*) 'data ...'
1199 !!$ write(unit_file,'(A,I7)') 'POINT_DATA ',mesh%np
1200 !!$ write(unit_file,'(A)') 'SCALARS scalars float 1'
1201 !!$ write(unit_file,'(A)') 'LOOKUP_TABLE default'
1202 !!$ DO i=1, mesh%np
1203 !!$ write(unit_file,'(e14.7,2x)') champ(i)
1204 !!$ ENDDO
1205 !!$ CLOSE(unit_file)
1206  END SUBROUTINE vtk_2d
1207 
1208  SUBROUTINE trace_profile(mesh, v, it, freq_plot, list_mode, nom_champ, num_dom)
1209 
1210  USE chaine_caractere
1211  USE def_type_mesh
1212 
1213  IMPLICIT NONE
1214  TYPE(mesh_type) :: mesh
1215  REAL(KIND=8) , DIMENSION(:,:,:), INTENT(INOUT):: v
1216  INTEGER, INTENT(IN) :: it, freq_plot
1217  INTEGER, DIMENSION(:), INTENT(IN) :: list_mode
1218  CHARACTER(len=3), INTENT(IN) :: nom_champ
1219  INTEGER, OPTIONAL, INTENT(IN) :: num_dom
1220 
1221  INTEGER :: l, lblank, rang_S
1222  CHARACTER(len=3) :: tit, tmode, tit_S
1223  CHARACTER(len=1) :: tcomp
1224  INTEGER :: i, k
1225 
1226  IF (PRESENT(num_dom)) THEN
1227  rang_s = num_dom
1228  ELSE
1229  rang_s = 0
1230  END IF
1231 
1232  WRITE(tit_s,'(i3)') rang_s
1233  lblank = eval_blank(3,tit_s)
1234  DO l = 1, lblank - 1
1235  tit_s(l:l) = '0'
1236  END DO
1237 
1238  WRITE(tit,'(i3)') it/freq_plot
1239  lblank = eval_blank(3,tit)
1240  DO l = 1, lblank - 1
1241  tit(l:l) = '0'
1242  END DO
1243 
1244  DO i= 1, SIZE(v,3)
1245 
1246  WRITE(tmode,'(i3)') list_mode(i)
1247  lblank = eval_blank(3,tmode)
1248  DO l = 1, lblank - 1
1249  tmode(l:l) = '0'
1250  END DO
1251 
1252  DO k= 1, size(v,2)
1253  WRITE(tcomp,'(i1)') k
1254  CALL plot_scalar_field(mesh%jj, mesh%rr, &
1255  v(:,k,i) , nom_champ//tcomp//'_m='//tmode//'_'//tit_s//'_'//tit//'.plt')
1256  ENDDO
1257 
1258  END DO
1259 
1260  END SUBROUTINE trace_profile
1261 
1262 END MODULE sub_plot
sub_plot::plot_arrow_label
subroutine plot_arrow_label(jj, rr, vv, file_name)
Definition: sub_plot.f90:468
sub_plot
Definition: sub_plot.f90:4
sub_plot::plot_stream
subroutine plot_stream(jj, rr, pp, i)
Definition: sub_plot.f90:665
sub_plot::plot_pressure_2d
subroutine plot_pressure_2d(jj, rr, uu)
Definition: sub_plot.f90:117
sub_plot::plot_ensight_vecteur
subroutine plot_ensight_vecteur(u8, vit)
Definition: sub_plot.f90:584
sub_plot::plot_scalar_field
subroutine plot_scalar_field(jj, rr, uu, file_name)
Definition: sub_plot.f90:704
def_type_mesh::mesh_type
Definition: def_type_mesh.f90:74
sub_plot::vtk_p1_2d
subroutine vtk_p1_2d(mesh, champ, unit_file, file_name)
Definition: sub_plot.f90:1062
chaine_caractere::eval_blank
integer function eval_blank(len_str, string)
Definition: chaine_caractere.f90:33
sub_plot::plot_pressure_p1_label
subroutine plot_pressure_p1_label(jj, rr, uu, file_name)
Definition: sub_plot.f90:298
sub_plot::plot_const_p1_label
subroutine plot_const_p1_label(jj, rr, uu, file_name)
Definition: sub_plot.f90:265
sub_plot::plot_two_scalar_field
subroutine plot_two_scalar_field(jj, rr, uu, jj2, rr2, uu2, file_name)
Definition: sub_plot.f90:854
sub_plot::plot_vorticity
subroutine plot_vorticity(jj, rr, zz, i)
Definition: sub_plot.f90:628
sub_plot::plot_vit_2d
subroutine plot_vit_2d(jj, rr, uu)
Definition: sub_plot.f90:9
sub_plot::vtk_2d
subroutine vtk_2d(mesh, champ, unit_file, file_name)
Definition: sub_plot.f90:1109
sub_plot::plot_arrow
subroutine plot_arrow(jj, rr, vv)
Definition: sub_plot.f90:73
sub_plot::plot_scalar_field_domain
subroutine plot_scalar_field_domain(jj, rr, id, index_dom, file_name)
Definition: sub_plot.f90:981
sub_plot::plot_pressure_label
subroutine plot_pressure_label(jj, rr, uu, file_name)
Definition: sub_plot.f90:433
sub_plot::plot_p1_cont_label
subroutine plot_p1_cont_label(jj, jjs, sides, list, rr, uu, file_name)
Definition: sub_plot.f90:332
sub_plot::plot_p1_matiere_label
subroutine plot_p1_matiere_label(jj, neigh, i_d, rr, uu, file_name)
Definition: sub_plot.f90:377
sub_plot::trace_profile
subroutine trace_profile(mesh, v, it, freq_plot, list_mode, nom_champ, num_dom)
Definition: sub_plot.f90:1209
sub_plot::plot_loc_rel_var
subroutine plot_loc_rel_var(jj, rr, vo, vv)
Definition: sub_plot.f90:183
chaine_caractere
Definition: chaine_caractere.f90:4
def_type_mesh
Definition: def_type_mesh.f90:4
sub_plot::plot_ensight_scalaire
subroutine plot_ensight_scalaire(p8, pres)
Definition: sub_plot.f90:607
sub_plot::plot_pressure
subroutine plot_pressure(jj, rr, uu)
Definition: sub_plot.f90:236
sub_plot::plot_pressure_p2_label
subroutine plot_pressure_p2_label(jj, rr, uu, file_name)
Definition: sub_plot.f90:522
form
section doc_intro_frame_work_num_app Numerical approximation subsection doc_intro_fram_work_num_app_Fourier_FEM Fourier Finite element representation The SFEMaNS code uses a hybrid Fourier Finite element formulation The Fourier decomposition allows to approximate the problem’s solutions for each Fourier mode modulo nonlinear terms that are made explicit The variables are then approximated on a meridian section of the domain with a finite element method The numerical approximation of a function f $f f is written in the following generic form
Definition: doc_intro.h:190