Changeset 9750 for NEMO/trunk/src/ICE/icethd_sal.F90

Timestamp:

2018-06-06T14:49:34+02:00 (6 years ago)

Author:

clem

Message:

make ice ponds working in debug mode

File:

• NEMO/trunk/src/ICE/icethd_sal.F90 (modified) (1 diff)

NEMO/trunk/src/ICE/icethd_sal.F90

@@ -77 +77 @@
             !---------------------------------------------------------
             IF( h_i_1d(ji) > 0._wp ) THEN
-               zs_sni   = sss_1d(ji) * ( rhoic - rhosn ) * r1_rhoic                                 ! Salinity of snow ice
-               zs_i_si = ( zs_sni      - s_i_1d(ji) ) *             dh_snowice(ji)  / h_i_1d(ji) ! snow-ice    
-               zs_i_bg = ( s_i_new(ji) - s_i_1d(ji) ) * MAX( 0._wp, dh_i_bott(ji) ) / h_i_1d(ji) ! bottom growth
+               zs_sni   = sss_1d(ji) * ( rhoic - rhosn ) * r1_rhoic                 ! Salinity of snow ice
+               zs_i_si = ( zs_sni      - s_i_1d(ji) ) * dh_snowice(ji) / h_i_1d(ji) ! snow-ice    
+               zs_i_bg = ( s_i_new(ji) - s_i_1d(ji) ) * dh_i_bog  (ji) / h_i_1d(ji) ! bottom growth
                ! Update salinity (nb: salt flux already included in icethd_dh)
                s_i_1d(ji) = s_i_1d(ji) + zs_i_bg + zs_i_si