/[lmdze]/trunk/dyn3d/laplacien.f
ViewVC logotype

Annotation of /trunk/dyn3d/laplacien.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 137 - (hide annotations)
Wed May 6 15:51:03 2015 UTC (9 years ago) by guez
Original Path: trunk/Sources/dyn3d/laplacien.f
File size: 832 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 guez 61 module laplacien_m
2 guez 57
3 guez 61 IMPLICIT NONE
4 guez 3
5 guez 61 contains
6 guez 3
7 guez 61 SUBROUTINE laplacien(klevel, teta)
8 guez 3
9 guez 61 ! 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 guez 3
13 guez 106 use dimens_m, only: iim, jjm
14 guez 65 use divergf_m, only: divergf
15 guez 137 use filtreg_scal_m, only: filtreg_scal
16 guez 61 use grad_m, only: grad
17     USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1
18 guez 3
19 guez 61 INTEGER, intent(in):: klevel
20 guez 106 REAL, intent(inout):: teta(iim + 1, jjm + 1, klevel)
21 guez 61
22     ! Variables locales:
23 guez 103 REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
24 guez 61
25     !-----------------------------------------------------------------
26    
27 guez 137 CALL filtreg_scal(teta, direct = .true., intensive = .true.)
28 guez 61 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