Program: CREATE
Purpose: Creates a SET (and subsets)
Category: UTILITY
File: create.shl
Author: K.G. Begeman
Keywords:
OUTSET= Name of set to create.
If the set already exists, the user is asked whether
the set should be deleted with DELETE=
DELETE= Delete this set ? [N]
This keyword is asked only when the set you want to
create already exists.
If you do not want to delete this set, you are asked
to give another set to create.
INSTRUME= Give type of instrument.
This item will be put in the header and will determine
the coordinate system
CTYPE1= Enter axis name of first axis [stop axis input]
For the n_th axis the keyword will be CTYPEn=, like
all the following keywords
NAXIS1= Size of axis in grids.
CRPIX1= Reference pixel of the axis.
CUNIT1= Physical units of axis [depends on CTYPE]
CRVAL1= Value at reference pixel in CUNITs.
CDELT1= Grid separation in CUNISs.
CROTA1= Rotation angle of axis in degrees [0.0].
This keyword is only asked in case of a DEC axis.
DUNIT1= Secondary units of axis [depends on CTYPE]
This keyword is only asked in case of a FREQ axis.
DRVAL1= Secondary reference value in DUNITs.
This keyword is only asked in case of a FREQ axis.
FREQ0= Give rest frequency in Hz.
This keyword is only asked in case of a FREQ axis.
BUNIT= Give data units: [W.U.]
FUNCTION= F(RA,...)= [empty]
Give an expression using the syntax in the
description. The number of parameters cannot
exceed the number of created axes. The names of
a parameter is the name of an axes. FUNCTION=
will prompt you with all available parameter names.
If carriage return is pressed, the image will be an
empty set.
Description: Create a set from scratch and fill it with values.
The values are calculated from a user given
expression. The indices of the parameters in that
expression correspond to the indices of the created
axes of the set and each grid is a parameter value.
Illegal operations, like dividing by zero or the
logarithm of a negative value, will cause the result
to be set to BLANK.
For complicated functions the user is advised to
make use of 'recall' files for easy input of the
function definition.
Available mathematics:
1) functions; syntax ff(..); where ff is one of the
following available functions:
abs(x) absolute value of x
sqrt(x) square root of x
sin(x) sine of x
asin(x) inverse sine of x
cos(x) cosine of x
acos(x) inverse cosine of x
tan(x) tangent of x
atan(x) inverse tan of x
atan2(x,y) inverse tan (mod 2pi) x = sin, y = cos
exp(x) exponential of x
ln(x) natural log of x
log(x) log (base 10) of x
sinh(x) hyperbolic sine of x
cosh(x) hyperbolic cosine of x
tanh(x) hyperbolic tangent of x
rad(x) convert x to radians
deg(x) convert x to degrees
erf(x) error function of x
erfc(x) 1-error function
max(x,y) maximum of x and y
min(x,y) minimum of x and y
sinc(x) sin(x)/x (sinc function)
sign(x) sign of x (-1,0,1)
mod(x,y) gives remainder of x/y
int(x) truncates to integer
nint(x) nearest integer
ranu(x,y) generates uniform noise between x and y
rang(x,y) generates gaussian noise with mean x
and dispersion y
ranp(x) generates poisson noise with mean x
ifeq(x,y,a,b) returns a if x == y, else b
ifne(x,y,a,b) returns a if x != y, else b
ifgt(x,y,a,b) returns a if x > y, else b
ifge(x,y,a,b) returns a if x >= y, else b
iflt(x,y,a,b) returns a if x < y, else b
ifle(x,y,a,b) returns a if x <= y, else b
ifblank(x,a,b) returns a if x == BLANK, else b
2) constants; syntax cc; where cc is one of the following
available constants:
PI 3.14159....
C speed of light (SI)
H Planck (SI)
K Boltzmann (SI)
G gravitation (SI)
S Stefan-Boltzman (SI)
M mass of sun (SI)
P parsec (SI)
BLANK Universal undefined value
Note: the Hubble constant is not included.
3) operators; syntax op; where op is one of the following
available operators:
+ addition
- subtraction
* multiplication
/ division
** power
4) parameters; syntax axisname; (CTYPEn=), max 32 names.
Updates: Jul 5, 1990: KGB, document created.
Jul 24, 1990: VOG, function evaluation.
Nov 14, 1991: VOG, $n replaced by axis name.
Oct 20, 1994: KGB, declared USERCHAR and FIEPAR.
Apr 2, 1998: VOG, Use usertext for instrument name.
Apr 27, 2011: JPT, Also check other gds_extend errors.