Changeset 8522 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90

Timestamp:

2017-09-14T17:52:02+02:00 (7 years ago)

Author:

clem

Message:

changes in style - part6 - ice diffusion (hope the scheme still converges)

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90 (modified) (2 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90

@@ -280 +280 @@
       !!-------------------------------------------------------------------
       INTEGER  ::   ji, jk   ! dummy loop indices
-      REAL(wp) ::   ztmelts, z1_2cp, zbbb, zccc  ! local scalar 
+      REAL(wp) ::   ztmelts, zbbb, zccc  ! local scalar 
       !!-------------------------------------------------------------------
       ! Recover ice temperature
-      z1_2cp = 1._wp / ( 2._wp * cpic )
       DO jk = 1, nlay_i
          DO ji = 1, nidx
@@ -290 +289 @@
             zbbb          = ( rcp - cpic ) * ztmelts + e_i_1d(ji,jk) * r1_rhoic - lfus
             zccc          = SQRT( MAX( zbbb * zbbb - 4._wp * cpic * lfus * ztmelts, 0._wp ) )
-            t_i_1d(ji,jk) = rt0 - ( zbbb + zccc ) * z1_2cp
+            t_i_1d(ji,jk) = rt0 - ( zbbb + zccc ) * 0.5_wp * r1_cpic
             
             ! mask temperature