38
39
40
44 USE create_h3d_output_per_part_mod
45
46
47
48#include "implicit_f.inc"
49
50
51
52#include "com04_c.inc"
53#include "scr17_c.inc"
54
55
56
57 TYPE (H3D_DATABASE) :: H3D_DATA
58 INTEGER ID_1D_SCALAR,ID_INPUT,STEXT1,SCOMMENT,IPT
59 INTEGER IPART(LIPART1,*)
60 CHARACTER(LEN=STEXT1) :: TEXT1
61 CHARACTER(LEN=NCHARLINE100) :: KEY3_GLOB
62 CHARACTER(LEN=SCOMMENT) :: COMMENT
63 TYPE(GROUP_) ,DIMENSION(NGRPART) ,INTENT(IN) :: IGRPART
64
65
66
67 INTEGER NVAR
68
69
70
71 CHARACTER(LEN=NCHARLINE100) :: CARTE,CARTE1
72 INTEGER I,J,L,N_H3D_PART
73
74 h3d_data%N_OUTP_H3D = h3d_data%N_OUTP_H3D + 1
75 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%OK = 1
76 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%KEYWORD = key3_glob
77 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%N_OUTP = h3d_data%N_OUTP_H3D + 3
78 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%ETYPE = 4
79 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%OUTP_TYPE = 1
80 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IPT = ipt
81 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IS_INFO1 = 0
82 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_STRING1 = stext1
83 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%STRING1(1:stext1) = text1(1:stext1)
84 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_STRING2 = 0
85 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_COMMENT = scomment
86 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%COMMENT(1:scomment) = comment(1:scomment)
87
88
89
90
91 n_h3d_part = h3d_data%INPUT_LIST(id_input)%NB_PART
92 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%N_H3D_PART_LIST = n_h3d_part
93
94 IF (n_h3d_part /= 0 )
95 . ALLOCATE(h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART_LIST(n_h3d_part))
96 ALLOCATE(h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(npart))
97
98
99
100
101
102
103 IF (n_h3d_part /= 0 ) THEN
104 DO i=1,npart
105 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(i) = 0
106 ENDDO
107 ELSE
108 DO i=1,npart
109 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(i) = 1
110 ENDDO
111 ENDIF
112
113
114
115
116
117 CALL create_h3d_output_per_part(n_h3d_part,h3d_data,id_input,lipart1,npart,ipart,ngrpart,igrpart)
118
119
120
121
122
123 RETURN
124 999 print *,'error lecture'
integer, parameter ncharline100