error while executing gtmodel during enrico_sed operation

[basuparth]

I was running enrico_sed to generate flux points and I ran into the following error. Prior to gtmodel all the other steps were executed successfully by enrico.

Error report

time -p gtmodel srcmaps=/home/parth/SINP/1ES1218+304/1ES1218+304_LAT_Analysis_PowerLaw2_srcMap.fits srcmdl=/home/parth/SINP/1ES1218+304/1ES1218+304_PowerLaw2_LAT_Analysis_out.xml outfile=/home/parth/SINP/1ES1218+304/1ES1218+304_LAT_Analysis_PowerLaw2_ModelMap.fits irfs="P8R3_SOURCE_V2" evtype=3 expcube=/home/parth/SINP/1ES1218+304/1ES1218+304_LAT_Analysis_ltCube.fits bexpmap=/home/parth/SINP/1ES1218+304/1ES1218+304_LAT_Analysis_BinnedMap.fits convol=yes resample=yes rfactor=2 outtype="CMAP" psfcorr=yes edisp=no chatter=2 clobber=yes debug=no gui=no mode="ql"

WARNING: version mismatch between CFITSIO header (v3.43) and linked library (v3.41).

Using evtype=3 (i.e., FRONT/BACK irfs)
real 61.48
user 59.56
sys 0.77
Traceback (most recent call last):
File "/home/parth/SINP/enrico/bin/enrico_sed", line 53, in
sed(config,infile)
File "/home/parth/SINP/enrico/bin/enrico_sed", line 13, in sed
run(infile)
File "/home/parth/SINP/enrico/enrico/RunGTlike.py", line 196, in run
FitRunner.ModelMap(outXml)
File "/home/parth/SINP/enrico/enrico/fitmaker.py", line 473, in ModelMap
self.obs.ModelMaps(xml)
File "/home/parth/SINP/enrico/enrico/gtfunction.py", line 424, in ModelMaps
utils.SubtractFits(self.cmapfile,self.ModelMap,self.Configuration)
File "/home/parth/SINP/enrico/enrico/utils.py", line 120, in SubtractFits
fits.writeto(abs_diff_file, data1 - data2, head)
File "/home/parth/anaconda2/envs/fermi/lib/python2.7/site-packages/astropy/utils/decorators.py", line 508, in wrapper
return function(*args, **kwargs)
File "/home/parth/anaconda2/envs/fermi/lib/python2.7/site-packages/astropy/io/fits/convenience.py", line 433, in writeto
checksum=checksum)
File "/home/parth/anaconda2/envs/fermi/lib/python2.7/site-packages/astropy/utils/decorators.py", line 508, in wrapper
return function(*args, **kwargs)
File "/home/parth/anaconda2/envs/fermi/lib/python2.7/site-packages/astropy/io/fits/hdu/base.py", line 377, in writeto
checksum=checksum)
File "/home/parth/anaconda2/envs/fermi/lib/python2.7/site-packages/astropy/utils/decorators.py", line 508, in wrapper
return function(*args, **kwargs)
File "/home/parth/anaconda2/envs/fermi/lib/python2.7/site-packages/astropy/io/fits/hdu/hdulist.py", line 872, in writeto
self.verify(option=output_verify)
File "/home/parth/anaconda2/envs/fermi/lib/python2.7/site-packages/astropy/io/fits/verify.py", line 121, in verify
raise VerifyError('\n' + '\n'.join(messages))
astropy.io.fits.verify.VerifyError:
Verification reported errors:
HDU 0:
NAXISj keyword out of range ('NAXIS3' when NAXIS == 2)
Note: astropy.io.fits uses zero-based indexing.

[davidsanchez]

Hi
I am aware of this and I have no solution yet
cheers

[mireianievas]

Is this already fixed? I see a head.remove('NAXIS3') statement in the code and I have never encountered this issue.