/[MITgcm]/MITgcm_contrib/jscott/igsm/src/chemtemco2.F
ViewVC logotype

Annotation of /MITgcm_contrib/jscott/igsm/src/chemtemco2.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (hide annotations) (download)
Fri Aug 11 19:35:30 2006 UTC (18 years, 11 months ago) by jscott
Branch: MAIN
CVS Tags: HEAD
atm2d package

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    

  ViewVC Help
Powered by ViewVC 1.1.22