For obtaining positions in a subset the function
GDSPOS is used
in GIPSY. It decodes the user input into grid positions, which are
returned in a double precision array.
Note that the positions returned by GDSPOS should be converted
to integers before using it for calculating coordinate words etc. In
this case the nearest integer must be taken, i.e. with the function NINT.