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
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