set_plot,'x' ;read radiosonde data from cart if (!d.name eq 'PS') then device,xsize=6.5,ysize=5.5,$ yoffset=1.,xoffset=1.,/inches,/portrait if (!d.name eq 'X') then window, 0, retain=2, ysize=400, xsize=600 read, prompt = 'Please Input Day on Interest (yymmdd): ', begindate ;begindate=980814 date=long(begindate) ;for i=0,3 do begin colour_hp !P.background=1 ;date=long(begindate+i) sdate=strcompress(string(date(0)),/remove_all) filename='../aeri/netcdf/'+sdate+'C1.CDF' filename2='../aeri/netcdf/'+sdate+'C2.CDF' fid = ncdf_open(filename) result=ncdf_varid(fid,'wnum') if result ge 0 then begin print,'ARM DOE NETCDF TYPE' ncdf_varget, fid, 'base_time', base_time ncdf_varget, fid, 'time_offset', time_offset ncdf_varget, fid, 'wnum', wave1 ncdf_varget, fid, 'mean_rad', rads1 ncdf_close,fid time_convert,base_time+time_offset,date,time=time time_minute=floor((time-floor(time))*60) time_hour=floor(time)*100 timehhmm=time_hour+time_minute filename2=strmid(filename,0,strpos(filename,'C'))+'C2.CDF' fid2 = ncdf_open(filename2) ncdf_varget, fid2, 'base_time', base_time2 ncdf_varget, fid2, 'time_offset', time_offset2 ncdf_varget, fid2, 'wnum', wave2 ncdf_varget, fid2, 'mean_rad', rads2 ncdf_close,fid2 time_convert,base_time2+time_offset2,date2,time=time2 endif else begin print,'Wisconsin NETCDF TYPE' fid = ncdf_open(filename) ncdf_varget, fid, 'base_time', base_time ncdf_varget, fid, 'time_offset', time_offset ncdf_varget, fid, 'wnum1', wave1 ncdf_varget, fid, 'mean_rad', rads1 ncdf_close,fid time_convert,base_time+time_offset,date,time=time time_minute=floor((time-floor(time))*60) time_hour=floor(time)*100 timehhmm=time_hour+time_minute filename2=strmid(filename,0,strpos(filename,'C'))+'C2.CDF' fid2 = ncdf_open(filename2) ncdf_varget, fid2, 'base_time', base_time2 ncdf_varget, fid2, 'time_offset', time_offset2 ncdf_varget, fid2, 'wnum1', wave2 ncdf_varget, fid2, 'mean_rad', rads2 ncdf_close,fid2 time_convert,base_time2+time_offset2,date2,time=time2 endelse btspec,wave1,rads1,bt btspec,wave2,rads2,bt2 !P.multi=[0,1,1] plot, time, bt(321,*),color=0,xrange=[0,24],yrange=[150,350],$ xstyle=1, ystyle=1,title='Andros AERI Radiance Summary '+sdate,ytitle='Brightness Temperature(K)',$ xtitle='Time (UTC)', thick=2,charsize=.8 oplot, time, bt2(1473,*), thick=2,color=2 oplot, time, bt(964,*),line=2, thick=2,color=4 oplot,[1,2],[200,200],color=0 xyouts,2.3,200,'675 cm!e-1!n (15 um near Surface Temp)',charsize=.8,color=0 oplot,[1,2],[190,190],color=2 xyouts,2.3,190,'2510 cm!e-1!n (4 um Shortwave Microwindow)',charsize=.8,color=2 oplot,[1,2],[180,180],color=4,line=2 xyouts,2.3,180,'985 cm!e-1!n (10 um Longwave Microwindow)',charsize=.8,color=4 spawn, 'mkdir '+sdate tvlct,r,g,b,/Get giffile=tvrd() write_gif, sdate+'/'+sdate+'_aeri_sum.gif', giffile, r, g, b ;command=strcompress('mv ~/goes8gps/idl.gif ~/goes8gps/goes8gps_'+strcompress(string(i),/remove_all)+'.gif') ;spawn,command ans='y' ;print,'Print out?' ;read,ans if (ans eq 'y' or ans eq 'Y') then begin print,' **PRINTING**...' set_plot,'ps' device,yoffset=1,ysize=9,/inches plot, time, bt(321,*),color=0,xrange=[0,24],yrange=[150,350],$ xstyle=1, ystyle=1,title='Andros AERI Channel 1 '+sdate,ytitle='Brightness Temperature(K)',$ xtitle='Time (UTC)', thick=2,charsize=.8 oplot, time, bt(373,*),color=0,line=1, thick=2 oplot, time, bt(964,*),color=0,line=2, thick=2 oplot,[1,2],[330,330],color=0 xyouts,2.3,330,'675 cm!e-1!n',charsize=.8,color=0 oplot,[8,9],[330,330],line=1,color=0 xyouts,9.3,330,'700 cm!e-1!n',charsize=.8,color=0 oplot,[15,16],[330,330],color=0,line=2 xyouts,16.3,330,'985 cm!e-1!n',charsize=.8,color=0 plot, time, bt2(1000,*),color=0,xrange=[0,24],yrange=[150,350],$ xstyle=1, ystyle=1,title='Andros AERI Channel 2 '+sdate,ytitle='Brightness Temperature(K)',$ xtitle='Time (UTC)', thick=2,charsize=.8 oplot, time, bt2(1027,*),color=0,line=1, thick=2 oplot, time, bt2(1473,*),color=0,line=2, thick=2 oplot,[1,2],[330,330],color=0 xyouts,2.3,330,'2282 cm!e-1!n',charsize=.8,color=0 oplot,[8,9],[330,330],line=1,color=0 xyouts,9.3,330,'2295 cm!e-1!n',charsize=.8,color=0 oplot,[15,16],[330,330],line=2,color=0 xyouts,16.3,330,'2510 cm!e-1!n',charsize=.8,color=0 device,/close spawn,'mv idl.ps '+sdate+'/'+sdate+'radsum.ps' set_plot,'x' endif ;endfor end