source: LMDZ4/trunk/libf/dyn3d/gr_int_dyn.F @ 5048

Last change on this file since 5048 was 524, checked in by lmdzadmin, 20 years ago

Initial revision
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
RevLine 
[524]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.