1 |
module inter_barxy_m |
module inter_barxy_m |
2 |
|
|
3 |
! From inter_barxy.F, version 1.1.1.1 2004/05/19 12:53:07 |
! From inter_barxy.F, version 1.1.1.1, 2004/05/19 12:53:07 |
4 |
|
|
5 |
implicit none |
implicit none |
6 |
|
|
11 |
! Author: P. Le Van |
! Author: P. Le Van |
12 |
|
|
13 |
use comgeom, only: aire_2d, apoln, apols |
use comgeom, only: aire_2d, apoln, apols |
14 |
use dimens_m, only: iim, jjm |
use dimensions, only: iim, jjm |
15 |
use inter_barx_m, only: inter_barx |
use inter_barx_m, only: inter_barx |
16 |
use inter_bary_m, only: inter_bary |
use inter_bary_m, only: inter_bary |
17 |
use nr_util, only: assert_eq, assert |
use nr_util, only: assert_eq, assert |
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 |
31 |
! jjm+1 latitudes rlatu du modele (latitudes des scalaires et de U) |
! jjm+1 latitudes rlatu du modele (latitudes des scalaires et de U) |
32 |
! Si taille de la seconde dim = jjm, on veut interpoler sur les |
! Si taille de la seconde dim = jjm, on veut interpoler sur les |
33 |
! jjm latitudes rlatv du modèle (latitudes de V) |
! jjm latitudes rlatv du mod\`ele (latitudes de V) |
34 |
|
|
35 |
! Local: |
! Local: |
36 |
|
|
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, & |