Hi!
We, at CCCma, are having a heat conservation issue in our coupling with out atmospheric model CanAM, and looking at the NEMO-SI3 side, we have a question:
In ice_thd_dh
, where hfx_sub
is computed, a comment says
! comment: not counted in mass/heat exchange in iceupdate.F90 since this is an exchange with atm. (not ocean)
hfx_sub_1d (ji) = hfx_sub_1d (ji) + ze_s(ji,jk) * zdum * a_i_1d(ji) * r1_Dt_ice ! Heat flux of snw that sublimates [W.m-2], < 0
but in ice_update_flx
the heat transfered to the ocean is computed using this hfx_sub
IF( ln_icethd ) THEN
qt_oce_ai(ji,jj) = qt_atm_oi(ji,jj) - hfx_sum(ji,jj) - hfx_bom(ji,jj) - hfx_bog(ji,jj) &
& - hfx_dif(ji,jj) - hfx_opw(ji,jj) - hfx_snw(ji,jj) &
& + hfx_thd(ji,jj) + hfx_dyn(ji,jj) + hfx_res(ji,jj) &
& + hfx_sub(ji,jj) - SUM( qevap_ice(ji,jj,:) * a_i_b(ji,jj,:) ) + hfx_spr(ji,jj)
ENDIF
Also, I was wondering if the latent heat of fusion of ice and snow is not already removed from Qns in the atmospheric side of the coupler as Qns includes the latent heat fluxes (evaporation over the ocean and sublimation over sea ice) as it is done in IPSL-CM6 (https://forge.ipsl.jussieu.fr/igcmg/wiki/IPSLCM6/Coupling)
Thanks in advance!