Return to stars.c CVS log

Up to [Research Unix] / researchv10no / cmd / sky

Annotation of researchv10no/cmd/sky/stars.c, revision 1.1.1.1

1.1       root        1: #include "sky.h"
                      2: #include <stdio.h>
                      3: 
                      4: stars()
                      5: {
                      6:        char starname[13];
                      7:        register i;
                      8:        double temp1;
                      9:        FILE *f;
                     10: 
                     11:        starname[12] = 0;
                     12:        object = starname;
                     13: 
                     14: 
                     15: /*
                     16:  *     read epoch of star table
                     17:  *     and convert to internal format
                     18:  *     epoch 1900.0 is J.D. 2415020.313
                     19:  */
                     20: 
                     21:        f = fopen(startab, "r");
                     22:        if(f == NULL) {
                     23:                printf("%s?\n", startab);
                     24:                return;
                     25:        }
                     26:        rline(f);
                     27:        epoch = atof(line);
                     28:        epoch = (epoch-1900.0) * 365.2422 + 0.313;
                     29: 
                     30: /*
                     31:  *     read mean places of stars at epoch of star table
                     32:  */
                     33: 
                     34: loop:
                     35:        if(rline(f)) {
                     36:                fclose(f);
                     37:                return;
                     38:        }
                     39:        for(i=0;i<12;i++){
                     40:                starname[i] = line[i];
                     41:        }
                     42:        agc = (long)atof(line+12);
                     43:        rah = atof(line+17);
                     44:        ram = atof(line+20);
                     45:        ras = atof(line+23);
                     46:        da = atof(line+30);
                     47:        dday = atof(line+37);
                     48:        dmin = atof(line+41);
                     49:        dsec = atof(line+44);
                     50:        dd = atof(line+50);
                     51:        px = atof(line+57);
                     52:        mag = atof(line+64);
                     53:        const = (int)atof(line+70);
                     54: 
                     55: /*
                     56:  *     convert rt ascension and declination to internal format
                     57:  */
                     58: 
                     59:        alpha = rah + ram/60. + ras/3600.;
                     60: 
                     61: 
                     62:        delta = fabs(dday) + dmin/60. + dsec/3600.;
                     63:        if(dday < 0.)
                     64:                delta = -delta;
                     65:        alpha *= 15.*radian;
                     66:        delta *= radian;
                     67: 
                     68:        star();
                     69:        output();
                     70: /*
                     71:  *     Check for occultation.
                     72:  */
                     73: 
                     74:        if(!((flags&GEO)||(flags&HELIO))){
                     75:                temp1 = sin(decl2)*sin(moonde) + cos(decl2)*cos(moonde)*cos(ra-moonra);
                     76:                temp1 = atan2(sqrt(1.-temp1*temp1),temp1)/radsec;
                     77:                temp1 = temp1/moonsd;
                     78:                if(temp1 <= 1.0){
                     79:                        printf("Occultation: Dist. = %.4f\n", temp1);
                     80:                }else if(temp1 <= 1.2){
                     81:                        printf("Near occultation: Dist. = %.4f\n", temp1);
                     82:                }
                     83:        }
                     84: 
                     85:        goto loop;
                     86: }

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