def big2littlePTS(self, fileName, foleName):
# ~~ Openning files ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
fle = open(fileName, "rb")
fole = open(foleName, "wb")
print " +> writing the volumes-file: ", foleName
# ~~ Read/Write dimensions ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
n1 = self.NPOIN3 + 1
minvol = self.minvol * np.ones(self.NPOIN3, dtype=np.float32)
# ~~ Read volumes ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
pbar = ProgressBar(maxval=self.HYDROIT).start()
for i in range(self.HYDROIT):
l, it = unpack(">ii", fle.read(4 + 4))
VOLUME = np.asarray(
unpack(">" + str(self.NPOIN3) + "f",
fle.read(4 * self.NPOIN3)))
VOLUME = np.maximum(VOLUME, minvol)
fle.seek(4, 1)
if it >= self.tfrom and it <= self.tstop:
pbar.write(" ~> read iteration: " + str(it), i)
fole.write(pack("<ii", 4 * n1, it - self.HYDRO00))
fole.write(pack("<" + str(self.NPOIN3) + "f", *(VOLUME)))
fole.write(pack("<i", 4 * n1))
else:
pbar.write(" ~> ignore iteration: " + str(it), i)
pbar.update(i)
pbar.finish()
fle.close()
fole.close()
def big2littleARE(self, fileName, foleName):
# ~~ Openning files ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
fle = open(fileName, "rb")
fole = open(foleName, "wb")
print " +> writing the areas-file: ", foleName
# ~~ Read/Write dimensions ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
NSEG2 = (3 * self.geo.NELEM3 + self.conlim.NPTFR) / 2
MBND2 = np.count_nonzero(self.conlim.BOR["lih"] != 2)
n3 = NSEG2 + MBND2
# n4 = 2*( self.geo.NPLAN-1 )*self.geo.NPOIN3
# ~~ Read volumes ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
pbar = ProgressBar(maxval=self.HYDROIT).start()
for i in range(self.HYDROIT):
l, it = unpack(">ii", fle.read(4 + 4))
AREAS = np.asarray(unpack(">" + str(n3) + "f", fle.read(4 * n3)))
fle.seek(4, 1)
if it >= self.tfrom and it <= self.tstop:
pbar.write(" ~> read iteration: " + str(it), i)
fole.write(pack("<ii", 4 * n3, it - self.HYDRO00))
fole.write(pack("<" + str(n3) + "f", *(AREAS)))
fole.write(pack("<i", 4 * n3))
else:
pbar.write(" ~> ignore iteration: " + str(it), i)
pbar.update(i)
pbar.finish()
fle.close()
fole.close()
def big2littleVFX(self, fileName, foleName):
# ~~ Openning files ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
fle = open(fileName, 'rb')
fole = open(foleName, 'wb')
print ' +> writing the flows-file: ', foleName
# ~~ Read/Write dimensions ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
NSEG2 = (3 * self.geo.NELEM3 + self.conlim.NPTFR) / 2
MBND2 = np.count_nonzero(self.conlim.BOR['lih'] != 2)
n3 = (NSEG2 + MBND2)
#n4 = 2*( self.geo.NPLAN-1 )*self.geo.NPOIN3
# ~~ Read volumes ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
pbar = ProgressBar(maxval=self.HYDROIT).start()
for i in range(self.HYDROIT):
l, it = unpack('>ii', fle.read(4 + 4))
VFLUXES = np.asarray(unpack('>' + str(n3) + 'f', fle.read(4 * n3)))
fle.seek(4, 1)
if it >= self.tfrom and it <= self.tstop:
pbar.write(' ~> read iteration: ' + str(it), i)
fole.write(pack('<ii', 4 * n3, it - self.HYDRO00))
fole.write(pack('<' + str(n3) + 'f', *(VFLUXES)))
fole.write(pack('<i', 4 * n3))
else:
pbar.write(' ~> ignore iteration: ' + str(it), i)
pbar.update(i)
pbar.finish()
fle.close()
fole.close()
