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

Annotation of /trunk/Sources/dyn3d/laplacien_rot.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
File size: 1126 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 3
2 guez 81 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rot.F,v 1.1.1.1
3     ! 2004/05/19 12:53:07 lmdzadmin Exp $
4 guez 3
5 guez 81 SUBROUTINE laplacien_rot(klevel, rotin, rotout, ghx, ghy)
6    
7     ! P. Le Van
8    
9     ! ************************************************************
10     ! ... calcul de ( rotat x nxgrad ) du rotationnel rotin .
11     ! ************************************************************
12    
13     ! klevel et rotin sont des arguments d'entree pour le s-prog
14     ! rotout est un argument de sortie pour le s-prog
15    
16     USE dimens_m
17     USE paramet_m
18     USE comgeom
19 guez 137 USE filtreg_v_m, ONLY: filtreg_v
20 guez 107 use rotatf_m, only: rotatf
21    
22 guez 81 IMPLICIT NONE
23    
24    
25    
26     ! .......... variables en arguments .............
27    
28     INTEGER, INTENT (IN) :: klevel
29 guez 106 REAL rotin(iim + 1, jjm, klevel), rotout(ip1jm, klevel)
30 guez 81
31     ! .......... variables locales ................
32    
33     REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
34     ! ........................................................
35    
36    
37 guez 137 CALL filtreg_v(rotin, intensive = .true.)
38 guez 81
39     CALL nxgrad(klevel, rotin, ghx, ghy)
40     CALL rotatf(klevel, ghx, ghy, rotout)
41    
42     RETURN
43     END SUBROUTINE laplacien_rot

  ViewVC Help
Powered by ViewVC 1.1.21