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 |
|
|
|
! This file is clean: there is neither C preprocessor directive, nor |
|
|
! include line. |
|
|
|
|
5 |
implicit none |
implicit none |
6 |
|
|
7 |
private |
private |
13 |
|
|
14 |
! Author: P. Le Van |
! Author: P. Le Van |
15 |
|
|
16 |
use nrutil, only: assert_eq, assert |
use nr_util, only: assert_eq, assert |
|
|
|
17 |
use dimens_m, only: iim, jjm |
use dimens_m, only: iim, jjm |
18 |
use comgeom, only: aire_2d, apoln, apols |
use comgeom, only: aire_2d, apoln, apols |
19 |
|
|
111 |
! ( L'indice 1 correspond a l'interface mailLE 1 / maille 2) |
! ( L'indice 1 correspond a l'interface mailLE 1 / maille 2) |
112 |
! ( Les abscisses sont exprimées en degres) |
! ( Les abscisses sont exprimées en degres) |
113 |
|
|
114 |
use nrutil, only: assert_eq |
use nr_util, only: assert_eq |
115 |
|
|
116 |
IMPLICIT NONE |
IMPLICIT NONE |
117 |
|
|
296 |
! Version unidimensionnelle, en latitude. |
! Version unidimensionnelle, en latitude. |
297 |
! L'indice 1 correspond à l'interface maille 1 -- maille 2. |
! L'indice 1 correspond à l'interface maille 1 -- maille 2. |
298 |
|
|
299 |
use nrutil, only: assert |
use nr_util, only: assert |
300 |
|
|
301 |
IMPLICIT NONE |
IMPLICIT NONE |
302 |
|
|
345 |
dym = 0. |
dym = 0. |
346 |
jmod = jmod + 1 |
jmod = jmod + 1 |
347 |
ELSE |
ELSE |
348 |
! yjmod(jmod) == yjdat(jdat) |
! {yjmod(jmod) == yjdat(jdat)} |
349 |
dy = yjmod(jmod) - y0 |
dy = yjmod(jmod) - y0 |
350 |
dym = dym + dy |
dym = dym + dy |
351 |
inter_bary(jmod) = (inter_bary(jmod) + dy * fdat(jdat)) / dym |
inter_bary(jmod) = (inter_bary(jmod) + dy * fdat(jdat)) / dym |
373 |
! reverses their order. |
! reverses their order. |
374 |
! Finally, the procedure adds 90° as the last value of the array. |
! Finally, the procedure adds 90° as the last value of the array. |
375 |
|
|
376 |
use nrutil, only: assert_eq |
use nr_util, only: assert_eq, pi |
|
use comconst, only: pi |
|
377 |
|
|
378 |
IMPLICIT NONE |
IMPLICIT NONE |
379 |
|
|
431 |
! This procedure converts to degrees, if necessary, and inverts the |
! This procedure converts to degrees, if necessary, and inverts the |
432 |
! order. |
! order. |
433 |
|
|
434 |
use comconst, only: pi |
use nr_util, only: pi |
435 |
|
|
436 |
IMPLICIT NONE |
IMPLICIT NONE |
437 |
|
|