| Rev | Line | |
|---|
| [1816] | 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 |
|---|
| [2341] | 10 | !-------------------------------------------------- |
|---|
| [1816] | 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 |
|---|
| [2341] | 16 | real :: qco2 |
|---|
| [1816] | 17 | integer:: ig,l |
|---|
| 18 | |
|---|
| [1932] | 19 | A=dlog(1.382d12) |
|---|
| [1816] | 20 | B=-3182.48 |
|---|
| [2341] | 21 | qco2=0. |
|---|
| [1816] | 22 | DO l=1,nlay |
|---|
| 23 | DO ig=1,ngrid |
|---|
| [2341] | 24 | ! qco2 = min(1., q(ig,l))!added by CM not sure, to be verified before |
|---|
| 25 | ! cleaning |
|---|
| 26 | ! pco2 = qco2 * (mmean(ig,l)/44.01) * p(ig,l) |
|---|
| [1816] | 27 | pco2 = q(ig,l) * (mmean(ig,l)/44.01) * p(ig,l) |
|---|
| 28 | tcond(ig,l)=B/(dlog(pco2)-A) |
|---|
| 29 | enddo |
|---|
| 30 | enddo |
|---|
| 31 | |
|---|
| 32 | end |
|---|
Note: See
TracBrowser
for help on using the repository browser.
