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 |