aboutsummaryrefslogtreecommitdiff
path: root/noao/imred/quadred/src/quad/quadtest/artobs.cl
blob: 292ed8c6033a482e24425254a3509b656522d9df (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
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
# ARTOBS -- Simulate observe command using artificial data.

procedure artobs ()

begin
	string	image, oim, ccdt
	int	picnum, nexps, i
	real	exptime
	string	imtitle

	# Get ccdtype
	ccdt = obspars.ccdtype

	# Get number of pictures to take
	nexps   = npics

	# Get next picture number
	if (obspars.autopicnum) {
	    picnum = obspars.picture.p_value
	} else {
	    picnum = obspars.picture
	}

	# Get exposure time
	if (ccdt != "zero") {
	    exptime = obspars.exposure
	} else {
	    exptime = 0.0
	}

	# Set filter
	if (obspars.setfilter != "none" && ccdt != "zero" && ccdt != "dark") {
	    if (instrpars.instrname != "")
		mkquad.filter = instrpars.filter1
	}

	# Get imtitle. This MUST always be the last interactive prompt!
	imtitle = title

	for (i = picnum; i < picnum+nexps; i = i+1) {

	    # Make image name
	    if (ccdt == "object") {
		printf ("obj%03d\n", i) | scan (image)
	    } else {
		printf ("%s%03d\n", ccdt, i) | scan (image)
	    }
	    if (access (image//".imh")) {
		oim = image
		image = mktemp (image//".")
		printf ("Output image %s already exists...  %s used \n", oim, 
		image)
	    }

	    if (ccdt == "dflat" || ccdt == "pflat")
		ccdt = "flat"

	    if (ccdt == "sflat" || ccdt == "comp" )
	       ccdt = "other"

	    # Call MKQUAD task
	    mkquad (image, exptime, ccdt)
	    hedit  (image, "i_title", imtitle, add+, ver-, show-)
	    obspars.picture.p_value = i + 1
	    printf ("Image %s written to disk\n", image, > "STDERR")
	}

end