source: codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/cp_neutral.F90

Last change on this file was 313, checked in by ymipsl, 10 years ago
  • implement splitting of XIOS file for lmdz physics
  • Termination is done properly in parallel by calling MPI_ABORT instead of abort or stop

YM
  • Property svn:executable set to *
File size: 648 bytes
Line 
1double precision function cp_neutral(T)
2
3  use gases_h
4
5  implicit none
6
7  ! inputs
8  double precision T
9
10
11  ! this function has been disabled in gradients_kcm.F90 because it doesnt
12  ! work if you have gaseous mixtures. need to decide whether to generalise
13  ! it or simply remove entirely...
14
15  ! Cp_n : cf CO2 dans abe&matsui (1988)
16  !cp_neutral = (22.26+5.981d-2*T-3.501d-5*T**2+7.469d-9*T**3)/m_n
17
18  if(trim(gnom(1)).eq.'N2_')then
19     cp_neutral = 1040.0
20  elseif(trim(gnom(1)).eq.'H2_')then
21     cp_neutral = 14310.0
22  else
23     print*,'Gas not recognised in cp_neutral!'
24     call abort_physiq
25  endif
26
27 
28end function cp_neutral
29
30
31
32
33
34
35
36
37
Note: See TracBrowser for help on using the repository browser.