4.5.5 Some examples for command procedures

In order to provide a starting point for your data reduction, we will list in the following a simple command procedure, and will explain thereafter how you can improve your reduction by adding more and more features. We recommend that you enter all commands by hand at least once, in order to learn about MOPSI. Note that values given in this example are only ...examples - you have to adjust them to your data.

Let us assume you have created a file called examlist.LIST with all maps of one source and all skydips surrounding these map measurements. It will probably look like:

7706th11x.001
7707th11x.001
...
7712th11x.001

A basic reduction file would then be the following. This example is already an enhanced version, since it creates the lists by itself (if the directory structure is okay) and uses some despiking:

!  Example procedure for OTF reduction
!  Version: 03/2000, Procedure: restore
!  Source: N2146
!
in-d "/home/smtobs/kont/raw2146"
sys "!test -f raw2146.LIST && /bin/rm raw2146.LIST"
sys "!cd ../raw2146;ls *th11x.001 > ../mopsi/raw2146.LIST;cd ../mopsi"
init in-list
init out-list
list raw2146
select type map skydip
select N2146* SKYDIP
sort list
show list
!
! Reduction Skydips
!
for i 1 to las
  read next skydip
  clea plot
  take rc 20
  skydip
next
sys "test -f n2146.sky && /bin/rm n2146.sky"
write taus n2146.sky
!
! Reduction Maps
!
out-d "/home/smtobs/kont/mopsi/n2146"
sys "!test -f ./n2146/*b2_res.btd && /bin/rm ./n2146/*b2_res.btd"
for i 1 to las
  read next map
  dele rc 20 21
  taus n2146.sky 10000
  corr ext
  base-t 2
  base-r az < -60 > 60 /pa 128 /proj 6 10 40.1 78 22 23.0
  base-s 2
  init spike_pointer
  low -13000
  high 13000
  base-r az < -60 > 60 /pa 128 /proj 6 10 40.1 78 22 23.0
  base-s 2
  mask spike 18000
  base-r az < -60 > 60 /pa 128 /proj 6 10 40.1 78 22 23.0
  base-s 2
  base-l 0
  store wei rms2 in base-r
  db_f -98 !
  restore -98 /nofilter !
  write .btd -th11x +b2_res
next
!
! Rotation to Az-El and Summation
!
in-d "/home/smtobs/kont/mopsi/n2146"
out-d "/home/smtobs/kont/mopsi"
read 7707b2_res.btd
convert /comm
  ext 700 700
  blc 350 -350
  proj 6 10 40.1 78 22 23.0
  run
read 7708b2_res.btd
  run
read 7710b2_res.btd
  run
  quit
sys "test -f n2146_res.nod && /bin/rm n2146_res.nod"
file out n2146_res.nod
write
close out
clea plo
plot
sys "test -f n2146.ps && /bin/rm n2146.ps"
hardc n2146.ps /dev eps col
!

This procedure will produce the following files:
1) A table of skydips, called n2146.sky
2) restored .btd data of your maps, in files named e.g. 7708b2_res.btd
3) a restored, converted, and co-added output-map in the file n2146_res.nod
4) a plot of the final map in n2146.ps

If you use the Shift-And-Add method, it looks very similar until the double-beam filtering (db_f -98 !), except the file names which should somehow indicate that the written .btd files are not restored. We give the second half of an appropriate file below:

  ...
  db_f -98 !
  write .btd -th11x +b2
next
!
! Rotation to Az-El and Summation
!
in-d "/home/smtobs/kont/mopsi/n2146"
out-d "/home/smtobs/kont/mopsi"
read 7707b2.btd
convert /comm
  ext 700 700
  blc 350 -350
  proj 6 10 40.1 78 22 23.0
  mask az > 10 /pa 128 /proj 6 10 40.1 78 22 23.0
  shift az 49
  run
read 7707b2.btd
  mask az < -10 /pa 128 /proj 6 10 40.1 78 22 23.0
  shift az -49
  mul -1.0
  run
read 7708b2.btd
  mask az > 10 /pa 128 /proj 6 10 40.1 78 22 23.0
  shift az 49
  run
read 7708b2.btd
  mask az < -10 /pa 128 /proj 6 10 40.1 78 22 23.0
  shift az -49
  mul -1
  run
read 7710b2.btd
  mask az > 10 /pa 128 /proj 6 10 40.1 78 22 23.0
  shift az 49
  run
read 7710b2.btd
  mask az < -10 /pa 128 /proj 6 10 40.1 78 22 23.0
  shift az -49
  mul -1
  run
  quit
sys "test -f n2146_saa.nod && /bin/rm n2146_saa.nod"
file out n2146_saa.nod
write
close out
!