source: LMDZ6/trunk/libf/dyn3dmem/rotatf_loc.f90 @ 5281

Last change on this file since 5281 was 5281, checked in by abarral, 4 days ago

Turn comgeom.h comgeom2.h into 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.8 KB
Line 
1SUBROUTINE rotatf_loc (klevel, x, y, rot )
2  !
3  ! Auteur : P.Le Van
4  !**************************************************************
5  !.  calcule le rotationnel
6  ! a tous les niveaux d'1 vecteur de comp. x et y ..
7  !   x  et  y etant des composantes  covariantes  ...
8  !********************************************************************
9  !   klevel, x  et y   sont des arguments d'entree pour le s-prog
10  !        rot          est  un argument  de sortie pour le s-prog
11  !
12  USE comgeom_mod_h
13  USE parallel_lmdz
14  USE mod_filtreg_p
15  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
16USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
17          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
18IMPLICIT NONE
19  !
20
21
22  !
23  !   .....  variables en arguments  ......
24  !
25  INTEGER :: klevel
26  REAL :: rot( ijb_v:ije_v,klevel )
27  REAL :: x( ijb_u:ije_u,klevel ), y( ijb_v:ije_v,klevel )
28  !
29  !  ...   variables  locales  ...
30  !
31  INTEGER :: l, ij
32  INTEGER :: ijb,ije,jjb,jje
33  !
34  !
35  ijb=ij_begin
36  ije=ij_end
37  if(pole_sud) ije=ij_end-iip1
38
39!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
40  DO l = 1,klevel
41  !
42    DO   ij = ijb, ije - 1
43     rot( ij,l )  =    y( ij+1 , l )  -  y( ij,l )   + &
44           x(ij +iip1, l )  -  x( ij,l )
45    ENDDO
46  !
47  !    .... correction pour rot( iip1,j,l)  ....
48  !    ....   rot(iip1,j,l)= rot(1,j,l) ...
49  !DIR$ IVDEP
50    DO  ij = ijb+iip1-1, ije, iip1
51     rot( ij,l ) = rot( ij -iim,l )
52    ENDDO
53  !
54  END DO
55!$OMP END DO NOWAIT
56    jjb=jj_begin
57    jje=jj_end
58    if (pole_sud) jje=jj_end-1
59    CALL filtreg_p( rot, jjb_v, jje_v,jjb,jje,jjm, &
60          klevel, 2, 2, .FALSE., 1 )
61
62!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
63    DO l = 1, klevel
64      DO ij = ijb, ije
65       rot(ij,l) = rot(ij,l) * unsairez(ij)
66      ENDDO
67    ENDDO
68!$OMP END DO NOWAIT
69  !
70  !
71  RETURN
72END SUBROUTINE rotatf_loc
Note: See TracBrowser for help on using the repository browser.