source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gr_int_dyn.f90 @ 5114

Last change on this file since 5114 was 5113, checked in by abarral, 4 months ago

Rename modules in misc from *_mod > lmdz_*
Put cbrt.f90, ch*.f90, pch*.f90 in new lmdz_libmath_pch.f90
  • 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.1 KB
Line 
1
2! $Header$
3
4SUBROUTINE gr_int_dyn(champin,champdyn,iim,jp1)
5  IMPLICIT NONE
6  !=======================================================================
7  !   passage d'un champ interpole a un champ sur grille scalaire
8  !=======================================================================
9  !-----------------------------------------------------------------------
10  !   declarations:
11  !   -------------
12
13  INTEGER :: iim
14  integer :: ip1, jp1
15  REAL :: champin(iim, jp1)
16  REAL :: champdyn(iim+1, jp1)
17
18  INTEGER :: i, j
19  real :: polenord, polesud
20
21  !-----------------------------------------------------------------------
22  !   calcul:
23  !   -------
24
25  ip1 = iim + 1
26  polenord = 0.
27  polesud = 0.
28  do i = 1, iim
29    polenord = polenord + champin (i, 1)
30    polesud = polesud + champin (i, jp1)
31  enddo
32  polenord = polenord / iim
33  polesud = polesud / iim
34  do j = 1, jp1
35    do i = 1, iim
36      if (j == 1) then
37        champdyn(i, j) = polenord
38      else if (j == jp1) then
39        champdyn(i, j) = polesud
40      else
41        champdyn(i, j) = champin (i, j)
42      endif
43    enddo
44    champdyn(ip1, j) = champdyn(1, j)
45  enddo
46
47  RETURN
48END SUBROUTINE gr_int_dyn
49
Note: See TracBrowser for help on using the repository browser.