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 bernoui_m |
2 | |
3 | IMPLICIT NONE |
4 | |
5 | contains |
6 | |
7 | function bernoui(phi, ecin) |
8 | |
9 | ! From LMDZ4/libf/dyn3d/bernoui.F, version 1.1.1.1 2004/05/19 12:53:06 |
10 | |
11 | ! Author: P. Le Van |
12 | |
13 | ! Objet : calcul de la fonction de Bernouilli aux niveaux s. |
14 | |
15 | USE dimens_m, ONLY: iim, jjm, llm |
16 | USE filtreg_scal_m, ONLY: filtreg_scal |
17 | |
18 | REAL, INTENT(IN):: phi(:, :, :), ecin (:, :, :) ! (iim + 1, jjm + 1, llm) |
19 | |
20 | REAL bernoui(iim + 1, jjm + 1, llm) |
21 | ! fonction de Bernouilli = filtre de (géopotentiel + énergie cinétique) |
22 | |
23 | !----------------------------------------------------------------------- |
24 | |
25 | bernoui = phi + ecin |
26 | CALL filtreg_scal(bernoui, direct = .true., intensive = .true.) |
27 | |
28 | END function bernoui |
29 | |
30 | end module bernoui_m |
ViewVC Help | |
Powered by ViewVC 1.1.21 |