trunk/dyn3d/gr_int_dyn_m.f

module gr_int_dyn_m

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