source: branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_brine.f @ 67

Last change on this file since 67 was 39, checked in by vancop, 7 years ago

add tank mass and salt balance
File size: 2.0 KB
Line 
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      INCLUDE 'tank.com'
23
24!==============================================================================!
25
26!
27!------------------------------------------------------------------------------!
28! X) Starting the routine
29!------------------------------------------------------------------------------!
30!
31
32      WRITE(numout,*)
33      WRITE(numout,*) ' *** ice_brine : '
34      WRITE(numout,*) ' ~~~~~~~~~~~~~~~~'
35      WRITE(numout,*)
36
37!
38!------------------------------------------------------------------------------!
39! X) Brine characteristics
40!------------------------------------------------------------------------------!
41!
42         
43      DO layer = 1, nlay_bio
44         
45         tc_bio(layer)    = t_i_bio(layer) - 273.15                     ! Temperature in celsius
46         zt2 = tc_bio(layer) * tc_bio(layer)
47         zt3 = tc_bio(layer) * zt2
48         
49         sbr_bio(layer)   = -21.4*tc_bio(layer) - 0.886*zt2 -0.0107*zt3 ! Brine 3rd order salinity
50         
51         rhobr_bio(layer) = rho0 + beta_ocs*( sbr_bio(layer)-s_w )  ! Brine density
52           
53      END DO
54
55!
56!------------------------------------------------------------------------------!
57! X) End of the routine
58!------------------------------------------------------------------------------!
59!
60      WRITE(numout,*)
61      WRITE(numout,*) ' End of ice_brine '
62      WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
63      WRITE(numout,*)
64     
65      RETURN
66
67      END
Note: See TracBrowser for help on using the repository browser.