X-Git-Url: https://gerrit.simantics.org/r/gitweb?a=blobdiff_plain;f=org.simantics.maps.server%2Fserver%2Fdist%2Fshare%2Fproj%2Ftestvarious;fp=org.simantics.maps.server%2Fserver%2Fdist%2Fshare%2Fproj%2Ftestvarious;h=ef1bfbee1f1c6c62df12bd1b2f30286b19860c45;hb=a018cc96644b634cec0dd999298cd374648e2980;hp=0000000000000000000000000000000000000000;hpb=4fb69ee0635d0a03b6a4253d5c1cd920c5eb5366;p=simantics%2Fdistrict.git

diff --git a/org.simantics.maps.server/server/dist/share/proj/testvarious b/org.simantics.maps.server/server/dist/share/proj/testvarious
new file mode 100644
index 00000000..ef1bfbee
--- /dev/null
+++ b/org.simantics.maps.server/server/dist/share/proj/testvarious
@@ -0,0 +1,904 @@
+:
+# Script to do some testing of various transformation that do not depend
+# on datum files.
+#
+#
+NAD_DIR=`dirname $0`
+EXE=$1
+
+usage()
+{
+    echo "Usage: ${0} <path to 'cs2cs' program>"
+    echo
+    exit 1
+}
+
+if test -z "${EXE}"; then
+    EXE=../src/cs2cs
+fi
+
+if test ! -x ${EXE}; then
+    echo "*** ERROR: Can not find '${EXE}' program!"
+    exit 1
+fi
+
+# Would be great to have a universale way of selecting a locale with
+# a decimal separator that is not '.'
+if command locale >/dev/null 2>/dev/null; then
+    if test `locale -a | grep fr_FR.utf8`; then
+        echo "Using locale with comma as decimal separator"
+        export LC_ALL=fr_FR.UTF-8
+        export PROJ_USE_ENV_LOCALE=1
+    fi
+fi
+
+echo "============================================"
+echo "Running ${0} using ${EXE}:"
+echo "============================================"
+
+OUT=tv_out
+#EXE=../src/cs2cs
+#
+echo "doing tests into file ${OUT}, please wait"
+rm -f ${OUT}
+#
+echo "##############################################################" >> ${OUT}
+echo Test raw ellipse to raw ellipse >> ${OUT}
+#
+$EXE +proj=latlong +ellps=clrk66 \
+ +to +proj=latlong +ellps=bessel \
+ -E >>${OUT} <<EOF
+79d58'00.000"W 37d02'00.000"N 0.0
+79d58'00.000"W 36d58'00.000"N 0.0
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo Test NAD27 to raw ellipse >> ${OUT}
+#
+$EXE +proj=latlong +datum=NAD27 \
+ +to +proj=latlong +ellps=bessel \
+ -E >>${OUT} <<EOF
+79d00'00.000"W 35d00'00.000"N 0.0
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo Between two 3parameter approximations on same ellipsoid >> ${OUT}
+#
+$EXE +proj=latlong +ellps=bessel +towgs84=5,0,0 \
+ +to +proj=latlong +ellps=bessel +towgs84=1,0,0 \
+ -E >>${OUT} <<EOF
+0d00'00.000"W 0d00'00.000"N 0.0
+79d00'00.000"W 45d00'00.000"N 0.0
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo 3param to raw ellipsoid on same ellipsoid >> ${OUT}
+#
+$EXE +proj=latlong +ellps=bessel +towgs84=5,0,0 \
+ +to +proj=latlong +ellps=bessel  \
+ -E >>${OUT} <<EOF
+0d00'00.000"W 0d00'00.000"N 0.0
+79d00'00.000"W 45d00'00.000"N 0.0
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo Test simple prime meridian handling. >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84 +pm=greenwich  \
+ +to +proj=latlong +datum=WGS84 +pm=1 \
+ -E >>${OUT} <<EOF
+0d00'00.000"W 0d00'00.000"N 0.0
+79d00'00.000"W 45d00'00.000"N 0.0
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo Test support for the lon_wrap switch. >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84  \
+ +to +proj=latlong +datum=WGS84 +lon_wrap=180 \
+ -E >>${OUT} <<EOF
+1d00'00.000"W 10d00'00.000"N 0.0
+0d00'00.000"W 10d00'00.000"N 0.0
+0d00'00.000"E 10d00'00.000"N 0.0
+1d00'00.000"E 45d00'00.000"N 0.0
+179d00'00.000"E 45d00'00.000"N 0.0
+181d00'00.000"E 45d00'00.000"N 0.0
+350d00'00.000"E 45d00'00.000"N 0.0
+370d00'00.000"E 45d00'00.000"N 0.0
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo Test simple prime meridian handling within a projection. >> ${OUT}
+#
+$EXE +proj=utm +zone=11 +datum=WGS84 +pm=3 \
+ +to +proj=latlong +datum=WGS84 +pm=1w \
+ -E >>${OUT} <<EOF
+500000 3000000
+EOF
+echo "##############################################################" >> ${OUT}
+echo Test geocentric x/y/z generation. >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84  \
+ +to +proj=geocent +datum=WGS84  \
+ -E >>${OUT} <<EOF
+0d00'00.001"W 0d00'00.001"N 0.0
+0d00'00.001"W 0d00'00.001"N 10.0
+79d00'00.000"W 45d00'00.000"N 0.0
+45d00'00.000"W 89d59'59.990"N 0.0
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo Test geocentric x/y/z consumption. >> ${OUT}
+#
+$EXE +proj=geocent +datum=WGS84  \
+ +to +proj=latlong +datum=WGS84  \
+ -E >>${OUT} <<EOF
+6378137.00      -0.00 0.00
+6378147.00      -0.00 0.00
+861996.98       -4434590.01 4487348.41
+0.00    -0.00 6356752.31
+EOF
+#
+#echo "#############################################################" >> ${OUT}
+#echo Test conversion between geocentric latlong and geodetic latlong >> ${OUT}
+#
+# The +geoc flag does not currently work with pj_transform() so this is
+# disabled.
+#
+#$EXE +proj=latlong +datum=WGS84 +geoc \
+# +to +proj=latlong +datum=WGS84  \
+# -E >>${OUT} <<EOF
+#0d00'00.000"W 0d00'00.000"N 0.0
+#79d00'00.000"W 45d00'00.000"N 0.0
+#12d00'00.000"W 45d00'00.000"N 0.0
+#0d00'00.000"W 90d00'00.000"N 0.0
+#EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test stere projection (re: win32 ticket 12)" >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84 \
+ +to +proj=stere +lat_0=90 +lon_0=0 +lat_ts=70 +datum=WGS84 \
+ -E >>${OUT} <<EOF
+105 40
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test stere without lat_ts (#147)" >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84 \
+ +to +proj=stere +lat_0=40 +lon_0=10  +datum=WGS84 \
+ -E >>${OUT} <<EOF
+20 45
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test sts projection (re: ticket 12)" >> ${OUT}
+#
+$EXE +proj=latlong +ellps=WGS84 \
+ +to +proj=kav5 +ellps=WGS84 +units=m \
+ -E >>${OUT} <<EOF
+4.897000 52.371000
+EOF
+$EXE +proj=kav5 +ellps=WGS84 +units=m \
+ +to +proj=latlong +ellps=WGS84 \
+ -E >>${OUT} <<EOF
+383646.088858 5997047.888175
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test RSO Borneo projection (re: ticket 62)" >> ${OUT}
+#
+$EXE +proj=latlong +a=6377298.556 +rf=300.8017 \
+ +to +proj=omerc +a=6377298.556 +rf=300.8017 +lat_0=4 +lonc=115 \
+                 +alpha=53d18\'56.9537 +gamma=53d7\'48.3685  +k_0=0.99984 \
+                 +x_0=590476.87 +y_0=442857.65 \
+ -E >>${OUT} <<EOF
+116d2'11.12630 5d54'19.90183
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test extended transverse mercator (#97)" >> ${OUT}
+#
+$EXE +proj=etmerc +k=0.998 +lon_0=-20 +datum=WGS84 +x_0=10000 +y_0=20000 \
+ +to +proj=latlong +datum=WGS84 \
+ -E >>${OUT} <<EOF
+10000 20000
+500000 2000000
+1000000 2000000
+2000000 2000000
+4000000 2000000
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test extended transverse mercator inverse (#97)" >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84 \
+ +to +proj=etmerc +k=0.998 +lon_0=-20 +datum=WGS84 +x_0=10000 +y_0=20000 \
+ -E >>${OUT} <<EOF
+0dN 0.000
+15d22'16.108"W	17d52'53.478"N 0.000
+10d40'55.532"W	17d42'48.526"N 0.000
+1d32'21.33"W	17d3'47.233"N 0.000
+15d4'42.357"E	14d48'56.372"N 0.000
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test transverse mercator (#97)" >> ${OUT}
+#
+$EXE +proj=tmerc +k=0.998 +lon_0=-20 +datum=WGS84 +x_0=10000 +y_0=20000 \
+ +to +proj=latlong +datum=WGS84 \
+ -E >>${OUT} <<EOF
+10000 20000
+500000 2000000
+1000000 2000000
+2000000 2000000
+4000000 2000000
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test transverse mercator inverse (#97)" >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84 \
+ +to +proj=tmerc +k=0.998 +lon_0=-20 +datum=WGS84 +x_0=10000 +y_0=20000 \
+ -E >>${OUT} <<EOF
+0dN 0.000
+15d22'16.108"W	17d52'53.478"N 0.000
+10d40'55.532"W	17d42'48.526"N 0.000
+1d32'21.33"W	17d3'47.233"N 0.000
+15d4'42.357"E	14d48'56.372"N 0.000
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test robinson projection (#113)" >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84 \
+ +to +proj=robin +datum=WGS84 \
+ -E >>${OUT} <<EOF
+-30 40
+-35 45
+20 40
+EOF
+$EXE +proj=robin +datum=WGS84 \
+ +to +proj=latlong +datum=WGS84 \
+ -E >>${OUT} <<EOF
+-2612095.95     4276351.58 0.00
+-2963455.42     4805073.65 0.00
+1741397.30      4276351.58 0.00
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test hammer projection (pull request #329)" >> ${OUT}
+$EXE +proj=latlong +datum=WGS84 \
+ +to +proj=hammer  +datum=WGS84 \
+ -E >>${OUT} <<EOF
+-30 40
+-35 45
+20 40
+EOF
+$EXE +proj=hammer  +datum=WGS84 \
+ +to +proj=latlong +datum=WGS84 \
+ -E >>${OUT} <<EOF
+-2711575.08	4395506.62 0.00
+-2964412.70	4929091.33 0.00
+1811748.54	4377349.50 0.00
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test healpix forward projection on sphere" >> ${OUT}
+$EXE +proj=latlong +a=1 +lon_0=0 \
+ +to +proj=healpix +a=1 +lon_0=0  -f '%.'5'f' \
+ -E >>${OUT} <<EOF
+0 41.81031
+-90 0
+EOF
+$EXE +proj=latlong +R=5 \
+ +to +proj=healpix +R=5 -f '%.'5'f' \
+ -E >>${OUT} <<EOF
+0 0
+0 41.810314895778596
+0 -41.810314895778596
+90.0 0
+-90.0 0
+-180 0
+-180 90.0
+-180 -90.0
+0 60.0
+0 -60.0
+EOF
+echo "Test healpix forward projection on ellipsoid" >> ${OUT}
+$EXE +proj=latlong +a=1 +lon_0=0 +ellps=WGS84 \
+ +to +proj=healpix +a=1 +lon_0=0 +ellps=WGS84 -f '%.'5'f' \
+ -E >>${OUT} <<EOF
+0 41.937853904844985
+-90 0
+EOF
+$EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
+ +to +proj=healpix +a=5 +e=0.8 +r_a=4.3220011711888882 -f '%.'5'f' \
+ -E >>${OUT} <<EOF
+0 0
+0 41.810314895778596
+0 -41.810314895778596
+90.0 0
+-90.0 0
+-180 0
+-180 90.0
+-180 -90.0
+0 60.0
+0 -60.0
+EOF
+echo "Test healpix inverse projection on ellipsoid" >> ${OUT}
+$EXE +proj=latlong +a=1 +lon_0=0 +ellps=WGS84 \
+ +to +proj=healpix +a=1 +lon_0=0 +ellps=WGS84 -f '%.'5'f' -I\
+ -E >>${OUT} <<EOF
+0 0.7853981633974483
+-1.5707963267948966 0
+EOF
+$EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
+ +to +proj=healpix +a=5 +e=0.8 +r_a=4.3220011711888882 -f '%.'5'f' -I\
+ -E >>${OUT} <<EOF
+0.0 0.0
+0.0 2.0547874222147415
+0.0 -2.0547874222147415
+6.788983564106746 0.0
+-6.788983564106746 0.0
+-13.577967128213492 0.0
+-10.183475346160119 6.788983564106746
+-10.183475346160119 -6.788983564106746
+0.0 3.351278550178025
+0.0 -3.351278550178025
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test rHEALPix forward projection on sphere north=0 south=0" >> ${OUT}
+$EXE +proj=latlong +R=5 \
+ +to +proj=rhealpix +R=5 +north_square=0 +south_square=0 -f '%.'5'f' \
+  -E >>${OUT} <<EOF
+-180 30.0
+-180 -25.714285714285715
+0 0
+60.0 41.809314895778598
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test rHEALPix forward projection on sphere north=1 south=1" >> ${OUT}
+$EXE +proj=latlong +R=5 \
+ +to +proj=rhealpix +R=5 +north_square=1 +south_square=1 -f '%.'5'f' \
+  -E >>${OUT} <<EOF
+-180 30.0
+-180 -25.714285714285715
+0 0
+60.0 41.809314895778598
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test rHEALPix inverse projection on sphere north=0 south=0" >> ${OUT}
+$EXE +proj=latlong +R=5 \
+ +to +proj=rhealpix +R=5 +north_square=0 +south_square=0 -f '%.'5'f' -I\
+  -E >>${OUT} <<EOF
+0.0 0.0
+0.0 3.9269908169872414
+0.0 -3.9269908169872414
+7.853981633974483 0.0
+-7.853981633974483 0.0
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test rHEALPix inverse projection on sphere north=1 south=1" >> ${OUT}
+$EXE +proj=latlong +R=5 \
+ +to +proj=rhealpix +R=5 +north_square=1 +south_square=1 -f '%.'5'f' -I\
+  -E >>${OUT} <<EOF
+0.0 0.0
+0.0 3.9269908169872414
+0.0 -3.9269908169872414
+7.853981633974483 0.0
+-7.853981633974483 0.0
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test rHEALPix forward projection on ellipsoid north=0 south=0" >> ${OUT}
+$EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
+ +to +proj=rhealpix +a=5 +e=0.8 +r_a=4.3220011711888882 +north_square=0 +south_square=0 -f '%.'5'f' \
+  -E >>${OUT} <<EOF
+0 0
+0 41.810314895778596
+0 -41.810314895778596
+90.0 0
+-90.0 0
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test rHEALPix forward projection on ellipsoid north=1 south=1" >> ${OUT}
+$EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
+ +to +proj=rhealpix +a=5 +e=0.8 +r_a=4.3220011711888882 +north_square=1 +south_square=1 -f '%.'5'f' \
+  -E >>${OUT} <<EOF
+0 0
+0 41.810314895778596
+0 -41.810314895778596
+90.0 0
+-90.0 0
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test rHEALPix inverse projection on ellipsoid north=0 south=0" >> ${OUT}
+$EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
+ +to +proj=rhealpix +a=5 -I +e=0.8 +r_a=4.3220011711888882 +north_square=0 +south_square=0 -f '%.'5'f'\
+  -E >>${OUT} <<EOF
+0.0 0.0
+0.0 2.0547874222147415
+0.0 -2.0547874222147415
+6.788983564106746 0.0
+-6.788983564106746 0.0
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test rHEALPix inverse projection on ellipsoid north=1 south=1" >> ${OUT}
+$EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
+ +to +proj=rhealpix +a=5 -I +e=0.8 +r_a=4.3220011711888882 +north_square=1 +south_square=1 -f '%.'5'f'\
+ -E >>${OUT} <<EOF
+0.0 0.0
+0.0 2.0547874222147415
+0.0 -2.0547874222147415
+6.788983564106746 0.0
+-6.788983564106746 0.0
+EOF
+
+echo "##############################################################" >> ${OUT}
+echo "Test geos projection" >> ${OUT}
+echo "Test geos on a sphere" >> ${OUT}
+$EXE +proj=latlong +ellps=sphere \
+ +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=sphere -E >>${OUT} <<EOF
+16d11'8" 58d35'31"
+-43d11'47" -22d54'30"
+18d25'26" -33d55'31"
+47d58'42" 29d22'11"
+EOF
+echo "Test geos on a ellipsoid" >> ${OUT}
+$EXE +proj=latlong +ellps=sphere \
+ +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=WGS84 -E >>${OUT} <<EOF
+16d11'8" 58d35'31"
+-43d11'47" -22d54'30"
+18d25'26" -33d55'31"
+47d58'42" 29d22'11"
+EOF
+echo "Test inv geos on a sphere" >> ${OUT}
+$EXE +proj=latlong +ellps=sphere \
+ +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=sphere -I -E >>${OUT} <<EOF
+849736.77 4960015.43
+-3780930.93 -2326595.36
+1608689.65 -3412115.56
+3825202.59 2885980.79
+EOF
+echo "Test inv geos on a ellipsoid" >> ${OUT}
+$EXE +proj=latlong +ellps=sphere \
+ +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=WGS84 -I -E >>${OUT} <<EOF
+852862.53 4945122.70
+-3787026.57 -2314765.32
+1612331.00 -3397031.37
+3832522.65 2872185.29
+EOF
+echo "Test geos on a sphere with alternate sweep" >> ${OUT}
+$EXE +proj=latlong +ellps=sphere \
+ +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=sphere +sweep=x -E >>${OUT} <<EOF
+16d11'8" 58d35'31"
+-43d11'47" -22d54'30"
+18d25'26" -33d55'31"
+47d58'42" 29d22'11"
+EOF
+echo "Test geos on a ellipsoid with alternate sweep" >> ${OUT}
+$EXE +proj=latlong +ellps=sphere \
+ +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=WGS84 +sweep=x -E >>${OUT} <<EOF
+16d11'8" 58d35'31"
+-43d11'47" -22d54'30"
+18d25'26" -33d55'31"
+47d58'42" 29d22'11"
+EOF
+echo "Test inv geos on a sphere with alternate sweep" >> ${OUT}
+$EXE +proj=latlong +ellps=sphere \
+ +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=sphere +sweep=x -I -E >>${OUT} <<EOF
+841586.28 4961396.21
+-3772913.22 -2339604.71
+1601377.77 -3415545.15
+3812722.89 2902474.62
+EOF
+echo "Test inv geos on a ellipsoid with alternate sweep" >> ${OUT}
+$EXE +proj=latlong +ellps=sphere \
+ +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=WGS84 +sweep=x -I -E >>${OUT} <<EOF
+844731.03 4946509.59
+-3779077.27 -2327750.87
+1605067.15 -3400461.47
+3820138.08 2888664.15
+EOF
+#
+echo "##############################################################" >> ${OUT}
+echo "Test the Natural Earth Projection" >> ${OUT}
+$EXE +proj=latlong +a=1 +lon_0=0 \
+ +to +proj=natearth +a=6371008.7714 +b=6371008.7714 -f '%.'7'f' \
+ -E >>${OUT} <<EOF
+0.0 0.0 0 0.0 0.0
+0.0 22.5 0 0.0 2525419.569383768
+0.0 45.0 0 0.0 5052537.389973222
+0.0 67.5 0 0.0 7400065.6562573705
+0.0 90.0 0 0.0 9062062.394736718
+45.0 0.0 0 4356790.016612169 0.0
+45.0 22.5 0 4253309.544984069 2525419.569383768
+45.0 45.0 0 3924521.5829515466 5052537.389973222
+45.0 67.5 0 3354937.47115583 7400065.6562573705
+45.0 90.0 0 2397978.2448443635 9062062.394736718
+90.0 0.0 0 8713580.033224339 0.0
+90.0 22.5 0 8506619.089968137 2525419.569383768
+90.0 45.0 0 7849043.165903093 5052537.389973222
+90.0 67.5 0 6709874.94231166 7400065.6562573705
+90.0 90.0 0 4795956.489688727 9062062.394736718
+135.0 0.0 0 1.3070370049836507E7 0.0
+135.0 22.5 0 1.2759928634952208E7 2525419.569383768
+135.0 45.0 0 1.177356474885464E7 5052537.389973222
+135.0 67.5 0 1.0064812413467491E7 7400065.6562573705
+135.0 90.0 0 7193934.734533091 9062062.394736718
+180.0 0.0 0 1.7427160066448677E7 0.0
+180.0 22.5 0 1.7013238179936275E7 2525419.569383768
+180.0 45.0 0 1.5698086331806187E7 5052537.389973222
+180.0 67.5 0 1.341974988462332E7 7400065.6562573705
+180.0 90.0 0 9591912.979377454 9062062.394736718
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test the Natural Earth II Projection" >> ${OUT}
+$EXE +proj=latlong +a=1 +lon_0=0 \
+ +to +proj=natearth2 +a=6371008.7714 +b=6371008.7714 -f '%.'7'f' \
+ -E >>${OUT} <<EOF
+0.0 0.0 0 0.00000000 0.00000000
+0.0 22.5 0 0.00000000 2531453.57080958
+0.0 45.0 0 0.00000000 5051471.50086845
+0.0 67.5 0 0.00000000 7395411.22478983
+0.0 90.0 0 0.00000000 9073776.52662810
+45.0 0.0 0 4239151.18200719 0.00000000
+45.0 22.5 0 4138348.61904244 2531453.57080958
+45.0 45.0 0 3830621.33880773 5051471.50086845
+45.0 67.5 0 3158326.32836996 7395411.22478983
+45.0 90.0 0 957973.37034235 9073776.52662810
+90.0 0.0 0 8478302.36401439 0.00000000
+90.0 22.5 0 8276697.23808488 2531453.57080958
+90.0 45.0 0 7661242.67761547 5051471.50086845
+90.0 67.5 0 6316652.65673992 7395411.22478983
+90.0 90.0 0 1915946.74068470 9073776.52662810
+135.0 0.0 0 12717453.54602160 0.00000000
+135.0 22.5 0 12415045.85712730 2531453.57080958
+135.0 45.0 0 11491864.01642320 5051471.50086845
+135.0 67.5 0 9474978.98510988 7395411.22478983
+135.0 90.0 0 2873920.11102705 9073776.52662810
+180.0 0.0 0 16956604.72802880 0.00000000
+180.0 22.5 0 16553394.47616980 2531453.57080958
+180.0 45.0 0 15322485.35523090 5051471.50086845
+180.0 67.5 0 12633305.31347990 7395411.22478983
+180.0 90.0 0 3831893.48136940 9073776.52662810
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test the Compact Miller projection" >> ${OUT}
+$EXE +proj=latlong +a=1 +lon_0=0 \
+ +to +proj=comill +a=6371008.7714 +b=6371008.7714 -f '%.'7'f' \
+ -E >>${OUT} <<EOF
+0.0 0.0 0 0.0 0.0
+0.0 22.5 0 0.0 2537439.6610749415
+0.0 45.0 0 0.0 5391682.432264133
+0.0 67.5 0 0.0 8661480.510260897
+0.0 90.0 0 0.0 12009484.264916677
+45.0 0.0 0 5003778.588046594 0.0
+45.0 22.5 0 5003778.588046594 2537439.6610749415
+45.0 45.0 0 5003778.588046594 5391682.432264133
+45.0 67.5 0 5003778.588046594 8661480.510260897
+45.0 90.0 0 5003778.588046594 12009484.264916677
+90.0 0.0 0 10007557.176093187 0.0
+90.0 22.5 0 10007557.176093187 2537439.6610749415
+90.0 45.0 0 10007557.176093187 5391682.432264133
+90.0 67.5 0 10007557.176093187 8661480.510260897
+90.0 90.0 0 10007557.176093187 12009484.264916677
+135.0 0.0 0 15011335.76413978 0.0
+135.0 22.5 0 15011335.76413978 2537439.6610749415
+135.0 45.0 0 15011335.76413978 5391682.432264133
+135.0 67.5 0 15011335.76413978 8661480.510260897
+135.0 90.0 0 15011335.76413978 12009484.264916677
+180.0 0.0 0 20015114.352186374 0.0
+180.0 22.5 0 20015114.352186374 2537439.6610749415
+180.0 45.0 0 20015114.352186374 5391682.432264133
+180.0 67.5 0 20015114.352186374 8661480.510260897
+180.0 90.0 0 20015114.352186374 12009484.264916677
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test pconic (#148)" >> ${OUT}
+#
+$EXE +proj=latlong +datum=WGS84 \
+ +to +proj=pconic  +units=m +lat_1=20n +lat_2=60n +lon_0=60W +datum=WGS84 \
+ -E >>${OUT} <<EOF
+-70.4 -23.65
+EOF
+$EXE +proj=pconic  +units=m +lat_1=20n +lat_2=60n +lon_0=60W +datum=WGS84 \
+ +to +proj=latlong +datum=WGS84 \
+ -E >>${OUT} <<EOF
+-2240096.40 -6940342.15
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test laea" >> ${OUT}
+#
+$EXE -f '%.12f' \
+     +proj=laea +lat_0=45 +lon_0=-100 +units=m +datum=WGS84 +no_defs \
+ +to +proj=latlong +datum=WGS84 \
+ -E >>${OUT} <<EOF
+-6086629.0 4488761.0
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test forward calcofi projection" >> ${OUT}
+$EXE +proj=latlong +ellps=clrk66 \
+ +to +proj=calcofi +ellps=clrk66 \
+ -E >>${OUT} <<EOF
+120d40'42.273"W	38d56'50.766"N
+121d9'W	34d9'N
+123d59'56.066"W	30d25'4.617"N
+EOF
+echo "Test inverse calcofi projection" >> ${OUT}
+$EXE +proj=calcofi +ellps=clrk66 \
+ +to +proj=longlat +ellps=clrk66 \
+ -E >>${OUT} <<EOF
+60 20
+80 60
+90 120
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Check inverse error handling with ob_tran (#225)" >> ${OUT}
+$EXE +proj=ob_tran \
+    +o_proj=moll +o_lon_p=LON_POLE +o_lat_p=LAT_POLE +lon_0=180 +ellps=WGS84 \
+ -E >>${OUT} <<EOF
+300000 400000
+20000000 30000000
+EOF
+echo "Test inverse handling" >> ${OUT}
+$EXE -I +proj=ob_tran \
+    +o_proj=moll +o_lon_p=LON_POLE +o_lat_p=LAT_POLE +lon_0=180 +ellps=WGS84 \
+ -E >>${OUT} <<EOF
+10 20
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test MGI datum gives expected results (#207)" >> ${OUT}
+#
+$EXE -f '%.7f' \
+     +proj=latlong +datum=WGS84 \
+ +to +init=epsg:31284 \
+ -E >>${OUT} <<EOF
+16.33 48.20
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test omerc sensitivity with locations 90d from origin(#114)" >> ${OUT}
+#
+$EXE -f '%.8f' \
+     +proj=latlong +ellps=WGS84 \
+ +to +proj=omerc +ellps=WGS84 +lon_1=62.581150 +lat_1=74.856102 \
+                 +lon_2=53.942810 +lat_2=74.905884 +units=km +no_rot +no_defs \
+ -E >>${OUT} <<EOF
+56.958381652832 72.8798
+56.9584 72.8798
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test omerc differences between poles (#190)" >> ${OUT}
+# First, north pole.
+$EXE -f '%.3f' \
+     +proj=latlong +ellps=WGS84 \
+ +to +proj=omerc +ellps=WGS84 +datum=WGS84 +no_rot \
+                 +lon_1=-27 +lat_1=70 +lon_2=-38 +lat_2=80 +lat_0=70 \
+ -E >>${OUT} <<EOF
+-27 70
+-27 80
+-27 89.9
+163 89.9
+163 80
+EOF
+# Again, for the south pole.
+$EXE -f '%.3f' \
+     +proj=latlong +ellps=WGS84 \
+ +to +proj=omerc +ellps=WGS84 +datum=WGS84 +no_rot \
+                 +lon_1=-27 +lat_1=-70 +lon_2=-38 +lat_2=-80 +lat_0=-70 \
+ -E >>${OUT} <<EOF
+-27 -70
+-27 -80
+-27 -89.9
+163 -89.9
+163 -80
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test qsc" >> ${OUT}
+#
+$EXE -f '%.7f' \
+     +proj=latlong +datum=WGS84 \
+ +to +proj=qsc +datum=WGS84 +no_defs \
+ -E >>${OUT} <<EOF
+13 -10
+EOF
+$EXE -f '%.13f' \
+     +proj=qsc +datum=WGS84 +no_defs \
+ +to +proj=latlong +datum=WGS84 \
+ -E >>${OUT} <<EOF
+2073986.94908809568733	-1680858.27222427958623
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test bug 229" >> ${OUT}
+#
+$EXE -f '%.13f' \
+     +init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0 \
+ +to +proj=latlong +datum=WGS84 +no_defs \
+ -E >>${OUT} <<EOF
+13 -10
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test bug 229 (2)" >> ${OUT}
+#
+$EXE -f '%.13f' \
+     +init=epsg:4326 +to +init=epsg:4326 \
+ -E >>${OUT} <<EOF
+13 -10
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test bug 244 " >> ${OUT}
+#
+$EXE -f '%.8f' \
+     +init=epsg:4326 \
+     +to +proj=aeqd +lon_0=130.0 +lat_0=40.0 +a=6378137 +b=6378137 +units=m +no_defs \
+ -E >>${OUT} <<EOF
+-140.100000 -87.000000
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test bug 244 (2)" >> ${OUT}
+#
+$EXE -f '%.11f' \
+     +proj=aeqd +lon_0=130.0 +lat_0=40.0 +a=6378137 +b=6378137 +units=m +no_defs \
+     +to +init=epsg:4326 \
+ -E >>${OUT} <<EOF
+987122.418330284 -14429896.539530909
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test bug 245 (use +datum=carthage)" >> ${OUT}
+#
+$EXE -f '%.7f' \
+     +proj=longlat +datum=WGS84 +no_defs +to +proj=utm +zone=32 +datum=carthage +no_defs \
+ -E >>${OUT} <<EOF
+10 34
+EOF
+echo "##############################################################" >> ${OUT}
+echo "Test bug 245 (use expansion of +datum=carthage)" >> ${OUT}
+#
+$EXE -f '%.7f' \
+     +proj=longlat +datum=WGS84 +no_defs +to +proj=utm +zone=32 +a=6378249.2 +b=6356515 +towgs84=-263.0,6.0,431.0 +no_defs \
+ -E >>${OUT} <<EOF
+10 34
+EOF
+echo  "##############################################################" >> ${OUT}
+echo  "Test SCH forward projection" >> ${OUT}
+#
+$EXE -f '%.7f' \
+    +proj=latlong +datum=WGS84 +to +proj=sch +datum=WGS84 +plat_0=30.0 +plon_0=45.0 \
+    +phdg_0=-12.0 +nodefs \
+ -E >> ${OUT} <<EOF
+0.0 0.0
+0.0 90.0
+45.0 45.0
+45.1 44.9
+44.9 45.1
+30.0 45.0
+EOF
+echo  "##############################################################" >> ${OUT}
+echo  "Test SCH inverse projection" >> ${OUT}
+#
+$EXE -f '%.7f' \
+    +proj=sch +datum=WGS84 +plat_0=30.0 +plon_0=45.0 +phdg_0=-12.0 +nodefs +to \
+    +proj=latlong +datum=WGS84 \
+ -E >> ${OUT} <<EOF
+0. 0.
+0. 1000.
+1000. 0.
+1000. 1000.
+EOF
+##############################################################################
+echo "##############################################################" >> ${OUT}
+echo "Test issue #316 (switch utm to use etmerc)" >> ${OUT}
+#
+$EXE -f '%.6f' \
+     +proj=latlong +datum=WGS84 +no_defs +to +proj=utm +zone=35 +datum=WGS84 +no_defs \
+ -E >>${OUT} <<EOF
+0 83
+EOF
+##############################################################################
+echo "##############################################################" >> ${OUT}
+echo "Test issue #316 (switch utm to use etmerc)" >> ${OUT}
+#
+$EXE -f '%.6f' \
+     +proj=latlong +datum=WGS84 +no_defs +to +proj=etmerc +datum=WGS84 +no_defs +k=0.9996 +lon_0=27 +x_0=500000 -f %.6f \
+ -E >>${OUT} <<EOF
+0 83
+EOF
+##############################################################################
+echo  "##############################################################" >> ${OUT}
+echo  "Test nzmg forward projection" >> ${OUT}
+#
+$EXE -f '%.7f' \
+    +proj=latlong +datum=WGS84 +to \
+    +proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000 +y_0=6023150 +ellps=WGS84 +units=m +no_defs \
+ -E >> ${OUT} <<EOF
+175. -40. 0.
+EOF
+##############################################################################
+echo  "##############################################################" >> ${OUT}
+echo  "Test nzmg inverse projection" >> ${OUT}
+#
+$EXE -f '%.7f' \
+    +proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000 +y_0=6023150 +ellps=WGS84 +units=m +no_defs +to \
+    +proj=latlong +datum=WGS84 \
+ -E >> ${OUT} <<EOF
+2680778.57267967 6132228.07645127 0.
+EOF
+##############################################################################
+echo  "##############################################################" >> ${OUT}
+echo  "Test misrsom forward projection" >> ${OUT}
+#
+$EXE -f '%.7f' \
+    +proj=latlong +datum=WGS84 +no_defs +to \
+    +proj=misrsom +path=41 +ellps=WGS84 +units=m \
+ -E >> ${OUT} <<EOF
+48.64966165540372 66.2263195368941 0.
+EOF
+##############################################################################
+echo  "##############################################################" >> ${OUT}
+echo  "Test misrsom inverse projection" >> ${OUT}
+#
+$EXE -f '%.7f' \
+    +proj=misrsom +path=41 +ellps=WGS84 +units=m +no_defs +to \
+    +proj=latlong +datum=WGS84 -f '%0.7f' \
+ -E >> ${OUT} <<EOF
+7461300.0 528000.0 0.0
+EOF
+echo  "##############################################################" >> ${OUT}
+echo  "Test patterson forward projection" >> ${OUT}
+#
+$EXE -f '%0.8f' \
+    +proj=latlong +datum=WGS84 \
+    +to +proj=patterson +a=6371008.7714 +b=6371008.7714 +units=m \
+ -E >> ${OUT} <<EOF
+-180 90
+-135 67.5
+-90 45
+-45 22.5
+0 0
+45 -22.5
+90 -45
+135 -67.5
+180 -90
+EOF
+echo  "##############################################################" >> ${OUT}
+echo  "Test patterson inverse projection" >> ${OUT}
+#
+$EXE -f '%0.3f' \
+    +proj=patterson +a=6371008.7714 +b=6371008.7714 +units=m \
+    +to +proj=latlong +datum=WGS84 \
+ -E >> ${OUT} <<EOF
+-20015114.352186374 11409566.822831295
+-15011335.76413978 8729502.054111844
+-10007557.176093187 5366413.421153781
+-5003778.588046594 2551415.729669344
+0.0 0.0
+5003778.588046594 -2551415.729669344
+10007557.176093187 -5366413.421153781
+15011335.76413978 -8729502.054111844
+20015114.352186374 -11409566.822831295
+EOF
+
+
+# Done!
+# do 'diff' with distribution results
+echo "diff ${OUT} with ${OUT}.dist"
+diff -b ${OUT} ${NAD_DIR}/${OUT}.dist
+if [ $? -ne 0 ] ; then
+	echo  ""
+	echo "PROBLEMS HAVE OCCURED"
+	echo "test file ${OUT} saved"
+    echo
+	exit 100
+else
+	echo "TEST OK"
+	echo "test file ${OUT} removed"
+    echo
+	/bin/rm -f ${OUT}
+	exit 0
+fi