source: LMDZ6/trunk/libf/dyn3d_common/gr_v_scal.f90 @ 5274

Last change on this file since 5274 was 5272, checked in by abarral, 8 days ago

Turn paramet.h into a module
  • 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.3 KB
RevLine 
[524]1!
2! $Header$
3!
[5246]4SUBROUTINE gr_v_scal(nx,x_v,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  !=======================================================================
[5271]27  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
[5272]28USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
29          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
[5271]30IMPLICIT NONE
[5246]31  !-----------------------------------------------------------------------
32  !   Declararations:
33  !   ---------------
[524]34
[5271]35
[5272]36
[5246]37  INCLUDE "comgeom.h"
[524]38
[5246]39  !   Arguments:
40  !   ----------
[524]41
[5246]42  INTEGER :: nx
43  REAL :: x_v(ip1jm,nx),x_scal(ip1jmp1,nx)
[524]44
[5246]45  !   Local:
46  !   ------
[524]47
[5246]48  INTEGER :: l,ij
[524]49
[5246]50  !-----------------------------------------------------------------------
[524]51
[5246]52  DO l=1,nx
53     DO ij=iip2,ip1jm
54        x_scal(ij,l)= &
55              (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l)) &
56              /(airev(ij-iip1)+airev(ij))
57     ENDDO
58     DO ij=1,iip1
59        x_scal(ij,l)=0.
60     ENDDO
61     DO ij=ip1jm+1,ip1jmp1
62        x_scal(ij,l)=0.
63     ENDDO
64  ENDDO
[524]65
[5246]66  RETURN
67END SUBROUTINE gr_v_scal
Note: See TracBrowser for help on using the repository browser.