trunk/dyn3d/gr_int_dyn_m.f90

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	guez	3	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	guez	49	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	guez	3
29			END function gr_int_dyn
30
31			end module gr_int_dyn_m