pro plot_eit_entrance, filters, hc=hc, noprint=noprint, $
notime=notime, filein=filein, logscale=logscale
;+
; NAME:
; plot_eit_entrance
; PURPOSE:
; Plot the comparison between computed EIT entrance filter
; and Jean-Pierre Delaboudiniere measured values.
; CALLING SEQUENCE:
; plot_eit_entrance
; plot_eit_entrance,/hc ; Send to the PostScript printer
; HISTORY:
; 10-Jun-95, J. R. Lemen (LPARL), Written.
;-
;
;From CATURA@SAG.SPACE.LOCKHEED.COM Tue Jun 6 14:20:44 1995
;Received: from sxt3.space.lockheed.com by sxt7.space.lockheed.com; (5.65/1.1.8.2/21Feb95-1254PM)
; id AA00366; Tue, 6 Jun 1995 14:20:43 -0700
;Received: by sxt3.space.lockheed.com (5.65/DEC-Ultrix/4.3)
; id AA25029; Tue, 6 Jun 1995 14:20:12 -0700
;Date: Tue, 6 Jun 1995 14:20:40 -0700 (PDT)
;From: CATURA@SAG.SPACE.LOCKHEED.COM
;To: LEMEN@SAG.SPACE.LOCKHEED.COM
;Cc: CATURA@SAG.SPACE.LOCKHEED.COM
;Message-Id: <950606142040.20203fb0@SAG.SPACE.LOCKHEED.COM>
;Subject: EIT Filters
;Status: RO
;
;Hi Jim
;
;I talked to Forbes Powell. The polyimide formula is C22 H10 N2 O4 and the
;density is 1.4 g/cm3. The filter is 40A Al2O3 1460 A Al and 700 A polyimide.
;Delaboudiniere gave me results of their measurements on the sample from
;NRL as :
;
; lambda T
;
; 171 .30-.35
;
; 304 .04-.05
;
; 500 .01-.02
;
;
;Dick
;
; -------------
; Later Dick told me that Forbes Powell called back and corrected
; composition of polyimide:
; polyimide = C14H6N2O4
; --------------
Wmeas = [171,304,500]
Tmeas0 = [.30,.04,.01]
Tmeas1 = [.35,.05,.02]
if n_elements(filters) eq 0 then begin
dum = eit_eff_area(0,filein=filein,area_str=Area_str) ; Read the era*genx file
filters = area_str.filter
endif
clearplot
if keyword_set(hc) then begin
set_plot,'PS'
device,/land,/inch,xsize=7.5,ysize=6.,xoff=1.25,yoff=11.-(11.-7.5)/2.
!x.thick=3 & !y.thick=3 & !p.thick=3 & !p.charthick=3.6
font = 0
endif else begin
set_plot,'x' & font = -1
endelse
;-------------
; Plot 1:
;-------------
xtitle = 'Wavelength (A)'
ytitle = 'Transmission'
title = 'EIT Entrance Filter'
if keyword_set(logscale) then $
plot_oo,filters(0).wave,filters(0).trans, $
xtitle=xtitle,title=title, $
font=font,charsize=1.5,thick=2 $
else $
plot_oi,filters(0).wave,filters(0).trans, $
xtitle=xtitle,title=title, $
font=font,charsize=1.5,thick=2
errplot,Wmeas,Tmeas0,Tmeas1,width=.03
ans = '' & if !d.name ne 'PS' then read,'* Enter <cr> to continue: ',ans
if !d.name eq 'PS' then begin
if not keyword_set(notime) then begin
stime = !stime
text = string('Plot done at: ',stime,format='(2a)')
if n_elements(filein) ne 0 then $
text = text+string(' File: ',filein,format='(2a)') else $
text = text+string(' Version: ',filters(0).Version, format='(2a)')
xyouts,0.5,-0.05,/norm,alig=.5,text
endif
endif
;-------------
; Plot 2:
;-------------
if keyword_set(logscale) then $
plot_oo,filters(0).wave,filters(0).trans, $
xtitle=xtitle,title=title, $
font=font,charsize=1.5,thick=2,xr=[100,1000] $
else $
plot_oi,filters(0).wave,filters(0).trans, $
xtitle=xtitle,title=title, $
font=font,charsize=1.5,thick=2,xr=[100,1000]
errplot,Wmeas,Tmeas0,Tmeas1,width=1.1
if !d.name eq 'PS' then begin
if not keyword_set(notime) then begin
text = string('Plot done at: ',stime,format='(2a)')
if n_elements(filein) ne 0 then $
text = text+string(' File: ',filein,format='(2a)') else $
text = text+string(' Version: ',filters(0).Version, format='(2a)')
xyouts,0.5,-0.05,/norm,alig=.5,text
endif
pprint,/reset,noprint=noprint
font = -1
endif
clearplot & !x.thick=0 & !y.thick=0 & !p.thick=0 & !p.charthick=0
end
Last revised: - Wed May 9 21:45:08 2007- F. Auchère