Changeset 4099 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

Timestamp:

2013-10-22T14:07:21+02:00 (11 years ago)

Author:

clem

Message:

 

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90 (modified) (3 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

@@ -84 +84 @@
       REAL(wp) ::   zdhnm, zhnnew, zhisn, zihic, zzc       !
       REAL(wp) ::   zfracs       ! fractionation coefficient for bottom salt entrapment
-      REAL(wp) ::   zds          ! increment of bottom ice salinity
       REAL(wp) ::   zcoeff       ! dummy argument for snowfall partitioning over ice and leads
       REAL(wp) ::   zsm_snowice  ! snow-ice salinity
@@ -429 +428 @@
                   zfracs = zswi1  * 0.12 + zswi12 * ( 0.8925 + 0.0568 * LOG( 100.0 * zgrr ) )   &
                      &                   + zswi2  * 0.26 / ( 0.26 + 0.74 * EXP ( - 724300.0 * zgrr ) ) 
-                  zds         = zfracs * sss_m(ii,ij) - s_i_new(ji)
+                  zfracs = MIN( 0.5 , zfracs )
                   s_i_new(ji) = zfracs * sss_m(ii,ij)
                ENDIF ! fc_bo_i
@@ -438 +437 @@
          DO ji = kideb, kiut
             IF( ( fc_bo_i(ji) + fbif_1d(ji) + qlbbq_1d(ji) ) .LT. 0.0  ) THEN
-               ! New ice salinity must not exceed 15 psu
+               ! New ice salinity must not exceed 20 psu
                s_i_new(ji) = MIN( s_i_new(ji), s_i_max )
                ! Metling point in K