Changeset 30

Timestamp:

10/13/16 09:18:27 (8 years ago)

Author:

vancop

Message:

cleaning

Location:

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

Files:

• ice.com (modified) (2 diffs)
• ice_phy_param.f (modified) (3 diffs)
• ice_sal_adv.f (modified) (4 diffs)

branches/2016/dev_v3.20_2016_platelet/SOURCES/source_3.20/ice.com

@@ -179 +179 @@
       COMMON/heateqcoe/lfus,cpg,cpw,gammac,betak,tmut,tpw,
      &              cpoc,deltah, betak1, betak2, visc_br, thdiff_br,
-     &              beta_ocs
+     &              beta_ocs, rho_br
 
       COMMON/barrowconf/sal_read(11), hi_read(11), hgins, hnins,
@@ -190 +190 @@
       COMMON/tuneforcing/sf_mult, tabq_ano
 
+      REAL(8) ::
+     &   alpha_GN 
+
       COMMON/fluidtpt/flu_beta, rad_io,
      &                frtr_si_phy, qsummer, d_br_mol, d_br_tur,
      &                ra_c, ra_smooth, e_thr_flu, delta_cw, 
+     &                alpha_GN,
      &                rn_e_newice, ini_swi, 
      &                nn_thcond,

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

@@ -112 +112 @@
       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
+      rho_br  = 1020.             ! brine density
       beta_ocs  = 0.81            ! regulates density changes due to changes in salinity
 
@@ -196 +197 @@
       READ(25,'(a3)') c_permeff   ! type of effective permeability
       READ(25,*)
-      READ(25,*) delta_cw         ! Cox and weeks gravity drainage parameter
+      READ(25,*) alpha_GN         ! Griewank and Notz gravity drainage parameter
+      READ(25,*)
+      READ(25,*) delta_CW         ! Cox and weeks gravity drainage parameter
       READ(25,*)
       READ(25,*) rn_e_newice      ! Liquid fraction in new ice
@@ -261 +264 @@
       WRITE(numout,*) ' ra_smooth  :', ra_smooth
       WRITE(numout,*) ' c_gravdr   :', c_gravdr
-      WRITE(numout,*) ' delta_cw   :', delta_cw
+      WRITE(numout,*) ' alpha_GN   :', alpha_GN
+      WRITE(numout,*) ' delta_CW   :', delta_CW
       WRITE(numout,*) ' rn_e_newice:', rn_e_newice
       WRITE(numout,*) ' nn_thcond  :', nn_thcond

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

@@ -6 +6 @@
       !! ** Purpose :
       !!        This routine computes new salinities in the ice
+      !!        Based on the method of Griewank and Notz (2013) &
+      !!        Rees Jones and Worster (2014)
+      !!        dS/dt = -w dSbr / dz
+      !!        with a trick to retrieve the lost salt
       !!
       !! ** Method  : Vertical salinity profile computation 
@@ -18 +22 @@
       !! ** External
       !!
-      !! ** References : Vancop. et al., 2008
+      !! ** References : Griewank and Notz JGR 2013
+      !!                 Rees Jones & Worster JGR 2014
       !!
       !! ** History  : 
-      !!    (06-2003) Martin Vancop. LIM1D
-      !!    (06-2008) Martin Vancop. BIO-LIM
-      !!    (09-2008) Martin Vancop. Explicit gravity drainage
+      !!    (03-2016) - Martin Vancoppenolle + Max Thomas + Ludivine Conte: initial version
+      !!    (10-2016) - Martin Vancoppenolle + Pat Wongpan : revision
       !!
       !!------------------------------------------------------------------
@@ -48 +52 @@
      &   zindtbis           ,    !:
      &   zdiagbis                !:
-
 
       REAL(8), DIMENSION(nlay_i,3) :: !: dummy factors for tracer equation
@@ -279 +282 @@
 
       alpha_GN  = 1.56e-3
-      rho_br_GN = 1020.
-      Rc_GN     = 1.01
-      !alpha_GN = 1.0e-3
-      !Rc_GN    = 1.0
 
       w_adv_br(:) = 0.0
       zRae(:) = 0.
       DO layer = 1, nlay_i
-         zRae(layer) = MAX( rayleigh(layer) - Rc_GN, 0.) ! correction from the orig scheme
+         zRae(layer) = MAX( rayleigh(layer) - ra_c, 0.) ! correction from the orig scheme
       END DO
 
       DO layer = 1, nlay_i
-         w_adv_br(layer) = - alpha_GN / rho_br_GN * 
+         w_adv_br(layer) = - alpha_GN / rho_br * 
      &      SUM ( zRae(1:layer)*deltaz_i_phy(1:layer) )
       END DO