From 647e12fed0c9c5b373be8b03061d14bda40baea3 Mon Sep 17 00:00:00 2001 From: Frank Warmerdam Date: Sun, 15 Dec 2002 00:13:30 +0000 Subject: lat_0 may now be set by user, but still defaults to 49d30N git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@1061 4e78687f-474d-0410-85f9-8d5e500ac6b2 --- src/PJ_krovak.c | 15 +++++++++++---- 1 file changed, 11 insertions(+), 4 deletions(-) (limited to 'src') diff --git a/src/PJ_krovak.c b/src/PJ_krovak.c index 4781d727..8dc9630c 100644 --- a/src/PJ_krovak.c +++ b/src/PJ_krovak.c @@ -30,6 +30,9 @@ ****************************************************************************** * * $Log$ + * Revision 1.3 2002/12/15 00:13:30 warmerda + * lat_0 may now be set by user, but still defaults to 49d30N + * * Revision 1.2 2002/12/14 19:35:21 warmerda * updated headers * @@ -54,13 +57,13 @@ FORWARD(s_forward); /* spheroid */ char tmp[16]; /* Constants, identical to inverse transform function */ - double s45, s90, fi0, e2, e, alfa, uq, u0, g, k, k1, n0, ro0, ad, a, s0, n; - double gfi, u, lon17, lon42, lamdd, deltav, s, d, eps, ro; + double s45, s90, e2, e, alfa, uq, u0, g, k, k1, n0, ro0, ad, a, s0, n; + double gfi, u, fi0, lon17, lon42, lamdd, deltav, s, d, eps, ro; s45 = 0.785398163397448; /* 45° */ s90 = 2 * s45; - fi0 = 0.863937979737193; /* Latitude of projection centre 49°30' */ + fi0 = P->phi0; /* Latitude of projection centre 49° 30' */ /* Ellipsoid is used as Parameter in for.c and inv.c, therefore a must be set to 1 here. @@ -137,7 +140,7 @@ INVERSE(s_inverse); /* spheroid */ s45 = 0.785398163397448; /* 45° */ s90 = 2 * s45; - fi0 = 0.863937979737193; /* Latitude of projection centre 49° 30' */ + fi0 = P->phi0; /* Latitude of projection centre 49° 30' */ /* Ellipsoid is used as Parameter in for.c and inv.c, therefore a must @@ -216,6 +219,10 @@ ENTRY0(krovak) P->a = 6377397.155; P->e = sqrt(P->es = 0.006674372230614); + /* if latitude of projection center is not set, use 49d30'N */ + if (!pj_param(P->params, "tlat_0").i) + P->phi0 = 0.863937979737193; + /* always the same */ P->inv = s_inverse; P->fwd = s_forward; -- cgit v1.2.3