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)-s_w ) ! 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 |
---|