FUNCTION UITMAG,source_in$,date$,ts,INFO$=info$ source$ = source_in$ source$ = STRSPLIT(source$,':',/EXTRACT) CASE source$[0] OF 'query': retData = '^uitmag' ELSE : BEGIN startdate = TIME_DOUBLE(date$) enddate = startdate + (ts * 60.) fileName$ = source$[1] PRINT,fileName$ savFName$ = fileName$ + '.sav' files = FINDFILE(savFName$, count=numFiles) ;numFiles=0 IF numFiles EQ 0 THEN BEGIN PRINT,'No previous saveset file found; creating saveset file from ASCII file.' ;Load data in from file. OPENR, unit, fileName$, /GET_LUN header$ = '' ; READF,unit,header$ ; READF,unit,year,month,day,hour,min,second ; IF year LE 90 THEN year = year + 2000 ELSE year = year + 1900 ; READF,unit,header$ READF,unit,header$ READF,unit,header$ READF,unit,header$ READF,unit,header$ READF,unit,header$ c = 0LL sec = 0D bx = 0D by = 0D bz = 0D REPEAT BEGIN line$='' READF,unit,line$ dataVec = STRSPLIT(line$,'/: ',/EXTRACT) ;DD/MM/YYYY HH:MM:SS Dec Horiz Vert Incl Total ;0 1 2 3 4 5 6 7 8 9 10 IF (c NE 0 ) THEN retData = [[retData], [dataVec]] ELSE retData = dataVec ++c ; print,c,' ',line$ ENDREP UNTIL EOF(unit) FREE_LUN,unit ;Transform time vector into epoch time. ;time0 = TIMEYMDHMSTOEPOCH(LONG(year),LONG(month),LONG(day),LONG(hour),LONG(min),LONG(second)) time$ = $ retData[2,*] + '-' + $ retData[1,*] + '-' + $ retData[0,*] + '/' + $ retData[3,*] + ':' + $ retData[4,*] + ':' + $ retData[5,*] HELP,time$ timeVec = TIME_DOUBLE(time$) ; Bx = retData[7,*] * COS(retData[6,*] * !PI/180.) By = retData[7,*] * SIN(retData[6,*] * !PI/180.) Bz = retData[8,*] ;retData = TRANSPOSE([[timeVec], [Bx], [By], [Bz]]) retData = TRANSPOSE([timeVec, Bx, By, Bz]) ; ;Create a 4 column array containing: ; ; [TimeVector Bx By Bz] ; retData[*,0] = retData[*,0] + time0 ;Create Information Vector type$ = 'mag' name$ = STRUPCASE(source$[1]) + ' UiTromsoe Mag (' + source$[2] + ')' unit$ = 'nT' shortID$ = STRUPCASE(source$[1]) + ':' + source$[2] info$ = [source_in$, type$, name$, unit$, shortID$ ] SAVE,file=savFName$,retData,info$ ENDIF ELSE BEGIN RESTORE,file=savFName$ ENDELSE ;Ensure that the array only contains data within the specified time range. upper = WHERE(retData[*,0] GE startdate) retData = retData[upper,*] lower = WHERE(retData[*,0] LE enddate) retData = retData[lower,*] ;Select Magnetometer Component comp$ = STRMID(source$[2],0,1,/REVERSE_OFFSET) retData = comp(retData,comp$) END ENDCASE RETURN, retData END