38
39
40
44 USE create_h3d_output_per_part_mod
45
46
47
48#include "implicit_f.inc"
49
50
51
52 TYPE (H3D_DATABASE) :: H3D_DATA
53 INTEGER ID_QUAD_TENSOR,ID_INPUT,STEXT1,SCOMMENT,LAYER,IR,IS,IT,IUVAR
54 INTEGER IPART(LIPART1,*)
55 CHARACTER(LEN=STEXT1) :: TEXT1
56 CHARACTER(LEN=NCHARLINE100) :: KEY3_GLOB
57 CHARACTER(LEN=SCOMMENT) :: COMMENT
58 TYPE(GROUP_) ,DIMENSION(NGRPART) ,INTENT(IN) :: IGRPART
59
60
61
62#include "com04_c.inc"
63#include "scr17_c.inc"
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 = 6
79 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%OUTP_TYPE = 3
80 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%LAYER = layer
81 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IR = ir
82 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IS = is
83 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IT = it
84 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IUVAR = iuvar
85 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IS_INFO1 = 0
86 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_STRING1 = stext1
87 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%STRING1(1:stext1) = text1(1:stext1)
88 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_STRING2 = 0
89 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_COMMENT = scomment
90 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%COMMENT(1:scomment) = comment(1:scomment)
91
92 n_h3d_part = h3d_data%INPUT_LIST(id_input)%NB_PART
93 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%N_H3D_PART_LIST = n_h3d_part
94
95 IF (n_h3d_part /= 0 )
96 . ALLOCATE(h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART_LIST(n_h3d_part))
97 ALLOCATE(h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(npart))
98
99 IF (n_h3d_part /= 0 ) THEN
100 DO i=1,npart
101 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(i) = 0
102 ENDDO
103 ELSE
104 DO i=1,npart
105 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(i) = 1
106 ENDDO
107 ENDIF
108
109 CALL create_h3d_output_per_part(n_h3d_part,h3d_data,id_input,lipart1,npart,ipart,ngrpart,igrpart)
110
111
112
113
114
115 RETURN
116 999 print *,'error lecture'
integer, parameter ncharline100