Changeset 10564


Ignore:
Timestamp:
2019-01-23T23:46:34+01:00 (18 months ago)
Author:
clem
Message:

correct a bug in sea-ice heat flux

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ICE/icedyn.F90

    r10538 r10564  
    240240                  zhs = v_s(ji,jj,jl) / a_i(ji,jj,jl) 
    241241                  IF( v_s(ji,jj,jl) > 0._wp .AND. zhs > phs_max(ji,jj,jl) .AND. a_i(ji,jj,jl) < 0.15 ) THEN 
    242                      zfra = a_i(ji,jj,jl) * phs_max(ji,jj,jl) / MAX( v_s(ji,jj,jl), epsi20 ) 
     242                     zfra = phs_max(ji,jj,jl) / MAX( zhs, epsi20 ) 
    243243                     ! 
    244244                     wfx_res(ji,jj) = wfx_res(ji,jj) + ( v_s(ji,jj,jl) - a_i(ji,jj,jl) * phs_max(ji,jj,jl) ) * rhos * r1_rdtice 
     
    255255                  zvs_excess = MAX( 0._wp, v_s(ji,jj,jl) - v_i(ji,jj,jl) * (rau0-rhoi) * r1_rhos ) 
    256256                  IF( zvs_excess > 0._wp ) THEN 
    257                      zfra = zvs_excess / MAX( v_s(ji,jj,jl), epsi20 ) 
     257                     zfra = ( v_s(ji,jj,jl) - zvs_excess ) / MAX( v_s(ji,jj,jl), epsi20 ) 
    258258                     wfx_res(ji,jj) = wfx_res(ji,jj) + zvs_excess * rhos * r1_rdtice 
    259259                     hfx_res(ji,jj) = hfx_res(ji,jj) - SUM( e_s(ji,jj,1:nlay_s,jl) ) * ( 1._wp - zfra ) * r1_rdtice ! W.m-2 <0 
Note: See TracChangeset for help on using the changeset viewer.