|
Last change
on this file since 2639 was
2456,
checked in by aslmd, 5 years ago
|
|
MESOSCALE Mars: fixing a couple of problems with double types related to introduction of CO2 model in r2362
|
|
File size:
711 bytes
|
| Line | |
|---|
| 1 | SUBROUTINE tcondco2(ngrid,nlay,p,q,tcond) |
|---|
| 2 | USE comcstfi_h |
|---|
| 3 | use conc_mod, only: mmean |
|---|
| 4 | |
|---|
| 5 | IMPLICIT NONE |
|---|
| 6 | |
|---|
| 7 | !--------------------------------------------------- |
|---|
| 8 | ! Condensation temperature for co2 ice; based on |
|---|
| 9 | ! the saturation in co2sat.F JA17 |
|---|
| 10 | !-------------------------------------------------- |
|---|
| 11 | |
|---|
| 12 | integer, intent(in) :: ngrid,nlay |
|---|
| 13 | real, intent(in), dimension(ngrid,nlay):: p,q |
|---|
| 14 | double precision, intent(out), dimension(ngrid,nlay):: tcond ! CO2 condensation temperature (atm) |
|---|
| 15 | double precision:: A,B,pco2 |
|---|
| 16 | real :: qco2 |
|---|
| 17 | integer:: ig,l |
|---|
| 18 | |
|---|
| 19 | A=log(1.382d12) |
|---|
| 20 | B=-3182.48 |
|---|
| 21 | qco2=0. |
|---|
| 22 | DO l=1,nlay |
|---|
| 23 | DO ig=1,ngrid |
|---|
| 24 | pco2 = q(ig,l) * (mmean(ig,l)/44.01) * p(ig,l) |
|---|
| 25 | tcond(ig,l)=B/(log(pco2)-A) |
|---|
| 26 | enddo |
|---|
| 27 | enddo |
|---|
| 28 | |
|---|
| 29 | end |
|---|
Note: See
TracBrowser
for help on using the repository browser.
