|
Last change
on this file since 135 was
135,
checked in by aslmd, 14 years ago
|
|
CHANGEMENT ARBORESCENCE ETAPE 2 -- NON COMPLET
|
|
File size:
1022 bytes
|
| Rev | Line | |
|---|
| [135] | 1 | subroutine calc_cpp3d(cppNI,rcpNI,t,p) |
|---|
| 2 | |
|---|
| 3 | !================================================================== |
|---|
| 4 | ! Purpose |
|---|
| 5 | ! ------- |
|---|
| 6 | ! Compute the atmospheric specific heat capacity as a |
|---|
| 7 | ! function of pressure and temperature (for CO2 gas only) |
|---|
| 8 | ! |
|---|
| 9 | ! Authors |
|---|
| 10 | ! ------- |
|---|
| 11 | ! Robin Wordsworth (2009) |
|---|
| 12 | ! |
|---|
| 13 | !================================================================== |
|---|
| 14 | |
|---|
| 15 | ! you still need to validate this equation vs something.... |
|---|
| 16 | |
|---|
| 17 | implicit none |
|---|
| 18 | |
|---|
| 19 | #include "comcstfi.h" |
|---|
| 20 | |
|---|
| 21 | real cp0, dB2dT2 |
|---|
| 22 | real cppNI ! specific heat capacity at const. pressure |
|---|
| 23 | real rcpNI ! R / specific heat capacity |
|---|
| 24 | real t |
|---|
| 25 | real p |
|---|
| 26 | real calmol |
|---|
| 27 | real bar |
|---|
| 28 | |
|---|
| 29 | calmol = 94.9784 ! 4.18/(mco2/1000) |
|---|
| 30 | bar = 100000 |
|---|
| 31 | |
|---|
| 32 | cp0 = (7.70+5.3e-3*t - 8.3e-7*t**2)*calmol |
|---|
| 33 | dB2dT2 = 2.69e-5*t - 0.0098 |
|---|
| 34 | cppNI = cp0 - t*(p/bar)*dB2dT2 |
|---|
| 35 | !cppNI = 1000*(t/460)**0.35 ! Sebastian's version |
|---|
| 36 | |
|---|
| 37 | rcpNI = R/cppNI |
|---|
| 38 | |
|---|
| 39 | return |
|---|
| 40 | end subroutine calc_cpp3d |
|---|
Note: See
TracBrowser
for help on using the repository browser.
