Initial commit

1 files changed, 168 insertions, 0 deletions

diff --git a/noao/onedspec/irsiids/batchred.cl b/noao/onedspec/irsiids/batchred.cl
new file mode 100644
index 00000000..1bd5a3a3
--- /dev/null
+++ b/noao/onedspec/irsiids/batchred.cl
@@ -0,0 +1,168 @@
+#{ BATCHRED -- Script file to generate another script file
+#  which runs several ONEDSPEC tasks in an automated fashion.
+#
+#  Currently the following procedures are automated:
+#    1. STANDARD
+#    2. SENSFUNC
+#    3. BSWITCH
+#    4. CALIBRATE
+#    5. ADDSETS
+#
+
+{
+# Say hello to the guy on the other side of the screen and check batch file.
+print ("\n----B A T C H   I I D S / I R S   F I L E   G E N E R A T O R----\n")
+
+s2 = "process.cl"			# Batch file to be created.
+if (access (s2)) {
+    print ("A batch file already exists - ")
+    if (query)
+        delete (s2, verify=no)
+}
+
+# Initialize
+rt = input				# Root name for spectra
+ot = output				# Output root name
+ttl = ">>&'" + ttylog + "')\n"		# Log file for tty output
+
+out = rt
+s1 = ""
+st = ""
+sns = ""
+stat = ""
+print ("i = ", start_rec, >>s2)
+
+if (standard) {				# STANDARD?
+    print ("\n#---STANDARD---\n")
+    print ("\n#---STANDARD---\n", >>s2)
+
+    st = std				# STD file
+    if (access (st)) {
+	print (st, " - already exists")
+	if (query)
+            delete (st, verify=no)
+    }
+
+    # Loop over all stars
+    b1 = yes
+    while (b1) {
+        # Check that the last entry was different - otherwise end input
+	records = ""
+    	s3 = records
+    	if (s3 == "")
+    	    b1 = no
+        else {
+    	    print ("standard (input='",rt,"',output='",st,"',",>>s2)
+    	    print ("\trecords='",s3,"',",>>s2)
+    	    print ("\tstar_name='",star_name,"',beam_switch=yes,",>>s2)
+	    print ("\tsamestar=yes,apertures='',bandwidth=INDEF,",>>s2)
+	    print ("\tbandsep=INDEF,interact=no,",ttl,>>s2)
+        }
+    }
+
+    print ("")
+}
+
+
+if (sensfunc) {				# SENSFUNC?
+    print ("\n#---SENSFUNC---\n")
+    print ("\n#---SENSFUNC---\n", >>s2)
+
+    if (st == "")
+	st = std			# STD file
+    sns = sensitivity			# Sensitivity image
+    stat = stats			# Statistics file
+
+    print ("\nsensfunc (standards='",st,"',sensitivity='",sns,"',",>>s2)
+    print ("\tlogfile='",stat,"',apertures='',ignoreaps=no,",>>s2)
+    print ("\tfunction='",function,"',order=",order,",",>>s2)
+    print ("\tinteract=no,",ttl,>>s2)
+
+    print ("")
+}
+
+
+if (bswitch) {				# BSWITCH?
+    print ("\n#---BSWITCH---\n")
+    print ("\n#---BSWITCH---\n", >>s2)
+
+    # Save starting output record number
+    in = out
+    out = "b" // ot
+    wt = weight				# Weighting?
+    if (stat == "")
+	stat = stats			# Statistics file
+
+    # Accumulate records
+    print ("next_rec = i", >>s2)
+    b1 = yes
+    while (b1) {
+	records = ""
+        s3 = records
+        if (s3 == "")
+    	    b1 = no
+        else {
+	    print ("j = next_rec", >> s2)
+    	    print ("bswitch (input='",in,"',output='",out,"',",>>s2)
+    	    print ("\trecords='",s3,"',stats='",stat,"',",>>s2)
+    	    print ("\tweighting=",wt,",subset=",subset,",",>>s2)
+	    print ("\tstart_rec=j,",>>s2)
+    	    print ("\twave1=",wave1,",wave2=",wave2,",",ttl,>>s2)
+        }
+    }
+
+    # Output records
+    print ("j = next_rec", >>s2)
+    s1 = "str (i) // '-' // str(j-1)"
+    print ("s1 = ", s1, >>s2)
+
+    print ("")
+}
+
+if (calibrate) {			# CALIBRATE?
+    print ("\n#---CALIBRATE---\n")
+    print ("\n#---CALIBRATE---\n", >>s2)
+
+    in = out
+    out = "c" // ot
+    if (sns == "")
+        sns = sensitivity		# Sensivity file name
+
+    if (s1 == "") {
+	records = ""
+	s1 = records
+        print ("s1 = '", s1, "'", >>s2)
+    }
+
+    print ("calibrate (input='",in,"',output='",out,"',records=s1,",>>s2)
+    print ("\tignoreaps=no,",>>s2)
+    print ("\textinct=no,flux=yes,",>>s2)
+    print ("\tsensitivity='",sns,"',fnu=",fnu,",",ttl,>>s2)
+
+    print ("")
+}
+
+if (addsets) {				# ADDSETS?
+    print ("\n#---ADDSETS---\n")
+    print ("\n#---ADDSETS---\n", >>s2)
+
+    in = out
+    out = "a" // ot
+    if (s1 == "") {
+	records = ""
+	s1 = records
+        print ("s1 = '", s1, "'", >>s2)
+    }
+
+    print ("addsets (input='",in,"',output='",out,"',records=s1,",>>s2)
+    print ("\tstart_rec=i,subset=2,",ttl,>>s2)
+}
+
+# All done with generator.  Ask whether to execute it.
+print ("File generation complete - filename=",s2)
+if (proceed == no)			# Execute batch file?
+    bye
+}
+
+# Execute generated batch file
+process &