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