xref: /illumos-gate/usr/src/test/util-tests/tests/awk/gnu/tweakfld.awk (revision e6d6c189fa3a95d7aa27bbe0aeacf7c1a6b57c8c)
1# To: bug-gnu-utils@prep.ai.mit.edu
2# Cc: arnold@gnu.ai.mit.edu
3# Date: Mon, 20 Nov 1995 11:39:29 -0500
4# From: "R. Hank Donnelly" <emory!head-cfa.harvard.edu!donnelly>
5#
6# Operating system: Linux1.2.13 (Slackware distrib)
7# GAWK version: 2.15 (?)
8# compiler: GCC (?)
9#
10# The following enclosed script does not want to fully process the input data
11# file. It correctly executes the operations on the first record, and then dies
12# on the second one. My true data file is much longer but this is
13# representative and it does fail on a file even as short as this one.
14# The failure appears to occur in the declared function add2output. Between the
15# steps of incrementing NF by one and setting $NF to the passed variable
16# the passed variable appears to vanish (i.e. NF does go from 68 to 69
17# and before incrementing it "variable" equals what it should but after
18# "variable" has no value at all.)
19#
20# The scripts have been developed using nawk on a Sun (where they run fine)
21# I have tried gawk there but get a different crash which I have not yet traced
22# down. Ideally I would like to keep the script the same so that it would run
23# on either gawk or nawk (that way I can step back and forth between laptop and
24# workstation.
25#
26# Any ideas why the laptop installation is having problems?
27# Hank
28#
29#
30# #!/usr/bin/gawk -f
31
32BEGIN {
33	# set a few values
34	FS = "\t"
35	OFS = "\t"
36	pi = atan2(0, -1)
37# distance from HRMA to focal plane in mm
38	fullradius = 10260.54
39
40	# set locations of parameters on input line
41	nf_nrg = 1
42	nf_order = 3
43	nf_item = 4
44	nf_suite = 5
45	nf_grating = 8
46	nf_shutter = 9
47	nf_type = 13
48	nf_src = 14
49	nf_target = 15
50	nf_voltage = 16
51	nf_flux = 17
52	nf_filt1 = 20
53	nf_filt1_th = 21
54	nf_filt2 = 22
55	nf_filt2_th = 23
56	nf_bnd = 24
57	nf_hrma_polar = 27
58	nf_hrma_az = 28
59	nf_detector = 30
60	nf_acis_read = 32
61	nf_acis_proc = 33
62	nf_acis_frame = 34
63	nf_hxda_aplist = 36
64	nf_hxda_y_range = 37
65	nf_hxda_z_range = 38
66	nf_hxda_y_step = 39
67	nf_hxda_z_step = 40
68	nf_sim_z = 41
69	nf_fam_polar = 43
70	nf_fam_az = 44
71	nf_fam_dither_type = 45
72	nf_mono_init = 51
73	nf_mono_range = 52
74	nf_mono_step = 53
75	nf_defocus = 54
76	nf_acis_temp = 55
77	nf_tight = 59
78	nf_offset_y = 64
79	nf_offset_z = 65
80
81	while( getline < "xrcf_mnemonics.dat" > 0 ) {
82		mnemonic[$1] = $2
83	}
84
85#	"date" | getline date_line
86# ADR: use a fixed date so that testing will work
87	date_line = "Sun Mar 10 23:00:27 EST 1996"
88        split(date_line, in_date, " ")
89        out_date = in_date[2] " " in_date[3] ", " in_date[6]
90}
91
92function add2output( variable ) {
93#print("hi1") >> "debug"
94	NF++
95#print("hi2") >> "debug"
96 	$NF = variable
97#print("hi3") >> "debug"
98}
99
100function error( ekey, message ) {
101	print "Error at input line " NR ", anode " ekey >> (ENVIRON["WORKDIR"] "/test.temp.errors.cleanup")
102	print "   " message "." >> (ENVIRON["WORKDIR"] "/test.temp.errors.cleanup")
103}
104
105function hxda_na() {
106	$nf_hxda_aplist = $nf_hxda_y_range = $nf_hxda_z_range = "N/A"
107	$nf_hxda_y_step = $nf_hxda_z_step = "N/A"
108}
109
110function acis_na() {
111	$nf_acis_read = $nf_acis_proc = $nf_acis_frame = $nf_acis_temp = "N/A"
112}
113
114function hrc_na() {
115#        print ("hi") >> "debug"
116}
117
118function fpsi_na() {
119	acis_na()
120	hrc_na()
121	$nf_sim_z = $nf_fam_polar = $nf_fam_az = $nf_fam_dither_type = "N/A"
122}
123
124function mono_na() {
125	$nf_mono_init = $nf_mono_range = $nf_mono_step = "N/A"
126}
127
128# this gives the pitch and yaw of the HRMA and FAM
129# positive pitch is facing the source "looking down"
130# positive yaw is looking left
131# 0 az is north 90 is up
132# this also adds in the FAM X,Y,Z positions
133
134function polaz2yawpitch(polar, az) {
135	theta = az * pi / 180
136	phi = polar * pi / 180 / 60
137
138
139	if( polar == 0 ) {
140		add2output( 0 )
141		add2output( 0 )
142	} else {
143		if(az == 0 || az == 180)
144			add2output( 0 )
145		else
146			add2output( - polar * sin(theta) )
147
148
149#			x = cos (phi)
150#			y = sin (phi) * cos (theta)
151#			add2output( atan2(y,x)*180 / pi * 60 )
152
153		if(az == 90 || az ==270 )
154			add2output( 0 )
155		else
156			add2output( - polar * cos(theta) )
157
158	}
159#			x = cos (phi)
160#			z= sin (phi) * sin (theta)
161#			add2output( atan2(z,x)*180 / pi * 60 )
162
163	if(config !~ /HXDA/) {
164# negative values of defocus move us farther from the source thus
165# increasing radius
166		radius = fullradius - defocus
167
168# FAM_x; FAM_y;  FAM_z
169	   	if((offset_y == 0) && (offset_z == 0)){
170			add2output( fullradius - radius * cos (phi) )
171
172			if (az == 90 || az ==270)
173				add2output( 0 )
174			else
175				add2output(  radius * sin (phi) * cos (theta) )
176
177			if (az == 0 || az == 180)
178				add2output( 0 )
179			else
180				add2output( - radius * sin (phi) * sin (theta) )
181	   	} else {
182# ******* THIS SEGMENT OF CODE IS NOT MATHEMATICALLY CORRECT FOR ****
183# OFF AXIS ANGLES AND IS SUPPLIED AS A WORKAROUND SINCE IT WILL
184# PROBABLY ONLY BE USED ON AXIS.
185			add2output( defocus )
186			add2output( offset_y )
187			add2output( offset_z )
188		}
189
190	} else {
191		add2output( "N/A" )
192		add2output( "N/A" )
193		add2output( "N/A" )
194	}
195}
196
197# set TIGHT/LOOSE to N/A if it is not one of the two allowed values
198function tight_na() {
199	if( $nf_tight !~ /TIGHT|LOOSE/ ) {
200		$nf_tight == "N/A"
201	}
202}
203
204# this entry is used to give certain entries names
205{
206	type = $nf_type
207	item = $nf_item
208	suite = $nf_suite
209	order = $nf_order
210	detector = $nf_detector
211	grating = $nf_grating
212	offset_y= $nf_offset_y
213	offset_z= $nf_offset_z
214	bnd = $nf_bnd
215	defocus = $nf_defocus
216}
217
218{
219	# make configuration parameter
220	# as well as setting configuration-dependent N/A values
221
222	if( $nf_bnd ~ "SCAN" ) {
223		# BND is scanning beam
224		config = "BND"
225		hxda_na()
226		fpsi_na()
227	} else {
228		if( grating == "NONE" ) {
229			config = "HRMA"
230		} else {
231			if( grating == "HETG" ) {
232				if( order != "Both" ) {
233				    $nf_shutter = order substr($nf_shutter, \
234					index($nf_shutter, ",") )
235				}
236			} else {
237				order = "N/A"
238			}
239			config = "HRMA/" grating
240		}
241
242		if( detector ~ /ACIS|HRC/ ) {
243			detsys = detector
244			nsub = sub("-", ",", detsys)
245			config = config "/" detsys
246			hxda_na()
247		} else {
248			config = config "/HXDA"
249			fpsi_na()
250			if( detector == "HSI" ) {
251				hxda_na()
252			}
253		}
254	}
255
256	add2output( config )
257
258	if( $nf_src ~ /EIPS|Penning/ ) mono_na()
259
260	if( $nf_src == "Penning" ) $nf_voltage = "N/A"
261
262	itm = sprintf("%03d", item)
263
264	if(config in mnemonic) {
265		if( type in mnemonic ) {
266		    ID = mnemonic[config] "-" mnemonic[type] "-" suite "." itm
267		    add2output( ID )
268		} else {
269			error(type, "measurement type not in list")
270		}
271	} else {
272		error(config, "measurement configuration not in list")
273	}
274
275	# add date to output line
276	add2output( out_date )
277
278	# Convert HRMA polar and azimuthal angles to yaw and pitch
279	polaz2yawpitch($nf_hrma_polar, $nf_hrma_az)
280
281	# set TIGHT/LOOSE to N/A if it is not one of the two allowed values
282	tight_na()
283
284	# compute number of HXDA apertures
285	if( config ~ /HXDA/ && $nf_hxda_aplist != "N/A")
286		add2output( split( $nf_hxda_aplist, dummy, "," ) )
287	else
288		add2output( "N/A" )
289
290	# make sure the BND value is properly set
291	if($nf_bnd == "FIXED" && detector ~ /ACIS/)
292		$nf_bnd =bnd"-SYNC"
293	else
294		$nf_bnd = bnd"-FREE"
295	print
296}
297