missing line in wl.f90

mishakb · Feb 17, 2020 · 5808bf4 · 5808bf4
1 parent 6e0f355
commit 5808bf4
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/source/wl.f90 b/source/wl.f90
@@ -3,7 +3,7 @@
     !AL 2018, following exactly the same approximations as in the DES papers
     !(can only use Weyl potential for lensing)
     ! MR 2019 update to use Weyl potential for galaxy-lensing cross
-	!CGQ 2019 add small changes to compute the cls and correlation functions for non-flat cases
+	!CGQ 2019 add small changes to obtain the cls and correlation functions for non-flat cases
 
     module wl
 
@@ -531,6 +531,9 @@ subroutine calc_theory(this,CMB,Theory,corrs, DataParams)
     omm = CMB%omdm+CMB%omb
 
     allocate(chis(this%num_z_p), dchis(this%num_z_p))
+	!>ISiTGR MOD START
+	allocate(term1(this%num_z_p,this%num_z_p), term2(this%num_z_p))
+	!<ISiTGR MOD END
     call this%Calculator%ComovingRadialDistanceArr(this%z_p, chis, this%num_z_p)
     dchis(1) = (chis(2) + chis(1))/2
     dchis(this%num_z_p) = chis(this%num_z_p) - chis(this%num_z_p-1)