igsm/src/chemtemco2.F


#include "ctrparam.h"

!       ==========================================================
!
!       CHEMTEMCO2.F:  Subroutine for deduce CO2 amount caused by
!                  land (TEM)  uptake of MIT Global Chemistry Model
!
!       ----------------------------------------------------------
!
!       Author:         Chien Wang
!                       MIT Joint Program on Science and Policy
!                               of Global Change
!
!       ----------------------------------------------------------
! 
!       Revision History:
!       
!       When    Who             What
!       ----    ----------      ------- 
!       101200                  developed from chemoceanco2
!
!       ==========================================================

!       ==============================
        subroutine chemtemco2 (fnep)
!       ==============================

#include "chem_para"
#include "chem_com"
#include "BD2G04.COM"

        dimension fnep(nlat)    !flux toward land of carbon in Gt

!       ------------------------------------------------------------

#if ( defined CPL_CHEM) && ( defined CPL_TEM )

        xconv = 44.e6           !from mole(C) to 10^-9 kg(CO2)
     &   *1.e15/12.            !from Gt (C) to mole(C)
c----
        i     = 1

c 052295:
        ktop  = 2
     
        do 1 j=3,nlat   ! throwout 2 south pole points
c
c  Calculate total mass of air and CO2:
c
          tmass    = 0.0                !Total air mass in kg
          xmass    = 0.0                !Total co2 mass in 10^-9 kg
        do k=1,ktop
          tmass    = tmass + airmass(i,j,k)
          xmass    = xmass
     &             + airmass(i,j,k)*zco2(i,j,k)
        enddo
          tmass1   = 1./tmass
c
c Convert flux from mole(C) to ppb(CO2) in mass:
c
          yamount = fnep(j)*xconv
c
c Map CO2 sink into two layer:
c
          xemi     = max(0.0, (xmass - yamount))
     &                            * tmass1     !ppbm

        do k=1,ktop
          zco2(i,j,k) = xemi
        enddo

1       continue

#endif

        return
         end

1	jscott	1.1
2			#include "ctrparam.h"
3
4			! ==========================================================
5			!
6			! CHEMTEMCO2.F: Subroutine for deduce CO2 amount caused by
7			! land (TEM) uptake of MIT Global Chemistry Model
8			!
9			! ----------------------------------------------------------
10			!
11			! Author: Chien Wang
12			! MIT Joint Program on Science and Policy
13			! of Global Change
14			!
15			! ----------------------------------------------------------
16			!
17			! Revision History:
18			!
19			! When Who What
20			! ---- ---------- -------
21			! 101200 developed from chemoceanco2
22			!
23			! ==========================================================
24
25			! ==============================
26			subroutine chemtemco2 (fnep)
27			! ==============================
28
29			#include "chem_para"
30			#include "chem_com"
31			#include "BD2G04.COM"
32
33			dimension fnep(nlat) !flux toward land of carbon in Gt
34
35			! ------------------------------------------------------------
36
37			#if ( defined CPL_CHEM) && ( defined CPL_TEM )
38
39			xconv = 44.e6 !from mole(C) to 10^-9 kg(CO2)
40			& *1.e15/12. !from Gt (C) to mole(C)
41			c----
42			i = 1
43
44			c 052295:
45			ktop = 2
46
47			do 1 j=3,nlat ! throwout 2 south pole points
48			c
49			c Calculate total mass of air and CO2:
50			c
51			tmass = 0.0 !Total air mass in kg
52			xmass = 0.0 !Total co2 mass in 10^-9 kg
53			do k=1,ktop
54			tmass = tmass + airmass(i,j,k)
55			xmass = xmass
56			& + airmass(i,j,k)*zco2(i,j,k)
57			enddo
58			tmass1 = 1./tmass
59			c
60			c Convert flux from mole(C) to ppb(CO2) in mass:
61			c
62			yamount = fnep(j)*xconv
63			c
64			c Map CO2 sink into two layer:
65			c
66			xemi = max(0.0, (xmass - yamount))
67			& * tmass1 !ppbm
68
69			do k=1,ktop
70			zco2(i,j,k) = xemi
71			enddo
72
73			1 continue
74
75			#endif
76
77			return
78			end
79