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