different ir_tol for 1 camera

[mohammad safarzadeh]

this is the output I get from this code:

no variation with that parameter seems something is wrong?

enzo_09=pyfits.open('enzo_broadband_0.001_ir_on_0.9_tol_npix_1000.fits')

enzo_05=pyfits.open('../npix_1000_ql_0.5_lores/enzo_broadband_0.001_ir_on_0.5_tol_npix_1000.fits')

enzo_03=pyfits.open('../npix_1000_ql_0.3_lores/enzo_broadband_0.001_ir_on_0.3_tol_npix_1000.fits')

L=enzo_09[5].data['lambda']

F09=enzo_09[5].data['L_lambda_ir1']

F05=enzo_05[5].data['L_lambda_ir1']

F03=enzo_03[5].data['L_lambda_ir1']

plt.figure()

ax=plt.gca()

ax.set_yscale('log')

ax.set_xscale('log')

ax.set_xlim(0,0.001)

plt.plot(L,L*F03)

plt.plot(L,L*F05)

plt.plot(L,L*F09)

plt.savefig('spectrum_ir_on_all_cameras_log_ldotlnu_cam1.png')

while this is the header file :

In [26]: enzo_09.info()

Filename: enzo_broadband_0.001_ir_on_0.9_tol_npix_1000.fits

No.    Name         Type      Cards   Dimensions   Format

0    PRIMARY     PrimaryHDU      24   ()           uint8   

1    BROADBAND   BinTableHDU     21   0R x 0C      []   

2    CAMERA0-PARAMETERS  ImageHDU        33   (1000, 1000)   float32   

3    CAMERA1-PARAMETERS  ImageHDU        33   (1000, 1000)   float32   

4    GADGET      BinTableHDU     93   0R x 0C      []   

5    INTEGRATED_QUANTITIES  BinTableHDU     59   120R x 12C   [D, D, D, D, D, D, D, D, D, D, D, D]   

6    MAKEGRID    BinTableHDU     11   0R x 0C      []   

7    MCRX        BinTableHDU     62   0R x 0C      []   

8    SCATTERING_LAMBDAS  BinTableHDU     35   0R x 5C      [J, D, L, D, J]   

9    SFRHIST     BinTableHDU     68   0R x 0C      []   

10   STELLARMODEL  BinTableHDU     60   0R x 0C      []   

11   FILTERS     BinTableHDU     41   22R x 12C    [30A, D, D, D, D, D, D, D, D, D, D, D]   

12   CAMERA0-BROADBAND  ImageHDU        49   (1000, 1000, 22)   float32   

13   CAMERA1-BROADBAND  ImageHDU        49   (1000, 1000, 22)   float32   

14   CAMERA0-BROADBAND-NONSCATTER  ImageHDU        49   (1000, 1000, 22)   float32   

15   CAMERA1-BROADBAND-NONSCATTER  ImageHDU        49   (1000, 1000, 22)   float32  

and 

In [37]: enzo_09[5].header

Out[37]: XTENSION= 'BINTABLE'           / binary table extension                         

BITPIX  =                    8 / 8-bit bytes                                    

NAXIS   =                    2 / 2-dimensional binary table                     

NAXIS1  =                   96 / width of table in bytes                        

NAXIS2  =                  120 / number of rows in table                        

PCOUNT  =                    0 / size of special data area                      

GCOUNT  =                    1 / one data group (required keyword)              

TFIELDS =                   12 / number of fields in each row                   

EXTNAME = 'INTEGRATED_QUANTITIES' / name of this binary table extension         

COMMENT This HDU contains spatially integrated SED's of the various quantities  

COMMENT For the camera quantities, the luminosity is calculated from the camera 

COMMENT flux by multiplying by 4*pi*d^2.  It's not a fundamental quantity.      

TTYPE1  = 'lambda  '           / label for field                                

TFORM1  = 'D       '           / format of field                                

TUNIT1  = 'm       '                                                            

TTYPE2  = 'lambda_intensity'   / label for field                                

TFORM2  = 'D       '           / format of field                                

TUNIT2  = 'm       '                                                            

TTYPE3  = 'L_lambda'           / label for field                                

TFORM3  = 'D       '           / format of field                                

TUNIT3  = 'W/m     '                                                            

TTYPE4  = 'L_lambda_absorbed'  / label for field                                

TFORM4  = 'D       '           / format of field                                

TUNIT4  = 'W/m     '                                                            

TTYPE5  = 'L_lambda_nonscatter0' / label for field                              

TFORM5  = 'D       '           / format of field                                

TUNIT5  = 'W/m     '                                                            

TTYPE6  = 'L_lambda_scatter0'  / label for field                                

TFORM6  = 'D       '           / format of field                                

TUNIT6  = 'W/m     '                                                            

TTYPE7  = 'L_lambda_ir0'       / label for field                                

TFORM7  = 'D       '           / format of field                                

TUNIT7  = 'W/m     '                                                            

TTYPE8  = 'L_lambda_nonscatter1' / label for field                              

TFORM8  = 'D       '           / format of field                                

TUNIT8  = 'W/m     '                                                            

TTYPE9  = 'L_lambda_scatter1'  / label for field                                

TFORM9  = 'D       '           / format of field                                

TUNIT9  = 'W/m     '                                                            

TTYPE10 = 'L_lambda_ir1'       / label for field                                

TFORM10 = 'D       '           / format of field                                

TUNIT10 = 'W/m     '                                                            

HIERARCH nlambda_intensity = 25                                                 

HIERARCH L_bol_grid = 3.51431616016778E+39 / [W] Bolometric source luminosity   

HIERARCH L_bol_absorbed = 3.07288730727783E+40 / [W] Bolometric luminosity absor

L_SCAT0 = 6.61819869361711E+38 / [W] Bolometric luminosity in  CAMERA0          

L_NS0   = 3.36079181957006E+39 / [W] Bolometric luminosity in  CAMERA0-NONSCATTE

L_IR0   = 2.48716201322247E+39 / [W] Bolometric luminosity in  CAMERA0-IR       

L_OUT0  = 3.10557232792972E+39 / [W] Total bolometric luminosity in camera 0, in

TTYPE11 = 'L_lambda_out0'      / label for field                                

TFORM11 = 'D       '           / format of field                                

TUNIT11 = 'W/m     '                                                            

L_SCAT1 = 5.62256168724911E+38 / [W] Bolometric luminosity in  CAMERA1          

L_NS1   =  3.2183044620205E+39 / [W] Bolometric luminosity in  CAMERA1-NONSCATTE

L_IR1   = 1.71525453525165E+39 / [W] Bolometric luminosity in  CAMERA1-IR       

L_OUT1  = 2.24649597925075E+39 / [W] Total bolometric luminosity in camera 1, in

TTYPE12 = 'L_lambda_out1'      / label for field                                

TFORM12 = 'D       '           / format of field                                

TUNIT12 = 'W/m     '                              

+ how come HIERARCH L_bol_absorbed  is greater than HIERARCH L_bol_grid ?