blob: 2a45abf4be357a46b185c6d0a0e198b72db3984d (
plain) (
blame)
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
|
C----------------------------------------------------------------------
subroutine ftuscc(input,np,scaled,scale,zero,output)
C unscale the array of complex numbers, prior to writing to the FITS file
C input r array of complex numbers (pairs of real/imaginay numbers)
C np i total number of values to scale (no. of pairs times 2)
C scaled l is the data scaled?
C scale d scale factor
C zero d offset
C output r output array
integer np,i,j
logical scaled
real input(np),output(np)
double precision scale,zero
j=1
if (scaled)then
do 10 i=1,np/2
output(j)=(input(j)-zero)/scale
j=j+1
C the imaginary part of the number is not offset!!
output(j)=input(j)/scale
j=j+1
10 continue
else
do 20 i=1,np
output(i)=input(i)
20 continue
end if
end
|
