Changeset 6470 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

Timestamp:

2016-04-13T17:41:02+02:00 (8 years ago)

Author:

clem

Message:

solve issue #1712

File:

• trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90 (modified) (1 diff)

trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

@@ -650 +650 @@
             sfx_dyn(ji,jj) = sfx_dyn(ji,jj) - smsw(ij) * rhoic * r1_rdtice
             wfx_dyn(ji,jj) = wfx_dyn(ji,jj) - vsw (ij) * rhoic * r1_rdtice   ! increase in ice volume due to seawater frozen in voids
+
+            ! virtual salt flux to keep salinity constant
+            IF( nn_icesal == 1 .OR. nn_icesal == 3 )  THEN
+               srdg2(ij)      = srdg2(ij) - vsw(ij) * ( sss_m(ji,jj) - sm_i(ji,jj,jl1) )           ! ridge salinity = sm_i
+               sfx_bri(ji,jj) = sfx_bri(ji,jj) + sss_m(ji,jj)    * vsw(ij) * rhoic * r1_rdtice  &  ! put back sss_m into the ocean
+                  &                            - sm_i(ji,jj,jl1) * vsw(ij) * rhoic * r1_rdtice     ! and get  sm_i  from the ocean 
+            ENDIF
             
             !------------------------------------------