DVO_Challenge – Pan-STARRS IPP

 ============================================================
 DVO tutorial with Gene Magnier
 03/10/08
 ============================================================

 http://kiawe.ifa.hawaii.edu/Elixir/Ohana/DVO/

 **ra values are in DEGREES**

 ssh host_with_dvo_installed

 Make sure you source the psconfig script that's appropriate for your shell,
 and set up your environment for using ipp, e.g.:
   psconfig ipp-2.5

 Mucking around with dvo

   > script dvo.logfile      make recording of session
   > dvo

   # Set up variables with catdir names
   $catdir_2mass =  /local/path/to/catdir.2mass                  2MASS catalog
   $catdir_merge =  /local/path/to/catdir.2mass/catdir.merge     2MASS+USNOB catalog

   $a = 0      set scalar
   echo $a     
   $a = 5*6    do math

   echo {5*6}           things in {} mean math operation

   create x 0 100 0.1   make an array full of values
   echo x[]
   echo x[0]
   echo x[999]

   set y = x^2          do an array operation
   set b = $a^2         need to handle scalars with "$", but not arrays  

   create x 0 100 0.1
   set x = y
   limits x y           set limits for plotting box
   box
   set y = x^2
   plot x y
   png -name plot.png
   ps -name plot.ps

   mcreate foo 100 100    make array of 100x100 ("m" = matrix)
   set xr = xramp(foo)    'xramp' is a function, makes 2d image as a slope in x-axis
   set yr = xramp(foo)    'xramp' is a function, makes 2d image as a slope in y-axis
   set zr = xr^2 + yr^2   
   stats xr               give image statistics

   tv xr 0 100            display the results with a range (must give a range)
   clear -image
   clear -graph

   buffers  		 list all arrays
   vectors		 list all images
   ??			 list all scalars

 Messing around with databases

    catdir $catdir_2mass     2MASS catalog
    catdir $catdir_merge     2MASS+USNOB catalog

    region 0 0 0.5       set center of region to (RA,Dec)=(0,0) with 0.5 deg half-box size
    region 0 0 0.5 sin   sine projection
    region 0 0 0.5 tan   tan projection
    region 0 0 90 ait ; cgrid     Aitoff projection of the whole sky & show celestial grid

    region 0 0 0.5 tan   tan projection
    pmeasure -all -m 10 15 -pt 7 -c blue    plot all mag=[10,15] detections (J, H, and K)
                                            using blue point symbol

    pmeasure -all -m 10 15 -pt 7 -c blue -photcode 2MASS_J   plot only J-band (bigger symbol = brighter)
    pmeasure -all -m 10 15 -pt 7 -c red -photcode 2MASS_K

    mextract -region ra, dec, photcode, mag     extract info from define plotting region
    region                                      show current value of region
    vectors                                     show results
    cplot ra dec        plot results on celestial sky

    pcat -c red         overlay box sizes for sub-regions

    skyregion 0 1 0 1   set search region (alternative to using 'region' command)
    mextract RA, DEC, MAG   now don't have to use "-region" flag
                            use capitals to define different input arrays
                            (but catalog fields are case-insensitive)

    czplot ra dec mag 10 15    plot {ra,dec} results with mag value encoded in symbol size
                               must be given the min & max value

    mextract -region RA, DEC, MAG where (MAG < 12.0)      restricted extraction
    mextract -region RA, DEC, MAG where (MAG < 12.0) && (MAG > 2.0) && (RA < 0.25)
    cplot RA DEC -c red

    delete ra dec photcode mag RA DEC MAG
    avextract -region ra dec 2MASS_J 2MASS_K    
    avextract -region ra dec 2MASS_J 2MASS_K 2MASS_J:err
    vectors           (note that resulting array names can't start with a number, so "v_" appended)
    clear
    czplot ra dec v_2MASS_J 15 10 -c blue -pt 7
    pmeasure -all -m 10 15 -pt 7 -c red -photcode 2MASS_J

    avextract ra, dec, 2MASS_J, glon, glat where (abs(glat) > 10)     extract based on galactic latitude

    limits -n 1 -0.5 3.5 17.0 8.0; clear; box
    plot {v_2MASS_J - v_2MASS_K} v_2MASS_K
         NOTE: now that window 1 is active, this will be the default range for "-region"
    dev -n 0       change back to device 0 
    avextract -region ra, dec, 2MASS_J, 2MASS_K where (2MASS_J-2MASS_K > 0) && (2MASS_J-2MASS_K < 0.5)
    plot -n 1 -c red {v_2MASS_J - v_2MASS_K} v_2MASS_K

    cursor     poke plotting window, click on mouse to return coordinates
               "q" to quit

    write -f "%10.6f %10.6f %7.3f" test.dat ra dec      dump to file
    !more test.dat         "!" is shell escape
    exec more test.dat     also shell escape

    list data -x "head test.dat"                        input files

    subset r = ra if (dec < 0.2)
    subset d = dec if (dec < 0.2)
    subset m = mag if (dec < 0.2)

    set bool = (dec < 0.2)       create boolean variable for selection
    setset r = ra if bool

    macro foo                    create a macro
      echo $1
      for i 0 10
      echo $i
      end
      end
    foo 5
    foo sdf sdfsd sdfdsf

    input filename               import file of commands
      in the file, can use "#" as comment character (though this doesn't work on the command line directly)
      cannot use "!" in the file as shell escape, must use "exec"

    vstat              get statistics for a vector (no "s" on the command name)
                         results are stored in variables $MEAN, $MIN, etc.
    stats              get statistics for a scalar

    sort ra dec mag    will re-sort the variables using first variable as the sorting field

 ----------------------------------------------------------------------

 the 'date' command is what you want here:

 dvo: date -h
 USAGE: date [-var variable] [-seconds] [-reftime seconds]

 * by itself, just echoes the date & time
 * to set a variable with the date, add -var name, where 'name' is the
 name of your variable.  by default, this is a human readable date, eg
 Wed Mar 12 08:30:11 2008
 * to calculate time differences, add the -seconds flag, and the value
 will be in seconds since 1970 Jan 01 00:00:00.  NOTE: There is a bug
 in that the variable is represented as a float, and the value is now >
 1.2e9 seconds.  The difference is thus not very accurate.  To get
 around this, use the -reftime feature to report the seconds since a
 more recent date.  A value of something like 1205340000 works at the
 moment.

 examples:

 dvo: date
 Wed Mar 12 08:30:11 2008

 dvo: date -var foo -seconds
 dvo: echo $foo
 1205346664

 dvo: date -var foo -seconds -reftime 1205340000
 dvo: echo $foo
 6700
 dvo: date -var bar -seconds -reftime 1205340000
 dvo: echo $bar
 6708
 dvo: echo {$bar - $foo}
 8

Context Navigation

FAQ / Suggestions for Changes

Questions

Change requests/suggestions

Examples

DVO tutorial

Sample dvo script

DVO tutorial (old)

Download in other formats: