Parent Directory | Revision Log
Removed unused argument missval in ma_fucoll_r[1-3]1, ma_fufill_r[1-3]1. Split filtreg into two procedures: filtreg_scal and filtreg_v. I did not like the test on the extent of the argument and there was no common code between the two cases: jjm and jjm + 1. Also, it is simpler now to just remove the argument "direct" from filtreg_v instead of allowing it and then stopping the program if it is false. Removed the computation of pkf in reanalyse2nat, was not used. As a consequence of the split of filtreg, had to extract the computation of pkf out of exner_hyb. This is clearer anyway because we want to be able to call exner_hyb with any size in the first two dimensions (as in test_disvert). But at the same time exner_hyb required particular sizes for the computation of pkf. It was awkward. The only computation of pkf is now in leapfrog.
1 | module laplacien_m |
2 | |
3 | IMPLICIT NONE |
4 | |
5 | contains |
6 | |
7 | SUBROUTINE laplacien(klevel, teta) |
8 | |
9 | ! From LMDZ4/libf/dyn3d/laplacien.F, version 1.1.1.1 2004/05/19 12:53:06 |
10 | ! P. Le Van |
11 | ! Calcul de div(grad) de teta. |
12 | |
13 | use dimens_m, only: iim, jjm |
14 | use divergf_m, only: divergf |
15 | use filtreg_scal_m, only: filtreg_scal |
16 | use grad_m, only: grad |
17 | USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1 |
18 | |
19 | INTEGER, intent(in):: klevel |
20 | REAL, intent(inout):: teta(iim + 1, jjm + 1, klevel) |
21 | |
22 | ! Variables locales: |
23 | REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel) |
24 | |
25 | !----------------------------------------------------------------- |
26 | |
27 | CALL filtreg_scal(teta, direct = .true., intensive = .true.) |
28 | CALL grad(klevel, teta, ghx, ghy) |
29 | CALL divergf(klevel, ghx, ghy, teta) |
30 | |
31 | END SUBROUTINE laplacien |
32 | |
33 | end module laplacien_m |
ViewVC Help | |
Powered by ViewVC 1.1.21 |