Parent Directory
| 
Revision Log
Split inter_barxy.f : one procedure per module, one module per file. Grouped the files into a directory. Split orbite.f. Value of raz_date read from the namelist is taken into account (resetting the step counter) even if annee_ref == anneeref and day_ref == dayref. raz_date is no longer modified by gcm main unit. (Following LMDZ.) Removed argument klon of interfsur_lim. Renamed arguments lmt_alb, lmt_rug to alb_new, z0_new (same name as corresponding actual arguments in interfsurf_hq). Removed argument klon of interfsurf_hq. Removed arguments qs and d_qs of diagetpq. Were always zero. Downgraded arguments d_qw, d_ql of diagetpq to local variables, they were not used in physiq. Removed all computations for solid water in diagetpq, was just zero. Downgraded arguments fs_bound, fq_bound of diagphy to local variables, they were not used in physiq. Encapsulated in a test on iprt all computations in diagphy. Removed parameter nbtr of module dimphy. Replaced it everywhere in the program by nqmx - 2. Removed parameter rnpb of procedure physiq. Kept the true case in physiq and phytrac. Could not work with false case anyway. Removed arguments klon, llm, airephy of qcheck. Removed argument ftsol of initrrnpb, was not used.
| 1 | module radiornpb_m |
| 2 | |
| 3 | IMPLICIT none |
| 4 | |
| 5 | contains |
| 6 | |
| 7 | function radiornpb(tr_seri, pdtphys, tautr) |
| 8 | |
| 9 | ! From phylmd/radiornpb.F, v 1.2 2005/05/25 13:10:09 |
| 10 | |
| 11 | ! Auteurs: AA + CG (LGGE/CNRS) Date 24-06-94 |
| 12 | ! Objet: Decroissance radioactive d'un traceur dans l'atmosphere |
| 13 | !G 24 06 94 : Pour un traceur, le radon |
| 14 | !G 16 12 94 : Plus un 2eme traceur, le 210Pb. Le radon decroit en plomb. |
| 15 | |
| 16 | ! Le pas de temps "pdtphys" est supposé beaucoup plus petit que la |
| 17 | ! constante de temps de décroissance. |
| 18 | |
| 19 | use dimens_m, only: llm, nqmx |
| 20 | use dimphy, only: klon |
| 21 | use nr_util, only: assert |
| 22 | |
| 23 | REAL, intent(in):: tr_seri(:, :, :), pdtphys, tautr(:) |
| 24 | real radiornpb(klon, llm, 2) |
| 25 | |
| 26 | ! Variable local to the procedure: |
| 27 | INTEGER it |
| 28 | |
| 29 | !----------------------------------------------- |
| 30 | |
| 31 | call assert(shape(tr_seri) == (/klon, llm, nqmx - 2/), "radiornpb tr_seri") |
| 32 | call assert(size(tautr) == nqmx - 2, "radiornpb tautr") |
| 33 | |
| 34 | DO it = 1, 2 |
| 35 | IF (tautr(it) > 0.) THEN |
| 36 | radiornpb(:, :, it) = - tr_seri(:, :, it) * pdtphys / tautr(it) |
| 37 | ELSE |
| 38 | radiornpb(:, :, it) = 0. |
| 39 | END IF |
| 40 | END DO |
| 41 | |
| 42 | !G161294 : Cas particulier radon 1 => plomb 2 |
| 43 | radiornpb(:, :, 2) = radiornpb(:, :, 2) - radiornpb(:, :, 1) |
| 44 | |
| 45 | END function radiornpb |
| 46 | |
| 47 | end module radiornpb_m |
| ViewVC Help | |
| Powered by ViewVC 1.1.21 |