def _storeData(self, s=None):
"""
stores the data. The strategy is to store it all, it has to be processed later.
The purpose of this function is to standardize the column "design", so the numerous data files are easily
accessed and generic scripts can be written for them.
s has to be given if this is the first time the method is called and an output should be created.
"""
e=self.excelOut
if not e: #start a new instance.
if not s:
raise Exception('simulation instance has to be provided at first call for _storeData')
e=ExcelOutput(template='sim/planting/template.xls', out=self.filename)
self.excelOut=e
e.changeSheet(0)
#set all the parameters
paramRow=3
type=copy.copy(s.m.type) #we need to shift 2aXh because of new directions...
e.modify(paramRow, 0, s.m.type) #machine type
e.modify(paramRow, 1, s.G.terrain.ttype) #terrain type
e.modify(paramRow, 2, s.G.terrain.stumpFile) #terrain type
e.modify(paramRow, 3, s.G.terrain.stumpsPerH) #stumps/ha
e.modify(paramRow, 4, s.G.terrain.groundModel)
e.modify(paramRow, 5, s.G.terrain.blockQuota*100) #percent
e.modify(paramRow, 6, s.G.terrain.boulderFreq)
e.modify(paramRow, 7, s.G.terrain.meanBoulderV*1000)
e.modify(paramRow, 9, s.m.craneMaxL)
if s.m.type[0:2]=='2a':
e.modify(paramRow, 10, s.m.craneIntersect)
e.modify(paramRow, 11, s.m.angleLim[0]*180/pi)
e.modify(paramRow, 12, s.m.angleLim[1]*180/pi)
if len(s.m.pDevs[0].plantHeads)==2: e.modify(paramRow, 13, s.m.pDevs[0].plantSepDist*100)
e.modify(paramRow, 14, s.m.pDevs[0].plantHeads[0].width*100)
e.modify(paramRow, 15, s.m.stockingRate)
e.modify(paramRow, 16, s.m.plantMinDist)
e.modify(paramRow, 17, s.m.times['switchFocus'])
e.modify(paramRow, 18, str(s.G.automatic['automove']))
e.modify(paramRow, 19, str(s.G.automatic['mound']))
e.modify(paramRow, 20, s.m.velocities['radial'])
e.modify(paramRow, 21, str(self.bestCase))
e.modify(paramRow, 22, str(self.worstCase))
e.changeSheet(1) #the time-data sheet
tdrow=2+self.sims
for col,val in enumerate(self.timeData):
print tdrow, col, val
e.modify(tdrow,col,round(val,2)) #writes the time data. Request from BTE to round off
e.save()
return e
def __init__(self,it=1): self.G=globalVar() self.G.areaPoly=[(0,0), (25,0), (25,40), (0,40)] self.G.terrain=Terrain(G=self.G) self.folder='outputFiles/NewThinning_2012' today=datetime.date.today() self.filename=self.folder+'/'+'ThinningWawoBundler_'+'%s'%(string.join([string.strip(a) for a in string.split(time.ctime(time.time()))], "_"))+'.xls' e=ExcelOutput(template='sim/thinning/template.xls', out=self.filename) self.Paramrow=1 for treeFile in self.G.terrain.thinningFiles: self.G.terrain.treeFile=treeFile self.G.terrain.readTrees() for head in ['BC','convAcc']: for nCranes in [1,2]: for bundler in [False, True]: for twigCrack in [False, True]: for simNumber in range(1,it+1): G=copy.deepcopy(self.G) self.s=ThinningSim(G=G, vis=False, anim=False, head=head, nCranes=nCranes, bundler=bundler, twigCrack=twigCrack, observer=True) self.s.stats['machineConfig']=self.getConfig(head,nCranes,twigCrack,bundler) self.s.stats['simNumber']=simNumber self.s.stats['treeFile']=treeFile self.s.stats['noHarvTrees']=sum([len(b.trees) for b in G.terrain.piles]) self.s.stats['noCraneCycles']=sum([pb.craneCycles for pb in G.terrain.piles]) self.s.stats['harvBiomass']=sum([b.weight for b in self.s.m.trees]) #total mass of the trees that were harvested, before tc? self.s.stats['harvStemMass']=sum([t.logWeight for t in self.s.m.trees])#total weight of the stems of the trees that are harvested. Not what's in the bundles or piles, but rather what was in terrain before chop. self.s.stats['harvStemVol']=sum([t.vol for t in self.s.m.trees])# total volume of the stems of the trees that are harvested. Not what's in the bundles or piles self.s.stats['noBundlesOrPiles']=len(G.terrain.piles) self.s.stats['minBunPileMass']=min([b.biomass for b in G.terrain.piles]) self.s.stats['maxBunPileMass']=max([b.biomass for b in G.terrain.piles]) self.s.stats['totBunPileMass']=sum([b.biomass for b in G.terrain.piles]) self.s.stats['minBunPileVol']=min([self.getVol(b) for b in G.terrain.piles]) self.s.stats['maxBunPileVol']=max([self.getVol(b) for b in G.terrain.piles]) self.s.stats['totBunPileVol']=sum([self.getVol(b) for b in G.terrain.piles]) self.s.stats['noMainStops']=len(self.s.m.positions)-1 #Here I assume what is meant is number of stops on the mainroad for harvesting. Not number of places with piles close by. Is this a good assumption. Gives seven all the time.. should be len(self.s.m.positions)-1 maybe? self.s.stats['totTimeConsumed']=self.s.now() if bundler==True: #self.s.stats['bundlingTime']=self.o.tstep*sum([c[1] for c in self.o.bundlerActiveMoni]) can be used instead, but does not give exact result as expected.see below self.s.stats['bundlingTime']=self.s.stats['noBundlesOrPiles']*self.s.m.bundler.timeBundle else: self.s.stats['bundlingTime']=0 self.s.stats['work time']#operator active time print self.s.stats['noBundlesOrPiles'], 'piles or bundles' print self.s.stats['noCraneCycles'], 'was the number of crane cycles' print self.s.stats['bundlingTime'], self.s.stats['totTimeConsumed'] print '---------------------------------------------------------------------' """This part here writes the data to the excel file""" e.modify(self.Paramrow,0,self.s.stats['machineConfig']) e.modify(self.Paramrow,1,self.s.stats['treeFile']) e.modify(self.Paramrow,2,self.s.stats['simNumber']) e.modify(self.Paramrow,3,self.s.stats['noHarvTrees']) e.modify(self.Paramrow,4,self.s.stats['noCraneCycles']) e.modify(self.Paramrow,5,self.s.stats['harvBiomass']) e.modify(self.Paramrow,6,self.s.stats['harvStemMass']) e.modify(self.Paramrow,7,self.s.stats['harvStemVol']) e.modify(self.Paramrow,8,self.s.stats['noBundlesOrPiles']) e.modify(self.Paramrow,9,self.s.stats['minBunPileMass'])# e.modify(self.Paramrow,10,self.s.stats['maxBunPileMass']) e.modify(self.Paramrow,11,self.s.stats['totBunPileMass']) e.modify(self.Paramrow,12,self.s.stats['minBunPileVol']) e.modify(self.Paramrow,13,self.s.stats['maxBunPileVol']) e.modify(self.Paramrow,14,self.s.stats['totBunPileVol']) e.modify(self.Paramrow,15,self.s.stats['noMainStops']) e.modify(self.Paramrow,16,self.s.stats['totTimeConsumed']) e.modify(self.Paramrow,17,self.s.stats['bundlingTime']) e.modify(self.Paramrow,18,self.s.stats['work time']) e.modify(self.Paramrow,19,self.s.stats['oneCraneWorkTime']) e.modify(self.Paramrow,20,self.s.stats['twoCranesWorkTime']) e.modify(self.Paramrow,21,self.s.stats['oneCraneWaitDriverTime']) e.modify(self.Paramrow,22,self.s.stats['twoCranesWaitDriverTime']) e.modify(self.Paramrow,23,self.s.stats['oneCraneWaitBundlerTime']) e.modify(self.Paramrow,24,self.s.stats['twoCranesWaitBundlerTime']) e.modify(self.Paramrow,25,self.s.stats['noCraneWaitings']) e.modify(self.Paramrow,26,self.s.stats['noCraneWaitingsTwo']) self.Paramrow+=1 e.save()#To be sure to save after each simulation, if something should go wrong print 'Congratulations, all your simulations has been run and the data has successfully been stored in the excel-file named ThinningWawoBundler_date and time.xls, to be found in tota/outputFiles/NewThinning_2012.'
def __init__(self,it=1, head='BC', nCranes=1, bundler=False, twigCrack=False, simNumber=1, rowNumber=1, treeFile=210):
if bundler==1: bundler=False
elif bundler==2: bundler=True
else: raise Exception('bundler should be 1 or 2 from bash')
if twigCrack==1: twigCrack=False
elif twigCrack==2: twigCrack=True
else: raise Exception('twigCrack should be 1 or 2 from bash')
self.G=globalVar()
#self.G.plotDelay=500 #for debug only
self.G.areaPoly=[(0,0), (25,0), (25,40), (0,40)]
self.G.terrain=Terrain(G=self.G)
self.folder='outputFiles/NewThinning_2012'
today=datetime.date.today()
self.filename=self.folder+'/'+'ThinningWawoBundler.xls'
e=ExcelOutput(template='outputFiles/NewThinning_2012/ThinningWawoBundler.xls', out=self.filename)
self.Paramrow=rowNumber
self.G.terrain.treeFile=treeFile
self.G.terrain.readTrees()
G=copy.deepcopy(self.G)
self.s=ThinningSim(G=G, vis=False, anim=False, head=head, nCranes=nCranes, bundler=bundler, twigCrack=twigCrack, observer=True)
self.s.stats['machineConfig']=self.getConfig(head,nCranes,twigCrack,bundler)
self.s.stats['simNumber']=simNumber
self.s.stats['treeFile']=treeFile
self.s.stats['noHarvTrees']=sum([len(b.trees) for b in G.terrain.piles])
self.s.stats['noCraneCycles']=sum([pb.craneCycles for pb in G.terrain.piles])
self.s.stats['harvBiomass']=sum([b.weight for b in self.s.m.trees]) #total mass of the trees that were harvested, before tc?
self.s.stats['harvStemMass']=sum([t.logWeight for t in self.s.m.trees])# total weight of the stems of the trees that are harvested. Not what's in the bundles or piles, but rather what was in terrain before chop.
self.s.stats['harvStemVol']=sum([t.vol for t in self.s.m.trees])# total volume of the stems of the trees that are harvested. Not what's in the bundles or piles
self.s.stats['noBundlesOrPiles']=len(G.terrain.piles)
self.s.stats['minBunPileMass']=min([b.biomass for b in G.terrain.piles])
self.s.stats['maxBunPileMass']=max([b.biomass for b in G.terrain.piles])
self.s.stats['totBunPileMass']=sum([b.biomass for b in G.terrain.piles])
self.s.stats['minBunPileVol']=min([self.getVol(b) for b in G.terrain.piles])
self.s.stats['maxBunPileVol']=max([self.getVol(b) for b in G.terrain.piles])
self.s.stats['totBunPileVol']=sum([self.getVol(b) for b in G.terrain.piles])
self.s.stats['noMainStops']=len(self.s.m.positions)-1 #Here I assume what is meant is number of stops on the mainroad for harvesting.
self.s.stats['totTimeConsumed']=self.s.now()
if bundler==True:
#self.s.stats['bundlingTime']=self.o.tstep*sum([c[1] for c in self.o.bundlerActiveMoni]) could be used instead, but does not give exact result as expected. Deviates by some 5%
self.s.stats['bundlingTime']=self.s.stats['noBundlesOrPiles']*self.s.m.bundler.timeBundle
else: self.s.stats['bundlingTime']=0
self.s.stats['work time']#operator active time
print self.s.stats['noBundlesOrPiles'], 'piles or bundles'
print self.s.stats['noCraneCycles'], 'was the number of crane cycles'
print self.s.stats['bundlingTime'], self.s.stats['totTimeConsumed']
print '---------------------------------------------------------------------'
"""This part here writes the data to the excel file"""
e.modify(self.Paramrow,0,self.s.stats['machineConfig'])
e.modify(self.Paramrow,1,self.s.stats['treeFile'])
e.modify(self.Paramrow,2,self.s.stats['simNumber'])
e.modify(self.Paramrow,3,self.s.stats['noHarvTrees'])
e.modify(self.Paramrow,4,self.s.stats['noCraneCycles'])
e.modify(self.Paramrow,5,self.s.stats['harvBiomass'])
e.modify(self.Paramrow,6,self.s.stats['harvStemMass'])
e.modify(self.Paramrow,7,self.s.stats['harvStemVol'])
e.modify(self.Paramrow,8,self.s.stats['noBundlesOrPiles'])
e.modify(self.Paramrow,9,self.s.stats['minBunPileMass'])#
e.modify(self.Paramrow,10,self.s.stats['maxBunPileMass'])
e.modify(self.Paramrow,11,self.s.stats['totBunPileMass'])
e.modify(self.Paramrow,12,self.s.stats['minBunPileVol'])
e.modify(self.Paramrow,13,self.s.stats['maxBunPileVol'])
e.modify(self.Paramrow,14,self.s.stats['totBunPileVol'])
e.modify(self.Paramrow,15,self.s.stats['noMainStops'])
e.modify(self.Paramrow,16,self.s.stats['totTimeConsumed'])
e.modify(self.Paramrow,17,self.s.stats['bundlingTime'])
e.modify(self.Paramrow,18,self.s.stats['work time'])
e.modify(self.Paramrow,19,self.s.stats['oneCraneWorkTime'])
e.modify(self.Paramrow,20,self.s.stats['twoCranesWorkTime'])
e.modify(self.Paramrow,21,self.s.stats['oneCraneWaitDriverTime'])
e.modify(self.Paramrow,22,self.s.stats['twoCranesWaitDriverTime'])
e.modify(self.Paramrow,23,self.s.stats['oneCraneWaitBundlerTime'])
e.modify(self.Paramrow,24,self.s.stats['twoCranesWaitBundlerTime'])
e.modify(self.Paramrow,25,self.s.stats['noCraneWaitings'])
e.modify(self.Paramrow,26,self.s.stats['noCraneWaitingsTwo'])
e.save()#To be sure to save after each simulation, if something should go wrong
def _storeData(self, s=None):
"""
Stores the data. The strategy is to store it all, it has to be processed later.
The purpose of this function is to standardize the column 'design', so the numerous data files are easily
accessed and generic scripts can be written for them.
s has to be given if this is the first time the method is called and an output should be created.
"""
e=self.excelOut
if not e: #start a new instance.
if not s:
raise Exception('simulation instance has to be provided at first call for _storeData')
e=ExcelOutput(template='sim/planting/template_art3.xls', out=self.filename)
self.excelOut=e
e.changeSheet(0)
#set all the parameters
paramRow=3
type=copy.copy(s.m.type) #we need to shift 2aXh because of new directions...
e.modify(paramRow, 0, s.m.type) #machine type
e.modify(paramRow, 1, s.G.terrain.ttype) #terrain type
e.modify(paramRow, 2, s.G.terrain.stumpFile) #terrain type
e.modify(paramRow, 3, s.G.terrain.stumpsPerH) #stumps/ha
e.modify(paramRow, 4, s.G.terrain.groundModel)
e.modify(paramRow, 5, s.G.terrain.blockQuota*100) #percent
e.modify(paramRow, 6, s.G.terrain.boulderFreq)
e.modify(paramRow, 7, s.G.terrain.meanBoulderV*1000)
e.modify(paramRow, 8, s.G.terrain.humusType) #humustype
e.modify(paramRow, 9, s.stats['noSurfBoulders'])#set in plantmsim
e.modify(paramRow, 10, round(s.stats['meanSurfBoulderDiam']*10,1))#changed to dm, set in plantmsim
e.modify(paramRow, 12, s.m.craneMaxL)
e.modify(paramRow, 13, s.m.plantMinDist*100)
e.modify(paramRow, 14, s.m.pDevs[0].plantHeads[0].width*100)
e.modify(paramRow, 15, s.G.simParam['TSR'])
if s.G.simParam['inverting']==True and s.G.simParam['ExcavatorInverting']==True:
noManRemInv = 5
cycleTimeInv = s.G.simParam['tInvExcavator']
elif s.G.simParam['inverting']==True and s.G.simParam['KOInverting']==True:
noManRemInv = 11
cycleTimeInv = s.G.simParam['tCWhenInvKO']
elif s.G.simParam['inverting']==False:
noManRemInv = 5
cycleTimeInv = 0
e.modify(paramRow, 16, cycleTimeInv)
e.modify(paramRow, 17, noManRemInv)
if s.G.simParam['rectangular']==False: scoopshape='SemiCyl'
elif s.G.simParam['rectangular']==True: scoopshape='Rect'
e.modify(paramRow, 18, scoopshape)
e.modify(paramRow, 19, s.G.simParam['multiplierFindMuSite'])
e.modify(paramRow, 20, s.G.simParam['moundRadius'])
if s.G.simParam['inverting']==True:
e.modify(paramRow, 21, s.G.simParam['SBMinvert'])
elif s.G.simParam['inverting']==False:
e.modify(paramRow, 21, s.G.simParam['SBMmound'])
e.changeSheet(1) #the time-data sheet
tdrow=2+self.sims
for col,val in enumerate(self.timeData):
e.modify(tdrow,col,round(val,2)) #writes the time data. Request from BTE to round off
e.save()
return e
