RunFileListTuples = [
              (['/AstroData/GALPROP/ics_isotropic_skymap_49_700405'],'ICS1-1000GeV.060529.fits'),
              (['/AstroData/GALPROP/pion_decay_skymap_49_700405'],'Pion1-1000GeV.060529.fits'),
              (['/AstroData/GALPROP/bremss_skymap_49_700405'],'Bremss1-1000GeV.060529.fits')
                  ]
ExposrFile = '/AstroData/exposr_g1-9_g004.fits'
ExposrHistory = ExposrFile
WantEBand = [1000.,1000000.]
WhichExp3rdBin = 1
WantExpRebin = [6,4]	# Usual exposre number of bins is 2x GALPROP
WantFlxRebin = [3,2]    # Rebinning to 120 x 90 prior to padding -- for speed
IntTimeLst = [806.,242.,73.,22.,7.]	# In "GALPROP" units (1500 ???).

#OutFileList = ['ModelFlux_SpherePadded.060521.fits', 'Exposure_SpherePadded.060521.fits', 'ModelData_SpherePadded.060521.fits', 'PoissDatons_Matched.060521.fits']

basic_string1a = 'ModelFluxRebinnd'
basic_string1b = 'ModelFluxPadded'

#from numarray import * ## Changed to numpy Aug 2007
from numpy import *
import pyfits
import random
from simpler_skymap_datons import *
from run_simpler_skymap_datons import *


#
# Beginning 1st module:
# 1/ Start accumulating flux model components:
EBandString = ''
for GalpropFileList,NameMeme in RunFileListTuples:
    OutModelFluxFitsFile = [basic_string1a+NameMeme,basic_string1b+NameMeme]
    model_flux_HDUpair = SkymapModelFluxHDU(GalpropFileList,WantEBand,WantFlxRebin,EBandString,OutModelFluxFitsFile)
#

# -------------------------------------------------------------------------------------