forked from ehsteve/microflare_dem
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathget_aia_file_list.pro
48 lines (37 loc) · 1.29 KB
/
get_aia_file_list.pro
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
FUNCTION get_aia_file_list, dir, WAVE_LIST = wave_list ,fileset = fileset
; PURPOSE: Given a directory get all of the aia cutout files in the directory.
;
; KEYWORDS:
; wave_list -
;
; WRITTEN: Steven Christe (8-Oct-2011)
default, wave_list, ['131','171','193','211','335','94']
default, fileset, 'ssw_cutout'
nwave =n_elements(wave_list)
nfiles = fltarr(nwave)
; search for all of the files
FOR i = 0, nwave-1 DO BEGIN
IF fileset EQ 'ssw_cutout' THEN files = file_search(dir + fileset + '*' + wave_list[i] + '_.fts')
IF fileset EQ 'AIA' THEN files = file_search(dir + fileset + '*0' + wave_list[i] + '.fits')
nfiles[i] = n_elements(files)
ENDFOR
; what to do if number of files is not the same?
; right now just get the minimum number of files
file_list = strarr(min(nfiles), nwave)
FOR i = 0, nwave-1 DO BEGIN
IF fileset EQ 'ssw_cutout' THEN BEGIN
files = file_search(dir + fileset + '*' + wave_list[i] + '_.fts')
nfiles[i] = n_elements(files)
times = anytim(aiacutout_to_time(files))
ENDIF
IF fileset EQ 'AIA' THEN BEGIN
files = file_search(dir + fileset + '*0' + wave_list[i] + '.fits')
nfiles[i] = n_elements(files)
times = anytim(aiaprep_to_time(files))
ENDIF
s = sort(times)
files = files[s]
file_list[*,i] = files[0:min(nfiles)-1]
ENDFOR
RETURN, file_list
END