[4] | 1 | SUBROUTINE ice_brine |
---|
| 2 | |
---|
| 3 | !------------------------------------------------------------------------------! |
---|
| 4 | ! *** ice_brine *** |
---|
| 5 | ! |
---|
| 6 | ! This routine computes brine density (kg/m3) and salinity (g/kg) |
---|
| 7 | ! from K-temperature |
---|
| 8 | ! in each vertical layer |
---|
| 9 | ! |
---|
| 10 | ! Refs: Notz (2005) and Vancoppenolle et al (TCD 2015) |
---|
| 11 | ! |
---|
| 12 | !------------------------------------------------------------------------------! |
---|
| 13 | |
---|
| 14 | USE lib_fortran |
---|
| 15 | |
---|
| 16 | INCLUDE 'type.com' |
---|
| 17 | INCLUDE 'para.com' |
---|
| 18 | INCLUDE 'const.com' |
---|
| 19 | INCLUDE 'ice.com' |
---|
| 20 | INCLUDE 'thermo.com' |
---|
| 21 | INCLUDE 'bio.com' |
---|
| 22 | |
---|
| 23 | !==============================================================================! |
---|
| 24 | |
---|
| 25 | ! |
---|
| 26 | !------------------------------------------------------------------------------! |
---|
| 27 | ! X) Starting the routine |
---|
| 28 | !------------------------------------------------------------------------------! |
---|
| 29 | ! |
---|
| 30 | |
---|
| 31 | WRITE(numout,*) |
---|
| 32 | WRITE(numout,*) ' *** ice_brine : ' |
---|
| 33 | WRITE(numout,*) ' ~~~~~~~~~~~~~~~~' |
---|
| 34 | WRITE(numout,*) |
---|
| 35 | |
---|
| 36 | ! |
---|
| 37 | !------------------------------------------------------------------------------! |
---|
| 38 | ! X) Brine characteristics |
---|
| 39 | !------------------------------------------------------------------------------! |
---|
| 40 | ! |
---|
| 41 | |
---|
| 42 | DO layer = 1, nlay_bio |
---|
| 43 | |
---|
| 44 | tc_bio(layer) = t_i_bio(layer) - 273.15 ! Temperature in celsius |
---|
| 45 | zt2 = tc_bio(layer) * tc_bio(layer) |
---|
| 46 | zt3 = tc_bio(layer) * zt2 |
---|
| 47 | |
---|
| 48 | sbr_bio(layer) = -21.4*tc_bio(layer) - 0.886*zt2 -0.0107*zt3 ! Brine 3rd order salinity |
---|
| 49 | |
---|
| 50 | rhobr_bio(layer) = rho0 + beta_ocs*( sbr_bio(layer)-oce_sal ) ! Brine density |
---|
| 51 | |
---|
| 52 | END DO |
---|
| 53 | |
---|
| 54 | ! |
---|
| 55 | !------------------------------------------------------------------------------! |
---|
| 56 | ! X) End of the routine |
---|
| 57 | !------------------------------------------------------------------------------! |
---|
| 58 | ! |
---|
| 59 | WRITE(numout,*) |
---|
| 60 | WRITE(numout,*) ' End of ice_brine ' |
---|
| 61 | WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' |
---|
| 62 | WRITE(numout,*) |
---|
| 63 | |
---|
| 64 | RETURN |
---|
| 65 | |
---|
| 66 | END |
---|