Index: /trunk/Ohana/src/addstar/src/dump.c
===================================================================
--- /trunk/Ohana/src/addstar/src/dump.c	(revision 3404)
+++ /trunk/Ohana/src/addstar/src/dump.c	(revision 3404)
@@ -0,0 +1,21 @@
+# include "addstar.h"
+
+int dump_rawstars (Stars *stars, int Nstars) {
+
+  int i;
+  FILE *f;
+
+  f = fopen ("stars.dat", "w");
+
+  for (i = 0; i < Nstars; i++) {
+    fprintf (f, "%4d  %10.6f %10.6f  %8.2f %8.2f  %7.2f %7.2f\n", 
+	     i, 
+	     stars[i].R, stars[i].D,
+	     stars[i].X, stars[i].Y,
+	     stars[i].M, stars[i].dM);
+  }
+
+  fclose (f);
+  exit (1);
+}
+
Index: /trunk/Ohana/src/mosastro/src/mosastro.c
===================================================================
--- /trunk/Ohana/src/mosastro/src/mosastro.c	(revision 3403)
+++ /trunk/Ohana/src/mosastro/src/mosastro.c	(revision 3404)
@@ -63,4 +63,12 @@
   Cerror   = GetScatter (&Nastro);
   fprintf (stderr, "scatter (chip) : %f for %d stars\n", Cerror, Nastro);
+
+  ClipOnFP (3.0);
+  Cerror   = GetScatter (&Nastro);
+  fprintf (stderr, "scatter (clip) : %f for %d stars\n", Cerror, Nastro);
+
+  FitChips (ChipOrder);
+  Cerror   = GetScatter (&Nastro);
+  fprintf (stderr, "scatter (chip) : %f for %d stars\n", Cerror, Nastro);
   if ((DUMP != NULL) && !strcmp (DUMP, "fitchips")) dump_match();
 
Index: /trunk/Ohana/src/mosastro/src/testcoords.c
===================================================================
--- /trunk/Ohana/src/mosastro/src/testcoords.c	(revision 3404)
+++ /trunk/Ohana/src/mosastro/src/testcoords.c	(revision 3404)
@@ -0,0 +1,66 @@
+# include "mosastro.h"
+
+void testcoords () {
+
+  int i;
+  FILE *f;
+  double R, D, P, Q, L, M, X, Y;
+  
+  fprintf (stderr, "starting test.dat\n");
+
+  /* version 1: transform with three separate stages */
+  /* X,Y -> L,M -> P,Q -> R,D */
+  f = fopen ("test.1.dat", "w");
+  if (f == (FILE *) NULL) {
+    fprintf (stderr, "ERROR: can't create output file \n");
+    exit (1);
+  }
+  for (i = 0; i < chip[0].Nmatch; i++) {
+    X = chip[0].raw[i].X;
+    Y = chip[0].raw[i].Y;
+    XY_to_RD (&L, &M, X, Y, &chip[0].map);
+    XY_to_RD (&P, &Q, L, M, &field.distort);
+    XY_to_RD (&R, &D, P, Q, &field.project);
+    fprintf (f, "%4d  %10.6f %10.6f  %8.2f %8.2f %8.2f %8.2f  %8.2f %8.2f\n", 
+	     i, R, D, P, Q, L, M, X, Y);
+  }
+  fclose (f);
+
+  /* version 2: transform with two separate stages */
+  /* X,Y -> L,M -> R,D */
+  f = fopen ("test.2.dat", "w");
+  if (f == (FILE *) NULL) {
+    fprintf (stderr, "ERROR: can't create output file \n");
+    exit (1);
+  }
+  field_combine ();
+  for (i = 0; i < chip[0].Nmatch; i++) {
+    X = chip[0].raw[i].X;
+    Y = chip[0].raw[i].Y;
+    XY_to_RD (&L, &M, X, Y, &chip[0].map);
+    XY_to_RD (&R, &D, L, M, &field.project);
+    fprintf (f, "%4d  %10.6f %10.6f   %8.2f %8.2f  %8.2f %8.2f\n", 
+	     i, R, D, L, M, X, Y);
+  }
+  fclose (f);
+
+
+  /* version 3: transform with automatic two stages in coordops */
+  /* X,Y -> R,D */
+  f = fopen ("test.3.dat", "w");
+  if (f == (FILE *) NULL) {
+    fprintf (stderr, "ERROR: can't create output file \n");
+    exit (1);
+  }
+  RegisterMosaic (&field.project);
+  strcpy (chip[0].map.ctype, "DEC--WRP");
+  for (i = 0; i < chip[0].Nmatch; i++) {
+    X = chip[0].raw[i].X;
+    Y = chip[0].raw[i].Y;
+    XY_to_RD (&R, &D, X, Y, &chip[0].map);
+    fprintf (f, "%4d  %10.6f %10.6f   %8.2f %8.2f\n", 
+	     i, R, D, X, Y);
+  }
+  fclose (f);
+  exit (1);
+}
Index: /trunk/Ohana/src/opihi/cmd.astro/drizzle.c
===================================================================
--- /trunk/Ohana/src/opihi/cmd.astro/drizzle.c	(revision 3403)
+++ /trunk/Ohana/src/opihi/cmd.astro/drizzle.c	(revision 3404)
@@ -40,4 +40,10 @@
   GetCoords (&coords_in, &in[0].header);
   GetCoords (&coords_out, &out[0].header);
+
+  /* for the moment, disable WRP / DIS */
+  if (!strcmp(&coords_in.ctype[4], "-WRP") || !strcmp(&coords_out.ctype[4], "-WRP")) {
+    fprintf (stderr, "WRP mode not implemented for astrom\n");
+    return (FALSE);
+  }
   
   scale_in = sqrt(fabs(coords_in.cdelt1*coords_in.cdelt2*(coords_in.pc1_1*coords_in.pc2_2 - coords_in.pc1_2*coords_in.pc2_1)));
Index: /trunk/Ohana/src/opihi/cmd.astro/transform.c
===================================================================
--- /trunk/Ohana/src/opihi/cmd.astro/transform.c	(revision 3403)
+++ /trunk/Ohana/src/opihi/cmd.astro/transform.c	(revision 3404)
@@ -23,4 +23,10 @@
   GetCoords (&coords_in, &in[0].header);
   GetCoords (&coords_out, &out[0].header);
+
+  /* for the moment, disable WRP / DIS */
+  if (!strcmp(&coords_in.ctype[4], "-WRP") || !strcmp(&coords_out.ctype[4], "-WRP")) {
+    fprintf (stderr, "WRP mode not implemented for astrom\n");
+    return (FALSE);
+  }
   
   scale_in = sqrt(fabs(coords_in.cdelt1*coords_in.cdelt2*(coords_in.pc1_1*coords_in.pc2_2 - coords_in.pc1_2*coords_in.pc2_1)));