Return to homing.c CVS log

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

Annotation of researchv8dc/cmd/map/libmap/homing.c, revision 1.1.1.1

1.1       root        1: #include "map.h"
                      2: 
                      3: struct coord stdpar;
                      4: 
                      5: Xmecca(place, x, y)
                      6: struct place *place;
                      7: float *x, *y;
                      8: {
                      9:        *x = - place->wlon.l;
                     10:        *y = -(place->nlat.c*stdpar.s -
                     11:                place->nlat.s*stdpar.c*place->wlon.c)/stdpar.c;
                     12:        if(fabs(place->wlon.l) > 0.01)
                     13:                *y *= *x/sin(*x);
                     14:        if(fabs(*y)>2.0)
                     15:                return(0);
                     16:        if(place->nlat.s*stdpar.s +
                     17:           place->nlat.c*stdpar.c*place->wlon.c <0)
                     18:                return(-1);
                     19:        return(1);
                     20: }
                     21: 
                     22: int (*mecca(par))()
                     23: float par;
                     24: {
                     25:        if(fabs(par)>80.)
                     26:                return(0);
                     27:        deg2rad(par,&stdpar);
                     28:        return(Xmecca);
                     29: }
                     30: 
                     31: Xhoming(place, x, y)
                     32: struct place *place;
                     33: float *x, *y;
                     34: {
                     35:        int i = Xmecca(place, x, y);
                     36:        float r, h;
                     37:        if(i != 1)
                     38:                return(i);
                     39:        r = acos(stdpar.s*place->nlat.s +
                     40:                stdpar.c*place->nlat.c*place->wlon.c);
                     41:        h = hypot(*x, *y);
                     42:        r = h<.001? r: r/h;
                     43:        *x *= r;
                     44:        *y *= r;
                     45:        return(1);
                     46: }
                     47: 
                     48: int (*homing(par))()
                     49: float par;
                     50: {
                     51:        if(mecca(par)==0)
                     52:                return(0);
                     53:        return(Xhoming);
                     54: }

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