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_SOLID_TENSOR,ID_INPUT,STEXT1,SCOMMENT,
55 . LAYER,,IR,IS,IT
56 INTEGER IPART(LIPART1,*)
57 CHARACTER(LEN=STEXT1) :: TEXT1
58 CHARACTER(LEN=NCHARLINE100) :: KEY3_GLOB
59 CHARACTER(LEN=SCOMMENT) :: COMMENT
60 TYPE(GROUP_) ,DIMENSION(NGRPART) ,INTENT(IN)
61
62
63
64# "com04_c.inc"
65#include "scr17_c.inc"
66
67
68
69 INTEGER NVAR
70
71
72
73 CHARACTER(LEN=NCHARLINE100) :: CARTE,CARTE1
74 INTEGER I,J,L,N_H3D_PART
75
76 h3d_data%N_OUTP_H3D = h3d_data%N_OUTP_H3D + 1
77 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%OK = 1
78 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%KEYWORD = key3_glob
79 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%N_OUTP = h3d_data%N_OUTP_H3D + 3
80 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%ETYPE = 7
81 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%OUTP_TYPE = 3
82 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%LAYER = layer
83 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IR = ir
84 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IS = is
85 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IT = it
86 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IUVAR = iuvar
87 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%IS_INFO1 = 0
88 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_STRING1
89 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%STRING1(1:stext1) = text1
90 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_STRING2 = 0
91 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%S_COMMENT = scomment
92 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%COMMENT(1:scomment) = comment(1:scomment)
93
94 n_h3d_part = h3d_data%INPUT_LIST(id_input)%NB_PART
95 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%N_H3D_PART_LIST = n_h3d_part
96
97 IF (n_h3d_part /= 0 )
98 . ALLOCATE(h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART_LIST(n_h3d_part))
99 ALLOCATE(h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(npart))
100
101 IF (n_h3d_part /= 0 ) THEN
102 DO i=1,npart
103 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(i) = 0
104 ENDDO
105 ELSE
106 DO i=1,npart
107 h3d_data%OUTPUT_LIST(h3d_data%N_OUTP_H3D)%PART(i) = 1
108 ENDDO
109 ENDIF
110
111 CALL create_h3d_output_per_part(n_h3d_part,h3d_data,id_input,lipart1,npart,ipart,ngrpart,igrpart)
112
113
114
115
116
117 RETURN
118 999 print *,'error lecture'
integer, parameter ncharline100