16 |
! Le pas de temps "pdtphys" est supposé beaucoup plus petit que la |
! Le pas de temps "pdtphys" est supposé beaucoup plus petit que la |
17 |
! constante de temps de décroissance. |
! constante de temps de décroissance. |
18 |
|
|
19 |
use dimens_m, only: llm |
use dimens_m, only: llm, nqmx |
20 |
use dimphy, only: klon, nbtr |
use dimphy, only: klon |
21 |
use numer_rec, only: assert |
use nr_util, only: assert |
22 |
|
|
23 |
REAL, intent(in):: tr_seri(:, :, :), pdtphys, tautr(:) |
REAL, intent(in):: tr_seri(:, :, :), pdtphys, tautr(:) |
24 |
real radiornpb(klon, llm, 2) |
real radiornpb(klon, llm, 2) |
28 |
|
|
29 |
!----------------------------------------------- |
!----------------------------------------------- |
30 |
|
|
31 |
call assert(shape(tr_seri) == (/klon, llm, nbtr/), "radiornpb tr_seri") |
call assert(shape(tr_seri) == (/klon, llm, nqmx - 2/), "radiornpb tr_seri") |
32 |
call assert(size(tautr) == nbtr, "radiornpb tautr") |
call assert(size(tautr) == nqmx - 2, "radiornpb tautr") |
33 |
|
|
34 |
DO it = 1, 2 |
DO it = 1, 2 |
35 |
IF (tautr(it) > 0.) THEN |
IF (tautr(it) > 0.) THEN |