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

Last change on this file since 38 was 27, checked in by vancop, 8 years ago

First physical fixes, new ice liquid fraction and mushy-layer thermal conductivity
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
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
Note: See TracBrowser for help on using the repository browser.