source: LMDZ.3.3/trunk/libf/dyn3d/gr_int_dyn.F @ 346

Last change on this file since 346 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
RevLine 
[2]1      subroutine gr_int_dyn(champin,champdyn,iim,jp1)
2      implicit none
3c=======================================================================
4c   passage d'un champ interpole a un champ sur grille scalaire
5c=======================================================================
6c-----------------------------------------------------------------------
7c   declarations:
8c   -------------
9
10      INTEGER iim
11      integer ip1, jp1
12      REAL champin(iim, jp1)
13      REAL champdyn(iim+1, jp1)
14
15      INTEGER i, j
16      real polenord, polesud
17
18c-----------------------------------------------------------------------
19c   calcul:
20c   -------
21
22      ip1 = iim + 1
23      polenord = 0.
24      polesud = 0.
25      do i = 1, iim
26        polenord = polenord + champin (i, 1)
27        polesud = polesud + champin (i, jp1)
28      enddo
29      polenord = polenord / iim
30      polesud = polesud / iim
31      do j = 1, jp1
32        do i = 1, iim
33          if (j .eq. 1) then
34            champdyn(i, j) = polenord
35          else if (j .eq. jp1) then
36            champdyn(i, j) = polesud
37          else
38            champdyn(i, j) = champin (i, j)
39          endif
40        enddo
41        champdyn(ip1, j) = champdyn(1, j)
42      enddo
43
44      RETURN
45      END
46
Note: See TracBrowser for help on using the repository browser.