Ignore:
Timestamp:
2017-09-14T17:52:02+02:00 (3 years ago)
Author:
clem
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90

    r8518 r8522  
    280280      !!------------------------------------------------------------------- 
    281281      INTEGER  ::   ji, jk   ! dummy loop indices 
    282       REAL(wp) ::   ztmelts, z1_2cp, zbbb, zccc  ! local scalar  
     282      REAL(wp) ::   ztmelts, zbbb, zccc  ! local scalar  
    283283      !!------------------------------------------------------------------- 
    284284      ! Recover ice temperature 
    285       z1_2cp = 1._wp / ( 2._wp * cpic ) 
    286285      DO jk = 1, nlay_i 
    287286         DO ji = 1, nidx 
     
    290289            zbbb          = ( rcp - cpic ) * ztmelts + e_i_1d(ji,jk) * r1_rhoic - lfus 
    291290            zccc          = SQRT( MAX( zbbb * zbbb - 4._wp * cpic * lfus * ztmelts, 0._wp ) ) 
    292             t_i_1d(ji,jk) = rt0 - ( zbbb + zccc ) * z1_2cp 
     291            t_i_1d(ji,jk) = rt0 - ( zbbb + zccc ) * 0.5_wp * r1_cpic 
    293292             
    294293            ! mask temperature 
Note: See TracChangeset for help on using the changeset viewer.