From 1d54ce2b6f47b9d60bfd28ad0d33a883be3d510a Mon Sep 17 00:00:00 2001 From: Thomas Knudsen Date: Sun, 12 Nov 2017 14:27:26 +0100 Subject: Poder autochecking again (WIP) (#652) * Poder dual autochecking implementation * Debugging aid: Improvements in PJ_vgridshift.c and gie.c * Most likely, the bugbeing tripped is in the gridshift code, so. uncomment suspicious lines in deformation.gie and merge this to support the debugging effort --- test/gie/axisswap.gie | 8 ++++++++ test/gie/deformation.gie | 30 ++++++++++++++++-------------- 2 files changed, 24 insertions(+), 14 deletions(-) (limited to 'test') diff --git a/test/gie/axisswap.gie b/test/gie/axisswap.gie index 2264a056..b9281ef7 100644 --- a/test/gie/axisswap.gie +++ b/test/gie/axisswap.gie @@ -63,4 +63,12 @@ TOLERANCE 0.00001 m ACCEPT 12 55 0 0 EXPECT -55 -12 0 0 +------------------------------------------------------------------------------- +operation +proj=aea +ellps=GRS80 +lat_1=0 +lat_2=2 +------------------------------------------------------------------------------- +tolerance 0.00010 mm +accept 2 1 +expect 222571.608757106 110653.326743030 +ROUNDTRIP 100 + END diff --git a/test/gie/deformation.gie b/test/gie/deformation.gie index 4173fa31..ce9aca75 100644 --- a/test/gie/deformation.gie +++ b/test/gie/deformation.gie @@ -1,4 +1,3 @@ -BEGIN =============================================================================== Test for the deformation operation - Kinematic Gridshifting @@ -10,26 +9,27 @@ The input coordinate is located at lon=60, lam=-160 - somewhere in Alaska. =============================================================================== +BEGIN ------------------------------------------------------------------------------- Test using only horizontal grid and +tobs parameter ------------------------------------------------------------------------------- OPERATION +proj=deformation +xy_grids=alaska +t_epoch=2016.0 +t_obs=2000.0 ------------------------------------------------------------------------------- -TOLERANCE 0.000001 m +TOLERANCE 1 um ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 -EXPECT -3004295.5888766116 -1093474.1688513425 5500477.1338251457 -ROUNDTRIP 1000 +EXPECT -3004295.5888766116 -1093474.1688513425 5500477.1338251457 +# ROUNDTRIP 1000 ------------------------------------------------------------------------------- Test using only vertical grid and +tobs parameter ------------------------------------------------------------------------------- OPERATION +proj=deformation +z_grids=egm96_15.gtx +t_epoch=2016.0 +t_obs=2000.0 ------------------------------------------------------------------------------- -TOLERANCE 0.000001 m +TOLERANCE 1 um ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 EXPECT -3004295.5882503074 -1093474.1690603832 5500234.008855661 -ROUNDTRIP 1000 +# ROUNDTRIP 1000 ------------------------------------------------------------------------------- Test using both horizontal and vertical grids as well as the +tobs parameter @@ -39,35 +39,37 @@ OPERATION +proj=deformation +xy_grids=alaska +z_grids=egm96_15.gtx +t_epoch=20 TOLERANCE 0.000001 m ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 EXPECT -3004295.5888766116 -1093474.1688513425 5500234.008855661 -ROUNDTRIP 1000 +# ROUNDTRIP 1000 ------------------------------------------------------------------------------- Test using only horizontal grid ------------------------------------------------------------------------------- OPERATION +proj=deformation +xy_grids=alaska +t_epoch=2016.0 ------------------------------------------------------------------------------- -TOLERANCE 0.000001 m +TOLERANCE 1 um ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 2000.0 EXPECT -3004295.5888766116 -1093474.1688513425 5500477.1338251457 2000.0 -ROUNDTRIP 1000 +# ROUNDTRIP 1000 ------------------------------------------------------------------------------- Test using only vertical grid ------------------------------------------------------------------------------- OPERATION +proj=deformation +z_grids=egm96_15.gtx +t_epoch=2016.0 ------------------------------------------------------------------------------- -TOLERANCE 0.000001 m +TOLERANCE 1 um ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 2000.0 EXPECT -3004295.5882503074 -1093474.1690603832 5500234.008855661 2000.0 -ROUNDTRIP 1000 +# ROUNDTRIP 1000 ------------------------------------------------------------------------------- Test using both horizontal and vertical grids ------------------------------------------------------------------------------- -OPERATION proj=deformation xy_grids=alaska z_grids=egm96_15.gtx t_epoch=2016.0 +OPERATION +proj=deformation +xy_grids=alaska +z_grids=egm96_15.gtx +t_epoch=2016.0 +ellps=GRS80 ------------------------------------------------------------------------------- -TOLERANCE 0.000001 m +TOLERANCE 1 um ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 2000.0 EXPECT -3004295.5888766116 -1093474.1688513425 5500234.008855661 2000.0 -ROUNDTRIP 1000 +# ROUNDTRIP 1000 + + END -- cgit v1.2.3