41#include "implicit_f.inc"
50 CHARACTER (LEN=255) :: VARNAME
52 DOUBLE PRECISION TEMP_DOUBLE
54 INTEGER :: NDETPS,NDETSG,NECRAN,NDETPL,NDETCORD
60 ndetps = detonators%N_DET_POINT
61 ndetsg = detonators%N_DET_LINE
62 necran = detonators%N_DET_WAVE_SHAPER
63 ndetpl = detonators%N_DET_PLANE
64 ndetcord = detonators%N_DET_CORD
70 WRITE(varname,
'(A,I0,A)')
'DFS_DETPOINT_',k ,
'__ISHADOW_'
71 temp_integer = detonators%POINT(k)%SHADOW
72 IF(temp_integer > 0)
CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
74 WRITE(varname,
'(A,I0,A)')
'DFS_DETPOINT_',k ,
'__MAT_'
75 temp_integer = detonators%POINT(k)%MAT
76 CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
78 WRITE(varname,
'(A,I0,A)')
'DFS_DETPOINT_',k ,
'__TDET_'
79 temp_double = detonators%POINT(k)%TDET
80 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
82 WRITE(varname,
'(A,I0,A)')
'DFS_DETPOINT_',k ,
'__XDET_'
83 temp_double = detonators%POINT(k)%XDET
84 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
86 WRITE(varname,
'(A,I0,A)')
'DFS_DETPOINT_',k ,
'__YDET_'
87 temp_double = detonators%POINT(k)%YDET
88 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
90 WRITE(varname,
'(A,I0,A)')
'DFS_DETPOINT_',k ,
'__ZDET_'
91 temp_double = detonators%POINT(k)%ZDET
92 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
94 WRITE(varname,
'(A,I0,A)')
'DFS_DETPOINT_',k ,
'__GRNOD_ID_'
95 temp_integer = detonators%POINT(k)%GRNOD_ID
96 IF(temp_integer > 0)
CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
102 WRITE(varname,
'(A,I0,A)')
'DFS_DETLINE_',k ,
'__MAT_'
103 temp_integer = detonators%LINE(k)%MAT
104 CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
106 WRITE(varname,
'(A,I0,A)')
'DFS_DETLINE_',k ,
'__TDET_'
107 temp_double = detonators%LINE(k)%TDET
108 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim
110 WRITE(varname,
'(A,I0,A)')
'DFS_DETLINE_',k ,
'__XDET_1_'
111 temp_double = detonators%LINE(k)%XDET_1
112 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
114 WRITE(varname,
'(A,I0,A)')
'DFS_DETLINE_',k ,
'__YDET_1_'
115 temp_double = detonators%LINE(k)%YDET_1
116 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
118 WRITE(varname,
'(A,I0,A)')
'DFS_DETLINE_',k ,
'__ZDET_1_'
119 temp_double = detonators%LINE(k)%ZDET_1
120 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
122 WRITE(varname,
'(A,I0,A)')
'DFS_DETLINE_',k ,
'__XDET_2_'
123 temp_double = detonators%LINE(k)%XDET_2
124 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
126 WRITE(varname,
'(A,I0,A)')
'DFS_DETLINE_',k ,
'__YDET_2_'
127 temp_double = detonators%LINE(k)%YDET_2
128 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
130 WRITE(varname,
'(A,I0,A)')
'DFS_DETLINE_',k ,
'__ZDET_2_'
131 temp_double = detonators%LINE(k)%ZDET_2
132 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
137 WRITE(varname,
'(A,I0,A)')
'DFS_WAVSHA_',k ,
'__MAT_'
138 temp_integer = detonators%WAVE_SHAPER(k)%MAT
139 CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
141 WRITE(varname,
'(A,I0,A)')
'DFS_WAVSHA_',k ,
'__NUMNOD_'
142 temp_integer = detonators%WAVE_SHAPER(k)%NUMNOD
143 CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
145 WRITE(varname,
'(A,I0,A)')
'DFS_WAVSHA_',k ,
'__VDET_'
146 temp_double = detonators%WAVE_SHAPER(k)%VDET
147 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
149 WRITE(varname,
'(A,I0,A)')
'DFS_WAVSHA_',k ,
'__TDET_'
150 temp_double = detonators%WAVE_SHAPER(k)%TDET
151 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
153 WRITE(varname,
'(A,I0,A)')
'DFS_WAVSHA_',k ,
'__XDET_'
154 temp_double = detonators%WAVE_SHAPER(k)%XDET
155 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
157 WRITE(varname,
'(A,I0,A)')
'DFS_WAVSHA_',k ,
'__YDET_'
158 temp_double = detonators%WAVE_SHAPER(k)%YDET
159 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
161 WRITE(varname,
'(A,I0,A)')
'DFS_WAVSHA_',k ,
'__ZDET_'
162 temp_double = detonators%WAVE_SHAPER(k)%ZDET
163 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
168 WRITE(varname,
'(A,I0,A)')
'DFS_DETPLANE_',k ,
'__MAT_'
169 temp_integer = detonators%PLANE(k)%MAT
170 CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
172 WRITE(varname,
'(A,I0,A)')
'DFS_DETPLANE_',k ,
'__TDET_'
173 temp_double = detonators%PLANE(k)%TDET
174 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
176 WRITE(varname,
'(A,I0,A)')
'DFS_DETPLANE_',k ,
'__XDET_'
177 temp_double = detonators%PLANE(k)%XDET
178 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
180 WRITE(varname,
'(A,I0,A)')
'DFS_DETPLANE_',k ,
'__YDET_'
181 temp_double = detonators%PLANE(k)%YDET
182 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
184 WRITE(varname,
'(A,I0,A)')
'DFS_DETPLANE_',k ,
'__ZDET_'
185 temp_double = detonators%PLANE(k)%ZDET
186 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
188 WRITE(varname,
'(A,I0,A)')
'DFS_DETPLANE_',k ,
'__NX_'
189 temp_double = detonators%PLANE(k)%NX
190 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
192 WRITE(varname,
'(A,I0,A)')
'DFS_DETPLANE_',k ,
'__NY_'
193 temp_double = detonators%PLANE(k)%NY
194 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
196 WRITE(varname,
'(A,I0,A)')
'DFS_DETPLANE_',k ,
'__NZ_'
197 temp_double = detonators%PLANE(k)%NZ
198 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
203 WRITE(varname,
'(A,I0,A)')
'DFS_DETCORD_',k ,
'__MAT_'
204 temp_integer = detonators%CORD(k)%MAT
205 CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
207 WRITE(varname,
'(A,I0,A)')
'DFS_DETCORD_',k ,
'__IOPT_'
208 temp_integer = detonators%CORD(k)%IOPT
209 CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
211 WRITE(varname,
'(A,I0,A)')
'DFS_DETCORD_',k ,
'__NUMNOD_'
212 temp_integer = detonators%CORD(k)%NUMNOD
213 CALL qaprint(varname(1:len_trim(varname)),temp_integer,0.0_8)
215 WRITE(varname,
'(A,I0,A)')
'DFS_DETCORD_',k ,
'__TDET_'
216 temp_double = detonators%CORD(k)%TDET
217 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
219 WRITE(varname,
'(A,I0,A)')
'DFS_DETCORD_',k ,
'__VDET_'
220 temp_double = detonators%CORD(k)%VDET
221 IF(temp_double/=zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)