42
42
P->rok = P->a / P->k0;
43
43
P->rtk = P->a * P->k0;
44
if ( pj_param(P->params, "talpha").i) {
45
alpha = pj_param(P->params, "ralpha").f;
46
lonz = pj_param(P->params, "rlonc").f;
44
if ( pj_param(P->ctx, P->params, "talpha").i) {
45
alpha = pj_param(P->ctx, P->params, "ralpha").f;
46
lonz = pj_param(P->ctx, P->params, "rlonc").f;
47
47
P->singam = atan(-cos(alpha)/(-sin(phi_0) * sin(alpha))) + lonz;
48
48
P->sinphi = asin(cos(phi_0) * sin(alpha));
50
phi_1 = pj_param(P->params, "rlat_1").f;
51
phi_2 = pj_param(P->params, "rlat_2").f;
52
lam_1 = pj_param(P->params, "rlon_1").f;
53
lam_2 = pj_param(P->params, "rlon_2").f;
50
phi_1 = pj_param(P->ctx, P->params, "rlat_1").f;
51
phi_2 = pj_param(P->ctx, P->params, "rlat_2").f;
52
lam_1 = pj_param(P->ctx, P->params, "rlon_1").f;
53
lam_2 = pj_param(P->ctx, P->params, "rlon_2").f;
54
54
P->singam = atan2(cos(phi_1) * sin(phi_2) * cos(lam_1) -
55
55
sin(phi_1) * cos(phi_2) * cos(lam_2),
56
56
sin(phi_1) * cos(phi_2) * sin(lam_2) -