def testUnits(self):
app1=self.app1
time,targetTime,timestepnumber=0,10,0
while True:
#determine critical time step
dt = app1.getCriticalTimeStep()
#update time
time = time+dt
if (time > targetTime):
#make sure we reach targetTime at the end
time = targetTime
timestepnumber = timestepnumber+1
# print ("Step: ", timestepnumber, time, dt)
# create a time step
istep = TimeStep.TimeStep(time, dt, timestepnumber)
#solve problem 1
app1.solveStep(istep)
#request Concentration property from app1
v = app1.getProperty(PropertyID.PID_Velocity, istep)
#Create a PhysicalQuantity object
V = PQ(v.getValue(), v.getUnits())
velocity = V.inBaseUnits()
self.assert_(v.getValue()==time)
#can be converted in km/s?
self.assert_(V.isCompatible('km/s'))
#can be converted in km?
self.assertFalse(V.isCompatible('km'))
# convert in km/h
V.convertToUnit('km/h')
#sys.stderr.write(str(V))
#give only the value
value = float(str(V).split()[0])
self.assertAlmostEqual(value,3.6*v.getValue())
if time==targetTime: break
time = time+dt
if (time > targetTime):
#make sure we reach targetTime at the end
time = targetTime
timestepnumber = timestepnumber+1
mupif.log.debug("Step: %g %g %g"%(timestepnumber,time,dt))
# create a time step
istep = TimeStep.TimeStep(time, dt, timestepnumber)
#solve problem 1
app1.solveStep(istep)
#request Concentration property from app1
v = app1.getProperty(PropertyID.PID_Velocity, istep)
#Create a PhysicalQuantity object
V = PQ(v.getValue(), v.getUnits())
velocity = V.inBaseUnits()
mupif.log.debug(velocity)
#can be converted in km/s?
mupif.log.debug(V.isCompatible('km/s'))
#can be converted in km?
mupif.log.debug(V.isCompatible('km'))
# convert in km/h
V.convertToUnit('km/h')
mupif.log.debug(V)
#give only the value