source: LMDZ6/branches/Amaury_dev/libf/phylmd/radio_decay.F90 @ 5116

Last change on this file since 5116 was 5103, checked in by abarral, 4 months ago

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.8 KB
RevLine 
[5099]1
[1191]2! $Id $
[5099]3
[1191]4SUBROUTINE radio_decay(radio,rnpb,dtime,tautr,tr,d_tr)
[5099]5
[1191]6! Caluclate radioactive decay for all tracers with radio(it)=true
[5099]7
[1191]8  USE dimphy
[5101]9  USE infotrac_phy, ONLY: nbtr
10  USE traclmdz_mod, ONLY: id_rn, id_pb
[1191]11  IMPLICIT NONE
12!-----------------------------------------------------------------------
13! Auteur(s): AA + CG (LGGE/CNRS) Date 24-06-94
14! Objet: Calcul de la tendance radioactive des traceurs type radioelements
[1409]15!        Cas particulier pour le couple radon-plomb : Le radon decroit en plomb
[1191]16!-----------------------------------------------------------------------
[5099]17
[1191]18! Entrees
[5099]19
[5103]20  LOGICAL,DIMENSION(nbtr),INTENT(IN)        :: radio ! .TRUE. = traceur radioactif
21  LOGICAL,INTENT(IN)                        :: rnpb  ! .TRUE. = decroissance RN = source PB
[1191]22  REAL,INTENT(IN)                           :: dtime ! Pas de temps physique (secondes)
23  REAL,DIMENSION(nbtr),INTENT(IN)           :: tautr ! Constante de decroissance radioactive
24  REAL,DIMENSION(klon,klev,nbtr),INTENT(IN) :: tr    ! Concentrations traceurs U/kgA
[5099]25
[1191]26! Sortie
[5099]27
[1191]28  REAL,DIMENSION(klon,klev,nbtr),INTENT(OUT) :: d_tr  ! Tendance de decroissance radioactive
[5099]29
[1191]30! Locales
[5099]31
[1191]32  INTEGER  :: i,k,it
[524]33
[1191]34
35  DO it = 1,nbtr
[1409]36     d_tr(:,:,it) = 0.
[1191]37     IF ( radio(it) ) THEN
[5082]38        IF (tautr(it) > 0.) THEN
[1191]39           DO k = 1,klev
40              DO i = 1,klon
41                 d_tr(i,k,it) = - tr(i,k,it) * dtime / tautr(it)
42              END DO
43           END DO
44        END IF
45     END IF
46  END DO
[1409]47
[1191]48!-------------------------------------------------------
[1409]49! Cas particulier radon (id_rn) => plomb (id_pb)
[1191]50!-------------------------------------------------------
51  IF ( rnpb ) THEN
52     DO k = 1,klev
[524]53        DO i = 1,klon
[1409]54           d_tr(i,k,id_pb) = d_tr(i,k,id_pb) - d_tr(i,k,id_rn)
[524]55        ENDDO
[1191]56     ENDDO
57  ENDIF
58
59END SUBROUTINE radio_decay
Note: See TracBrowser for help on using the repository browser.