/[lmdze]/trunk/Sources/dyn3d/bernoui.f
ViewVC logotype

Contents of /trunk/Sources/dyn3d/bernoui.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 137 - (show annotations)
Wed May 6 15:51:03 2015 UTC (9 years ago) by guez
File size: 742 byte(s)
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