Changeset 5202 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

Timestamp:

2015-04-07T17:40:16+02:00 (9 years ago)

Author:

clem

Message:

LIM3: important bug fix to avoid crashes. See ticket #1508

File:

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

trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

@@ -819 +819 @@
       DO jk = 1, nlay_i             ! Sea ice energy of melting
          DO ji = kideb, kiut
-            ztmelts      = - tmut  * s_i_1d(ji,jk) + rt0 
-            rswitch      = MAX( 0._wp , SIGN( 1._wp , -(t_i_1d(ji,jk) - rt0) - epsi20 ) )
-            q_i_1d(ji,jk) = rhoic * ( cpic * ( ztmelts - t_i_1d(ji,jk) )                                                &
-               &                   + lfus * ( 1.0 - rswitch * ( ztmelts-rt0 ) / MIN( t_i_1d(ji,jk) - rt0, -epsi20 ) )   &
-               &                   - rcp  *                   ( ztmelts-rt0 )  ) 
+            ztmelts      = - tmut  * s_i_1d(ji,jk) + rt0
+            t_i_1d(ji,jk) = MIN( t_i_1d(ji,jk), ztmelts ) ! Force t_i_1d to be lower than melting point
+                                                          !   (sometimes dif scheme produces abnormally high temperatures)   
+            q_i_1d(ji,jk) = rhoic * ( cpic * ( ztmelts - t_i_1d(ji,jk) )                           &
+               &                    + lfus * ( 1.0 - ( ztmelts-rt0 ) / ( t_i_1d(ji,jk) - rt0 ) )   &
+               &                    - rcp  *         ( ztmelts-rt0 )  ) 
          END DO
       END DO