source: LMDZ6/trunk/libf/dyn3dmem/gr_u_scal_loc.f90 @ 5435

Last change on this file since 5435 was 5285, checked in by abarral, 2 months 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.4 KB
RevLine 
[1632]1!
2! $Header$
3!
[5246]4SUBROUTINE gr_u_scal_loc(nx,x_u,x_scal)
5  !%W%    %G%
6  !=======================================================================
7  !
8  !   Author:    Frederic Hourdin      original: 11/11/92
9  !   -------
10  !
11  !   Subject:
12  !   ------
13  !
14  !   Method:
15  !   --------
16  !
17  !   Interface:
18  !   ----------
19  !
20  !  Input:
21  !  ------
22  !
23  !  Output:
24  !  -------
25  !
26  !=======================================================================
[5281]27  USE comgeom_mod_h
[5246]28  USE parallel_lmdz
[5271]29  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
[5285]30USE paramet_mod_h
[5271]31IMPLICIT NONE
[5246]32  !-----------------------------------------------------------------------
33  !   Declararations:
34  !   ---------------
[1632]35
[5271]36
[5272]37
[1632]38
[5246]39  !   Arguments:
40  !   ----------
[1632]41
[5246]42  INTEGER :: nx
43  REAL :: x_u(ijb_u:ije_u,nx),x_scal(ijb_u:ije_u,nx)
[1632]44
[5246]45  !   Local:
46  !   ------
[1632]47
[5246]48  INTEGER :: l,ij
49  INTEGER :: ijb,ije
[1632]50
[5246]51  !-----------------------------------------------------------------------
52  ijb=ij_begin
53  ije=ij_end
[1632]54
[5246]55!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
56  DO l=1,nx
57     DO ij=ijb+1,ije
58        x_scal(ij,l)= &
59              (aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l)) &
60              /(aireu(ij)+aireu(ij-1))
61     ENDDO
62  ENDDO
[1632]63!$OMP ENDDO NOWAIT
64
[5246]65  ijb=ij_begin
66  ije=ij_end
[1632]67
[5246]68!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
69  DO l=1,nx
70     DO ij=ijb,ije-iip1+1,iip1
71       x_scal(ij,l)=x_scal(ij+iip1-1,l)
72     ENDDO
73  ENDDO
[1632]74!$OMP ENDDO NOWAIT
[5246]75  RETURN
76
77END SUBROUTINE gr_u_scal_loc
Note: See TracBrowser for help on using the repository browser.