Context Navigation

ice_gas_column.f @ 57

Last change on this file since 57 was 57, checked in by vancop, 7 years ago
initial import of source_3.5
File size: 1.7 KB

Line
1	SUBROUTINE ice_gas_column( ntra, zc_i, zdeltaz, ln_write,
2	& zmt_i )
3
4	! gas conservation
5	INCLUDE 'type.com'
6	INCLUDE 'para.com'
7
8	INTEGER ::
9	& jk ! : index for ice layers
10	& jn ! : index for tracers
11
12	REAL(8), DIMENSION(ntra_bio_max) ::
13	& zmt_i ! Total gas content (./m2)
14	REAL(8), DIMENSION(ntra_bio_max,nlay_bio) ::
15	& zc_i ! Input total concentrations (./m3)
16	REAL(8), DIMENSION(nlay_bio) ::
17	& zdeltaz ! Input layer thicknesses
18	LOGICAL ::
19	& ln_write ! Write or not
20
21	!=============================================================================!
22
23	IF ( ln_write ) THEN
24
25	WRITE(numout,*)
26	WRITE(numout,) ' ** ice_gas_column : '
27	WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~ '
28	WRITE(numout,*)
29
30	WRITE(numout,*) ' nlay_bio : ', nlay_bio
31	WRITE(numout,*) ' ntra : ', ntra
32	!DO jn = 1, ntra
33	! WRITE(numout,*) ' zc_i : ', zc_i(jn,1:nlay_bio)
34	!END DO
35
36	ENDIF
37
38	zmt_i(:) = 0.0
39
40	DO jn = 1, ntra
41	DO jk = 1, nlay_bio
42	zmt_i(jn) = zmt_i(jn) + zc_i(jn,jk) *
43	& zdeltaz(jk)
44	END DO
45
46	! IF ( ln_write ) THEN
47	! WRITE(numout,*) ' jn : ', jn
48	! WRITE(numout,*) ' zmt_i : ', zmt_i(jn)
49	! ENDIF
50
51	END DO
52
53	RETURN
54
55	!=============================================================================!
56	!-- End of ice_gas_column --
57
58	END

Note: See TracBrowser for help on using the repository browser.