igsm/src_chem/chemstrat.F


#include "ctrparam.h"

!       ============================================================
!
!       CHEMEMISSION.F:  Subroutine for calculating stratospheric
!                         decuction of MIT Global Chemistry Model
!
!       ------------------------------------------------------------
!
!       Author:         Chien Wang
!                       MIT Joint Program on Science and Policy
!                               of Global Change
!
!       ----------------------------------------------------------
! 
!       Revision History:
!       
!       When    Who             What
!       ----    ----------      -------
!       102694  Chien Wang      rev.    
!       080200  Chien Wang      repack based on CliChem3 & add cpp
!
!       ==========================================================

      subroutine chemstrat (dtb)

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

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

#if ( defined CPL_CHEM )

c------------------------
c  Deduction of CFC11
c    in stratosphere:
c
      dta = dtb *3.0
      i=1
      k=nlev

c 4/19/95:

      tf11=1./3.1536e8*dta      !10 yr
c      tf12=1./1.173139e9*dta
c      tn2o=1./1.110067e9*dta

      tf12=tf11
      tn2o=tf11

      do j=1,nlat
c        cfc11sd     = cfc11sd 
c     &              + airmass0(i,j,k)*p(i,j)
c     &              * tf11 * cfc11(i,j,k)
        cfc11(i,j,k)= cfc11(i,j,k)
     &              * (1.-tf11)

c        cfc12sd     = cfc12sd 
c     &              + airmass0(i,j,k)*p(i,j)
c     &              * tf12 * cfc12(i,j,k)
        cfc12(i,j,k)= cfc12(i,j,k)
     &              * (1.-tf12)

c        xn2osd      = xn2osd 
c     &              + airmass0(i,j,k)*p(i,j)
c     &              * tn2o * xn2o(i,j,k)
        xn2o (i,j,k)= xn2o(i,j,k)
     &              * (1.-tn2o)
      enddo

#endif

      return
       end
1	jscott	1.1
2			#include "ctrparam.h"
3
4			! ============================================================
5			!
6			! CHEMEMISSION.F: Subroutine for calculating stratospheric
7			! decuction 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			! 102694 Chien Wang rev.
22			! 080200 Chien Wang repack based on CliChem3 & add cpp
23			!
24			! ==========================================================
25
26			subroutine chemstrat (dtb)
27
28			#include "chem_para"
29			#include "chem_com"
30			#include "BD2G04.COM"
31
32			! ----------------------------------------------------------
33
34			#if ( defined CPL_CHEM )
35
36			c------------------------
37			c Deduction of CFC11
38			c in stratosphere:
39			c
40			dta = dtb *3.0
41			i=1
42			k=nlev
43
44			c 4/19/95:
45
46			tf11=1./3.1536e8*dta !10 yr
47			c tf12=1./1.173139e9*dta
48			c tn2o=1./1.110067e9*dta
49
50			tf12=tf11
51			tn2o=tf11
52
53			do j=1,nlat
54			c cfc11sd = cfc11sd
55			c & + airmass0(i,j,k)*p(i,j)
56			c & * tf11 * cfc11(i,j,k)
57			cfc11(i,j,k)= cfc11(i,j,k)
58			& * (1.-tf11)
59
60			c cfc12sd = cfc12sd
61			c & + airmass0(i,j,k)*p(i,j)
62			c & * tf12 * cfc12(i,j,k)
63			cfc12(i,j,k)= cfc12(i,j,k)
64			& * (1.-tf12)
65
66			c xn2osd = xn2osd
67			c & + airmass0(i,j,k)*p(i,j)
68			c & * tn2o * xn2o(i,j,k)
69			xn2o (i,j,k)= xn2o(i,j,k)
70			& * (1.-tn2o)
71			enddo
72
73			#endif
74
75			return
76			end