/[lmdze]/trunk/filtrez/inifgn.f
ViewVC logotype

Annotation of /trunk/filtrez/inifgn.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 167 - (hide annotations)
Mon Aug 24 16:30:33 2015 UTC (8 years, 8 months ago) by guez
Original Path: trunk/Sources/filtrez/inifgn.f
File size: 2162 byte(s) 
Added program test_inifilr.

Encapsulated ppm3d into a module and added implicit none. Removed
unused argument dum.

Encountered a problem in procedure invert_zoom_x. With grossismx=2.9,
DZOOMX=0.3, taux=5, for xuv = -0.25, for i = 1, rtsafe fails because
fval is about 1e-16 instead of 0 at xval = pi. So distinguished the
cases abs_y = 0 or pi. Needed then to add argument beta to
invert_zoom_x.

Moved the output of eignvalues of differentiation matrix from inifilr
to inifgn, where they are computed.

Simpler definition of j1 in inifilr.
1 guez 113 module inifgn_m
2 guez 3
3 guez 140 use dimens_m, only: iim
4    
5 guez 113 IMPLICIT NONE
6 guez 3
7 guez 140 private iim
8    
9 guez 154 real sddu(iim), sddv(iim)
10     ! sdd[uv] = sqrt(2 pi / iim * (derivative of the longitudinal zoom
11     ! function)(rlon[uv]))
12    
13 guez 140 real unsddu(iim), unsddv(iim)
14    
15 guez 113 contains
16 guez 3
17 guez 154 SUBROUTINE inifgn(eignval_v, eignfnu, eignfnv)
18 guez 3
19 guez 113 ! From LMDZ4/libf/filtrez/inifgn.F, v 1.1.1.1 2004/05/19 12:53:09
20 guez 3
21 guez 154 ! Authors: H. Upadyaya, O. Sharma
22 guez 3
23 guez 154 ! Computes the eigenvalues and eigenvectors of the discrete analog
24     ! of the second derivative with respect to longitude.
25    
26 guez 143 use acc_m, only: acc
27 guez 113 USE dimens_m, ONLY: iim
28 guez 139 USE dynetat0_m, ONLY: xprimu, xprimv
29 guez 167 use jumble, only: new_unit
30 guez 143 use numer_rec_95, only: jacobi, eigsrt
31 guez 3
32 guez 154 real, intent(out):: eignval_v(:) ! (iim)
33     ! eigenvalues sorted in descending order
34 guez 81
35 guez 154 real, intent(out):: eignfnu(:, :), eignfnv(:, :) ! (iim, iim) eigenvectors
36    
37 guez 113 ! Local:
38 guez 154
39 guez 161 REAL delta(iim, iim) ! second derivative, symmetric, elements are angle^{-2}
40 guez 154
41     REAL deriv_u(iim, iim), deriv_v(iim, iim)
42     ! first derivative at u and v longitudes, elements are angle^{-1}
43    
44     REAL eignval_u(iim)
45 guez 167 INTEGER i, unit
46 guez 81
47 guez 113 !----------------------------------------------------------------
48 guez 81
49 guez 151 print *, "Call sequence information: inifgn"
50    
51 guez 140 sddv = sqrt(xprimv(:iim))
52     sddu = sqrt(xprimu(:iim))
53     unsddu = 1. / sddu
54     unsddv = 1. / sddv
55 guez 81
56 guez 154 deriv_u = 0.
57     deriv_u(iim, 1) = unsddu(iim) * unsddv(1)
58     forall (i = 1:iim) deriv_u(i, i) = - unsddu(i) * unsddv(i)
59     forall (i = 1:iim - 1) deriv_u(i, i + 1) = unsddu(i) * unsddv(i + 1)
60 guez 81
61 guez 154 deriv_v = - transpose(deriv_u)
62 guez 140
63 guez 161 delta = matmul(deriv_v, deriv_u) ! second derivative at v longitudes
64     CALL jacobi(delta, eignval_v, eignfnv)
65 guez 143 CALL acc(eignfnv)
66 guez 154 CALL eigsrt(eignval_v, eignfnv)
67 guez 81
68 guez 161 delta = matmul(deriv_u, deriv_v) ! second derivative at u longitudes
69     CALL jacobi(delta, eignval_u, eignfnu)
70 guez 143 CALL acc(eignfnu)
71 guez 154 CALL eigsrt(eignval_u, eignfnu)
72 guez 81
73 guez 167 call new_unit(unit)
74     open(unit, file = "inifgn_out.txt", status = "replace", action = "write")
75     write(unit, fmt = *) '"eignval_v"', eignval_v
76     close(unit)
77    
78 guez 113 END SUBROUTINE inifgn
79 guez 81
80 guez 113 end module inifgn_m
  ViewVC Help
Powered by ViewVC 1.1.21