trunk/filtrez/inifgn.f

module inifgn_m

  IMPLICIT NONE

contains

  SUBROUTINE inifgn(dv)

    ! 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 dv(iim)

    ! Local:
    REAL vec(iim, iim), vec1(iim, 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.)

    DO i = 1, iim
       sddv(i) = sqrt(xprimv(i))
       sddu(i) = sqrt(xprimu(i))
       unsddu(i) = 1./sddu(i)
       unsddv(i) = 1./sddv(i)
    END DO

    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)

  END SUBROUTINE inifgn

end module inifgn_m
1	guez	113	module inifgn_m
2	guez	3
3	guez	113	IMPLICIT NONE
4	guez	3
5	guez	113	contains
6	guez	3
7	guez	113	SUBROUTINE inifgn(dv)
8	guez	3
9	guez	113	! From LMDZ4/libf/filtrez/inifgn.F, v 1.1.1.1 2004/05/19 12:53:09
10	guez	3
11	guez	113	! H.Upadyaya, O.Sharma
12	guez	3
13	guez	113	USE dimens_m, ONLY: iim
14			USE comgeom, ONLY: xprimu, xprimv
15			USE coefils, ONLY: eignfnu, eignfnv, sddu, sddv, unsddu, unsddv
16	guez	3
17	guez	113	real dv(iim)
18	guez	81
19	guez	113	! Local:
20			REAL vec(iim, iim), vec1(iim, iim)
21			REAL du(iim)
22			real d(iim)
23			REAL pi
24			INTEGER i, j, k, imm1, nrot
25	guez	81
26	guez	113	EXTERNAL acc, jacobi
27	guez	81
28	guez	113	!----------------------------------------------------------------
29	guez	81
30	guez	113	imm1 = iim - 1
31			pi = 2.*asin(1.)
32	guez	81
33	guez	113	DO i = 1, iim
34	guez	121	sddv(i) = sqrt(xprimv(i))
35			sddu(i) = sqrt(xprimu(i))
36	guez	113	unsddu(i) = 1./sddu(i)
37			unsddv(i) = 1./sddv(i)
38	guez	81	END DO
39
40	guez	113	DO j = 1, iim
41			DO i = 1, iim
42			vec(i, j) = 0.
43			vec1(i, j) = 0.
44			eignfnv(i, j) = 0.
45			eignfnu(i, j) = 0.
46			END DO
47			END DO
48	guez	81
49	guez	113	eignfnv(1, 1) = -1.
50			eignfnv(iim, 1) = 1.
51			DO i = 1, imm1
52			eignfnv(i+1, i+1) = -1.
53			eignfnv(i, i+1) = 1.
54	guez	81	END DO
55	guez	113	DO j = 1, iim
56			DO i = 1, iim
57			eignfnv(i, j) = eignfnv(i, j)/(sddu(i)*sddv(j))
58			END DO
59	guez	81	END DO
60	guez	113	DO j = 1, iim
61			DO i = 1, iim
62			eignfnu(i, j) = -eignfnv(j, i)
63			END DO
64			END DO
65	guez	81
66	guez	113	DO j = 1, iim
67			DO i = 1, iim
68			vec(i, j) = 0.0
69			vec1(i, j) = 0.0
70			DO k = 1, iim
71			vec(i, j) = vec(i, j) + eignfnu(i, k)*eignfnv(k, j)
72			vec1(i, j) = vec1(i, j) + eignfnv(i, k)*eignfnu(k, j)
73			END DO
74			END DO
75	guez	81	END DO
76
77	guez	113	CALL jacobi(vec, iim, iim, dv, eignfnv, nrot)
78			CALL acc(eignfnv, d, iim)
79			CALL eigen_sort(dv, eignfnv, iim, iim)
80	guez	81
81	guez	113	CALL jacobi(vec1, iim, iim, du, eignfnu, nrot)
82			CALL acc(eignfnu, d, iim)
83			CALL eigen_sort(du, eignfnu, iim, iim)
84	guez	81
85	guez	113	END SUBROUTINE inifgn
86	guez	81
87	guez	113	end module inifgn_m