/[MITgcm]/MITgcm_contrib/jscott/igsm/src_chem/chemtmp.F
ViewVC logotype

Contents of /MITgcm_contrib/jscott/igsm/src_chem/chemtmp.F

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


Revision 1.1 - (show annotations) (download)
Thu Sep 17 17:40:33 2009 UTC (15 years, 10 months ago) by jscott
Branch: MAIN
CVS Tags: HEAD
chem module archive

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

  ViewVC Help
Powered by ViewVC 1.1.22