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