6 |
|
|
7 |
function radiornpb(tr_seri, pdtphys, tautr) |
function radiornpb(tr_seri, pdtphys, tautr) |
8 |
|
|
9 |
! From phylmd/radiornpb.F, v 1.2 2005/05/25 13:10:09 |
! From phylmd/radiornpb.F, version 1.2, 2005/05/25 13:10:09 |
10 |
|
|
11 |
! Auteurs: AA + CG (LGGE/CNRS) Date 24-06-94 |
! Auteurs : AA + Christophe Genthon (LGGE/CNRS) |
12 |
! Objet: Decroissance radioactive d'un traceur dans l'atmosphere |
! Date: June 1994 |
13 |
!G 24 06 94 : Pour un traceur, le radon |
! Objet: d\'ecroissance radioactive d'un traceur dans l'atmosph\`ere |
14 |
!G 16 12 94 : Plus un 2eme traceur, le 210Pb. Le radon decroit en plomb. |
|
15 |
|
! Pour un traceur, le radon. Plus un deuxi\`eme traceur, le 210 |
16 |
! Le pas de temps "pdtphys" est supposé beaucoup plus petit que la |
! Pb. Le radon d\'ecro\^it en plomb. Le pas de temps "pdtphys" est |
17 |
! constante de temps de décroissance. |
! suppos\'e beaucoup plus petit que la constante de temps de |
18 |
|
! d\'ecroissance. |
19 |
|
|
20 |
use dimensions, only: llm, nqmx |
use dimensions, only: llm, nqmx |
21 |
use dimphy, only: klon |
use dimphy, only: klon |
22 |
use nr_util, only: assert |
use nr_util, only: assert |
23 |
|
|
24 |
REAL, intent(in):: tr_seri(:, :, :), pdtphys, tautr(:) |
REAL, intent(in):: tr_seri(:, :, :) ! (klon, llm, nqmx - 2) |
25 |
|
REAL, intent(in):: pdtphys |
26 |
|
REAL, intent(in):: tautr(:) ! (nqmx - 2) |
27 |
real radiornpb(klon, llm, 2) |
real radiornpb(klon, llm, 2) |
28 |
|
|
29 |
! Variable local to the procedure: |
! Local: |
30 |
INTEGER it |
INTEGER it |
31 |
|
|
32 |
!----------------------------------------------- |
!----------------------------------------------- |
33 |
|
|
34 |
call assert(shape(tr_seri) == (/klon, llm, nqmx - 2/), "radiornpb tr_seri") |
call assert(shape(tr_seri) == [klon, llm, nqmx - 2], "radiornpb tr_seri") |
35 |
call assert(size(tautr) == nqmx - 2, "radiornpb tautr") |
call assert(size(tautr) == nqmx - 2, "radiornpb tautr") |
36 |
|
|
37 |
DO it = 1, 2 |
DO it = 1, 2 |
42 |
END IF |
END IF |
43 |
END DO |
END DO |
44 |
|
|
45 |
!G161294 : Cas particulier radon 1 => plomb 2 |
! Cas particulier radon 1 => plomb 2 |
46 |
radiornpb(:, :, 2) = radiornpb(:, :, 2) - radiornpb(:, :, 1) |
radiornpb(:, :, 2) = radiornpb(:, :, 2) - radiornpb(:, :, 1) |
47 |
|
|
48 |
END function radiornpb |
END function radiornpb |