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

Last change on this file since 5274 was 5246, checked in by abarral, 11 days ago

Convert fixed-form to free-form sources .F -> .{f,F}90
(WIP: some .F remain, will be handled in subsequent commits)
  • 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
RevLine 
[524]1!
2! $Header$
3!
[5246]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  !   -------------
[524]12
[5246]13  INTEGER :: iim
14  integer :: ip1, jp1
15  REAL :: champin(iim, jp1)
16  REAL :: champdyn(iim+1, jp1)
[524]17
[5246]18  INTEGER :: i, j
19  real :: polenord, polesud
[524]20
[5246]21  !-----------------------------------------------------------------------
22  !   calcul:
23  !   -------
[524]24
[5246]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
[524]46
[5246]47  RETURN
48END SUBROUTINE gr_int_dyn
[524]49
Note: See TracBrowser for help on using the repository browser.