diff --git a/Makefile b/Makefile
index 8d26b02..aa7bfd2 100644
--- a/Makefile
+++ b/Makefile
@@ -15,7 +15,7 @@ URL_BIN_REL="$(URL_BIN)/releases"
 URL_R_REL="$(URL_R)/releases"
 VERSION1=1
 VERSION2=0
-VERSION3=6
+VERSION3=7
 LICENSE=AGPL-3
 LICENSE_FULL="GNU Affero General Public License, Version 3"
 LICENSE_URL="https://www.gnu.org/licenses/agpl-3.0"
diff --git a/UPDATES b/UPDATES
index d6dd91f..1b08471 100644
--- a/UPDATES
+++ b/UPDATES
@@ -1,3 +1,5 @@
+1.0.7:
+	Added compatibility warning for interface changes in future versions.
 1.0.6:
 	Corrected a bug that may produce biased output in pij_gassist, pij_gassist_trad, and pijs_gassist when nodiag is set or when memlimit is small so computation is split into chunks.
 	Now setting histogram bounds based on the maximum of all LLRs (as opposed to the maximum of the chunk when memlimit is small) in pij_gassist, pij_gassist_trad, and pijs_gassist. This ensures the output is independent of memlimit (related to question from sritchie73@github).
diff --git a/doc.pdf b/doc.pdf
index 22c7dfa..baa4a9b 100644
Binary files a/doc.pdf and b/doc.pdf differ
diff --git a/pij/cassist/cassist.c b/pij/cassist/cassist.c
index c9026b6..047430a 100644
--- a/pij/cassist/cassist.c
+++ b/pij/cassist/cassist.c
@@ -76,6 +76,7 @@ int pijs_cassist_pv(const MATRIXF* g,const MATRIXF* t,const MATRIXF* t2,VECTORF*
 
 	//Step 1: Supernormalization
 	LOG(9,"Supernormalizing...")
+	LOG(6,"Compatibility: From version 2.x.x, Continuous data will not go through supernormalization by default. A separate function interface will be provided for supernormalization. Simple normalization to zero mean and unit variance will still proceed as before.")
 	MATRIXFF(memcpy)(gnew,g);
 	ret=supernormalizea_byrow(gnew);
 	MATRIXFF(memcpy)(tnew,t);
@@ -150,6 +151,7 @@ int pijs_cassist(const MATRIXF* g,const MATRIXF* t,const MATRIXF* t2,VECTORF* p1
 
 	//Step 1: Supernormalization
 	LOG(9,"Supernormalizing...")
+	LOG(6,"Compatibility: From version 2.x.x, Continuous data will not go through supernormalization by default. A separate function interface will be provided for supernormalization. Simple normalization to zero mean and unit variance will still proceed as before.")
 	MATRIXFF(memcpy)(gnew,g);
 	ret=supernormalizea_byrow(gnew);
 	MATRIXFF(memcpy)(tnew,t);
diff --git a/pij/gassist/gassist.c b/pij/gassist/gassist.c
index 081f531..4ef9d19 100644
--- a/pij/gassist/gassist.c
+++ b/pij/gassist/gassist.c
@@ -81,6 +81,7 @@ int pijs_gassist_pv(const MATRIXG* g,const MATRIXF* t,const MATRIXF* t2,VECTORF*
 
 	//Step 1: Supernormalization
 	LOG(9,"Supernormalizing...")
+	LOG(6,"Compatibility: From version 2.x.x, Continuous data will not go through supernormalization by default. A separate function interface will be provided for supernormalization. Simple normalization to zero mean and unit variance will still proceed as before.")
 	MATRIXFF(memcpy)(tnew,t);
 	ret=supernormalizea_byrow(tnew);
 	MATRIXFF(memcpy)(tnew2,t2);
@@ -169,6 +170,7 @@ int pijs_gassist(const MATRIXG* g,const MATRIXF* t,const MATRIXF* t2,VECTORF* p1
 
 	//Step 1: Supernormalization
 	LOG(9,"Supernormalizing...")
+	LOG(6,"Compatibility: From version 2.x.x, Continuous data will not go through supernormalization by default. A separate function interface will be provided for supernormalization. Simple normalization to zero mean and unit variance will still proceed as before.")
 	MATRIXFF(memcpy)(tnew,t);
 	ret=supernormalizea_byrow(tnew);
 	MATRIXFF(memcpy)(tnew2,t2);
diff --git a/pij/rank.c b/pij/rank.c
index 3892272..487fd84 100644
--- a/pij/rank.c
+++ b/pij/rank.c
@@ -125,6 +125,7 @@ int pij_rank_pv(const MATRIXF* t,const MATRIXF* t2,MATRIXF* p,size_t memlimit)
 
 	//Step 1: Supernormalization
 	LOG(9,"Supernormalizing...")
+	LOG(6,"Compatibility: From version 2.x.x, Continuous data will not go through supernormalization by default. A separate function interface will be provided for supernormalization. Simple normalization to zero mean and unit variance will still proceed as before.")
 	MATRIXFF(memcpy)(tnew,t);
 	ret=supernormalizea_byrow(tnew);
 	MATRIXFF(memcpy)(tnew2,t2);
@@ -207,6 +208,7 @@ int pij_rank(const MATRIXF* t,const MATRIXF* t2,MATRIXF* p,char nodiag,size_t me
 
 	//Step 1: Supernormalization
 	LOG(9,"Supernormalizing...")
+	LOG(6,"Compatibility: From version 2.x.x, Continuous data will not go through supernormalization by default. A separate function interface will be provided for supernormalization. Simple normalization to zero mean and unit variance will still proceed as before.")
 	MATRIXFF(memcpy)(tnew,t);
 	ret=supernormalizea_byrow(tnew);
 	MATRIXFF(memcpy)(tnew2,t2);