########set a WISE filter model
from sedfitter.filter import Filter
from astropy import units as u
from sedfitter import fit
from sedfitter.extinction import Extinction
import numpy as np
from sedfitter import plot
from sedfitter import filter_output
from sedfitter import write_parameters, write_parameter_ranges
from astropy import constants as const
f1 = Filter()
f1.name = "ws1" ## reference name
f1.central_wavelength = 3.3526 * u.micron ## effective wavelength
f1.nu = const.c.value/ np.loadtxt('RSR-W1.txt',usecols=[0])/1e-6*u.Hz ## observed bandpass frequency
f1.response = np.loadtxt('RSR-W1.txt',usecols=[1]) ## filter response curve
f1.normalize() ## normalization
from sedfitter.convolve import convolve_model_dir
print('f1 started')
convolve_model_dir('/Users/ygong/soft/models_r06', [f1]) ## convolve all available models