Changeset 9888 for NEMO/trunk/src/ICE

Timestamp:

2018-07-06T12:33:01+02:00 (6 years ago)

Author:

clem

Message:

improve and debug BDY with sea ice. With this commit there should not be anymore problems in regional configurations (hopefully).

Location:

NEMO/trunk/src/ICE

Files:

• icedyn_adv.F90 (modified) (1 diff)
• icevar.F90 (modified) (1 diff)

NEMO/trunk/src/ICE/icedyn_adv.F90

@@ -105 +105 @@
       ! ==> 1) remove negative ice areas and volumes (conservation is ensure)
       CALL ice_var_zapsmall 
-      ! ==> 2) remove remaining negative advected fields (conservation is not preserved)
+      ! ==> 2) remove remaining negative advected fields (conservation is not preserved) => conservation issue
       WHERE( v_s (:,:,:)   < 0._wp )   v_s (:,:,:)   = 0._wp
       WHERE( sv_i(:,:,:)   < 0._wp )   sv_i(:,:,:)   = 0._wp

NEMO/trunk/src/ICE/icevar.F90

@@ -853 +853 @@
          DO ji = 1, npti
             ztmelts      = - tmut  * sz_i_1d(ji,jk)
-            t_i_1d(ji,jk) = MIN( t_i_1d(ji,jk), ztmelts + rt0 ) ! Force t_i_1d to be lower than melting point
+            t_i_1d(ji,jk) = MIN( t_i_1d(ji,jk), ztmelts + rt0 ) ! Force t_i_1d to be lower than melting point => likely conservation issue
                                                                 !   (sometimes zdf scheme produces abnormally high temperatures)   
             e_i_1d(ji,jk) = rhoic * ( cpic * ( ztmelts - ( t_i_1d(ji,jk) - rt0 ) )           &