diff --git a/src/include/Commons.h b/src/include/Commons.h
index db9a134de6f7521d9f60c8a8c4b8b82b0ceb252f..4aeb5ef024cfd61b3d6c038d2b02a8b358629d41 100644
--- a/src/include/Commons.h
+++ b/src/include/Commons.h
@@ -214,6 +214,8 @@ public:
 	//! \brief QUESTION: definition?
 	std::complex<double> *vj0;
 	//! \brief QUESTION: definition?
+	std::complex<double> *vj;
+	//! \brief QUESTION: definition?
 	std::complex<double> *vyhj;
 	//! \brief QUESTION: definition?
 	std::complex<double> *vyj0;
diff --git a/src/libnptm/Commons.cpp b/src/libnptm/Commons.cpp
index 004ee59e73c2c06753ea3313f0b3cea2c04ecbd1..342381f8eec4ac6a21944bc3e0fb5b6104f935c4 100644
--- a/src/libnptm/Commons.cpp
+++ b/src/libnptm/Commons.cpp
@@ -99,6 +99,7 @@ C1_AddOns::C1_AddOns(C4 *c4) {
 	nlemt = 2 * c4->nlem;
 	vh = new complex<double>[(nsph * nsph - 1) * c4->litpo];
 	vj0 = new complex<double>[nsph * c4->lmtpo];
+	vj = new complex<double>[1]; // QUESTION: is 1 really enough for a general case?
 	vyhj = new complex<double>[nsph * (nsph - 1) * c4->litpo];
 	vyj0 = new complex<double>[nsph * c4->lmtpo];
 	am0m = new complex<double>*[nlemt];