|
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
|
| Line | |
|---|
| 1 | subroutine blackl(blalong,blat,blae) |
|---|
| 2 | |
|---|
| 3 | implicit double precision (a-h,o-z) |
|---|
| 4 | |
|---|
| 5 | ! physical constants |
|---|
| 6 | sigma=5.670374D-8 |
|---|
| 7 | pi=datan(1.d0)*4.d0 |
|---|
| 8 | c0=2.997925d+08 |
|---|
| 9 | h=6.62607d-34 |
|---|
| 10 | cbol=1.380649d-23 |
|---|
| 11 | rind=1.d0 |
|---|
| 12 | c=c0/rind |
|---|
| 13 | c1=h*(c**2) |
|---|
| 14 | c2=h*c/cbol |
|---|
| 15 | |
|---|
| 16 | |
|---|
| 17 | blae=2.d0*pi*c1/blalong**5/(dexp(c2/blalong/blat)-1.d0) |
|---|
| 18 | |
|---|
| 19 | |
|---|
| 20 | return |
|---|
| 21 | end |
|---|
| 22 | |
|---|
| 23 | subroutine blackn(blalong,blat,blae) |
|---|
| 24 | |
|---|
| 25 | implicit double precision (a-h,o-z) |
|---|
| 26 | |
|---|
| 27 | ! physical constants |
|---|
| 28 | sigma=5.670374D-8 |
|---|
| 29 | pi=datan(1.d0)*4.d0 |
|---|
| 30 | c0=2.997925d+08 |
|---|
| 31 | h=6.62607d-34 |
|---|
| 32 | cbol=1.380649d-23 |
|---|
| 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.
