source:
trunk/LMDZ.GENERIC/libf/phystd/blackl.F
@
3728
Last change on this file since 3728 was 3663, checked in by , 4 months ago | |
---|---|
File size: 775 bytes |
Rev | Line | |
---|---|---|
[135] | 1 | subroutine blackl(blalong,blat,blae) |
[253] | 2 | |
[135] | 3 | implicit double precision (a-h,o-z) |
[253] | 4 | |
5 | ! physical constants | |
[3663] | 6 | sigma=5.670374D-8 |
[135] | 7 | pi=datan(1.d0)*4.d0 |
[3663] | 8 | c0=2.997925d+08 |
9 | h=6.62607d-34 | |
10 | cbol=1.380649d-23 | |
[135] | 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 | |
[3654] | 22 | |
23 | subroutine blackn(blalong,blat,blae) | |
24 | ||
25 | implicit double precision (a-h,o-z) | |
26 | ||
27 | ! physical constants | |
[3663] | 28 | sigma=5.670374D-8 |
[3654] | 29 | pi=datan(1.d0)*4.d0 |
[3663] | 30 | c0=2.997925d+08 |
31 | h=6.62607d-34 | |
32 | cbol=1.380649d-23 | |
[3654] | 33 | rind=1.d0 |
34 | c=c0/rind | |
35 | c1=h*(c**2) | |
36 | c2=h*c/cbol | |
37 | ||
38 | ||
39 | blae=2.d0*pi*c1*blalong**3/(dexp(c2*blalong/blat)-1.d0) | |
40 | ||
41 | ||
42 | return | |
43 | end |
Note: See TracBrowser
for help on using the repository browser.