~ubuntu-branches/ubuntu/vivid/proj/vivid

« back to all changes in this revision

Viewing changes to src/geod_set.c

Committer: Bazaar Package Importer
Author(s): Peter S Galbraith
Date: 2002-01-11 10:27:12 UTC
Revision ID: james.westby@ubuntu.com-20020111102712-ayi18r8y2eesv0y9

Tags: upstream-4.4.5

Import upstream version 4.4.5

files added:

AUTHORS

COPYING

ChangeLog

INSTALL

Makefile.am

Makefile.in

NEWS

README

aclocal.m4

config.guess

config.sub

configure

configure.in

install-sh

ltconfig

ltmain.sh

man/Makefile.am

man/Makefile.in

man/man1

man/man1/Makefile.am

man/man1/Makefile.in

man/man1/cs2cs.1

man/man1/geod.1

man/man1/nad2nad.1

man/man1/proj.1

man/man3

man/man3/Makefile.am

man/man3/Makefile.in

man/man3/pj_init.3

missing

mkinstalldirs

nad/GL27

nad/Makefile.am

nad/Makefile.in

nad/README

nad/epsg

nad/nad.lst

nad/nad27

nad/nad83

nad/pj_out27.dist

nad/pj_out83.dist

nad/proj_def.dat

nad/test27

nad/test83

nad/world

src/Makefile.am

src/Makefile.in

src/PJ_aea.c

src/PJ_aeqd.c

src/PJ_airy.c

src/PJ_aitoff.c

src/PJ_august.c

src/PJ_bacon.c

src/PJ_bipc.c

src/PJ_boggs.c

src/PJ_bonne.c

src/PJ_cass.c

src/PJ_cc.c

src/PJ_cea.c

src/PJ_chamb.c

src/PJ_collg.c

src/PJ_crast.c

src/PJ_denoy.c

src/PJ_eck1.c

src/PJ_eck2.c

src/PJ_eck3.c

src/PJ_eck4.c

src/PJ_eck5.c

src/PJ_eqc.c

src/PJ_eqdc.c

src/PJ_fahey.c

src/PJ_fouc_s.c

src/PJ_gall.c

src/PJ_gins8.c

src/PJ_gn_sinu.c

src/PJ_gnom.c

src/PJ_goode.c

src/PJ_hammer.c

src/PJ_hatano.c

src/PJ_imw_p.c

src/PJ_labrd.c

src/PJ_laea.c

src/PJ_lagrng.c

src/PJ_larr.c

src/PJ_lask.c

src/PJ_lcc.c

src/PJ_loxim.c

src/PJ_lsat.c

src/PJ_mbt_fps.c

src/PJ_mbtfpp.c

src/PJ_mbtfpq.c

src/PJ_merc.c

src/PJ_mill.c

src/PJ_mod_ster.c

src/PJ_moll.c

src/PJ_mpoly.c

src/PJ_nell.c

src/PJ_nell_h.c

src/PJ_nocol.c

src/PJ_nsper.c

src/PJ_nzmg.c

src/PJ_ob_tran.c

src/PJ_ocea.c

src/PJ_oea.c

src/PJ_omerc.c

src/PJ_ortho.c

src/PJ_poly.c

src/PJ_putp2.c

src/PJ_putp3.c

src/PJ_putp4p.c

src/PJ_putp5.c

src/PJ_putp6.c

src/PJ_robin.c

src/PJ_rpoly.c

src/PJ_sconics.c

src/PJ_somerc.c

src/PJ_stere.c

src/PJ_sts.c

src/PJ_tcc.c

src/PJ_tcea.c

src/PJ_tmerc.c

src/PJ_tpeqd.c

src/PJ_urm5.c

src/PJ_urmfps.c

src/PJ_vandg.c

src/PJ_vandg2.c

src/PJ_vandg4.c

src/PJ_wag2.c

src/PJ_wag3.c

src/PJ_wag7.c

src/PJ_wink1.c

src/PJ_wink2.c

src/aasincos.c

src/adjlon.c

src/bch2bps.c

src/bchgen.c

src/biveval.c

src/cs2cs.c

src/dmstor.c

src/emess.c

src/emess.h

src/gen_cheb.c

src/geocent.c

src/geocent.h

src/geod.c

src/geod_for.c

src/geod_inv.c

src/geod_set.c

src/geodesic.h

src/makefile.vc

src/mk_cheby.c

src/nad2bin.c

src/nad2nad.c

src/nad_cvt.c

src/nad_init.c

src/nad_intr.c

src/nad_list.h

src/p_series.c

src/pj_apply_gridshift.c

src/pj_auth.c

src/pj_datum_set.c

src/pj_datums.c

src/pj_deriv.c

src/pj_ell_set.c

src/pj_ellps.c

src/pj_errno.c

src/pj_factors.c

src/pj_fwd.c

src/pj_init.c

src/pj_inv.c

src/pj_latlong.c

src/pj_list.c

src/pj_list.h

src/pj_malloc.c

src/pj_mlfn.c

src/pj_msfn.c

src/pj_open_lib.c

src/pj_param.c

src/pj_phi2.c

src/pj_pr_list.c

src/pj_qsfn.c

src/pj_release.c

src/pj_strerrno.c

src/pj_transform.c

src/pj_tsfn.c

src/pj_units.c

src/pj_utils.c

src/pj_zpoly1.c

src/proj.c

src/proj.def

src/proj_api.h

src/proj_config.h.in

src/projects.h

src/rtodms.c

src/stamp-h.in

src/strtod.c

src/vector1.c

Show diffs side-by-side

added added

removed removed

src/geod_set.c

#ifndef lint

static const char SCCSID[]="@(#)geod_set.c 4.8 95/09/23 GIE REL";

#endif

#include "projects.h"

#include "geodesic.h"

void

geod_set(int argc, char **argv) {

paralist *start = 0, *curr;

double es;

char *name;

int i;

/* put arguments into internal linked list */

if (argc <= 0)

emess(1, "no arguments in initialization list");

for (i = 0; i < argc; ++i)

if (i)

curr = curr->next = pj_mkparam(argv[i]);

else

start = curr = pj_mkparam(argv[i]);

/* set elliptical parameters */

if (pj_ell_set(start, &a, &es)) emess(1,"ellipse setup failure");

/* set units */

if (name = pj_param(start, "sunits").s) {

char *s;

for (i = 0; (s = pj_units[i].id) && strcmp(name, s) ; ++i) ;

if (!s)

emess(1,"%s unknown unit conversion id", name);

fr_meter = 1. / (to_meter = atof(pj_units[i].to_meter));

} else

to_meter = fr_meter = 1.;

if (ellipse = es != 0.) {

onef = sqrt(1. - es);

f = 1 - onef;

f2 = f/2;

f4 = f/4;

f64 = f*f/64;

} else {

onef = 1.;

f = f2 = f4 = f64 = 0.;

}

/* check if line or arc mode */

if (pj_param(start, "tlat_1").i) {

double del_S;

#undef f

phi1 = pj_param(start, "rlat_1").f;

lam1 = pj_param(start, "rlon_1").f;

if (pj_param(start, "tlat_2").i) {

phi2 = pj_param(start, "rlat_2").f;

lam2 = pj_param(start, "rlon_2").f;

geod_inv();

geod_pre();

} else if (S = pj_param(start, "dS").f) {

al12 = pj_param(start, "rA").f;

geod_pre();

geod_for();

} else emess(1,"incomplete geodesic/arc info");

if ((n_alpha = pj_param(start, "in_A").i) > 0) {

if (!(del_alpha = pj_param(start, "rdel_A").f))

emess(1,"del azimuth == 0");

} else if (del_S = fabs(pj_param(start, "ddel_S").f)) {

n_S = S / del_S + .5;

} else if ((n_S = pj_param(start, "in_S").i) <= 0)

emess(1,"no interval divisor selected");

}

/* free up linked list */

for ( ; start; start = curr) {

curr = start->next;

pj_dalloc(start);

}

Older »