1 | SUBROUTINE ice_sal_column( kideb,kiut,zmt_i,zs_i,zdeltaz,nlay_i, |
---|
2 | & ln_write ) |
---|
3 | |
---|
4 | !-----------------------------------------------------------------------! |
---|
5 | ! This routine sums the mass of salt over the whole sea ice column |
---|
6 | ! (c) Martin Vancoppenolle, September 2008 |
---|
7 | !-----------------------------------------------------------------------! |
---|
8 | INCLUDE 'type.com' |
---|
9 | INCLUDE 'para.com' |
---|
10 | |
---|
11 | INTEGER :: |
---|
12 | & ji , ! : index for space |
---|
13 | & jk ! : index for ice layers |
---|
14 | |
---|
15 | REAL(8) :: |
---|
16 | & zmt_i ! : total mass of salt |
---|
17 | REAL(8), DIMENSION(maxnlay) :: |
---|
18 | & zs_i ! : bulk salinity |
---|
19 | REAL(8), DIMENSION(maxnlay) :: |
---|
20 | & zdeltaz ! : thickness of the layers |
---|
21 | |
---|
22 | LOGICAL :: |
---|
23 | & ln_write |
---|
24 | |
---|
25 | IF ( ln_write ) THEN |
---|
26 | WRITE(numout,*) |
---|
27 | WRITE(numout,*) ' *** ice_sal_column : ' |
---|
28 | WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~ ' |
---|
29 | WRITE(numout,*) ' kideb : ', kideb |
---|
30 | WRITE(numout,*) ' kiut : ', kiut |
---|
31 | WRITE(numout,*) ' zmt_i : ', zmt_i |
---|
32 | WRITE(numout,*) ' zs_i : ', ( zs_i(jk), jk = 1, nlay_i ) |
---|
33 | WRITE(numout,*) ' zdeltaz:', ( zdeltaz(jk), jk = 1, nlay_i ) |
---|
34 | ENDIF |
---|
35 | |
---|
36 | zmt_i = 0.0 |
---|
37 | DO jk = 1, nlay_i |
---|
38 | zmt_i = zmt_i + zs_i(jk) * zdeltaz(jk) |
---|
39 | END DO |
---|
40 | |
---|
41 | IF ( ln_write ) THEN |
---|
42 | WRITE(numout,*) ' zmt_i : ', zmt_i |
---|
43 | ENDIF |
---|
44 | |
---|
45 | |
---|
46 | !=============================================================================! |
---|
47 | !-- End of ice_sal_column -- |
---|
48 | RETURN |
---|
49 | |
---|
50 | END |
---|