Changeset 29 for branches/2016/dev_v3.20_2016_platelet/SOURCES/source_3.20/ice_phy_param.f

Timestamp:

10/12/16 18:08:11 (8 years ago)

Author:

vancop

Message:

Add ice_sal_adv routine for GN13

File:

• branches/2016/dev_v3.20_2016_platelet/SOURCES/source_3.20/ice_phy_param.f (modified) (3 diffs)

branches/2016/dev_v3.20_2016_platelet/SOURCES/source_3.20/ice_phy_param.f

@@ -110 +110 @@
       s_w     = 34.0              ! 34. is the control
 
-      visc_br = 1.79e-3           ! dynamic viscosity of water at 0C 
-      beta_ocs= 0.81              ! regulates density changes due to changes in salinity
+      visc_br = 1.9e-3            ! dynamic viscosity of water at 0C (kg/(m.s)) - Sharqawy et al, 2010, eq. 22
+      thdiff_br = 1.2e-7          ! thermal diffusivity of brine (m2/s) - Notz and Worster, 2008
+      beta_ocs  = 0.81            ! regulates density changes due to changes in salinity
 
       !-----------------
@@ -187 +188 @@
       READ(25,*) ra_smooth        ! coefficient to smooth the hyperbolic tangential for ra_c
       READ(25,*)
-      READ(25,'(a2)') gravdr      ! type of gravity drainage
+      READ(25,'(a2)') c_gravdr    ! type of gravity drainage
+      READ(25,*)
+      READ(25,'(a3)') c_sbr       ! type of brine salinity for salt advection
+      READ(25,*)
+      READ(25,'(a3)') c_perm      ! type of permeability 
+      READ(25,*)
+      READ(25,'(a3)') c_permeff   ! type of effective permeability
       READ(25,*)
       READ(25,*) delta_cw         ! Cox and weeks gravity drainage parameter
@@ -253 +260 @@
       WRITE(numout,*) ' ra_c       :', ra_c
       WRITE(numout,*) ' ra_smooth  :', ra_smooth
-      WRITE(numout,*) ' gravdr     :', gravdr
+      WRITE(numout,*) ' c_gravdr   :', c_gravdr
       WRITE(numout,*) ' delta_cw   :', delta_cw
       WRITE(numout,*) ' rn_e_newice:', rn_e_newice