Changeset 20 for branches/2016/dev_v3.20_2016_gravity_drainage/SOURCES/source_3.20/ice_sal_adv.f

Timestamp:

10/05/16 14:25:40 (8 years ago)

Author:

vancop

Message:

Ludivine source files

File:

• branches/2016/dev_v3.20_2016_gravity_drainage/SOURCES/source_3.20/ice_sal_adv.f (modified) (8 diffs)

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

@@ -43 +43 @@
      &   zb                 ,    !: gravity drainage
      &   zc                 ,    !: upward advective flow
-     &   w_adv_br           ,    !: brine velocity
      &   zRae               ,    !: effective Ra
      &   ze                 ,    !: downward advective flow
@@ -49 +48 @@
      &   zindtbis           ,    !:
      &   zdiagbis                !:
+
 
       REAL(8), DIMENSION(nlay_i,3) :: !: dummy factors for tracer equation
@@ -75 +75 @@
      &   zlevel             ,    !: height of the water column
      &   zthdiff                 !: thermal diffusivity
-
+ 
       REAL(8), DIMENSION(nlay_i+1) ::  
      &   z_sbr_int               !: brine salinity at layer interfaces
@@ -300 +300 @@
       rho_br_GN = 1020.
       Rc_GN     = 1.01
+      !alpha_GN = 1.0e-3
+      !Rc_GN    = 1.0
 
       w_adv_br(:) = 0.0
@@ -322 +324 @@
 !------------------------------------------------------------------------------|
 !
-
+  
       DO layer = 1, nlay_i
          za(layer) = w_adv_br(layer) * ddtb / deltaz_i_phy(layer)
@@ -330 +332 @@
       sn_i_b(1) = z_sbr_i(1) * ( e_i_b(1) + za(1) ) + 
      &            z_sbr_i(2) * ( - za(1) ) 
+
 
       ! inner layers
@@ -342 +345 @@
      &                   z_sbr_i(nlay_i)   * ( e_i_b(nlay_i) + 
      &                   za(nlay_i)/2. ) - za(nlay_i) * oce_sal
+ 
 
       IF ( ln_write ) THEN
          WRITE(numout,*) 
-         WRITE(numout,*)
          WRITE(numout,*) ' sn_i_b   : ', ( sn_i_b(layer) , 
      &                   layer = 1, nlay_i )
@@ -384 +387 @@
 
       ENDIF ! ln_sal
+
+
 !
 !------------------------------------------------------------------------------|