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

Last change on this file since 5281 was 5281, checked in by abarral, 3 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.4 KB
Line 
1SUBROUTINE covcont_loc (klevel,ucov, vcov, ucont, vcont )
2  USE comgeom_mod_h
3  USE parallel_lmdz
4  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
5  USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
6          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
7IMPLICIT NONE
8
9  !=======================================================================
10  !
11  !   Auteur:  P. Le Van
12  !   -------
13  !
14  !   Objet:
15  !   ------
16  !
17  !  *********************************************************************
18  !    calcul des compos. contravariantes a partir des comp.covariantes
19  !  ********************************************************************
20  !
21  !=======================================================================
22
23  INTEGER :: klevel
24  REAL :: ucov( ijb_u:ije_u,klevel ),  vcov( ijb_v:ije_v,klevel )
25  REAL :: ucont( ijb_u:ije_u,klevel ), vcont( ijb_v:ije_v,klevel )
26  INTEGER :: l,ij
27  INTEGER :: ijbu,ijbv,ijeu,ijev
28
29
30  ijbu=ij_begin-iip1
31  ijbv=ij_begin-iip1
32  ijeu=ij_end+iip1
33  ijev=ij_end+iip1
34
35  if (pole_nord) then
36    ijbu=ij_begin+iip1
37    ijbv=ij_begin
38  endif
39
40  if (pole_sud) then
41    ijeu=ij_end-iip1
42    ijev=ij_end-iip1
43  endif
44
45!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
46  DO l = 1,klevel
47
48  DO  ij = ijb_u,ije_u
49  ucont( ij,l ) = ucov( ij,l ) * unscu2( ij )
50  END DO
51
52  DO ij = ijb_v,ije_v
53  vcont( ij,l ) = vcov( ij,l ) * unscv2( ij )
54  END DO
55
56  END DO
57!$OMP END DO NOWAIT
58  RETURN
59END SUBROUTINE covcont_loc
Note: See TracBrowser for help on using the repository browser.