1 |
|
2 |
#include "ctrparam.h" |
3 |
|
4 |
! ============================================================ |
5 |
! |
6 |
! CHEMTMP1.F: Subroutine for calculating some temperary |
7 |
! parameters for monthly averaging |
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 |
! 103194 Chien Wang rev. |
22 |
! 080200 Chien Wang repack based on CliChem3 & add cpp |
23 |
! |
24 |
! ========================================================== |
25 |
|
26 |
subroutine chemtmp1 (dth,airmass0,p,pvv,xx1) |
27 |
|
28 |
#include "chem_para" |
29 |
#include "chem_tmp" |
30 |
#include "BD2G04.COM" |
31 |
|
32 |
dimension airmass0(nlon,nlat,nlev) |
33 |
dimension pvv (nlon,nlat,nlev) |
34 |
dimension xx1 (nlon,nlat,nlev) |
35 |
|
36 |
! ---------------------------------------------------------- |
37 |
|
38 |
#if ( defined CPL_CHEM ) |
39 |
|
40 |
c |
41 |
c CFC11 N-S transport amount in 10^-9 kg |
42 |
c |
43 |
i = 1 |
44 |
do 1 k=1,nlev |
45 |
cfcnsf(k) = cfcnsf(k) |
46 |
& +( max(0.0,pvv(i,13,k)) |
47 |
& *airmass0(i,12,k)*p(i,12)*xx1(i,12,k) |
48 |
& + min(0.0,pvv(i,13,k)) |
49 |
& *airmass0(i,13,k)*p(i,13)*xx1(i,13,k) ) |
50 |
& /dyv(13)*dth *1.e-9 !kg |
51 |
c & +( (pvv(i,14,k)+pvv(i,13,k))/dyv(13) |
52 |
c & *airmass0(i,13,k)*p(i,13)*xx1(i,13,k) |
53 |
c & -(pvv(i,13,k)+pvv(i,12,k))/dyv(12) |
54 |
c & *airmass0(i,12,k)*p(i,12)*xx1(i,12,k) ) |
55 |
c & *dth*0.5e-9 !kg |
56 |
1 continue |
57 |
|
58 |
#endif |
59 |
|
60 |
return |
61 |
end |
62 |
|
63 |
c==================================================== |
64 |
|
65 |
subroutine chemtmp2 |
66 |
|
67 |
c==================================================================c |
68 |
c c |
69 |
c CHEMTMP1.F: Subroutine for calculating some temperary c |
70 |
c parameters for monthly averaging c |
71 |
c ------------------------------------------------- c |
72 |
c Author: Chien Wang c |
73 |
c MIT Joint Program on Science and Policy c |
74 |
c of Global Change c |
75 |
c Last Revised on: October 31, 1994 c |
76 |
c c |
77 |
c==================================================================c |
78 |
|
79 |
#include "chem_para" |
80 |
#include "chem_tmp" |
81 |
#include "BD2G04.COM" |
82 |
|
83 |
#if ( defined CPL_CHEM ) |
84 |
|
85 |
write(168)cfcnsf |
86 |
|
87 |
do k=1,nlev |
88 |
cfcnsf(k) = 0.0 |
89 |
end do |
90 |
|
91 |
#endif |
92 |
|
93 |
return |
94 |
end |
95 |
|
96 |
|