--- trunk/Sources/filtrez/inifilr.f	2015/06/23 15:14:20	151
+++ trunk/Sources/filtrez/inifilr.f	2015/07/16 17:39:10	156
@@ -7,12 +7,14 @@
   ! jfilts index of the first line filtered in SH
 
   ! North:
+
   real, allocatable:: matriceun(:, :, :), matrinvn(:, :, :) 
   ! (iim, iim, 2:jfiltnu)
 
   real, allocatable:: matricevn(:, :, :) ! (iim, iim, jfiltnv)
 
   ! South:
+
   real, allocatable:: matriceus(:, :, :), matrinvs(:, :, :)
   ! (iim, iim, jfiltsu:jjm)
 
@@ -25,13 +27,13 @@
     ! From filtrez/inifilr.F, version 1.1.1.1 2004/05/19 12:53:09
     ! H. Upadhyaya, O. Sharma
 
-    ! This routine computes the eigenfunctions of the laplacian on the
+    ! This routine computes the eigenvectors of the laplacian on the
     ! stretched grid, and the filtering coefficients. The modes are
     ! filtered from modfrst to iim.
 
     USE dimens_m, ONLY : iim, jjm
     USE dynetat0_m, ONLY : rlatu, rlatv, xprimu, grossismx
-    use inifgn_m, only: inifgn, eignfnu, eignfnv
+    use inifgn_m, only: inifgn
     use jumble, only: new_unit
     use nr_util, only: pi
 
@@ -44,6 +46,9 @@
     REAL colat0 ! > 0
     REAL eignft(iim, iim), coff
 
+    real eignfnu(iim, iim), eignfnv(iim, iim)
+    ! eigenvectors of the discrete laplacian
+
     ! Filtering coefficients (lamda_max * cos(rlat) / lamda):
     real coefilu(iim, jjm), coefilv(iim, jjm)
     real coefilu2(iim, jjm), coefilv2(iim, jjm)
@@ -56,9 +61,9 @@
 
     print *, "Call sequence information: inifilr"
 
-    CALL inifgn(eignvl)
+    CALL inifgn(eignvl, eignfnu, eignfnv)
 
-    ! compute eigenvalues and eigenfunctions
+    ! compute eigenvalues and eigenvectors
     ! compute the filtering coefficients for scalar lines and
     ! meridional wind v-lines
     ! we filter all those latitude lines where coefil < 1