Index: /trunk/psastro/src/psastroWCS.c
===================================================================
--- /trunk/psastro/src/psastroWCS.c	(revision 9642)
+++ /trunk/psastro/src/psastroWCS.c	(revision 9643)
@@ -178,4 +178,7 @@
 	chip->fromFPA = p_psPlaneTransformLinearInvert(toFPA);
 
+	while (fpa->toSky->R <        0) fpa->toSky->R += 2.0*M_PI;
+	while (fpa->toSky->R > 2.0*M_PI) fpa->toSky->R -= 2.0*M_PI;
+
 	// remove the correction to the common plate scale
 	if (isMosaic) {
@@ -186,18 +189,17 @@
 	}
 
-	fprintf (stderr, "toFPA: %f %f  (%f,%f),(%f,%f)\n", 
+	psTrace ("psastro", 5, "toFPA: %f %f  (%f,%f),(%f,%f)\n", 
 		 chip->toFPA->x->coeff[0][0], chip->toFPA->y->coeff[0][0], 
 		 chip->toFPA->x->coeff[1][0], chip->toFPA->x->coeff[0][1], 
 		 chip->toFPA->y->coeff[1][0], chip->toFPA->y->coeff[0][1]);
 
-	fprintf (stderr, "frFPA: %f %f  (%f,%f),(%f,%f)\n", 
+	psTrace ("psastro", 5, "frFPA: %f %f  (%f,%f),(%f,%f)\n", 
 		 chip->fromFPA->x->coeff[0][0], chip->fromFPA->y->coeff[0][0], 
 		 chip->fromFPA->x->coeff[1][0], chip->fromFPA->x->coeff[0][1], 
 		 chip->fromFPA->y->coeff[1][0], chip->fromFPA->y->coeff[0][1]);
 
-	fprintf (stderr, "field: %f,%f  %f,%f\n", 
-		 DEG_RAD*fpa->toSky->R, DEG_RAD*fpa->toSky->D, 
-		 3600*DEG_RAD*fpa->toSky->Xs, 3600*DEG_RAD*fpa->toSky->Ys);
-
+	psLogMsg ("psastro", 3, "field center: %f,%f, plate scale: %f,%f (arcsec/pixel)\n", 
+		  DEG_RAD*fpa->toSky->R, DEG_RAD*fpa->toSky->D, 
+		  3600*DEG_RAD*fpa->toSky->Xs, 3600*DEG_RAD*fpa->toSky->Ys);
     }
     return true;