source: LMDZ6/trunk/libf/dyn3dmem/laplacien_rotgam_loc.f90 @ 5407

Last change on this file since 5407 was 5285, checked in by abarral, 7 weeks ago

As discussed internally, remove generic ONLY: ... for new _mod_h modules

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.3 KB
Line 
1SUBROUTINE laplacien_rotgam_loc ( klevel, rotin, rotout )
2  !
3  ! P. Le Van
4  !
5  !   ************************************************************
6  !   ... calcul de  (rotat x nxgrad)_gam  du rotationnel rotin ..
7  !   ************************************************************
8  ! klevel et teta  sont des arguments  d'entree pour le s-prog
9  !  divgra     est  un argument  de sortie pour le s-prog
10  !
11  USE comgeom_mod_h
12  USE parallel_lmdz
13  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
14USE paramet_mod_h
15IMPLICIT NONE
16  !
17
18
19
20  !
21  !    .............   variables  en  arguments    ...........
22  !
23  INTEGER :: klevel
24  REAL :: rotin( ijb_v:ije_v,klevel ), rotout( ijb_v:ije_v,klevel )
25  !
26  !   ............     variables   locales     ...............
27  !
28  INTEGER :: l, ij
29  REAL :: ghy(ijb_v:ije_v,llm), ghx(ijb_u:ije_u,llm)
30  !   ........................................................
31  !
32  INTEGER :: ijb,ije
33
34  !
35
36  CALL   nxgrad_gam_loc ( klevel, rotin,   ghx ,   ghy  )
37  CALL   rotat_nfil_loc ( klevel, ghx  ,   ghy , rotout )
38  !
39  ijb=ij_begin
40  ije=ij_end
41  if(pole_sud) ije=ij_end-iip1
42!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
43  DO l = 1, klevel
44    DO ij = ijb, ije
45     rotout(ij,l) = rotout(ij,l) * unsairz_gam(ij)
46    ENDDO
47  ENDDO
48!$OMP END DO NOWAIT
49  RETURN
50END SUBROUTINE laplacien_rotgam_loc
Note: See TracBrowser for help on using the repository browser.