libf/dyn3d/gr_int_dyn_m.f90

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	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	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
31	END function gr_int_dyn
32
33	end module gr_int_dyn_m