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

Last change on this file since 862 was 776, checked in by Laurent Fairhead, 17 years ago

Suite du merge entre la version et la HEAD: quelques modifications
de Yann sur le

LF

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