Opened 8 years ago

Closed 8 years ago

Last modified 3 years ago

#997 closed Bug (fixed)

Heat flux arrays not initialised in LIM3

Reported by: ufla Owned by: vancop
Priority: normal Milestone:
Component: LIM3 Version: release-3.4
Severity: Keywords: LIM*
Cc: ufla

Description

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)...)
#endif

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

Commit History (2)

ChangesetAuthorTimeChangeLog
3808gm2013-02-13T07:31:26+01:00

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

3807gm2013-02-13T07:29:43+01:00

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 (6)

comment:1 Changed 8 years ago by ufla

  • Cc ufla added

comment:2 Changed 8 years ago by vichi

  • Owner changed from NEMO team to vancop
  • Version changed from nemo_v3_3_1 to nemo_v3_4

comment:3 Changed 8 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 3 years ago by nemo

  • Keywords LIM* added

comment:5 Changed 3 years ago by nemo

  • Keywords release-3.4* added

comment:6 Changed 3 years ago by nemo

  • Keywords release-3.4* removed
Note: See TracTickets for help on using tickets.