New URL for NEMO forge!

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#997 (Heat flux arrays not initialised in LIM3) – NEMO

Opened 11 years ago

Closed 10 years ago

Last modified 16 months ago

#997 closed Bug (fixed)

Heat flux arrays not initialised in LIM3

Reported by: ufla Owned by: vancop
Priority: normal Milestone:
Component: LIM3 Version: v3.4
Severity: Keywords: LIM* v3.4
Cc: ufla


In lim_thd_dif, latent heat flux is updated as

qla_ice_1d(ji) = qla_ice_1d(ji) + dqla_ice_1d(ji) * (...)

However, the qla_ice_1d and dqla_ice_1d arrays are initialised in lim_thd from their 2d counterpart only in uncoupled mode:

#if ! defined key_coupled
  CALL tab_2d_1d(nbpb,qla_ice_1d(1:nbpb),qla_ice(:,:,jl)... )
  CALL tab_2d_1d(nbpb,dqla_ice_1d(1:nbpb),dqla_ice(:,:,jl)...)

which means that they are not initialised (but used) for coupled runs.

Commit History (2)


dev_MERGE_2012: #997 and #898: in LIM3, no use of qla_ice in coupled mode and its capping in forced mode


dev_v3_4_STABLE_2012: #997 and #898: in LIM3, no use of qla_ice in coupled mode and its capping in forced mode

Change History (7)

comment:1 Changed 11 years ago by ufla

  • Cc ufla added

comment:2 Changed 11 years ago by vichi

  • Owner changed from NEMO team to vancop
  • Version changed from v3.3 to v3.4

comment:3 Changed 10 years ago by gm

  • Resolution set to fixed
  • Status changed from new to closed

Thanks for pointing out this issue.

The current version of NEMO (and the forthcoming v3.4.1 stable and v3.5 alpha) does not include a coupled interface with LIM3. In particular, in sbcice_lim.F90, there is no specification of ice-atmosphere fluxes (case ksbc=5) as it is done in LIM2.

In order to just solve the qla_ice_1d issue reported here (i.e. we do not introduced a coupled interface with LIM3), the use of qla_ice_1d in both lim_thd_dif AND lim_thd_dh has been surrounded by "#if ! defined key_coupled" ... "#endif"

This mirror what has been done in LIM2: no update of sublimation inside LIM2 in coupled mode

Note that the specification of ice-atmosphere fluxes (call to sbc_cpl_ice_flx and set sprecip to minus emp_ice) is still missing in sbcice_lim. 

see Changeset 3808 for v3.5 alpha branch  and Changeset 3807 for the v3.4.1 stable branch

N.B. this changeset also include the Latent heat flux capping (see #898)

Gurvan and Martin

comment:4 Changed 6 years ago by nemo

  • Keywords LIM* added

comment:5 Changed 6 years ago by nemo

  • Keywords release-3.4* added

comment:6 Changed 6 years ago by nemo

  • Keywords release-3.4* removed

comment:7 Changed 16 months ago by nemo

  • Keywords v3.4 added
Note: See TracTickets for help on using tickets.