source: LMDZ6/branches/SETHET_DECOUPLE/libf/dyn3d_common/gr_int_dyn.F @ 5440

Last change on this file since 5440 was 1945, checked in by lguez, 11 years ago

Duplicated files moved to dyn3d_common.
  • 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.2 KB
Line 
1!
2! $Header$
3!
4      subroutine gr_int_dyn(champin,champdyn,iim,jp1)
5      implicit none
6c=======================================================================
7c   passage d'un champ interpole a un champ sur grille scalaire
8c=======================================================================
9c-----------------------------------------------------------------------
10c   declarations:
11c   -------------
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
21c-----------------------------------------------------------------------
22c   calcul:
23c   -------
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 .eq. 1) then
37            champdyn(i, j) = polenord
38          else if (j .eq. 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
48      END
49
Note: See TracBrowser for help on using the repository browser.