 
 	Tenv = himan::metutil::VirtualTemperature_(Tenv, Penv * 100);
 	Tparcel = himan::metutil::VirtualTemperature_(Tparcel, Penv * 100);
+
+	prevTenv = himan::metutil::VirtualTemperature_(prevTenv, prevPenv * 100);
+	prevTparcel = himan::metutil::VirtualTemperature_(prevTparcel, prevPenv * 100);
 
 	if (Tparcel >= Tenv && prevTparcel >= Tenv)
 	{

	Tenv = himan::metutil::VirtualTemperature_(Tenv, Penv * 100);
	Tparcel = himan::metutil::VirtualTemperature_(Tparcel, Penv * 100);

	double coldColderLimit = 233.15;
	double coldWarmerLimit = 263.15;


	Tenv = himan::metutil::VirtualTemperature_(Tenv, Penv * 100);
	Tparcel = himan::metutil::VirtualTemperature_(Tparcel, Penv * 100);

	prevTenv = himan::metutil::VirtualTemperature_(prevTenv, prevPenv * 100);
	prevTparcel = himan::metutil::VirtualTemperature_(prevTparcel, prevPenv * 100);

	if (Tparcel > Tenv)
	{
		// Have buoyancy at current height
