def run(self):
if not self._kill:
self.dataobj = casio()
self.dataobj.readinp(self.parent._filename)
self.dataobj.readcax()
self.emit(SIGNAL('progressbar_update(int)'),90)
if not self._kill:
self.dataobj.calcbtf()
if not self._kill:
self.parent.dataobj = self.dataobj
def load_pickle(self,filename):
self.statusBar().showMessage('Importing data from %s' % filename, 2000)
self.dataobj = casio()
self.dataobj.loadpic(filename)
self.init_pinobjects()
#fuetype = 'SVEA-96'
#self.dataobj.btf = btf(self.dataobj,fuetype)
#self.setpincoords()
#self.draw_fuelmap()
#self.set_pinvalues()
# Update case number list box
ncases = len(self.dataobj.cases)
for i in range(1,ncases+1):
self.case_cbox.addItem(str(i))
self.connect(self.case_cbox, SIGNAL('currentIndexChanged(int)'), self.fig_update)
self.fig_update()
SK = min([SK1,SK2])
# Sum weighting factors
SWJ = (WP[i-1,j] + WP[i+1,j] + WP[i,j-1] + WP[i,j+1])*WJ # Side rods
SWK = (WP[i-1,j-1] + WP[i+1,j-1] + WP[i-1,j+1] + WP[i+1,j+1])*WK # Diagonal rods
DOW[i-1,j-1] = (RP[i,j] + SJ + SK)/(1.0 + SWJ + SWK)*np.sqrt(Ntotrods/float(Nhotrods)) + AC[i-1,j-1]
return DOW
def calc_btf(fuetype,POW3):
naxial_nodes = 25
DOW = np.zeros((naxial_nodes,POW3.shape[1],POW3.shape[2]))
# Calculate lattice K-factors for all nodes
for z in range(naxial_nodes):
DOW[z,:,:] = rfact_axial(fuetype,POW3[z,:,:])
# Integrate along z-direction to get pinwise R-factors
DOX = sum(DOW,0)/naxial_nodes
return DOX
if __name__ == '__main__':
casobj = casio()
casobj.loadpic('caxfiles_a10xm.p')
POW3 = casobj.pow3(casobj)
