Return to bonne.c CVS log

Up to [Research Unix] / researchv10no / cmd / map / libmap

Annotation of researchv10no/cmd/map/libmap/bonne.c, revision 1.1.1.1

1.1       root        1: #include "map.h"
                      2: 
                      3: static struct coord stdpar;
                      4: static double r0;
                      5: 
                      6: static
                      7: Xbonne(struct place *place, double *x, double *y)
                      8: {
                      9:        double r, alpha;
                     10:        r = r0 - place->nlat.l;
                     11:        if(r<.001)
                     12:                if(fabs(stdpar.c)<1e-10)
                     13:                        alpha = place->wlon.l;
                     14:                else if(fabs(place->nlat.c)==0)
                     15:                        alpha = 0;
                     16:                else 
                     17:                        alpha = place->wlon.l/(1+
                     18:                                stdpar.c*stdpar.c*stdpar.c/place->nlat.c/3);
                     19:        else
                     20:                alpha = place->wlon.l * place->nlat.c / r;
                     21:        *x = - r*sin(alpha);
                     22:        *y = - r*cos(alpha);
                     23:        return(1);
                     24: }
                     25: 
                     26: proj
                     27: bonne(double par)
                     28: {
                     29:        if(fabs(par*RAD) < .01)
                     30:                return(Xsinusoidal);
                     31:        deg2rad(par, &stdpar);
                     32:        r0 = stdpar.c/stdpar.s + stdpar.l;
                     33:        return(Xbonne);
                     34: }

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.