Psrav

psrav is script to process new arrival time data, and combine it with existing data.

Mods: make seperate 'date'av.tpo and 'date'av.prf files from each 'date'av.bin file. This should shorten processing time, as only new data will need to be run through psrtime to produce the toas and profiles, rather than the entire psrav.bin. Nov 2000 CAJ

Separate out 42ft data into _2 directory.

Add command a time options to auto file

It goes roughly as follows:-

Usage : psrav <<1 or 2(optional, default 1)>> <pulsarname>  <<date (optional)>> 

• Go to directory /psrdata/timing/$psrname or /psrdata/timing/$psrname/_2
• Check for file auto (exit if not found, or file contains veto)
• Check for correct.dat. (warn if not found, but carry on).
• Split auto file into auto.tmp containing dates and com.tmp containing any lines prefixed by the key word COMMAND or TIME
• Call script sortbinfile which creates a new version of datefile containing the names of all the data files 7 to F??.bin and 7 to F?_2.bin, in ascending order.
• Call script checkdate which removes the names of any .bin files named in auto.tmp from datefile.
• Call script checkav to create a list of dates which have a .bin file, but either nonexistant av.bin or one which has been created before the <date>.bin file was last modified. lso create a file of dates for which a .tpo and/or .prf does not exist or is earlier thena the <date>.bin file.
• .... or make "lists" containing only the date specified by the command line instructions.
• Shift any existing .com and psrav.* files into a temporary directory.
• remove any *av.bin, *.tpo and *.prf files that are going to be updated.
• Run mkpsravcomprogram to create psrav.com file which goes

  		obey psravedit.dat
  
  		set in1f 98a_2.bin
  		set out1f 98a_2av.bin
  		do read in 500 avch a4 write a4
  		avch a4 write a4 clear in
  		close out1f
  

  for all required .bin files.

  (Note: if the auto contains any lines with the key word COMMAND, a symbol autocom will be defined as the contents of all these lines. If the auto contains any lines with the key word TIME the command read in 500 will be replaced with read in xxxx xxx xxx
  where xxx xxx xxx is a legitimate qualifier for read.
  For example psrav.com may read:

  		obey psravedit.dat
                  set autocom = "reject a4 rms 2 ignore 0.1 periods save"
  
  		set in1f 98a_2.bin
  		set out1f 98a_2av.bin
  		do read in 7200 mins -1 autocom avch a4 write a4
  		autocom avch a4 write a4 clear in
  		close out1f
  

  If a file called psravedit.dat is not found, it is created.

  		set autoedit on
  		edit clear
  		ob correct.dat
  		edit add
  		avch in 2 1
  		setvar in nch 1
  		fold in
  		setvar in psrnam 1133+16                                                                 
          
  		changetbin in -400
  		correct in
  		window in start 0.0 1.5 periods
  		setbases
  		fitbases in
  		edit end
  
  

  then the psrav2.com script file is written.

  		 edit set
  		 setvar in pol 1 1 'I' 
  		 setvar in psrnam 1133+16 
                                                                  
  		 window in start 0.0 2.0 periods
  		 setbases fitbases in
  		 edit end
   
  		 symbol wprf="set outch 2 perc write toaprf dummy set outchan 1 perc"
   
  		 set in1f 98a_2av.bin
  		 set out1f 98a_2av.tpo
  		 set out2f 98a_2av.prf
  		 do read in findtemp in getarrtime wprf write line tempo
                   close out1f
                   close out2f
  
  

  for all required *av.bin files.

  Note: If the call to mkpsravcom is followed by the flag freq (case insensitive) then the line

   do read in findtemp in getarrtime wprf write line tempo
  

  will be replaced by

   do read in findtemp in changetbin in -1024 getarrtime freq wprf write line tempo
  

• Run psrprof obeying psrav.com
• Call crunch script to concatenate all the *av.bin files into a new psrav.bin. (Concatenate is brought to you at the special request of Mark Roberts!)
• Run psrprof, obeying psrav2.com to produce a series of *av.tpo and *av.prf files.
• Call crunchtpo script to concatenate all the *av.tpo files into one psrav.tpo, and all the *av.prf files into one psrav.prf file.
• Run mkrdavtpocom to create a file rdavtpo.com

  		 read ephem from psrav
  		 obey 1133+16.sat
  		 read tempo from psrav
  		 sort sats
  		 jplbat
  		 set source 1133+16 
  		 set start -1
  		 set finish -1
  		 write bats to psrav
  
  

• run psrtime obeying rdavtpo.com
• Run mkpsravdat to create psrav.dat

  		 set source 1133+16 
  		 read bats from psrav
  		 sort bats
  		 read ephem from psrav
  		 read del from psrav
  		 set start -1
  		 set finish -1
  

  Chris Jordan
  Aug. 2002