igsm/src_chem/chemoceanco2.F


#include "ctrparam.h"

!       ==========================================================
!
!       CHEMOCEANCO2.F:  Subroutine for deduce CO2 amount caused by
!                          ocean uptake of MIT Global Chemistry Model
!
!       ----------------------------------------------------------
!
!       Author:         Chien Wang
!                       MIT Joint Program on Science and Policy
!                               of Global Change
!
!       ----------------------------------------------------------
! 
!       Revision History:
!       
!       When    Who             What
!       ----    ----------      ------- 
!       102596  Chien Wang      rev.
!       080100  Chien Wang      repack based on CliChem3 & add cpp
!
!       ==========================================================

!       ==============================
        subroutine chemoceanco2 (fco2)
!       ==============================

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

        dimension fco2(nlat)    !flux toward ocean of co2 in mole(C)

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

#if ( defined CPL_CHEM) && ( defined CPL_OCEANCO2 )

        xconv = 44.e6         !from mole(C) to 10^-9 kg(CO2)
c----
        i     = 1

c 052295:
        ktop  = 2
     
#if ( N_LAT == 24 )
        do  j=3,nlat    ! throwout 2 south pole points
#endif
#if ( N_LAT == 46 )
        do  j=4,nlat    ! throwout 3 south pole points
#endif
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 = fco2(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

        enddo

#endif

        return
         end

1
2	#include "ctrparam.h"
3
4	! ==========================================================
5	!
6	! CHEMOCEANCO2.F: Subroutine for deduce CO2 amount caused by
7	! ocean 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	! 102596 Chien Wang rev.
22	! 080100 Chien Wang repack based on CliChem3 & add cpp
23	!
24	! ==========================================================
25
26	! ==============================
27	subroutine chemoceanco2 (fco2)
28	! ==============================
29
30	#include "chem_para"
31	#include "chem_com"
32	#include "BD2G04.COM"
33
34	dimension fco2(nlat) !flux toward ocean of co2 in mole(C)
35
36	! ------------------------------------------------------------
37
38	#if ( defined CPL_CHEM) && ( defined CPL_OCEANCO2 )
39
40	xconv = 44.e6 !from mole(C) to 10^-9 kg(CO2)
41	c----
42	i = 1
43
44	c 052295:
45	ktop = 2
46
47	#if ( N_LAT == 24 )
48	do j=3,nlat ! throwout 2 south pole points
49	#endif
50	#if ( N_LAT == 46 )
51	do j=4,nlat ! throwout 3 south pole points
52	#endif
53	c
54	c Calculate total mass of air and CO2:
55	c
56	tmass = 0.0 !Total air mass in kg
57	xmass = 0.0 !Total co2 mass in 10^-9 kg
58	do k=1,ktop
59	tmass = tmass + airmass(i,j,k)
60	xmass = xmass
61	& + airmass(i,j,k)*zco2(i,j,k)
62	enddo
63	tmass1 = 1./tmass
64	c
65	c Convert flux from mole(C) to ppb(CO2) in mass:
66	c
67	yamount = fco2(j)*xconv
68	c
69	c Map CO2 sink into two layer:
70	c
71	xemi = max(0.0, (xmass - yamount))
72	& * tmass1 !ppbm
73
74	do k=1,ktop
75	zco2(i,j,k) = xemi
76	enddo
77
78	enddo
79
80	#endif
81
82	return
83	end
84