/[lmdze]/trunk/phylmd/Radlwsw/radlwsw.f90
ViewVC logotype

Diff of /trunk/phylmd/Radlwsw/radlwsw.f90

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

revision 336 by guez, Thu Jun 13 14:40:06 2019 UTC revision 337 by guez, Mon Sep 16 16:54:50 2019 UTC
# Line 5  module radlwsw_m Line 5  module radlwsw_m
5  contains  contains
6    
7    SUBROUTINE radlwsw(dist, mu0, fract, paprs, play, tsol, albedo, t, q, wo, &    SUBROUTINE radlwsw(dist, mu0, fract, paprs, play, tsol, albedo, t, q, wo, &
8         cldfra, cldemi, cldtaupd, heat, heat0, cool, cool0, radsol, albpla, &         cldfra, cldemi, cldtau, heat, heat0, cool, cool0, radsol, albpla, &
9         topsw, toplw, solsw, sollw, sollwdown, topsw0, toplw0, solsw0, sollw0, &         topsw, toplw, solsw, sollw, sollwdown, topsw0, toplw0, solsw0, sollw0, &
10         lwdn0, lwdn, lwup0, lwup, swdn0, swdn, swup0, swup, ok_ade, topswad, &         lwdn0, lwdn, lwup0, lwup, swdn0, swdn, swup0, swup, ok_ade, topswad, &
11         solswad)         solswad)
# Line 47  contains Line 47  contains
47      real, intent(in):: cldemi(klon, klev)      real, intent(in):: cldemi(klon, klev)
48      ! emissivite des nuages dans l'IR (entre 0 et 1)      ! emissivite des nuages dans l'IR (entre 0 et 1)
49    
50      real, intent(in):: cldtaupd(klon, klev)      real, intent(in):: cldtau(klon, klev)
51      ! epaisseur optique des nuages dans le visible (present-day value)      ! \'epaisseur optique des nuages dans le visible (present-day value)
52    
53      real, intent(out):: heat(klon, klev)      real, intent(out):: heat(klon, klev)
54      ! échauffement atmosphérique (visible) (K/jour)      ! échauffement atmosphérique (visible) (K/jour)
# Line 187  contains Line 187  contains
187               PCLDLD(i, k) = cldfra(iof+i, k)*cldemi(iof+i, k)               PCLDLD(i, k) = cldfra(iof+i, k)*cldemi(iof+i, k)
188               PCLDLU(i, k) = cldfra(iof+i, k)*cldemi(iof+i, k)               PCLDLU(i, k) = cldfra(iof+i, k)*cldemi(iof+i, k)
189               PCLDSW(i, k) = cldfra(iof+i, k)               PCLDSW(i, k) = cldfra(iof+i, k)
190               PTAU(i, 1, k) = MAX(cldtaupd(iof+i, k), 1e-05)               PTAU(i, 1, k) = MAX(cldtau(iof+i, k), 1e-05)
191               ! (1e-12 serait instable)               ! (1e-12 serait instable)
192               PTAU(i, 2, k) = MAX(cldtaupd(iof+i, k), 1e-05)               PTAU(i, 2, k) = MAX(cldtau(iof+i, k), 1e-05)
193               ! (pour 32-bit machines)               ! (pour 32-bit machines)
194               POMEGA(i, 1, k) = 0.9999 - 5e-04 * EXP(-0.5 * PTAU(i, 1, k))               POMEGA(i, 1, k) = 0.9999 - 5e-04 * EXP(-0.5 * PTAU(i, 1, k))
195               POMEGA(i, 2, k) = 0.9988 - 2.5e-03 * EXP(-0.05 * PTAU(i, 2, k))               POMEGA(i, 2, k) = 0.9988 - 2.5e-03 * EXP(-0.05 * PTAU(i, 2, k))
Legend:
Removed from v.336  
changed lines
  Added in v.337
  ViewVC Help
Powered by ViewVC 1.1.21