/[lmdze]/trunk/phylmd/radiornpb.f
ViewVC logotype

Diff of /trunk/phylmd/radiornpb.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 265 by guez, Tue Mar 20 09:35:59 2018 UTC revision 266 by guez, Thu Apr 19 17:54:55 2018 UTC
# Line 6  contains Line 6  contains
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
# Line 39  contains Line 42  contains
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

Legend:
Removed from v.265  
changed lines
  Added in v.266

  ViewVC Help
Powered by ViewVC 1.1.21