--- trunk/filtrez/inifgn.f	2014/09/18 13:36:51	112
+++ trunk/filtrez/inifgn.f	2014/09/18 19:56:46	113
@@ -1,106 +1,93 @@
+module inifgn_m
 
-! $Header: /home/cvsroot/LMDZ4/libf/filtrez/inifgn.F,v 1.1.1.1 2004/05/19
-! 12:53:09 lmdzadmin Exp $
+  IMPLICIT NONE
 
-SUBROUTINE inifgn(dv)
+contains
 
-  ! ...  H.Upadyaya , O.Sharma  ...
+  SUBROUTINE inifgn(dv)
 
-  USE dimens_m
-  USE paramet_m
-  USE comgeom
-  USE serre
-  USE coefils
-  IMPLICIT NONE
+    ! From LMDZ4/libf/filtrez/inifgn.F, v 1.1.1.1 2004/05/19 12:53:09
 
+    ! H.Upadyaya, O.Sharma 
 
+    USE dimens_m, ONLY: iim
+    USE comgeom, ONLY: xprimu, xprimv
+    USE coefils, ONLY: eignfnu, eignfnv, sddu, sddv, unsddu, unsddv
 
-  REAL vec(iim, iim), vec1(iim, iim)
-  REAL dlonu(iim), dlonv(iim)
-  REAL du(iim), dv(iim), d(iim)
-  REAL pi
-  INTEGER i, j, k, imm1, nrot
-
-
-  EXTERNAL ssum, acc, jacobi
-  ! C      EXTERNAL eigen
-  REAL ssum
-
-
-  imm1 = iim - 1
-  pi = 2.*asin(1.)
-
-  DO i = 1, iim
-    dlonu(i) = xprimu(i)
-    dlonv(i) = xprimv(i)
-  END DO
-
-  DO i = 1, iim
-    sddv(i) = sqrt(dlonv(i))
-    sddu(i) = sqrt(dlonu(i))
-    unsddu(i) = 1./sddu(i)
-    unsddv(i) = 1./sddv(i)
-  END DO
+    real dv(iim)
 
-  DO j = 1, iim
-    DO i = 1, iim
-      vec(i, j) = 0.
-      vec1(i, j) = 0.
-      eignfnv(i, j) = 0.
-      eignfnu(i, j) = 0.
-    END DO
-  END DO
+    ! Local:
+    REAL vec(iim, iim), vec1(iim, iim)
+    REAL dlonu(iim), dlonv(iim)
+    REAL du(iim)
+    real d(iim)
+    REAL pi
+    INTEGER i, j, k, imm1, nrot
+
+    EXTERNAL acc, jacobi
 
+    !----------------------------------------------------------------
+
+    imm1 = iim - 1
+    pi = 2.*asin(1.)
 
-  eignfnv(1, 1) = -1.
-  eignfnv(iim, 1) = 1.
-  DO i = 1, imm1
-    eignfnv(i+1, i+1) = -1.
-    eignfnv(i, i+1) = 1.
-  END DO
-  DO j = 1, iim
     DO i = 1, iim
-      eignfnv(i, j) = eignfnv(i, j)/(sddu(i)*sddv(j))
+       dlonu(i) = xprimu(i)
+       dlonv(i) = xprimv(i)
     END DO
-  END DO
-  DO j = 1, iim
+
     DO i = 1, iim
-      eignfnu(i, j) = -eignfnv(j, i)
+       sddv(i) = sqrt(dlonv(i))
+       sddu(i) = sqrt(dlonu(i))
+       unsddu(i) = 1./sddu(i)
+       unsddv(i) = 1./sddv(i)
     END DO
-  END DO
 
-  DO j = 1, iim
-    DO i = 1, iim
-      vec(i, j) = 0.0
-      vec1(i, j) = 0.0
-      DO k = 1, iim
-        vec(i, j) = vec(i, j) + eignfnu(i, k)*eignfnv(k, j)
-        vec1(i, j) = vec1(i, j) + eignfnv(i, k)*eignfnu(k, j)
-      END DO
-    END DO
-  END DO
-
-
-  CALL jacobi(vec, iim, iim, dv, eignfnv, nrot)
-  CALL acc(eignfnv, d, iim)
-  CALL eigen_sort(dv, eignfnv, iim, iim)
-
-  CALL jacobi(vec1, iim, iim, du, eignfnu, nrot)
-  CALL acc(eignfnu, d, iim)
-  CALL eigen_sort(du, eignfnu, iim, iim)
-
-  ! c   ancienne version avec appels IMSL
-
-  ! CALL MXM(eignfnu,iim,eignfnv,iim,vec,iim)
-  ! CALL MXM(eignfnv,iim,eignfnu,iim,vec1,iim)
-  ! CALL EVCSF(iim,vec,iim,dv,eignfnv,iim)
-  ! CALL acc(eignfnv,d,iim)
-  ! CALL eigen(eignfnv,dv)
-
-  ! CALL EVCSF(iim,vec1,iim,du,eignfnu,iim)
-  ! CALL acc(eignfnu,d,iim)
-  ! CALL eigen(eignfnu,du)
+    DO j = 1, iim
+       DO i = 1, iim
+          vec(i, j) = 0.
+          vec1(i, j) = 0.
+          eignfnv(i, j) = 0.
+          eignfnu(i, j) = 0.
+       END DO
+    END DO
+
+    eignfnv(1, 1) = -1.
+    eignfnv(iim, 1) = 1.
+    DO i = 1, imm1
+       eignfnv(i+1, i+1) = -1.
+       eignfnv(i, i+1) = 1.
+    END DO
+    DO j = 1, iim
+       DO i = 1, iim
+          eignfnv(i, j) = eignfnv(i, j)/(sddu(i)*sddv(j))
+       END DO
+    END DO
+    DO j = 1, iim
+       DO i = 1, iim
+          eignfnu(i, j) = -eignfnv(j, i)
+       END DO
+    END DO
+
+    DO j = 1, iim
+       DO i = 1, iim
+          vec(i, j) = 0.0
+          vec1(i, j) = 0.0
+          DO k = 1, iim
+             vec(i, j) = vec(i, j) + eignfnu(i, k)*eignfnv(k, j)
+             vec1(i, j) = vec1(i, j) + eignfnv(i, k)*eignfnu(k, j)
+          END DO
+       END DO
+    END DO
+
+    CALL jacobi(vec, iim, iim, dv, eignfnv, nrot)
+    CALL acc(eignfnv, d, iim)
+    CALL eigen_sort(dv, eignfnv, iim, iim)
+
+    CALL jacobi(vec1, iim, iim, du, eignfnu, nrot)
+    CALL acc(eignfnu, d, iim)
+    CALL eigen_sort(du, eignfnu, iim, iim)
 
-  RETURN
-END SUBROUTINE inifgn
+  END SUBROUTINE inifgn
 
+end module inifgn_m