;PRO GET_MASS ; ;PURPOSE: This script will open up the MASS files containing the CN2 values ; and dig out the values function return_mass,min1,hr1,date1,min2,hr2,date2 Mass_Dir='/nfs/slac/g/ki/ki08/lsst/CPanalysis/MASS/' Mass_files=['orig_turb_data0', 'orig_turb_data1', 'orig_turb_data2',$ 'orig_turb_data4', 'orig_turb_data8','orig_turb_data16'] CN2=dblarr(13,1553) for i=0, N_elements(Mass_files)-1 do begin readcol,Mass_Dir+Mass_files(i),date,time,CN2_temp,format='a,a,f' for j=0, N_elements(time)-1 do begin time_arr=strsplit(time(j),':',/extract) date_arr=strsplit(date(j),'/',/extract) CN2(i,j)=CN2_temp(j) CN2(6,j)=date_arr(0) ;Year CN2(7,j)=date_arr(1) ;Month CN2(8,j)=date_arr(2) ;Date CN2(9,j)=time_arr(0) ;Hour CN2(10,j)=time_arr(1) ;Minute CN2(11,j)=time_arr(2) ;Second juldate,[CN2(6,j),CN2(7,j),CN2(8,j),CN2(9,j),CN2(10,j),CN2(11,j)],jd CN2(12,j)=double(jd) endfor endfor if hr1 ne hr2 then begin region=where(((CN2(10,*) gt min1) and (CN2(9,*) eq hr1) and $ (CN2(8,*) eq date1))$ or ((CN2(10,*) lt min2) and (CN2(9,*) eq hr2) and $ (CN2(8,*) eq date2))) endif else begin region=where(((CN2(10,*) gt min1) and (CN2(9,*) eq hr1) and $ (CN2(8,*) eq date1))$ and ((CN2(10,*) lt min2) and (CN2(9,*) eq hr2) and $ (CN2(8,*) eq date2))) endelse if region(0) ne -1 then begin startreg=region(0) stopreg=region(n_elements(region)-1) return,double(CN2(*,startreg:stopreg)) endif else begin b=fltarr(13,1) return,b endelse end