#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 |
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)
| Changeset | Author | Time | ChangeLog |
|---|---|---|---|
| 3808 | gm | 2013-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 |
| 3807 | gm | 2013-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 (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
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.
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