19 |
use ord_coordm_m, only: ord_coordm |
use ord_coordm_m, only: ord_coordm |
20 |
|
|
21 |
REAL, intent(in):: dlonid(:) |
REAL, intent(in):: dlonid(:) |
22 |
! (longitude from input file, in rad, from -pi to pi) |
! longitude from input file, in rad, from -pi to pi |
23 |
|
|
24 |
REAL, intent(in):: dlatid(:), champ(:, :), rlonimod(:) |
REAL, intent(in):: dlatid(:), champ(:, :), rlonimod(:) |
25 |
|
|
26 |
REAL, intent(in):: rlatimod(:) |
REAL, intent(in):: rlatimod(:) |
27 |
! (latitude angle, in degrees or rad, in strictly decreasing order) |
! latitude angle, in degrees or rad, in strictly decreasing order |
28 |
|
|
29 |
real, intent(out):: champint(:, :) |
real, intent(out):: champint(:, :) |
30 |
! Si taille de la seconde dim = jjm + 1, on veut interpoler sur les |
! Si taille de la seconde dim = jjm + 1, on veut interpoler sur les |
45 |
|
|
46 |
!----------------------------------- |
!----------------------------------- |
47 |
|
|
48 |
jnterfd = assert_eq(size(champ, 2) - 1, size(dlatid), & |
jnterfd = assert_eq(size(champ, 2) - 1, size(dlatid), "inter_barxy jnterfd") |
|
"inter_barxy jnterfd") |
|
49 |
jmods = size(champint, 2) |
jmods = size(champint, 2) |
50 |
call assert(size(champ, 1) == size(dlonid), "inter_barxy size(champ, 1)") |
call assert(size(champ, 1) == size(dlonid), "inter_barxy size(champ, 1)") |
51 |
call assert((/size(rlonimod), size(champint, 1)/) == iim, & |
call assert((/size(rlonimod), size(champint, 1)/) == iim, & |