Changeset 115 for trunk/examples

Timestamp:

12/01/04 12:30:51 (19 years ago)

Author:

mar637

Message:

Aligned with asap v0.1a

File:

• trunk/examples/test.py (modified) (6 diffs)

trunk/examples/test.py

@@ -3 +3 @@
 #import asap module
 from asap import *
-# create a reader
-r = sdreader()
-# open an RPFITS file
-r.open('/u/mmarquar/zorro/singledish/data/2001-09-01_0332_P363.rpf') 
+# create a readersa and open an RPFITS file
+r = reader('/u/mmarquar/zorro/singledish/data/2001-09-01_0332_P363.rpf')
 # create a vector with numbers [0..109]
 integrations = range(110)
-r.read(integrations)
 # get the data out of  the reader
-scans = r.getdata()
+scans = r.read(integrations)
 # close the reader
 r = None
-# Test sdwriter.
-print 'Writing ASAP table to disk as /tmp/test.tbl...'
-scans.makepersistent('/tmp/test.tbl')
-print "removing /tmp/test.tbl ..."
-os.remove('/tmp/test_SDWriter.sdfits')
-print 'Begin sdwriter tests...'
-# Create an MS2 writer.
-w = sdwriter('MS2')
-# Change to SDFITS output format (the default).
-w.setformat()
-# Write out the spectra.
-w.write(scans, '/tmp/test_SDWriter.sdfits')
+print 'Begin export test...'
+# Write out the spectra in SDFITS format
+scans.save('/tmp/test_SDWriter.sdfits','SDFITS')
 # clean up
 print "removing test_SDWriter.sdfits ..."
@@ -31 +19 @@
 # print a short summary of the data
 scans.summary()
-# get the scan with the number '1'
-scan = scans.getscan(0)
-# get the scan with the number '1'
-ref = scans.getscan(1)
+# get the scan with the number '0'
+scan = scans.get_scan(0)
+# get the scan with the name 'ref_R'
+ref = scans.get_scan('ref_R')
 # close the data table
 scans = None
 # open the math server
 # average the "on" scan
-scanav = average(scan)
+scanav = average_time(scan)
 # get rid of the original scan
 scan = None
@@ -45 +33 @@
 scanav.summary()
 # average the "off" scan
-refav = average(ref)
+refav = average_time(ref)
 # get rid of the original scan
 ref = None
@@ -53 +41 @@
 quot = quotient(scanav,refav)
 # set the cursor to polarisation 0
-quot.setpol(0)
+quot.set_selection(thepol=0)
 # get the spectrum for polarisation 0
 v0 = quot.getspectrum()
@@ -59 +47 @@
 print v0[0:10]
 # set the cursor to polarisation 1
-quot.setpol(1)
+quot.set_selection(thepol=1)
 # get the spectrum for polarisation 1
 v1 = quot.getspectrum()
@@ -65 +53 @@
 print v1[0:10]
 # write it to disk for further use
-quot.makepersistent('/tmp/myfirstquotient.tbl')
+quot.save('/tmp/myfirstquotient.asap')
 # cleanup
-print "removing /tmp/myfirstquotient.tbl ..."
-os.system('rm -rf /tmp/myfirstquotient.tbl')
+print "removing /tmp/myfirstquotient.asap ..."
+os.system('rm -rf /tmp/myfirstquotient.asap')
 print "Test successful."