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