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

Last change on this file since 3728 was 3663, checked in by gmilcareck, 4 months ago

Generic PCM:
Cleaning up the code. The gas constant and Avogadro number values
was not the same between the files in the model.
We choose the value recommended by the 2019 revision of the SI.
R=8.314463 J.K-1.mol-1 and NA=6.022141e23 mol-1.
GM

File size: 775 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
[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.