1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <syserr.h>
include <fset.h>
include <imhdr.h>
include <imio.h>
include "idbc.h"
# IMADDFI -- Insert a user field in the image header after the specified
# keyword. It is an error if the named field already exists.
#procedure imaddfi (im, key, pkey, datatype, baf)
procedure iminfi (im, key, pkey, datatype, baf)
pointer im #I image descriptor
char key[ARB] #I name of the new parameter
char pkey[ARB] #I 'key' will be inserted bef/after pkey
char datatype[ARB] #I string permits generalization to domains
int baf # I Insert BEFORE or AFTER
pointer rp, sp, keyname, ua, ip
int fd, max_lenuserarea, curlen, buflen, nchars, piv
int idb_kwlookup(), idb_findrecord()
int strlen(), idb_filstr(), nowhite()
char card[IDB_RECLEN+1]
errchk syserrs, sprintf, pargstr, pargi
begin
call smark (sp)
call salloc (keyname, SZ_FNAME, TY_CHAR)
nchars = idb_filstr (key, Memc[keyname], IDB_SZFITSKEY)
nchars = nowhite (Memc[keyname], Memc[keyname], IDB_SZFITSKEY)
call strupr (Memc[keyname])
# Check for a redefinition.
if ((idb_kwlookup (key) > 0) || (idb_findrecord (im, key, rp) > 0))
call syserrs (SYS_IDBREDEF, key)
# Open the user area string for appending. 'buflen' is the malloc-ed
# buffer length in struct units; IMU is the struct offset to the user
# area, i.e., the size of that part of the image descriptor preceding
# the user area.
ua = IM_USERAREA(im)
curlen = strlen (Memc[ua])
buflen = LEN_IMDES + IM_LENHDRMEM(im)
max_lenuserarea = (buflen - IMU) * SZ_STRUCT - 1
if (curlen+81 >= max_lenuserarea) {
IM_HDRLEN(im) = LEN_IMHDR +
(curlen + 10*36*81 + SZ_STRUCT-1) / SZ_STRUCT
IM_LENHDRMEM(im) = IM_HDRLEN(im) + (SZ_UAPAD / SZ_STRUCT)
call realloc (im, IM_LENHDRMEM(im) + LEN_IMDES, TY_STRUCT)
buflen = LEN_IMDES + IM_LENHDRMEM(im)
max_lenuserarea = (buflen - IMU) * SZ_STRUCT - 1
}
# If the user area is not empty the last character must be the newline
# record delimiter, else the new record we add will be invalid.
if (curlen > 0 && Memc[ua+curlen-1] != '\n')
if (curlen >= max_lenuserarea) {
call syserrs (SYS_IDBOVFL, key)
} else {
Memc[ua+curlen] = '\n'
curlen = curlen + 1
Memc[ua+curlen] = EOS
}
# Find keyw_after
if (idb_findrecord (im, pkey, rp) == 0) {
# Keyw not found. Append the new keyword.
rp = ua+curlen
baf = BEFORE
} else {
# Shift cards after pivot or before pivot
if (baf == AFTER)
piv = rp
else
piv = rp - IDB_RECLEN - 1
for (ip= ua+curlen-IDB_RECLEN-1; ip>=piv; ip=ip-IDB_RECLEN-1) {
call amovc (Memc[ip], Memc[ip+IDB_RECLEN+1], IDB_RECLEN)
}
}
Memc[ua+curlen+IDB_RECLEN]='\n'
Memc[ua+curlen+IDB_RECLEN+1]=EOS
# Form a card with keyword name and placeholder for value.
call sprintf (card, IDB_RECLEN+10, "%-8s= %s%*t\n")
call pargstr (Memc[keyname])
if (datatype[1] == 'c') {
call pargstr ("' '")
call pargi (IDB_LENSTRINGRECORD + 1)
} else {
call pargstr ("")
call pargi (IDB_LENNUMERICRECORD + 1)
}
# Replace keyword at the position rp+81.
if (baf == AFTER)
call amovc (card, Memc[rp+IDB_RECLEN+1], IDB_RECLEN)
else
call amovc (card, Memc[rp], IDB_RECLEN)
#for (ip=1; ip<5; ip=ip+1) {
#call eprintf("<%40.40s>\n")
# call pargstr(Memc[rp-(2-ip)*(IDB_RECLEN+1)])
#}
call sfree (sp)
end
|
