.ec _
.ll 76
.ad b
.in 4
.ti -4
1 Language
.nf
COMBINE       : Combine Images with automatic resampling
EXPLORE       : Make a plot of spectra around a 2-D map
KEPLER        : Optimal line detection for Keplerian disk
SELFCAL       : Perform a self calibration
.fi
  
.ti -4
1 COMBINE
.ti +4
[BUNDLES_\]COMBINE OutCube CODE In1 In2 [Off] [/FACTOR A1 A2]
  [/THRESHOLD T1 T2] [/BLANKING Bval]
  
This is a superset of the MAP__COMBINE command, that automatically 
resamples and reproject the In1 and/or In2 data sets so that they match
spatially (according to In1) and spectrally. 

  See HELP MAP__COMBINE for details on the combination codes,
factors, thresholds and blanking.

.ti -4
2 CODE
.ti +4
[BUNDLES_\]COMBINE OutCube CODE In1 In2 [Off] [/FACTOR A1 A2]
  [/THRESHOLD T1 T2] [/BLANKING Bval]

  CODE is the operation code. Allowed values are
.nf
ADD           OutCube = A1*In1 + A2*In2 + Off  where In1>T1 and In2>T2
DIVIDE        OutCube = A1*In1 / A2*In2 + Off    
MULTIPLY      OutCube = A1*In1 * A2*In2 + Off    
OPACITY       OutCube = -Log( A1*In1 / A2*In2 + Off)  
INDEX         OutCube = Log( A1*In1 / A2*In2) / Log(Nu1/Nu2)  
.fi  
where Nu1 is the Frequency of In1, and Nu2 that of In2.
Off is 0 if not specified.
  
.ti -4
2 /BLANKING
.ti +4
[BUNDLES_\]COMBINE OutCube CODE In1 In2 [/FACTOR A1 A2]
  [/THRESHOLD T1 T2] /BLANKING Bval

Specify the Blanking value to be used in the OutCube. If not
specified, the Blanking from In1 is used instead, and
that of In2 if In1 has no Blanking.

.ti -4
2 /FACTOR 
.ti +4
[BUNDLES_\]COMBINE OutCube CODE In1 In2 /FACTOR A1 A2
  [/THRESHOLD T1 T2] [/BLANKING Bval]

Specify the factors to apply to In1 and In2. Default is 1.0

.ti -4
2 /THRESHOLD 
.ti +4
[BUNDLES_\]COMBINE OutCube CODE In1 In2 [/FACTOR A1 A2]
  /THRESHOLD T1 T] [/BLANKING Bval]

Specify the thresholds above which the computation is
valid for In1 and In2.  For pixels below this threshold,
the OutCube is blanked.

Default is no threshold (-huge(0.)).

.!======================================================================  
.ti -4
1 EXPLORE
.ti +4
[BUNDLES_\]EXPLORE  [DataCube] [/ADD Cx Cy [Box]] [/NOPAUSE]

Create a plot with a 2-D map in the middle, and up to 8 spectra
around it. 

DataCube is the name of an internal 3-D buffer, SIC variable or
data file. The 2-D map will be either a channel of this 3-D data set,
or a velocity integrated map obtained from it. It will displayed
at the center of the plot. If not present, argument DataCube keeps its
previous value.

Spectra can be displayed around, with arrows showing their position
in the map, either interactively (when the /NOPAUSE option is not
present), or manually through the /ADD option. In interactive mode,
the user can press S to extract and display the spectrum from
any place in the 2-D map.

EXPLORE ?  will list the control variables for this command. The
usual SIZE and CENTER control the 2-D map size, RANGE controls the
velocity range, and SCALE the spectral flux range.

.ti -4
2 /ADD
.ti +4
[BUNDLES_\]EXPLORE  /ADD Cx Cy [Box]

Add the spectrum coming from positin Cx Cy (in arcsec) in the specified
Box number. The Box number indicates the placement as in a numeric
keypad layout. Box 5 (the center) is reserved for the 2-D map.

If Box is not present, the spectrum is placed in the most plausible
position according to the coordinates (Cx,Cy).

.ti -4
2 /NOPAUSE
.ti +4
[BUNDLES_\]EXPLORE  [DataCube] /NOPAUSE

Just display the 2-D map and any already existing spectra, but do
not call the cursor.

.ti -4
1 KEPLER
.ti +4
[BUNDLES_\]KEPLER [?|DataCube] 

Re-align spectra from a Keplerian disk according to the projected
rotation velocity at any point in the disk, and computes the
combined integrated spectrum and brightness radial distribution. 
The projection assumes a thin disk (no flaring).

The command is controlled by a set of KEPLER__* variables, whose
values can be verified by KEPLER ?  and modified by the user.
In addition to the built-in KEPLER__* variables, KEPLER__RMAX
indicates the radius

DataCube is the name of the data cube to be used. It can a SIC 3-D
variable or the name of a GILDAS (or FITS) data file. The default is
CLEAN.

.ti -4
2 Variables:
.ti +4

  Most input variables are pre-defined by the code. The exceptions
are KEPLER__VDISK and KEPLER__RMAX, which can be defined by the user.
If these are not present, they are derived from the data (for the
disk velocity) or from KEPLER__ROUT (for KEPLER__RMAX).

.ti -4
2 KEPLER__X0
.ti +4
  X offset of disk center (in arcsec)
.ti -4
2 KEPLER__Y0
.ti +4
  Y offset of disk center (in arcsec)
.ti -4
2 KEPLER__ROTA
.ti +4
  Position angle of the projection of the disk rotation axis 
(in degree, East from North)
.ti -4
2 KEPLER__INCLI
.ti +4
  Disk inclination (in degree)
.ti -4
2 KEPLER__DIST
.ti +4
  Disk distance (in pc).
.ti -4
2 KEPLER__VMASS
.ti +4
  Disk rotation velocity at 100 au in km/s. Should be equal
  to 2.98 * sqrt(M*/Msun) from Kepler laws.
.ti -4
2 KEPLER__RINT
.ti +4
  Inner radius in au.
.ti -4
2 KEPLER__ROUT
.ti +4
  Outer radius in au. The integrated profile is computed by summing up
  from KEPLER__RINT to KEPLER__ROUT.
.ti -4
2 KEPLER__STEP
.ti +4
  Sampling step in au. Must be compatible with the angular resolution.
  Note that the spectral resolution allows some gain in spatial
  resolution compared to the beam size.
.ti -4
2 KEPLER__RMAX
.ti +4
  KEPLER__RMAX is an optional variable that specified up to which
maximum radius (in au) the radial profile is computed.
By default (i.e. if the variable has not been created by the
user), the maximum radius is taken as KEPLER__ROUT + 5*KEPLER__STEP.
.ti -4
2 KEPLER__THETA
.ti +4
  Maximum angle (in degree) from the disk projected major axis beyond 
which the spectra are not considered in the averaging process.  

This angle is needed because  the projected rotation velocity is degenerate
to zero along the minor axis, which results in a lower effective
angular resolution in this direction. This effect depends on the
disk inclination: a larger angle can be used for less inclined disks.

The default is 60 degree. 
.ti -4
2 KEPLER__VDISK
.ti +4
  KEPLER__VDISK is an optional variable that defines the 
velocity (in km/s) at which the radial profile will be computed. If this
variable is not defined by the user, the velocity will be taken
from the peak of the integrated, velocity aligned, spectrum.

.ti -4
2 Results:
.ti +4
  The KEPLER command produces 3 different SIC Tables: KEPLER__SPECTRUM,
KEPLER__SPECTRA, and KEPLER__PROFILE, available as SIC variables
for further writing or plotting.

.ti -4
2 KEPLER__SPECTRUM
.ti +4
  Velocity-aligned, disk-integrated spectrum, over the region defined
by KEPLER__RINT, KEPLER__ROUT and KEPLER__THETA. This is a 3 column
table containing the velocities (km/s) in column 1, the flux in
column 2, and an estimate of the flux error in column 3.
  
.ti -4
2 KEPLER__SPECTRA
.ti +4
  Velocity-aligned spectra as a function of radius. This is a Nrad column
table containing the spectra (averaged brightness temperature) for
each of the Nrad radii defined by the sampling KEPLER__RINT, KEPLER__RMAX
and KEPLER__STEP.
  The velocities are defined in KEPLER__SPECTRUM[1].

.ti -4
2 KEPLER__PROFILE
.ti +4
  Radial profile of the peak brightness temperature. The peak brightness
temperature is that appearing at KEPLER__VDISK if defined, or at the
velocity defined by the maximum of the integrated spectrum. This is
a 3 column  table containing the radii (in au), the brightness (in K)
and an estimate of its error.

.!======================================================================  
  
.ti -4
1 SELFCAL

.so imager-help-selfcal.nroff
 
.ti -4
1 ENDOFHELP


