35 |
|
|
36 |
! Local: |
! Local: |
37 |
|
|
38 |
REAL a(iim, iim) ! second derivative, symmetric, elements are angle^{-2} |
REAL delta(iim, iim) ! second derivative, symmetric, elements are angle^{-2} |
39 |
|
|
40 |
REAL deriv_u(iim, iim), deriv_v(iim, iim) |
REAL deriv_u(iim, iim), deriv_v(iim, iim) |
41 |
! first derivative at u and v longitudes, elements are angle^{-1} |
! first derivative at u and v longitudes, elements are angle^{-1} |
59 |
|
|
60 |
deriv_v = - transpose(deriv_u) |
deriv_v = - transpose(deriv_u) |
61 |
|
|
62 |
a = matmul(deriv_v, deriv_u) ! second derivative at v longitudes |
delta = matmul(deriv_v, deriv_u) ! second derivative at v longitudes |
63 |
CALL jacobi(a, eignval_v, eignfnv) |
CALL jacobi(delta, eignval_v, eignfnv) |
64 |
CALL acc(eignfnv) |
CALL acc(eignfnv) |
65 |
CALL eigsrt(eignval_v, eignfnv) |
CALL eigsrt(eignval_v, eignfnv) |
66 |
|
|
67 |
a = matmul(deriv_u, deriv_v) ! second derivative at u longitudes |
delta = matmul(deriv_u, deriv_v) ! second derivative at u longitudes |
68 |
CALL jacobi(a, eignval_u, eignfnu) |
CALL jacobi(delta, eignval_u, eignfnu) |
69 |
CALL acc(eignfnu) |
CALL acc(eignfnu) |
70 |
CALL eigsrt(eignval_u, eignfnu) |
CALL eigsrt(eignval_u, eignfnu) |
71 |
|
|