source: LMDZ4/trunk/libf/phylmd/radiornpb.F @ 593

Last change on this file since 593 was 524, checked in by lmdzadmin, 20 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.5 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE radiornpb(tr,dtime,tautr,d_tr)
5      IMPLICIT none
6c======================================================================
7c Auteur(s): AA + CG (LGGE/CNRS) Date 24-06-94
8c Objet: Decroissance radioactive d'un traceur dans l'atmosphere
9CG240694 : Pour un traceur, le radon
10CG161294 : Plus un 2eme traceur, le 210Pb. Le radon decroit en plomb.
11c======================================================================
12c Arguments:
13c======================================================================
14#include "dimensions.h"
15#include "dimphy.h"
16c======================================================================
17C
18      INTEGER i , k , it
19      REAL tr(klon,klev,nbtr) , d_tr(klon,klev,nbtr)
20      REAL dtime
21      REAL tautr(nbtr)
22C
23      WRITE(*,'(''PASSAGE radiornpb ... '',$)')
24C Attention, pour un pas de temps beaucoup plus petit que la decroissance!!!
25
26      DO it = 1,2
27           IF ( tautr(it) .GT. 0. ) THEN
28                DO k = 1,klev
29                DO i = 1,klon
30                d_tr(i,k,it) = - tr(i,k,it) * dtime / tautr(it)
31                END DO
32                END DO
33           ELSE
34                DO k = 1,klev
35                DO i = 1,klon
36                d_tr(i,k,it) = 0.
37                END DO
38                END DO
39           END IF
40      END DO
41C
42CG161294 : Cas particulier radon 1 => plomb 2
43c
44      DO k = 1,klev
45        DO i = 1,klon
46          d_tr(i,k,2) = d_tr(i,k,2) - d_tr(i,k,1)
47        ENDDO
48      ENDDO
49c
50      WRITE(*,*) ' radiornpb OK'
51c
52      RETURN
53      END
Note: See TracBrowser for help on using the repository browser.