%N colin %F $MIR/src/subs/co.for
%D Generate a linearized approximation to a nonlinear system.
%P rjs
%: coordinates
%B
      subroutine coLin(lu,in,x1,n,ctypel,crpixl,crvall,cdeltl)

      integer lu,n
      character in*(*),ctypel(n)*(*)
      double precision x1(*),crpixl(n),crvall(n),cdeltl(n)
  ---------------------------------------------------------------------
  Generate a linearised approximation to a non-linear coordinate
  system.

  Input:
    lu         Handle of the input coordinate system.
    x1         Input coordinate.  The linearised approximation
               is performed around this point.
    in         String describing the format of x1 coordinate, as
               per coCvt.
    n          Number of axes to output.
  Output:
    ctypel )
    crpixl )   Output coordinate descriptions.
    crvall )
    cdeltl )
