|
Last change
on this file since 374 was
305,
checked in by rwordsworth, 13 years ago
|
|
Several new files added as part of the climate evolution model
(main program kcm.F90). Some general cleanup in physiq.F90 and
callcorrk.F90. Bugs in dust radiative transfer and H2 Rayleigh
scattering corrected.
|
-
Property svn:executable set to
*
|
|
File size:
446 bytes
|
| Line | |
|---|
| 1 | double precision function cp_neutral(T) |
|---|
| 2 | |
|---|
| 3 | implicit none |
|---|
| 4 | |
|---|
| 5 | #include "gases.h" |
|---|
| 6 | |
|---|
| 7 | ! inputs |
|---|
| 8 | double precision T |
|---|
| 9 | |
|---|
| 10 | |
|---|
| 11 | ! Cp_n : cf CO2 dans abe&matsui (1988) |
|---|
| 12 | !cp_neutral = (22.26+5.981d-2*T-3.501d-5*T**2+7.469d-9*T**3)/m_n |
|---|
| 13 | |
|---|
| 14 | if(gnom(1).eq.'N2_')then |
|---|
| 15 | cp_neutral = 1040.0 |
|---|
| 16 | elseif(gnom(1).eq.'H2_')then |
|---|
| 17 | cp_neutral = 14310.0 |
|---|
| 18 | else |
|---|
| 19 | print*,'Gas not recognised in cp_neutral!' |
|---|
| 20 | call abort |
|---|
| 21 | endif |
|---|
| 22 | |
|---|
| 23 | |
|---|
| 24 | end function cp_neutral |
|---|
| 25 | |
|---|
| 26 | |
|---|
| 27 | |
|---|
| 28 | |
|---|
| 29 | |
|---|
| 30 | |
|---|
| 31 | |
|---|
| 32 | |
|---|
| 33 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
