trunk/dyn3d/gr_int_dyn_m.f

module gr_int_dyn_m

  ! Clean: no C preprocessor directive, no include line

  implicit none

contains

  function gr_int_dyn(champin)

    ! From dyn3d/gr_int_dyn.F,v 1.1.1.1 2004/05/19 12:53:07

    ! Passage d'un champ interpolé à un champ sur grille scalaire

    REAL, intent(in):: champin(:, :)
    REAL gr_int_dyn(size(champin, 1) + 1, size(champin, 2))

    ! Variables local to the procedure:
    integer iim, jp1

    !-----------------------------------------------------------------------

    iim = size(champin, 1)
    jp1 = size(champin, 2)

    gr_int_dyn(:, 1) = sum(champin(:, 1)) / iim ! north pole
    gr_int_dyn(:, jp1) = sum(champin(:, jp1)) / iim ! south pole
    gr_int_dyn(: iim, 2: jp1 - 1) = champin(:, 2: jp1 - 1)
    gr_int_dyn(iim + 1, 2: jp1 - 1) = gr_int_dyn(1, 2: jp1 - 1)

  END function gr_int_dyn

end module gr_int_dyn_m
1	guez	3	module gr_int_dyn_m
2
3			! Clean: no C preprocessor directive, no include line
4
5			implicit none
6
7			contains
8
9			function gr_int_dyn(champin)
10
11			! From dyn3d/gr_int_dyn.F,v 1.1.1.1 2004/05/19 12:53:07
12
13			! Passage d'un champ interpolé à un champ sur grille scalaire
14
15			REAL, intent(in):: champin(:, :)
16			REAL gr_int_dyn(size(champin, 1) + 1, size(champin, 2))
17
18			! Variables local to the procedure:
19			integer iim, jp1
20
21			!-----------------------------------------------------------------------
22
23			iim = size(champin, 1)
24			jp1 = size(champin, 2)
25
26			gr_int_dyn(:, 1) = sum(champin(:, 1)) / iim ! north pole
27			gr_int_dyn(:, jp1) = sum(champin(:, jp1)) / iim ! south pole
28	guez	49	gr_int_dyn(: iim, 2: jp1 - 1) = champin(:, 2: jp1 - 1)
29			gr_int_dyn(iim + 1, 2: jp1 - 1) = gr_int_dyn(1, 2: jp1 - 1)
30	guez	3
31			END function gr_int_dyn
32
33			end module gr_int_dyn_m