[4] | 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 |
---|