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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include "imwcs.h"
# IW_GBIGFITS -- Get a FITS string valued parameter of arbitrary size.
# Since a FITS string stored as a single parameter is limited to at most
# 69 characters, multiple FITS cards must be used to store longer strings.
# At the time that this routine is called, IW_RFITS has already been called
# to scan the FITS header and build up a list of WCS oriented cards,
# including card types and pointers to the card data. Our job is merely
# to take these cards in order and concatenate the string values into one
# large string, returning a pointer to the string as the function value.
# The caller must later make a MFREE call to free this buffer.
pointer procedure iw_gbigfits (iw, ctype, axis)
pointer iw #I pointer to IMWCS descriptor
int ctype #I card type
int axis #I axis to which card refers
int ncards, i, j, ch
pointer cp, bp, ip, op, rp
define put_ 10
begin
# How much space do we need?
ncards = 0
do i = 1, IW_NCARDS(iw) {
cp = IW_CARD(iw,i)
if (C_AXIS(cp) == axis && C_TYPE(cp) == ctype)
ncards = ncards + 1
}
# Allocate the space.
call calloc (bp, ncards * MAX_FITSCOLS, TY_CHAR)
# For successive cards 1, 2, 3, etc...
op = bp
do j = 1, ncards {
# Find the card.
rp = NULL
do i = 1, IW_NCARDS(iw) {
cp = IW_CARD(iw,i)
if (C_AXIS(cp) != axis)
next
if (C_INDEX(cp) != j)
next
if (C_TYPE(cp) != ctype)
next
rp = C_RP(cp)
break
}
# Append to the string buffer.
if (rp != NULL) {
#call amovc (Memc[rp+IDB_STARTVALUE+1], Memc[op], MAX_FITSCOLS)
#op = op + MAX_FITSCOLS
do i = 1, MAX_FITSCOLS {
ip = rp + IDB_STARTVALUE + i
ch = Memc[ip]
if (ch == EOS || ch == '\n') {
break
} else if (ch == '\'') {
if (Memc[ip+1] == '\'') {
goto put_
} else if (Memc[ip-1] == '\'') {
;
} else if (i > 1 && i <= MAX_FITSCOLS) {
# If we're not at the end of the card, we have a
# complete string, but add a space for appending
# so we don't concatenate.
Memc[op] = ' '
op = op + 1
break
} else
break
} else {
put_ Memc[op] = ch
op = op + 1
}
}
}
}
Memc[op] = EOS
return (bp)
end
|
