[GRASS-SVN] r45423 - grass/branches/develbranch_6/ps/ps.map

[svn_grass at osgeo.org]

Author: hamish
Date: 2011-02-17 22:48:55 -0800 (Thu, 17 Feb 2011)
New Revision: 45423

Modified:
   grass/branches/develbranch_6/ps/ps.map/do_geogrid.c
Log:
comment on common semi-bug

Modified: grass/branches/develbranch_6/ps/ps.map/do_geogrid.c
===================================================================
--- grass/branches/develbranch_6/ps/ps.map/do_geogrid.c	2011-02-17 19:45:09 UTC (rev 45422)
+++ grass/branches/develbranch_6/ps/ps.map/do_geogrid.c	2011-02-18 06:48:55 UTC (rev 45423)
@@ -373,10 +373,13 @@
 }
 
 
-/******************************************************
+/********************************************************************
  * check projected coords to make sure they do not 
- * go outside region -- if so re-project 
-********************************************************/
+ *   go outside region -- if so re-project 
+ * This gets a bit wobbly if the projection does not reverse-project
+ *   cleanly, e.g. if the eastern edge of a UTM region is many
+ *   degrees east of the end of the zone where things get distorted.
+ *******************************************************************/
 void check_coords(double e, double n, double *lon, double *lat, int par)
 {
     double x, y;
@@ -409,8 +412,6 @@
 	    G_fatal_error(_("Error in pj_do_proj"));
 
 	if (par == 1) {
-/* FIXME: something isn't quite right here- stuck at NE corner's longitude */
-	    /* lines of latitude -- const. northing */
 	    /* convert correct UTM to ll */
 	    if (pj_do_proj(&x, &y, &info_out, &info_in) < 0)
 		G_fatal_error(_("Error in pj_do_proj"));