Changeset 7976


Ignore:
Timestamp:
03/27/19 09:20:11 (3 weeks ago)
Author:
stefan
Message:

read arc file

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/libsrc/util/cosmos0putaway.c

    r7974 r7976  
    146146        char    mn[3]; 
    147147        char    sec[7]; 
    148         char    cos_date[35]; 
     148        char    cos_date[35]; /* Start date of data requested */ 
    149149        char    record_id[35]; 
    150150        char    str[COSMOSLINELEN]; 
     
    189189 
    190190        /* cos_date "06/17/2018, 18:34:38.004 GMT (Q=5) "<-35 chars; Q stands for Quality, we have no way to know that */ 
    191         sprintf(cos_date, "%s/%s/%s, %s:%s:%s  GMT", mo, dy, year, hr, mn, sec); 
     191        sprintf(cos_date, "%s/%s/%s, %s:%s:%s  UTC", mo, dy, year, hr, mn, sec); 
    192192 
    193193        RecursiveCreateDir(OutDir); 
     
    255255        } 
    256256        timefunc = localtime; 
    257  
     257        ot = (time_t)(arcmsg.sum.ot - GSEC1970);  
     258        timeinfo = timefunc(&ot); 
    258259        /* Copy the library file to the putaway file*/ 
    259260/*      c = fgetc(COSMOS0Libfp); 
     
    270271                                        records may use “Test Record of”, etc.). 
    271272                                        41-80 Earthquake date and time (including time zone). */ 
    272                             /* Maybe one day we'll read Hypoinverse output to detect events, but for now, this is blank */ 
    273                         fprintf(COSMOS0fp, "Record of                 Earthquake of                                         \n"); 
     273                        strftime(timestr, 80, "Record of                 Earthquake of %B %d, %Y, %H:%M:%S   UTC", timeinfo); 
     274                        fprintf(COSMOS0fp, "%s\n", timestr); 
     275                        break; 
     276                case 3: /*3 12-46 Hypocenter: 12-20 Latitude (positive North); 22-31 Longitude (positive East); and 
     277                                        35-37 Depth (km); 40-46 Source agency for hypocenter information. 
     278                                        47-80 Magnitude(s), including source agency for values (may be entered manually, so 
     279                                        spacing may be variable).**/ 
     280                        sprintf(tempstr, "Hypocenter: %9.4f %9.4f H=%3.0fkm M=%3.1f, Md=%3.1f", 
     281                                arcmsg.sum.lat, arcmsg.sum.lon, arcmsg.sum.z, arcmsg.sum.Mpref, arcmsg.sum.Md); 
     282                        fprintf(COSMOS0fp, "%s\n", tempstr); 
    274283                        break; 
    275284                case 4: /*      4 9-42 Earthquake origin date and time, in UTC (i.e., GMT), with source agency; 43-80, 
    276285                                        For use by preparing agency; may include note, e.g., “Preliminary Processing”, etc.*/ 
    277                         ot = (time_t)(arcmsg.sum.ot - GSEC1970);  
    278                         timeinfo = timefunc(&ot); 
    279                         strftime(timestr, 80, "%Y.%m.%d %H:%M:%S", timeinfo); // Prepare origin time 
     286 
     287                        /* Prepare origin time. Second can have tenths place, so leaving two spaces here to indicate  
     288                           we're not going to that level of precision, the timestruct doesn't have tenths. Probably worth revisiting.*/ 
     289                        strftime(timestr, 80, "%m/%d/%Y, %H:%M:%S  ", timeinfo);  
    280290                        sprintf(tempstr, "Origin: %s UTC", timestr); 
    281291                        fprintf(COSMOS0fp, "%s\n", tempstr); 
     
    306316                        time(&rawtime); 
    307317                        timeinfo = gmtime(&rawtime); 
    308                         sprintf(tempstr, "Processed: %d/%d/%d %d:%d GMT                                       ",  
     318                        sprintf(tempstr, "Processed: %d/%d/%d %d:%d UTC                                       ",  
    309319                                (timeinfo->tm_mon + 1), 
    310320                                timeinfo->tm_mday, 
     
    314324                        ); 
    315325                        fprintf(COSMOS0fp, "%s\n", tempstr); 
     326                        break; 
     327                case 47: /*AQMS normally writes this comment, but we're not using AQMS*/ 
     328                        /* Example: | RcrdId: NC.72282711.NC.C031.HNE.01 */ 
     329                        sprintf(tempstr, "| RcrdId: %s.72282711.%s.%s.%s.%s ", /* FIX: If anyone knows what I can put here instead of 72282711, I'll do it!*/ 
     330                                ptrReq->net, ptrReq->net, ptrReq->sta, ptrReq->chan, ptrReq->loc); 
     331                        fprintf(COSMOS0fp, "%s\n", tempstr); 
     332                        break; 
     333                case 48: /*AQMS normally writes this comment, but we're not using AQMS*/ 
     334                        /*Example:|<SCNL>C031.HNE.NC.01    <AUTH> 2015/03/01 16:46:25.000*/ 
     335                        sprintf(tempstr, "|<SCNL> %s.%s.%s.%s    <AUTH> %d/%d/%d %d:%d:%d.000", 
     336                                ptrReq->sta, ptrReq->chan, ptrReq->net, ptrReq->loc, 
     337                                (timeinfo->tm_year + 1900),     (timeinfo->tm_mon + 1), timeinfo->tm_mday, 
     338                                timeinfo->tm_hour, timeinfo->tm_min, timeinfo->tm_sec); 
     339                        fprintf(COSMOS0fp, "%s", tempstr); 
    316340                        break; 
    317341                default:  
Note: See TracChangeset for help on using the changeset viewer.