Permalink
Browse files

removed a bunch of little bugs. all working now.

  • Loading branch information...
1 parent db2c946 commit 763158b972175deff2b3aa4df43945af349a8d98 @skendrew committed Mar 30, 2010
View
2 .project
@@ -1,6 +1,6 @@
<?xml version="1.0" encoding="UTF-8"?>
<projectDescription>
- <name>METISpublic</name>
+ <name>midIR_sensitivity</name>
<comment></comment>
<projects>
</projects>
View
28 README
@@ -3,7 +3,7 @@ MIDIR_SENSITIVITY REPOSITORY
What is it?
-----------
-This repository contains idl code to calculate the sensitivity of a ground-based mid-infrared instrument for astronomy. The code was written during 2008-2009 for the Phase A study of the instrument METIS, the Mid-Infrared E-ELT Imager and Spectrograph, for the 42-m European Extremely Large Telescope. The code and all input parameters are highly tailored for the particular design parameters of the E-ELT and METIS, however, the programme is structured in such a way that the parameters can easily be adjusted for a different system, or alternative input files used.
+This repository contains idl code to calculate the sensitivity of a ground-based mid-infrared instrument for astronomy. The code was written during 2008-2009 for the Phase A study of the instrument METIS (http://www.strw.leidenuniv.nl/metis), the Mid-Infrared E-ELT Imager and Spectrograph, for the 42-m European Extremely Large Telescope. The code and all input parameters are highly tailored for the particular design parameters of the E-ELT and METIS, however, the programme is structured in such a way that the parameters can easily be adjusted for a different system, or alternative input files used.
Basic features of the E-ELT and METIS are:
- a 42-m aperture
@@ -32,13 +32,23 @@ The package contains 3 key files for sensitivity calculations in METIS' 3 main m
3. metis_ifu_sn.pro -> High-resolution integral field spectroscopy (R~100,000 in LM, 50,000 in N band; image slicer design)
- IF spectrosocpy specific parameter file: metis_sens_ifu_param.pro
-More parameter files:
+
+Dependencies
+------------
+
+Parameter files:
- metis_sens_param.pro: general input parameters related to telescope and atmosphere. Contains parameters for setting file locations to be adjusted for your own system
- metis_pre_param.pro: input parameters related to the pre-optics of METIS
-- atm/paranal_*.dat: atmospheric profile files for a Paranal-like site, containing both transmission and radiance info. Profiles have been pre-convolved to resolutions appropriate for the different modes of observation: 100 for imaging, 3000 for slit spectroscopy, 50000 for N band IFU spectrosocpy and 100000 for LM-band IFU spectroscopy
-- atm/high_and_dry_*.dat: dito Paranal for a high & dry site at 5000 m
-- ee_files/: this folder contains 2 lookup tables for AO-corrected encircled energies. (i) EnsquaredEnergy_AO_LMN_V10.dat is for single conjugate AO and is calculated with AO simulation tool PAOLA. Please refer to the tech note by Kendrew for detailed input details; (ii) atlas_eeav.dat is the same for the laser tomography AO mode, provided by the team designing the LTAO module for E-ELT at ONERA, France. These values are encircled energies instead of ensquared, so take care to compare. See tech note for more info.
-- profiles/: contains a number of input files loosely grouped for being "profiles" of sorts.
+
+Atmospheric profiles:
+- atm/paranal_*.dat: atmospheric profile files for a Paranal-like site, containing both transmission and radiance info. Profiles have been pre-convolved to resolutions appropriate for the different modes of observation: 100 for imaging, 3000 for slit spectroscopy, 50000 for N band IFU spectroscopy and 100000 for LM-band IFU spectroscopy
+- atm/high_and_dry_*.dat: dito Paranal for a high & dry site at 5000 m
+
+Encircled energy lookup tables:
+- ee_files/: this folder contains 2 lookup tables for AO-corrected encircled energies. (i) EnsquaredEnergy_AO_LMN_V10.dat is for single conjugate AO and is calculated with AO simulation tool PAOLA. Please refer to the tech note by Kendrew for detailed input details; (ii) atlas_eeav.dat is the same for the laser tomography AO mode, provided by the team designing the LTAO module for E-ELT at ONERA, France. These values are encircled energies instead of ensquared, so take care to compare. See tech note for more info.
+
+Profiles:
+The profiles/ folder contains a number of input files loosely grouped for being "profiles" of sorts.
-> bandic_prof.dat: input for calculation of the transmission profile of the dichroic that splits LM and N band radiation
-> filters.dat: filter transmission details for imager mode
-> grisms.dat: grisms info for low-resolution spectroscopy
@@ -57,8 +67,10 @@ Auxiliary functions:
Requirements
------------
+This code was written using IDL 7.0 under Windows XP SP3 and should be compatible with most reasonably current versions of IDL. I've also used it regularly under linux, which requires some minor modifications to the directory structure. The input directories are specified in the metis_sens_param.pro file.
-This code was written using IDL 7.0 under Windows XP and should be compatible with most reasonably current versions of IDL. I've also used it regularly under linux, which requires some minor modifications to the directory structure. The input directories are specified in the metis_sens_param.pro file.
+The IDL Astronomy Library is referenced throughout - you can download it here:
+http://idlastro.gsfc.nasa.gov/ftp/
Who wrote it?
@@ -68,7 +80,7 @@ This code was written by Sarah Kendrew, Postdoc in Astronomy Instrumentation at
Acknowledgment
--------------
If you use this code, either for METIS-specific calculations or for your own project, please cite:
-Kendrew et al., Proc. SPIE, in preparation(2010)
+Kendrew et al., Proc. SPIE, in preparation (2010)
-- this paper will be presented at SPIE San Diego 2010 so I'll update the reference asap.
If you use the code extensively, I'd like to hear about it and would appreciate being added as a co-author to any resulting publications.
View
2 metis_ifu_param.pro
@@ -1,5 +1,5 @@
-;pro metis_sens_ifu_param
+;pro metis_ifu_param
;parameters file for sensitivity calculations
;
;provisional temperatures
View
BIN metis_ifu_sensitivity.zip
Binary file not shown.
View
1 metis_ifu_sn.pro
@@ -163,7 +163,6 @@ endif
print, 'Atmosphere data read'
;perform some unit conversions on the emissivity values (from photons/s/um/m^2/arcsec^2 to W/cm^2/sr/um)
-;em=em*(h*c/l)*4.25d10. UNITS WOMAN, UNITS! PAY ATTENTION!
em_l=em*(h*c/(l*1d-6))*1d-4*4.25d10
; clip these arrays to wavelength range of interest
View
2 metis_img_param.pro
@@ -1,5 +1,5 @@
-;pro metis_sens_img_param
+;pro metis_img_param
;parameters file for sensitivity calculations for the IMAGER
;
View
BIN metis_img_sensitivity.zip
Binary file not shown.
View
2 metis_img_sn.pro
@@ -247,8 +247,6 @@ print, 'mean throughput of instrument = ', mean(thru)
dqe=metis_dqe(l, band)
print, 'mean DQE in filter band = ', mean(dqe)
-plot, l, dqe
-
; calculate efficiency
eff=mean(thru)*tel_thru*mean(dqe)*pcg
print, 'total efficiency including DQE, pcg, telescope, and instrument = ', eff
View
BIN metis_lres_sensitivity.zip
Binary file not shown.
View
2 metis_mres_sn.pro
@@ -323,7 +323,7 @@ if keyword_set(extended) then begin
endelse
endif else begin
print, 'sensitivity at s/n 10 in 1 hr = ', minsig[line]*1000., ' microJy'
- print, 'line sensitivity at s/n 10 in 1 hr = ', minsig_line[line]*1d19, '*1d-19 W/m2/arcsec^2'
+ print, 'line sensitivity at s/n 10 in 1 hr = ', minsig_line[line]*1d19, '*1d-19 W/m2'
endelse
View
23 metis_sens_ifu_param.pro
@@ -1,23 +0,0 @@
-
-;pro metis_sens_ifu_param
-;parameters file for sensitivity calculations
-
-;provisional temperatures
-hspec_temp_lm=65. ; Kelvin
-hspec_temp_n=25.
-
-;SPO path
-spo_no_surf=3 ;7 ;# of "normal" reflective optics, e.g. re-imaging & fold mirrors
-spo_ref_surf=0.99
-;spodic_trans=0.85
-slicer_ref=0.99 ;to be multiplied by slicing loss, calculated in code
-slices_lm=24.
-slices_n=15.
-
-;SMO path
-smo_no_surf=7 ;6 ;# "normal" surfaces
-smo_ref_surf=0.99
-no_grating=1 ;2
-grating_eff=0.8^no_grating
-pix_spat=2.d ;no. pixels per spatial res element
-pix_spec=2.5 ;dito spectral
View
22 metis_sens_img_param.pro
@@ -1,22 +0,0 @@
-
-;pro metis_sens_img_param
-;parameters file for sensitivity calculations for the IMAGER
-
-
-
-no_surf=15 ;number of cold surfaces. 4 pre-optics+3 de-rotator+3 collimator+2 folding mirrors+3 camera
-ref_surf=0.99 ;reflectance per cold surface
-window_trans=0.95
-no_dichroic=1
-dichroic_trans=0.8^no_dichroic
-;pix_spat=3.d ;no. pixels per spatial res element. assume nyquist sampled.
-thru=ref_surf^no_surf*window_trans*dichroic_trans
-
-pix_scale_lm=8.6 ; mas/pixel
-pix_scale_nq=17.2 ; mas/pixel
-;pix_spec=2. ; for testing - REMOVE
-
-print, 'imager throughput = ', thru
-
-;end
-
View
2 metis_sens_param.pro
@@ -4,7 +4,7 @@
;instrument parameters moved to metis_sen_img_param and metis_sens_ifu_param.pro
; file locations
-def_dir='c:/metis/code_public/'
+def_dir='c:/metis/midIR_sensitivity/'
atm_dir=def_dir+'atm/' ; location of atmospehric profiles
ee_dir=def_dir+'ee_files/' ; location of encircled energy tables
profiles_dir=def_dir+'profiles/' ; location of miscellaneous input data, mainly optical profiles e.g. throughputs, resolutions etc
View
10 metis_thru.pro
@@ -16,19 +16,15 @@ function metis_thru, mode, band, x, POL=pol
;************************************
;*CONSTANTS *
;************************************
-;pre_surf=9 ;no. of reflective optics in the pre-optics
-;pre_ref=0.99 ;reflectivity per pre-optics mirror
-;pre_thru=pre_ref^pre_surf
-;ar_factor=1.3
-;thru=1.
+
;wollaston_thru=0.9 ; throughput of the wollaston prism in METIS
; CONSTANTS ALL READ IN VIA PARAMETER FILES
@metis_sens_param
@metis_pre_param
-@metis_thermal_ifu_param
-@metis_thermal_img_param
+@metis_ifu_param
+@metis_img_param
;**********************************************************************
;* WINDOW TRANSMISSION CURVE also used for prism transmission for IFU*
View
5 metis_winbgr.pro
@@ -5,8 +5,9 @@ function metis_winbgr, x, temp
n=2.4028 ;real part of refractive index @ 10.6 um
r=(n-1)/(n+1)
win_thick=1.2 ;window thickness [cm]
- ;read in window absortion coefficient and calculate emissivity
-readcol, 'c:/metis/thermal/metis_thermal/znse_abs.dat', format='d,d', abs_l, alpha, skipline=4, /silent
+;read in window absortion coefficient and calculate emissivity
+absfile=profiles_dir+'znse_abs.dat'
+readcol, absfile, format='d,d', abs_l, alpha, skipline=4, /silent
win_em=((1-r^2)*(1-exp(-alpha*win_thick)))/(1-(r^2*exp(-alpha*win_thick))) ;r and d defined in metis_thermal_pre_param
win_em=interpol(win_em, abs_l, x)
;win_em[where(l gt 20.)]=1. ;set the emissivity to level off after 20 micron at 100%

0 comments on commit 763158b

Please sign in to comment.