source: trunk/LMDZ.GENERIC/libf/phystd/blackl.F @ 1150

Last change on this file since 1150 was 959, checked in by jleconte, 12 years ago

15/05/2013 == JL

  • correction in radiative scheme to enforce double precision
File size: 381 bytes
RevLine 
[135]1      subroutine blackl(blalong,blat,blae)
[253]2
[135]3      implicit double precision (a-h,o-z)
[253]4
5      ! physical constants
[959]6      sigma=5.67032D-8
[135]7      pi=datan(1.d0)*4.d0
8      c0=2.9979d+08
9      h=6.6262d-34
10      cbol=1.3806d-23
11      rind=1.d0
12      c=c0/rind
13      c1=h*(c**2)
14      c2=h*c/cbol
[253]15
16
[135]17      blae=2.d0*pi*c1/blalong**5/(dexp(c2/blalong/blat)-1.d0)
[253]18
19
[135]20      return
21      end
Note: See TracBrowser for help on using the repository browser.